# generic commands for twographbox output size
load "twographbox.in"



set data style lines

set xlabel "Probability of state 1"
set ylabel "Utility" 
set xrange [0:1]
#set yrange [0:1]
set nokey

f(x,u1,u2) = (1-x)*u1 + x*u2

set output "alpha-vectors1.eps"

set label "[{/Times-Italic Stay}]" at 0.6,1.6
set label "[{/Times-Italic Go}]" at 0.15,1.15
set yrange [0:3]

plot \
  f(x,0.1,1.9)  title "(A2)" lw 2 lt 1, \
  f(x,0.9,1.1)  title "(A1)" lw 2 lt 1, \
  "alpha-vectors1.data" lw 6 lt 1, \
  "alpha-vectors1-verticals.data" lw 1 lt 2

set output "alpha-vectors2-all.eps"
set nolabel
set yrange [0:3]

plot \
   f(x,1.72,1.28) lw 3 lt 1, \
   f(x,1.48,1.68) lw 3 lt 1, \
   f(x,1.32,1.52) lw 2 lt 2, \
   f(x,1.08,1.92) lw 2 lt 2, \
   f(x,0.28,2.72) lw 3 lt 1, \
   f(x,0.68,2.48) lw 3 lt 1, \
   f(x,0.52,2.32) lw 2 lt 2, \
   f(x,0.92,2.08) lw 2 lt 2
 
set output "alpha-vectors2.eps"
set yrange [0:3]

plot \
  f(x,0.28,2.72)  lw 2 lt 1, \
  f(x,0.68,2.48)  lw 2 lt 1, \
  f(x,1.72,1.28)  lw 2 lt 1, \
  f(x,1.48,1.68)  lw 2 lt 1, \
  "alpha-vectors2.data" lw 6 lt 1, \
  "alpha-vectors2-verticals.data" lw 1 lt 2
 
set output "alpha-vectors8.eps"
set yrange [4.5:7.5]

plot \
  "alpha-vectors8.data" lw 6 lt 1
#  "alpha-vectors8-verticals.data" lw 1 lt 2
 
# with a1 0.3 0.7 (70% chance of changing state)
# ((0.59999937 1) (A2 (1 A2 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 3.28 5.72)
# ((0.4833333 0.59999937) (A2 (1 A1 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 3.568 5.528)
# ((0 0.39999974) (A1 (1 A2 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 4.36 4.64)
# ((0.39999974 0.4833333) (A1 (1 A2 (1 A2) (2 A2)) (2 A1 (1 A2) (2 A2))) 4.264 4.784)

# with a1 0.1 0.9 (90% chance of changing state)
# ((0.5 1) (A2 (1 A2) (2 A2)) 2.1 3.9)
# ((0 0.5) (A1 (1 A2) (2 A2)) 2.9 3.1)

# ((0.6 1)   (A2 (1 A2 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 3.28 5.72)
# ((0.5 0.6) (A2 (1 A1 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 3.664 5.464)
# ((0 0.4)   (A1 (1 A2 (1 A2) (2 A2)) (2 A2 (1 A2) (2 A2))) 4.72 4.28)
# ((0.4 0.5  (A1 (1 A2 (1 A2) (2 A2)) (2 A1 (1 A2) (2 A2))) 4.464 4.664))


# with 0,1 instead of 1,2
# ((0.5 1) (A2 (0 A2) (1 A2)) 0.1 1.9)
# ((0 0.5) (A1 (0 A2) (1 A2)) 0.9 1.1)
#
# ((0.6 1) (A2 (0 A2 (0 A2) (1 A2)) (1 A2 (0 A2) (1 A2))) 0.28 2.72 )
# ((0.5 0.6) (A2 (0 A1 (0 A2) (1 A2)) (1 A2 (0 A2) (1 A2))) 0.664 2.46)
# ((0 0.4) (A1 (0 A2 (0 A2) (1 A2)) (1 A2 (0 A2) (1 A2)))   1.72 1.28)
# ((0.4 0.5) (A1 (0 A2 (0 A2) (1 A2)) (1 A1 (0 A2) (1 A2))) 1.464 1.664)