【正文】
for j=1:n+1 for i=1:n+1 for k=1:m+1 A(j,i)=A(j,i)+x(k)^(j+i2) end endend。,39。)title(39。for j=1:n+1 for i=1:m+1 B(j)=B(j)+y(i)*x(i)^(j1) endendB=B39。y軸39。擬合圖39。a=inv(A)*B。title 39。x=::。*39。n=3。plot(x,z) legend(39。plot(x,y,39。hold onm=6。z=a(1)+a(2)*x+a(3)*x.^2+a(4)*