solving systems of linear equations with symmetric band matrix using cholesky decomposition