public double Det() // Function for determinant { if (L == null) MakeLU(); double det = detOfP; for (int i = 0; i < rows; i++) det *= U[i, i]; return det; }