ţ

来自「现在插值函数成为许多人的需要」· 代码 · 共 40 行

   以下是牛顿插值的源代码 ,希望能给大家提供参考，程序不大却很有用，因为在C语言中一般找不到这个程序，都是用MATLAB编写的          
         
    function newton_interpolation(X,Y,m)
               n=2;
               f=1;
               F=0;
               x=sym('x');
               s=1;
               S=0;
               for i=1:m
                      for j=1:i
        for k=1:i
            if k~=j
                f=f*(X(j)-X(k));
                    end
                      end
                       F=F+Y(j)/f;
                             f=1;
                             if j~=1
                        s=s*(x-X(j-1));
                          end
                         end
                       S=S+F*s;
                          s=1;
                          F=0;
                           end
                            S
                          pretty(S)
                          collect(S)
                          subs(S,x,0.13)
                             subs(S,x,0.2)