% initialize for testing phase2  (file is init.m)
% enter:   init , reg
pt=6  % pause time between iterations, if pt > 20, then a stop occurs
       %    and user hits return to continue
a=[1 2 10 16 1 0 0
   1.5 2 4 5 0 1 0
   .5 .6 1 2 0 0 1 ]
b=[800 1000 340]'
c=-[8 14 30 50 0 0 0]
bas=[5 6 7]
bas0=bas;c0=c;b0=b;a0=a;
