x=[0.03 0.19 0.25 0.40 0.44 0.52 0.69 0.78 0.87]
y=[2.87 2.44 2.25 1.82 1.68 1.44 0.980 0.630 0.43]
s=columns(x);
sx=0,sy=0,sxx=0,sxy=0,syy=0;
for n=1:s
sx+=x(n);
sy+=y(n);
sxx+=x(n)^2;
sxy+=x(n)*y(n);
syy+=y(n)^2;
endfor
delta=s*sxx-sx^2;
a=(s*sxy-sx*sy)/delta;
b=(sxx*sy-sx*sxy)/delta;
err=0;
printf("\nC1: %d, C0: %d\n \n", a,b);
for n=1:s
printf("%i. %d\n",n,a*x(n)+b-y(n));
err+=(y(n)-a*x(n)-b)^2;
endfor
err/=s;
printf("\nErr: %d\n", err);