牛顿插值法.m

来自「本程序是用MATLAB编写」· M 代码 · 共 23 行

function gjy00
clc
X=[1 1.5 0 2];
Y=[3 13/4 3 5/3];
[m,n]=size(X) ;
F(:,1)=Y'
temp = 1;
g = F(1,1);
syms x;
for i=2:n  
    for j=2:i    
        F(i,j)=(F(i-1,j-1)-F(i,j-1))/(X(i-j+1)-X(i))   
    end
end

for i = 1:n 
  
    temp = temp*(x-X(i));
    g = g+F(i,i)*temp
end