code: def dephasing_kraus_map(p=0.10):
      
      return [np.sqrt(1 - p) * np.eye(2), np.sqrt(p) * np.diag([1, -1])]
docs: Generate the Kraus operators corresponding to a dephasing channel.
  
      :params float p: The one-step dephasing probability.
      :return: A list [k1, k2] of the Kraus operators that parametrize the map.
      :rtype: list