code: def relaxation_operators(p):
      
      k0 = np.array([[1.0, 0.0], [0.0, np.sqrt(1 - p)]])
      k1 = np.array([[0.0, np.sqrt(p)], [0.0, 0.0]])
      return k0, k1
docs: Return the amplitude damping Kraus operators