function dtdy = ODEinv1 (t, y)
dtdy = t / (t^3 - 2*y) ;