try staring at this for an hour and you will redefine borenessfunction P=pivoteer(M,pivotrow,pivotcolumn)
element=M(pivotrow,pivotcolumn)
for i=1:length(M(pivotrow,:))
M(pivotrow,i) = M(pivotrow,i)/element
end
for i=1:length(M(:,1))
fact=M(i,pivotcolumn)
for j=1:length(M(1,:))
if i<>pivotrow
M(i,j)=M(i,j)-(fact*M(pivotrow,j))
end
end
end
P=M
endfunction