horner.m

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

function  p = horner(a, x)
% Horner's rule to compute
% p = a_1*x^n + ... + a_n*x + a_(n+1)

% Version 27.01.2004.  INCBOX

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