%Calculo de los coeficientes para la reconstruccion IC.
%   coefi=coeficientesipmcfrac(nl,r)
%   Entradas:
%       nl  = numero de nodos a la izquierda de f^{k-1}_{i} (vale para valores negativos)
%       nr  = numero de nodos a la derecha de f^{k-1}_{i}
%   Salidas:
%       ciefi  = vector de coeficientes buscados ( multiplicara a fr(k-1,i-nl:i+nl))

function coeficientesipmcfrac(nl,nr)
r=nl+nr+1;
nl=-nl;

%ojo nl viene ya negatvo
for j=nl:nr
   for l=2:r+1
      A(j-nl+1,l-1)=(j+1)^(l-1)-j^(l-1);
   end
end 

for k=1:r
   b(k)=1/2^(k-1);
end

coefi=b*A^(-1);

format rational
  coefi
format
sprintf('Comienza en f^{k-1}_{i+(%d)} y acaba en f^{k-1}_(i+(%d))',nl,nr)