set term svg size 800,600 dynamic enhanced

set termoption dashed

set style line 1 lt 1 lc rgb "black" lw 1
set style line 2 lt 2 lc rgb "black" lw 1
set style line 3 lt 3 lc rgb "black" lw 1
set style line 4 lt 4 lc rgb "black" lw 1
set style line 5 lt 5 lc rgb "black" lw 1
set style line 6 lt 1 lc rgb "black" lw 2

ismin(x,y) = (x>y)?y:x
ismax(x,y) = (x>y)?x:y

set xlabel 'Water Content / Water Content at 1.5MPa'
set ylabel 'UW*ts (m/s)'

# data from hagen 1996 (taken from graph of dynamic, static increase = (dynamic/0.8)-0.26
#0	0
#0.21129	0.112349
#0.222581	0.075608
#0.256452	0.179454
#0.280645	0.222175
#0.28871	0.179331
#0.309677	0.252645
#0.314516	0.090547
#0.354839	0.133207
#0.467742	0.154182
#0.470968	0.252029
#0.51129	0.257991
#0.548387	0.211977
#0.559677	0.392362
#0.562903	0.370943
#0.716129	0.312254

wucwts(wet_rat) = 0.58 * (exp(wet_rat) - 1.0 - 0.7*wet_rat*wet_rat)

#set log x
#set log y
set xrange [0:1]
set yrange [0:]
set key bottom right

# using decimal in numbers to ensure real division, not integer divistion
plot wucwts(x) ls 1 with lines title 'predicted' \
   , '-' ls 5 with points title 'Saleh and Fryear, 1995'
0	0
0.21129	0.112349
0.222581	0.075608
0.256452	0.179454
0.280645	0.222175
0.28871	0.179331
0.309677	0.252645
0.314516	0.090547
0.354839	0.133207
0.467742	0.154182
0.470968	0.252029
0.51129	0.257991
0.548387	0.211977
0.559677	0.392362
0.562903	0.370943
0.716129	0.312254
EOF

#pause -1