clc; clear all; close all; %% =========linprog()========= % min(z) =cX %st. AX <=b % Aeq <=beq; %%========================= c= [3, -1, -1]; A = [1,-2,1;4,-1,-1]; b = [1,-3]; Aeq = [-2,0,1]; beq = 1; x = linprog(c,A,b,Aeq,beq); %%================= c =[1,2,3,4,1,2,3,4]; Aeq=[1, -1,-1, 1, -1, 1, 1, -1; 1, -1, 1, -3, -1, 1, -1, 3; 1, -1,- 2, 3, -1, 1,-2, -3] beq = [0; 1; -1/2]; lb =zeros(8,1); uv0 = ones(8,1); options = optimset('LargeScale', 'off', ... 'Simplex', 'on','Diagnostics','on','MaxIter',1000); [uv,fval,exitflag,output,lambda]=linprog(c,[],[],Aeq,beq,lb,[],uv0,options);