clear;clc;
A=[1 2 1 0; 1 2/3 0 1];
b1=4;
b2=3;
B=[b1 b2];
b=[8;4];
c=[-1;-1;0;0];
Ab=A(:,B);

APr=((inv(Ab))*A);
bPr=(inv(Ab)*b);
TransponC=(transpose(c));
Cpomoc=TransponC(:,B);
CPr=(TransponC-(Cpomoc*(inv(Ab))*A));
% h=cbt b s pruhem
HPr=(Cpomoc*bPr);
matice=[APr bPr;CPr HPr];

disp(matice);

s=find(min(CPr));
disp(s);    clear;clc;
A=[1 2 1 0; 1 2/3 0 1];
b1=4;
b2=3;
B=[b1 b2];
b=[8;4];
c=[-1;-1;0;0];
Ab=A(:,B);

APr=((inv(Ab))*A);
bPr=(inv(Ab)*b);
TransponC=(transpose(c));
Cpomoc=TransponC(:,B);
CPr=(TransponC-(Cpomoc*(inv(Ab))*A));
% h=cbt b s pruhem
HPr=(Cpomoc*bPr);
matice=[APr bPr;CPr HPr];

disp(matice);

s=find(min(CPr));
disp(s);