function [p,q]=horner(pol,pto)
% *********************************************************************
% Esta función evalua un polinomio y su derivada usando el método de 
% Horner. El polinomio debe venir dado en forma vectorial.
%
% Las variables de entrada son:
% pol: Polinomio [an,...,a1,a0]
% pto: Punto donde evaluamos
% **********************************************************************

n=size(pol,2);
p=pol(1);
q=pol(1);
for k=2:n-1
   p=pol(k)+pto*p;
   q=p+pto*q;
end;
p=pol(n)+pto*p;