h = .02;
t(1) = 1;
y(1) = 2;
f(1) = 3;
N = 100;
for n=1:N
	t(n+1)=t(n) + h;
	y(n+1)= y(n) + f(n)*h;
	f(n+1)= y(n+1) + t(n+1);
end
plot(t,y,'o')
hold on
plot(t,y,'--')
stop
s = 1 : .01 : 3;
z = 4*exp(s-1)-s-1;
plot(s,z)