
% Linear inequalities
A=[];
b=[];
% Linear equalities
Aeq=[1 1 1];
beq=[850];
x0=[500 200 10];
% lower and upper bound
lb=[];
ub=[];
lb=[150 100 50];
ub=[600 400 200];
[xopt, fopt] = fmincon(@ELD,x0,A,b,Aeq,beq,lb,ub);
%[xopt,fopt,flag,somethinguseful,lagrange] = fmincon(@ELD,x0,A,b,Aeq,beq,lb,ub);