code: def DoubleExcitation(phi):
      r"""Double excitation rotation.
  
      Args:
          phi (float): rotation angle
      Returns:
          array: the four-qubit Givens rotation describing the double excitation
      """
  
      c = math.cos(phi / 2)
      s = math.sin(phi / 2)
  
      U = math.eye(16)
      U[3, 3] = c  # 3 (dec) = 0011 (bin)
      U[3, 12] = -s  # 12 (dec) = 1100 (bin)
      U[12, 3] = s
      U[12, 12] = c
  
      return U
  
imports: import pennylane
  from pennylane import numpy
  from pennylane import math
label_desc: Use image node_id 26 for calling a global function with example usage: DoubleExcitation(phi) and returns: U