function rhs = f(x,u)
rhs =  [u(2); -sin(u(1))];

