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 'Flat Biomass Soil Cover Fraction'
set ylabel 'UC*ts (m/s)'

# data from hagen 1996 (taken from graph of dynamic, static increase = (dynamic/0.8)-0.26
#-0.00031125
#0.03643375
#0.13655625
#0.310985
#0.4498325

wucsts(bffcv,sfcv) = (1 - exp(-1.2*bffcv))*(exp(-0.3*sfcv))

#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 wucsts(x,0.0) ls 1 with lines title 'predicted' \
   , '-' ls 5 with points title 'flat wheat straw'
0.0 0.0
0.0619637 0.03643375
0.140919 0.13655625
0.307202 0.310985
0.597938 0.4498325
EOF

#pause -1