public double[][] cholesky(double[][] matriz) {
        double temp[][] = new double[linhas][colunas];
        for (int i = 0; i < linhas; i++)  {
            for (int j = 0; j <= i; j++) {
                double soma = 0.0;
                for (int k = 0; k < j; k++) {
                    soma += temp[i][k] * temp[j][k];
                }
                if (i == j) temp[i][i] = Math.sqrt(matriz[i][i] - soma);
                else        temp[i][j] = 1.0 / temp[j][j] * (matriz[i][j] - soma);
            }
            if (temp[i][i] <= 0) {
                System.err.println("Diagonal Principal não é positiva");
            }
        }
        return temp;
    }