Begin

m,n,x,y

i=1 шаг 1 до n

j=1 шаг 1 до m

F(i,j)=x(i)^(j-1)

a=inv(F’*F)*F’*y

yr=F*a; dy=y-yr

r=dy’*dy

a,dy,r

close plot(x,y,’ko’) hold on plot(xt,yt,’k-’)

End

12