function f=L(x) f=3*x(1)^2+x(2)^4 -x(1)*x(2);