function y=phase2(x) % %Author: David Ivers [m,n]=size(x); m=m-1; y(:,1)=x(2:m,2); k=1; for j=3:n if x(m+1,j) < 0.5 k=k+1; y(1:m-1,k)=x(2:m,j); y(m,k)=j-2; end end t=displaytableau(y);