f=inline('y+t','t','y');
sol=inline('-t-1+2*exp(t)','t');
h=0.5;
[teuler,yeuler]=eul(f,[0,3],1,h)
eulerror=max(abs(sol(teuler)-yeuler))
ezplot(sol,[0,3])
hold on
plot(teuler,yeuler,'o')
hold off