function rhs = simple_rhs(t,y)
rhs = t^2 + y;