function rhs = f(t,u);
rhs = -2*t*u^2;