function [ p, pp ] = horner ( c, x )
%
%  HORNER evaluates a polynomial and its derivative at a point X.
%
n = size ( c, 2 );

pp = 0.0;
p = c(1);
for i = 2: n
  pp = pp * x + p;
  p = p * x + c(i);
end