hw1p1b

Contents

data
solve
display solution

% solve LP:
%       min 3*x1 + 2*x2
%        st   x1 -   x2 <= 3
%             x1 - 3*x2 <= 1
%          -4*x1 -   x2 <= 0
%                    x2 >= 0

data

c = [3; 2];
b = [3; 1; 0];
A = [1 -1;
     1 -3;
    -4 -1];
Lb = [-inf; 0];
Ub = [];

solve

x = linprog(c,A,b,[],[],Lb,Ub);

Optimization terminated.

display solution

fprintf('obj = %6.2f, solution: x = (%6.2f, %6.2f)\n',c'*x,x)

obj =   0.00, solution: x = (  0.00,   0.00)