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)
}