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 'Biomass Flat Cover Fraction'
set ylabel 'Emission Ratio'

# data from hagen 1996 (Figure 2)
# 0,1
#0.0600001,0.715376
#0.139554,0.425865
#0.308185,0.174216
#0.599777,0.0962463

renv(bffcv) = 0.075 + 0.934*exp(-bffcv/0.149)

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

# using decimal in numbers to ensure real division, not integer divistion
plot renv(x) ls 1 with lines title 'predicted' \
   , '-' ls 5 with points title 'Hagen, 1996'
0 1
0.0600001 0.715376
0.139554 0.425865
0.308185 0.174216
0.599777 0.0962463
EOF

#pause -1