newtonchaoz.m

来自「一些数值计算用的MATLAB代码,适合初学MATLAB和数字计算的人使用,根据M」· M 代码 · 共 26 行

n=input('次数：');
x=input('x[]=');
f=input('y[]=');
z=f;
i=2;
while i<=n
    for j=i:n
        f(j)=(z(j)-z(j-1))/(x(j)-x(j-i+1));
    end
    for k=i:n
        z(k)=f(k);
    end
    i=i+1;
end
u=input('要求值：');
i=2;
while i<=n
   for h=i:n
       f(h)=(u-x(i-1))*f(h);
   end
 i=i+1;
end
l=sum(f)