Prob9 para quem nao conseguiu adaptar a solução

16 janeiro 2009, 16:18 Célia Maria Santos Cardoso de Jesus

%prob9.m

x0 = [0 0 0 0 0];

[x,fval,exitflag,output,lambda,grad,hessian] = ...
fmincon(@objfun9,x0,[],[],[],[],[],[],@confuneq9)
[c,ceq] = confuneq9(x) % Check the constraint values at x
L=x(1)+x(2)+x(3)-210

function f = objfun9(x)

p1 = x(1); p2 = x(2); p3 = x(3);c1 = 3000+270*p1+0.3*p1^2;
c2 = 8000+270*p2+0.2*p2^2;
c3 = 6000+220*p3+0.25*p3^2;f = c1+c2+c3;