horner1.m

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

function  [p, dp] = horner1(a, x)
% Horner's rule used to compute p and p'

% Version 11.12.2003.  INCBOX

p = a(1);  dp = 0;
for  i = 2 : length(a)
  dp = dp*x + p;
  p = p*x + a(i);
end