eulersteps:=proc(f,x0,y0,h,N,T)
	local xold,yold,xnew,ynew,n:
	T[0]:=[x0,y0]:
	for n from 1 to N do
	    xold:=T[n-1][1]: yold:=T[n-1][2]:
	    xnew:=xold+h:
	    ynew:=yold+h*f(xold,yold):
	    T[n]:=[xnew,ynew]:
	end do:
	end:

#f:=proc(x,y) y: end:
#p:=plot(eulersteps(f,0,1,.1,10),color=blue);

#Q:=dsolve({diff(y(x),x)=y(x),y(0)=1},type=numeric,range=0..1);
#q:=odeplot(Q);

#plotsetup(x11);
#with(plots);
#display(p,q);