%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Université de Lille 1
%MAO 2004/2005
%
%         Représentation graphique d'un polynome d'interpolation
%                          sur l'intervalle [a,b]
%                               Animation
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function interpanim(NMAX,a,b)
%************* Points d'interpolation *********************


for N=2:NMAX
%N=nb de points d'interpolation
dxu=(b-a)/N; %pas d'espace

xu=a:dxu:b;%pts régulièrement répartis


dx=0.1*dxu;
x=a:dx:b;%pts pour la représentation des fonctions


%******************************************************************* 
%         construction des polynomes d'interpolation 
%                 pour la fonction f=exp(sin(pi*x))*cos(3*pi*x) 
%******************************************************************


            %+++++++++++ Polynome d'interpolation en les xu ++++++++++++++++
            
            fxu=exp(sin(pi*xu)).*cos(3*pi*xu);
            %calcul des coefficients du polynome d'interpolation polyfit
            cxu=polyfit(xu,fxu,length(xu)-1);
            %evaluation aux points d'interpolation x
            pxu=polyval(cxu,x);
            
            
            
            
            %++++ évaluation de la fonction aux pts de discrétisation x ++++
            
            f=exp(sin(pi*x)).*cos(3*pi*x);
            
            
%******************************************************************* 
%                    Représentations graphiques
%                
%*******************************************************************
            figure(1)
            plot(x,pxu,xu,fxu,'+',x,f)
            title('Interpolation animéé')
            xlabel('x')
            legend('Poly. Pts unif.','pts interp.','la fonction f',-2)
            drawnow
            pause(1)
        end
            
 