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])
