$\newcommand{\ket}[1]{|{#1}\rangle}\newcommand{\bra}[1]{\langle{#1}|}$Utilizando un entorno de programación, genere los estados $\delta_m^{BB}$, $m=0,1,2,3$:
\begin{align*} \delta_{0}^{BB} = \dfrac{1}{2} \left[\begin{smallmatrix} 1 & 0 & 0 & 1\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ 1 & 0 & 0 & 1 \end{smallmatrix}\right],\quad \delta_{1}^{BB} = \dfrac{1}{2} \left[\begin{smallmatrix} 0 & 0 & 0 & 0\\ 0 & 1 & 1 & 0\\ 0 & 1 & 1 & 0\\ 0 & 0 & 0 & 0 \end{smallmatrix}\right],\quad \delta_{2}^{BB} = \dfrac{1}{2} \left[\begin{smallmatrix} 0 & 0 & 0 & 0\\ 0 & 1 &-1 & 0\\ 0 &-1 & 1 & 0\\ 0 & 0 & 0 & 0 \end{smallmatrix}\right],\quad \delta_{3}^{BB} = \dfrac{1}{2} \left[\begin{smallmatrix} 1 & 0 & 0 &-1\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ -1 & 0 & 0 & 1 \end{smallmatrix}\right]. \end{align*}# Generamos los estados:
# Comprobamos que todos los productos cruzados son 0:
# Definimos el POVM:
# Aplicamos la medida: