function [t,u]=eulero(f,tI,tF,yI,N)
h=(tF-tI)/N;
t=[tI:h:tF];
u(1)=yI;
for n=1:N
    K=feval(f,t(n),u(n));
    u(n+1)=u(n)+h*K;
end