test_cycle

clear;
load cycle;
disp(A); disp(b'); disp(c')

[m,n] = size(A);
nbs = 1:n-m;                     % indices for nonbasics
bas = n-m+1:n;                   % indices for basics
opts.freq = 6;
opts.maxit = 60;
x = yzsimplex(A,b,c,bas,nbs,opts);

    0.5000   -5.5000   -2.5000    9.0000    1.0000         0         0
    0.5000   -1.5000   -0.5000    1.0000         0    1.0000         0
    1.0000         0         0         0         0         0    1.0000

     0     0     1

   -10    57     9    24     0     0     0

Iter    6:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   12:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   18:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   24:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   30:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   36:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   42:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   48:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   54:  Obj = 0.00000000e+00   basis: 5,4,7,
Iter   60:  Obj = 0.00000000e+00   basis: 5,4,7,
maxiter reached