Init with probabilities in a thermal state

One can set the initial state to

$$\rho_{ini}=\begin{pmatrix}p_1&0\\ 0&1-p_1\end{pmatrix}\otimes \cdots \otimes \begin{pmatrix}p_n&0\\ 0&1-p_n\end{pmatrix}$$ (4)

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

init-mixed $p_1$ $p_2$ ... $p_n$;

Example in the read-line mode:
  
        $>$init-mixed 0.7 0.6;
        $>$show;
       0.42$\vert$0$>$$<$0$\vert$+ 0.28$\vert$1$>$$<$1$\vert$+ 0.18$\vert$2$>$$<$2$\vert$+ 0.12$\vert$3$>$$<$3$\vert$
        $>$CN(1,2);
        $>$show;
       0.42$\vert$0$>$$<$0$\vert$+ 0.28$\vert$1$>$$<$1$\vert$+ 0.12$\vert$2$>$$<$2$\vert$+ 0.18$\vert$3$>$$<$3$\vert$