function f=fEXS710(t,U)
f=[U(2);
    -2*t*U(2)-U(1)^2];