function t = objfun3(invec)

x = invec(1);
y = invec(2);
z = invec(3);

t = -(x^2*y + 3*x*y*z + 4*y^2*sqrt(z));