function deriv = f(x,u);
deriv = [ u(2) 4*u(1)-4*x u(4) 4*u(3) ];
