线性规划大M法

clc;
A=[111100100;-120010010;-10200-1001];
b=[6;4;2];
c=[-2;1;-1;0;0;0;99999;99999;99999];
[mn]=size（A）;
x0=[0;0;0;0;0;0;6;4;2];
v=[789123456];
x0=x0（v:）;

while 1
    m_B=A（:v（1:m））;
    m_N=A（:v（m+1:n））;
    m_I=eye（n-mn-m）;
    d=[-inv（m_B）*m_N;m_I];
    cp=c（v:）;
    r=cp‘*d;
    if all（r>=0）
        display（sprintf（‘x is optimal‘））;
        cp=c（v:）;
        z=cp‘*x0;
        return;
    end
    [rqq]=min（r）