Fixed bug when [L U P Q] = lusol(A) was called as matrices P and Q were formed