【正文】
r i=1:n+1 for k=1:m+1 A(j,i)=A(j,i)+x(k)^(j+i2) end endend。B=[0 0 0 0]。for j=1:n+1 for i=1:m+1 B(j)=B(j)+y(i)*x(i)^(j1) endendB=B39。a=inv(A)*B。x=[::]。z=a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^3。plot(x,z) legend(39。離散點(diǎn)39。,39。y=a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^339。)title(39。擬合圖39。)