function [x,sol]=transport_ani(a,b,N,T,M)
close all;
clc;
k=T/M;
j=1;
h=(b-a)/(N+1);
la=k/h;
x=linspace(a,b,N+2);
x=x';
plot(x,g(x));
 title(['lambda=',num2str(la,2)]);
 axis([-10,10.1,-1,2]);
pause;
while j<M
    t=j*k;
    [x,sol]=upwind(a,b,N,T,M,t);
    plot(x,sol,x,g(x-t));
   % legend('Forw Euler','Back Euler','uexact');
    %plot(x,sol,'--',x,sol2,'-.');
    %legend('Forw Euler','Back Euler');
    title(['lambda=',num2str(la,2), '  t=',num2str(t,'%3.2f')])
    axis([-10,10,-1,2]);
     pause(0.1);
     j=j+1;
end
