function [T,Y]=euleri(fnimi,Tspan,ya,n)
a=Tspan(1); b=Tspan(2);
h=(b-a)/n;
Y=zeros(1,n+1);
T=a:h:b;
Y(1)=ya;
for j=1:n
   Y(j+1)=Y(j)+h*feval(fnimi,T(j),Y(j));
end;
