newcz.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 18 行

function [y,R]=newcz(X,Y,x,M)
n=length(X);m=length(Y);
for t=1;m
    z=x(t);A=zeros(n,n);A(:,1)=Y';
    s=0.0;p=1.0;q1=1.0;c1=1.0;
    for j=2:n
        for i=j:n
            A(i,j)=(A(i,j-1)-A(i-1,j-1))/(x(i)-x(i-j+1));
        end
        q1=abs(q1*(z-x(j-1)));c1=c1*j;
    end
    c=A(n,n);q1=abs(q1*(z-X(n)));
    for k=(n-1):-1:1
        c=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k);
    end
    Y(k)=polyval(C,z);
end
R=M*q1/c1;