def METHOD_NAME(self, t):
    psi, theta, phi = self.rotation_angles(t)
    c = qibo.models.Circuit(self.nqubits, density_matrix=True)
    # H gate
    c.add(qibo.gates.RZ(q=self.q, theta=np.pi / 2, trainable=False))
    c.add(qibo.gates.RX(q=self.q, theta=np.pi / 2, trainable=False))
    c.add(qibo.gates.RZ(q=self.q, theta=np.pi / 2, trainable=False))
    # RZ(psi)
    c.add(qibo.gates.RZ(q=self.q, theta=psi))
    # RX(theta)
    c.add(qibo.gates.RZ(q=self.q, theta=np.pi / 2, trainable=False))
    c.add(qibo.gates.RX(q=self.q, theta=-np.pi / 2, trainable=False))
    c.add(qibo.gates.RZ(q=self.q, theta=-theta))
    c.add(qibo.gates.RX(q=self.q, theta=np.pi / 2, trainable=False))
    # RZ(phi)
    c.add(qibo.gates.RZ(q=self.q, theta=phi))
    c.add(qibo.gates.M(self.q))
    return c