load ex5-2-6
N = [-1 -2 -3 4 -5]; B = [6 7];
v = lb; v(4) = ub(4);
d= sign(A*v-b);
A = [A -diag(d)]; w = [zeros(5,1); ones(2,1)]; 
lb = [lb; zeros(2,1)]; ub = [ub; inf*ones(2,1)];
[x,B,N] = rsmbdd(A,b,w,lb,ub,B,N)

x =

     2
     1
     0
     5
     0
     0
     0


B =

     1     2


N =

    -6    -7    -3     4    -5

w'*x

ans =

     0

ub(6:7) = zeros(2,1); p(6:7) = zeros(2,1);
[x,Bnew,Nnew] = rsmbdd(A,b,p,lb,ub,B,N)

x =

         0
    5.0000
    0.8571
    4.4286
         0
         0
         0


Bnew =

     3     4


Nnew =

    -6    -7    -1     2    -5

p'*x

ans =

   -28


lb = lb(1:5); ub = ub(1:5); A = A(:,1:5); p = p(1:5); N = N(3:5);
[x,Bnew,Nnew] = rsmbdd(A,b,p,lb,ub,B,N)

x =

         0
    5.0000
    0.8571
    4.4286
         0


Bnew =

     3     4


Nnew =

    -1     2    -5

p'*x

ans =

   -28