% Calculo de coeficientes para IC.
% En este caso contreto 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,nr+1}(1)-q_{-1,nr+1}(1/2)
% suponiendo que extendemos a partir de datos a la derecha
%   preliipmccoefisr(nr)
%   Entradas
%       nr = número de nodos que utilizamos a la derecha de f^{k}_{i+1}
%   Salidas
%       Se muestren por pantalla los coeficientes buscados

function preliipmccoefisrder(nr)
nr=nr+1; %porque nos desplazamos uno a la derecha  (es como usar nr nodos pero a la derecha de f^{k-1}_{i+1}
         % se hace esto para que en la gui se mantenga el orden correcto
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   

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