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 'Ln(Crushing Energy (J/kg))'
set ylabel 'Abrasion Coefficient (1/m)'

canag(seags) = exp(-2.07-0.077*seags**2.5-0.119*log(seags))

#set log x
#set log y
#set xrange [0.0005:5]
#set yrange [0:3100]
set key top right

# using decimal in numbers to ensure real division, not integer divistion
plot canag(x) ls 1 with lines title 'Predicted' \
   , '-' using 1:2 ls 4 title 'Crust' with points \
   , '-' using 1:2 ls 5 title 'Aggregates' with points
0.022015 0.16892
0.106933 0.146202
0.142053 0.151862
0.227426 0.171808
0.306315 0.153778
0.309271 0.131089
0.638005 0.136155
0.717516 0.161774
EOF
1.57767 0.0788634
1.64349 0.0642908
1.68755 0.082278
1.69354 0.0773438
1.72613 0.068208
1.80729 0.0634944
1.83524 0.0671839
1.92357 0.0644407
1.97311 0.0745346
1.99301 0.0710752
2.11558 0.0688134
2.30644 0.066528
2.38866 0.0679793
2.46333 0.0652408
2.55676 0.0523848
2.57388 0.0326502
2.57637 0.0471993
2.61245 0.0585308
2.69324 0.0515979
2.69931 0.0471568
2.72516 0.0385166
2.74198 0.0570065
2.74805 0.0525654
2.77896 0.0335658
2.83478 0.0404516
2.89772 0.0490612
2.91678 0.0406699
2.92176 0.0298174
2.94254 0.0315365
2.98483 0.0391668
3.00975 0.0251014
3.05124 0.0280464
3.0528 0.0371704
3.09264 0.0304982
3.26226 0.0237812
3.30969 0.0215454
3.43255 0.0210097
3.44486 0.013114
3.4869 0.0192647
3.55326 0.00789776
3.58096 0.0101077
3.61026 0.0216882
3.67001 0.0115566
3.70382 0.00957205
3.73717 0.00487496
3.84013 0.00779872
3.89443 0.0058071
3.96299 0.00726309
3.99655 0.00379897
4.05077 0.00131416
4.09251 0.00573871
4.21537 0.00520308
4.31122 0.0054166
4.33089 0.00171067
4.39312 0.00612814
4.41973 0.0019266
4.50176 0.00189828
4.63151 0.0021001
4.75452 0.00205763
4.79128 0.00251075
4.87065 0.00251075
4.92516 0.0017521
EOF

#pause -1

#Hagen-Skidmore-Saleh_1992_crust
#ln(crushing energy J/kg), Abrasion Coefficient (1/m)
#0.022015,0.16892
#0.106933,0.146202
#0.142053,0.151862
#0.227426,0.171808
#0.306315,0.153778
#0.309271,0.131089
#0.638005,0.136155
#0.717516,0.161774

#Hagen-Skidmore-Saleh_1992_aggregates
#ln(crushing energy J/kg), Abrasion Coefficient (1/m)
#1.57767,0.0788634
#1.64349,0.0642908
#1.68755,0.082278
#1.69354,0.0773438
#1.72613,0.068208
#1.80729,0.0634944
#1.83524,0.0671839
#1.92357,0.0644407
#1.97311,0.0745346
#1.99301,0.0710752
#2.11558,0.0688134
#2.30644,0.066528
#2.38866,0.0679793
#2.46333,0.0652408
#2.55676,0.0523848
#2.57388,0.0326502
#2.57637,0.0471993
#2.61245,0.0585308
#2.69324,0.0515979
#2.69931,0.0471568
#2.72516,0.0385166
#2.74198,0.0570065
#2.74805,0.0525654
#2.77896,0.0335658
#2.83478,0.0404516
#2.89772,0.0490612
#2.91678,0.0406699
#2.92176,0.0298174
#2.94254,0.0315365
#2.98483,0.0391668
#3.00975,0.0251014
#3.05124,0.0280464
#3.0528,0.0371704
#3.09264,0.0304982
#3.26226,0.0237812
#3.30969,0.0215454
#3.43255,0.0210097
#3.44486,0.013114
#3.4869,0.0192647
#3.55326,0.00789776
#3.58096,0.0101077
#3.61026,0.0216882
#3.67001,0.0115566
#3.70382,0.00957205
#3.73717,0.00487496
#3.84013,0.00779872
#3.89443,0.0058071
#3.96299,0.00726309
#3.99655,0.00379897
#4.05077,0.00131416
#4.09251,0.00573871
#4.21537,0.00520308
#4.31122,0.0054166
#4.33089,0.00171067
#4.39312,0.00612814
#4.41973,0.0019266
#4.50176,0.00189828
#4.63151,0.0021001
#4.75452,0.00205763
#4.79128,0.00251075
#4.87065,0.00251075
#4.92516,0.0017521