function [temps,Y] = EULEREXPLICITE(t0,y0,T,N,f) 
h = T/N ; temps = [t0]; Y = [y0]; 
for k = 1:N 
  y      = Y(k) + h*f(temps(k),Y(k)); 
  temps  = [temps,temps(k)+h]; Y  = [Y,y];
end    
endfunction   
