func showLU(a *mat.DenseMatrix) { fmt.Printf("\na:\n%v\n", a) l, u, p := a.LU() fmt.Printf("l:\n%v\n", l) fmt.Printf("u:\n%v\n", u) fmt.Printf("p:\n%v\n", p) }