====== CalculateRho ======

###
//CalculateRho($H$)// computes the density matrix for a given Tight Binding model.

###

===== Input =====

  * $H$ : Tight Binding Object
  * Possible options are:
    * "Nk" List of integers the format {N$_x$,N$_y$,N$_z$} giving the number of Crystal momentum states used to solve the model. If the reciprocal lattice is zero, then Nk is always set to {0,0,0}, independent of the option given (Standard value {40,40,40})
    * "mu" Real number giving the chemical potential $\mu$. (Standard value 0)
    * "kbT" or "kBT" Positive real number indicating the value of $\beta=k_BT$ for the calculation. (Standard value DBL_EPSILON)
   

===== Output =====

  * $\rho$ : Matrix object with size N$_x$*N$_y$*N$_z$*N$_{orbs}$, where N$_{orbs}$ is the number of defined orbitals per unit cell in the Tight Binding input.  

===== Example =====

###
description text
###

==== Input ====
<code Quanty Example.Quanty>
-- some example code
</code>

==== Result ====
<file Quanty_Output>
text produced as output
</file>

===== Table of contents =====
{{indexmenu>.#1}}