function [u]=euler1(neqn,t0,tf,u0,npts)
h=(tf-t0)/npts;
for i=1:npts
[ut0]=derv(neqn,t0,u0);
u=u0+ut0*h;
t=t0+h;
u0=u;
t0=t;
end
