% Filtros para la reconstruccion IAP en el caso de estar realizando SR y querer extender a partir de datos a la izquierda del intervalo [x_{i-1},x_i]
%     solu=capiaiavpcoefiesquerrasr(nl,nr,r)
%     ENTRADAS:
%           nl+nr = tamaño del stencil
%           r     = grado del polinomio IA
%     SALIDAS:
%           solu  = coeficientes que multiplican a f^{k-1} para obtenener f^{k}_{2i-1} (también se muestran por pantalla)
%
% Autor: José Jaime Noguera.  +INFO: www.uv.es/anims/noguera

function solu=capiaiavpcoefiesquerrasr(nl,nr,r)

A=matriuiap(nl+nr-1,1,r);
C=(A'*A)^(-1)*A';
b0=(nl+nr-1:-1:2)';
b1=(-nl-nr+2:-1)';

for l=3:r+1
   coef1(l)=(3/2)^(l-1);
end 

multfi_1=(3/2)+(3/2)*ones(1,r-1)*C*b1-coef1(3:r+1)*C*b1;
multfi_2=1-(3/2)+(3/2)*ones(1,r-1)*C*b0-coef1(3:r+1)*C*b0;
multfi2_3fin=-(3/2)*ones(1,r-1)*C+coef1(3:r+1)*C;

solu=[multfi2_3fin , multfi_2, multfi_1];

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