x=polyregress(:,1);
y=polyregress(:,2);

n=length(x);
sumx=sum(x);
sumx2=sum(x.^2);
sumx3=sum(x.^3);
sumx4=sum(x.^4);
sumy=sum(y);
sumxy=sum(x.*y);
sumx2y=sum((x.^2).*y);

A=[n sumx sumx2;sumx sumx2 sumx3;sumx2 sumx3 sumx4]
B=[sumy sumxy sumx2y]'
ai=inv(A)*B