% Calculo de coeficientes para IC
% En este caso contreto es para los coeficientes
% de f^{k}_{2i-1} en el caso SR. Lo que calcula el programa es f^{k}_{2i} fruto de q_{-1,3}(1)-q_{-1,3}(1/2)
% suponiendo que extendemos a partir de datos a la derecha. Despues lo
% modificamos para que muestre la salida para la izquierda, ya que se
% pueden deducir unos de otros.
%   solu=preliipmccoefisrGUIder(nr)
%   Entradas
%       nr = número de nodos que utilizamos a la izquierda de f^{k-1}_{i-1}
%   Salidas
%      Coeficientes buscados

function preliipmccoefisrizq(nl)
nr=nl+1; %truco para no reprogramar % a parte cogemos una mas porque es como coger nl a la izquierda de f^{k-1}_{i-1}
nl=-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=1:nl+nr+1
   b(i)=2-(1/2)^(i-1);
end   

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



