% Calculo de coeficientes para la función G en la reconstruccion IC.
%   preliipmcfunciong(r,punt)
%   Entradas
%       r = grado reconstruccion para F (r-1 para f)
%       punt = punto donde se sutituye la funcion G (0 para x^{k-1}_{i-1}  ,  0.5 para x^{k}_{2i-1},  1 para x^{k-1}_{i})
%   Salidas
%      Coeficientes buscados

function preliipmcfunciong(r,punt)
nl=-1;
nr=r+1; %para que nl+nr=1;
for j=-nl:nr 
   for l=2:nl+nr+2
      A(j+nl+1,l-1)=((j+1)^(l-1)-j^(l-1));   
   end   
end   

for i=2:nl+nr+2
   b(i-1)=-1+(punt)^(i-1);
end   

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

nl=r+1; %para que nl+nr=1;
nr=-1;
for j=-nl:nr 
 
   for l=2:nl+nr+2
      A(j+nl+1,l-1)=((j+1)^(l-1)-j^(l-1));   
   end
   
end   


for i=2:nl+nr+2
   b(i-1)=-(punt)^(i-1);
end


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

sprintf('Hay que sumar f^{k-1}_i')
