intpval.m

来自「用matlab写的一些数值算法」· M 代码 · 共 10 行

function  p = intpval(x,c,t)
% Value p = P(t) of interpolating polynomial

% Version 11.12.2003.  INCBOX

m = length(x);  % number of interpolation points and coeffs
p = c(m) * ones(size(t));
for  k = m-1 : -1 : 1
  p = p .* (t - x(k)) + c(k);
end