function E = rosenbrock(u)
%x=[x,y]
Du(x)=gradient(u(x))
E = sum(int(f(x)).*u,int(g(x)).* abs(Du(x)))