clear all; close all; clc; %1 %a-b c=[40; 70]; A=[ 1, 1; 10, 50]; b=[100; 4000]; Aeq=[]; beq=[]; [x,V]=linprog(-c,A,b,Aeq,beq,[0;0]) [y,V2]=linprog(b,-A',-c,Aeq,beq,[0 0]) %c-d N=10001; M=1; for i=M:100:N b2=[100; i]; [x2,V3]=linprog(-c,A,b2,Aeq,beq,[0;0]); g(((i-1)/100)+1)=(V-V3)/(i-4000); end figure(1) plot(M:100:N,g) %f c=[40; 70]; A=[ 1, 1; 10, 50]; b=[100; 5001]; Aeq=[]; beq=[]; [x,V]=linprog(-c,A,b,Aeq,beq,[0;0]); [y,V2]=linprog(b,-A',-c,Aeq,beq,[0 0]); N=10001; M=1; for i=M:100:N b2=[100; i]; [x2,V3]=linprog(-c,A,b2,Aeq,beq,[0;0]); g(((i-1)/100)+1)=(V-V3)/(i-5000); end figure(2) plot(M:100:N,g) %2 %a-b c=[2; 0; 1]; A=[ -5, 2, 0; 2, -3, 2; 1, 1, -3; 1, 0, 0; 0, 1, 0; 0, 0, 1]; b=[0; 0; 0; 1; 1; 1]; Aeq=[]; beq=[]; [x,V]=linprog(-c,A,b,Aeq,beq,[0;0;0]) [y,V2]=linprog(b,-A',-c,Aeq,beq,[0 0 0 0 0 0])