Fixed a minor glitch with calculation of P, Q. Now adj(P).P gives a factor of identity matrix