%state limit function for problem # 2
%y = g(x,u)
function y = g(x,u)
v(1) = 400*u(2) + 2000;
v(2) = 250*u(3) + 1000;
v(3) = 100*u(1) + 500;
v(4) = exp(0.7386 + 0.0998*u(4));
y = 1 - 4*v(1)/(x(1)*x(2)^2*v(4)) - 4*v(2)/(x(1)^2*x(2)*v(4))...
                                  - (v(3)/(x(1)*x(2)*v(4)))^2;

y = -y;