LU decomposition

Any square matrix A (generically, up to a permutation P) factorises as A = PLU where L is unit lower-triangular and U is upper-triangular. Computed by Gaussian elimination in O(n³) operations; reduces solving Ax = b to two triangular…