function [coc,rest]=horner(p,pto)

% ************************************************
% Esta función calcula el cociente y 
% y el resto obtenidos al dividir un polinomio 
% por (x-pto) usando el método de 
% Horner. 
% El polinomio debe venir dado en forma vectorial.
%
% Las variables de entrada son:
% p: Polinomio [an,...,a1,a0]
% pto: Punto donde evaluamos
% Las variables de salida son: 
% coc : cociente [bn-1,...,b1,b0]
% rest: resto
% *************************************************

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

coc=q(1:n)
rest=q(n+1)