%s Sound speed vs solid fraction due to Atkinson, %Hampton, Urick. 
xmin=0;
delx=0.65;
ymin=1400;
dely=400;
xwidth=6.42;		%inches
ywidth=6.37;		%inches
xmax=xmin+delx;
ymax=ymin+dely;
V=[xmin xmax ymin ymax];
axis(V);
xdot=[.41
.9875
1.28
1.506
1.861
2.276
2.689
3.324
3.677
3.779
4.152
4.642
4.71
5.219
5.386
5.61
5.953
5.995];
ydot=[1.581
1.679
1.77
1.76
1.837
1.95
2.12
2.465
2.121
2.08
2.361
2.573
2.3
2.617
3.348
3.727
4.277
4.4];
xstar=[.898
.915
1.818
2.046
2.078
2.84
3.01
3.11
3.66
3.746
3.943
4.008
4.757
4.93
5.05
6.12
6.12
6.12];
ystar=[1.735
1.797
2.1
2.281
1.952
1.952
2.28
2.396
2.206
2.466
3.138
3.097
3.594
2.354
3.636
5.457
5.51
5.61];
xtriang=[.913
.967
1.992
2.932
3.655
3.786
4.69
5.054
6.12];
ytriang=[1.589
1.797
1.95
2.216
2.47
2.85
3.255
3.33
4.9];
xurick=[.3035
.646
.648
1.12
1.13
1.789
1.798
2.362
2.597
2.67
3.02
3.43];
yurick=[1.415
1.285
1.26
1.07
1.08
.92
.9
.969
.989
1.013
1.132
1.273];
xhampton=[.285
.4275
.652
.947
1.052
1.229
1.462
1.555
1.654
1.93
2.146
2.35
2.464
2.53
2.583
2.72
2.985
3.305];
yhampton=[1.47
1.384
1.31
1.158
1.086
1.042
.994
.948
.897
.848
.898
.824
.8675
.821
.872
.89
.9565
1.006];
xdot=xmin+(delx/xwidth)*xdot; 
ydot=ymin+(dely/ywidth)*ydot;
xstar=xmin+(delx/xwidth)*xstar; 
ystar=ymin+(dely/ywidth)*ystar;
xtriang=xmin+(delx/xwidth)*xtriang; 
ytriang=ymin+(dely/ywidth)*ytriang;
xurick=xmin+(delx/xwidth)*xurick; 
yurick=ymin+(dely/ywidth)*yurick;
xhampton=xmin+(delx/xwidth)*xhampton; 
yhampton=ymin+(dely/ywidth)*yhampton;
plot(xtriang,ytriang,'x',xdot,ydot,'o',xstar,ystar,'+',xurick,yurick,'*',xhampton,yhampton,'+r')