invint.m

来自「Interpolation routines in matlab」· M 代码 · 共 29 行

%
% Inverse Polynomial Interpolation
%

clear;
n = 5;
x = -1:2/(n-1):1;
for i = 1:n
  f(i) = 1-exp(x(i)^2);
end
z = [ 0 ];


p = AitkenNeville( f, x, z );


clear;
a = 1.6;
b = 2.1;
for n = 2:2:6
    h = (b-a)/(n-1);
    x = (a:h:b);
    f = cos(x).*cosh(x)+ones(size(x));
    z = [ 0 ];

    p = AitkenNeville( f, x, z );

end