Init with biases in a thermal state

One can set the initial state to

$$\rho_{ini}=\begin{pmatrix}\frac{1+\varepsilon_1}{2}&0\\ 0&\frac{1... ...{pmatrix}\frac{1+\varepsilon_n}{2}&0\\ 0&\frac{1-\varepsilon_n}{2}\end{pmatrix}$$ (5)

( $0\le \varepsilon_i\le 1$) with the command

init-thermal-biases $\varepsilon_1$ $\varepsilon_2$ ... $\varepsilon_n$;