function y=Horner(d,c,x,b)
% Usando la regla de Horner se evalúa el polinomio de grado d,
% definido en el vector c, en el punto x (puede ser vector).
% Entrada: c(d+1), coeficientes del polinomio: c1, c2,...;
% x, punto (vector) en el que se evalúa el polinomio;
% b, puntos base, si es el caso.
% Salida: y, valor del polinomio en x
 if nargin<4, b=zeros(d,1); end 
 y=c(d+1);
 for i=d:-1:1
    y = y.*(x-b(i))+c(i);
 end
end
