【正文】
*x.^2+a(4)*x.^339。y=[,]。x軸39。散點(diǎn)圖39。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。x=[::]。,39。)。)title(39。plot(x,z) legend(39。for j=1:n+1 for i=1:m+1 B(j)=B(j)+y(i)*x(i)^(j1) endendB=B39。n=3。y軸39。*39。擬合圖39。離散點(diǎn)39。a=inv(A)*B。A=zeros(n+1