function [Y,t] = euleroesp(lam,fun,a,b,y0,N)

%  [Y,t] = euleroesp(lam,fun,a,b,y0,N)
%
%  Applica il metodo di EULERO ESPLICITO per 
%  approssimare la soluzione del problema ai 
%  valori iniziali
%
%  y'(t) = lam*y(t) + fun(t,lam),        
%  y(a) = y0,    a<=t<=b
%
%  con fun(t,lam) = -lam*g(t) - g'(t).
%
%  N rappresenta il numero di punti della
%  discretizzazione dell'intervallo di integrazione.


Y   = y0;

h   = (b-a)/N;

t  = a + [0:N]*h;

for i =1:N
    
    fi = feval(fun,t(i),lam);
    
    k = lam*y0 +fi;
    
    y0 = y0 + h*k;
    
    Y(i+1)=y0;
    
end    