function dydt = ODEexp1 (t, y)
dydt = (t^3 - 2*y) / t;