function [] = nonAutonomos()
clf
hold on

tspan = [-3/2 3/2];

plotSolutionNA

y0 = -1.4;
h = 1/4;
for k = 1:4
	[Y,T] = expEuler(@f,tspan,y0,h);
	h = h/2
	plot(T,Y,'ro-','LineWidth',2)
	axis([-1.5 1.5 -1.5 1.5])
end
shg

function [dy] = f(t,y)
dy = t^2+y^2;