func demo(m *mat.DenseMatrix) {
	fmt.Println("A:")
	fmt.Println(m)
	l, err := m.Cholesky()
	if err != nil {
		fmt.Println(err)
		return
	}
	fmt.Println("L:")
	fmt.Println(l)
}