function [b, c] = Division_Sintetica_Derivadas(a, r, s)

n = length(a);

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