code: def max_entangled_circuit(wires, shots=10000, interface="autograd"):
      """maximally entangled state preparation circuit"""
      dev = qml.device("default.qubit", wires=wires, shots=shots)
  
      @qml.qnode(dev, interface=interface)
      def circuit():
          qml.Hadamard(wires=0)
          for i in range(1, wires):
              qml.CNOT(wires=[0, i])
          return qml.classical_shadow(wires=range(wires))
  
      return circuit
  
imports: import pytest
  import numpy
  import pennylane
  import pennylane.numpy
  from pennylane.shadows import ClassicalShadow, median_of_means, pauli_expval
label_desc: Use image node_id 3 for calling a global function with example usage: max_entangled_circuit(wires, shots, interface) and returns: circuit, qml