function [c, r] = Division_Sintetica2(a, r, s)
% Calcula la division Sintetica de un polinomio de grado dos
% f(x) = a^T x y g(x) = x^2 + rx + s. 
% Recibe como paramentros los coeficientes a, r y s y regresa
% el polinomio resultante c y el residuo r

n = length(a);

b=zeros(n,1);
 
b(n) = a(n);
b(n-1) = a(n-1) - r*b(n);

for k=n-2: -1:1
      b(k) = a(k) - r*b(k+1) - s*b(k+2);
end;

fprintf(1,'Residuo %6.4f x + %6.4f\n', b(2), b(1)+r*b(2)); 
b
c = b(3:n);
r = [b(2); b(1)+r*b(2)];
end

