code: def depolarizing_operators(p):
      
      k0 = np.sqrt(1.0 - p) * I
      k1 = np.sqrt(p / 3.0) * X
      k2 = np.sqrt(p / 3.0) * Y
      k3 = np.sqrt(p / 3.0) * Z
      return k0, k1, k2, k3
docs: Return the phase damping Kraus operators