% Filtros para la reconstruccion IAP en el caso de estar realizando SR y querer extender a partir de datos a la derecha 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=capiaiavpcoefidretasr(nl,nr,r)

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

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

multfi=1+1/2-(1/2)*ones(1,r-1)*C*b0-coef1(3:r+1)*C*b0;
multfi1=-1/2-1/2*ones(1,r-1)*C*b1-coef1(3:r+1)*C*b1;
multfi2fin=(1/2)*ones(1,r-1)*C+coef1(3:r+1)*C;

solu=[multfi , multfi1, multfi2fin];

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