code: def mixed_density(qubits: Union[int, Qubits]) -> Density:
      
      N, qubits = qubits_count_tuple(qubits)
      matrix = np.eye(2**N) / 2**N
      return Density(matrix, qubits)
docs: Returns the completely mixed density matrix