function [x,sol]=feuler_ani(a,b,N,T,M)
close all;
clc;
k=T/M;
j=1;
h=(b-a)/(N+1);
la=k/h^2;
x=linspace(a,b,N+2);
x=x';
plot(x,g(x));
 title(['lambda=',num2str(la,2)]);
 axis([-0.1,1.1,-1.1,1.1]);
pause(0.1);
while j<M
    t=j*k;
    [x,sol]=feuler(a,b,N,T,M,t);
    [x,sol2]=beuler(a,b,N,T,M,t);
    plot(x,sol,'--',x,sol2,'-.',x,uexact(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([-0.1,1.1,-1,1]);
     pause(0.1);
     j=j+1;
end
