Z=csvread('Cereals no alpha.csv');
sodium=Z(:,4);
fiber=Z(:,5);
s=10.3492;
pptdist=1.9926;
[db1,db2,dypbar,dypran] = inference2D(sodium,fiber,s,pptdist);
disp([db1,db2,dypbar,dypran])
b0=45.2616;
bs=-0.0606;
bf=3.2923;
yp=b0+200*bs+5*bf;
disp(yp)
