x=experiment(:,1);
y=experiment(:,2);
sumx=sum(x);
sumy=sum(y);
sumsquaredx=sumx^2;
xy=x.*y;
sumxy=sum(xy);
xsquared=x.^2;
sumxsquared=sum(xsquared);
averagex=mean(x);
averagey=mean(y);
n=length(x);
a1=(n*sumxy-sumx*sumy)/(n*sumxsquared-sumsquaredx)
a0=averagey-a1*averagex
ymodel=a0+a1*x