function y=pivot(ix,jx,x) % %Author: David Ivers if x(ix,jx) == 0 error('Pivot element is zero.') end; y(ix,:) = x(ix,:)/x(ix,jx); [m,n]=size(x); m=m-1; for k=1:ix-1 y(k,:) = x(k,:) - x(k,jx)*y(ix,:); end; for k=ix+1:m y(k,:) = x(k,:) - x(k,jx)*y(ix,:); end; y(:,jx) = zeros(m,1); y(ix,jx) = 1; y(m+1,:) = x(m+1,:); t=displaytableau(y);