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 'Soil Clay Fraction'
set ylabel 'Suspension Fraction from Clod Abrasion'

sfssan(sfcla) = 0.4 - 4.83*sfcla + 27.18*sfcla**2 - 53.7*sfcla**3 + 42.25*sfcla**4 - 10.7*sfcla**5

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

# using decimal in numbers to ensure real division, not integer divistion
plot sfssan(x) ls 1 with lines title 'Predicted' \
   , '-' using 1:2 ls 4 title 'Data' with points
0.055 0.2303
0.087 0.1378
0.227 0.1602
0.284 0.2704
EOF

#pause -1