function f = fcurl(x,u)
f = [u(2) 0.5*cos(x).*sqrt(1+u(2).^2).^3]';