densityMatrixSolver Calculates system expectation values from the Density Matrix of a given Hamiltonian, including superconducting Order Parameter