function r = rhs(t,y)

r = y-t^2+1;
