$\newcommand{\ket}[1]{|{#1}\rangle}\newcommand{\bra}[1]{\langle{#1}|}$Considere el estado dado por la matriz de densidad de probabilidad $\delta_1=\ket{00}\bra{00}$.
import numpy as np
# Definimos las puertas identidad, H y CNOT:
# y el estado ket00: