function [t,u]=assolutastab(lambda,h,met,T)
N=ceil(T/h);
t=zeros(1,N+1);
u=zeros(1,N+1);
t(1)=0;
u(1)=1;
switch met
    case 1
        for n=1:N    
           t(n+1)=t(n)+h;    
           u(n+1)=u(n)*(1+h*lambda);
        end
    case 2
        for n=1:N
            t(n+1)=t(n)+h;
            u(n+1)=u(n)*(2+h*lambda)/(2-h*lambda);
        end
   case 3
       t(2)=t(1)+h;
       u(2)=exp(lambda*h);
       for n=2:N
            t(n+1)=t(n)+h;
            u(n+1)=u(n-1)+2*h*lambda*u(n);
       end
end

