mynewton.m

来自「自己编写的matlab程序 数值分析里的程序」· M 代码 · 共 22 行

function X=myNewton(A,B,C)
syms x;
m=length(A);
n=length(B);
if m~=n
    'error'
end
for i=2:n
    for j=i:n
        B(i,j)=(B(i-1,j)-B(i-1,j-1))/(A(j)-A(j-i+1));
    end
end
Y=B(1,1);
X=sym('0');
t=sym('1');
for i=2:n
    t=t*(x-A(i-1)); 
    X=X+B(i,i)*t;
end
X=simple(X);
X=X+B(1,1);