function twoD = twoders(t,y)

twoD = [y-t^2+1 y-t^2+1-2*t];
