function?[CRBANS]=AZ（abc）
ST=zeros（1length（b））;%⊙seita值
b=b（:）;%基解
times=0;%debug用的，檢測循環過幾次
flag=1;%用于標記是否繼續迭代的值
ZH=0;%換入基變量列和換出基變量行的相交值
H=0;%換入基變量所在的航
L=0;%換入基變量所在的列
A=[a?eye（length（b））];%基變量和非基變量系數
C=[c?zeros（1length（b））];%沒什么好說的
g=zeros（length（C）1）;%A和F的積
f=C;%單純性法中的-z檢驗數
F=zeros（length（b）1）;%單純性的CR有的地方叫CB我叫他CR
RES=0;%MAXZ結果值
while?（flag）%進行無限迭代（后面有寫跳出條件）
times=times+1;%每次迭代times+1
%%%%%%%確定-Z%%%%%%%%%
for?i=1:length（C）
????for?ii=1:length（F）
g（i）=g（i）+A（iii）*F（ii）;
????end??
f（i）=C（i）-g（i）;
end
g=zeros（length（C）1）;%把臨時變量g歸零
%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%找出A中換入基變量的列數%%%%%
CHU=max（f）;%找到-Z中最大的數
for?i=1:length（C）
????for?ii=1:length（F）
g（i）=g（i）+A（iii）*F（ii）;
????end??
f（i）=C（i）-g（i）;
????if（f（i）==CHU）
????L=i;
????end
end
g=zeros（length（C）1）;%把臨時變量g歸零
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%確定換出基變量%%%%%%%%%
for?i=1:length（F）
????ST（i）=b（i）/A（iL）;
????if（ST（i）<0）
????????ST（i）=inf;
????end
end
RU=min（ST）;%找到seit