code: def OrbitalRotation(phi):
      r"""Quantum number preserving four-qubit one-parameter gate.
  
      Args:
          phi (float): rotation angle
      Returns:
          array: the four-qubit matrix describing the operation
      """
  
      c = math.cos(phi / 2)
      s = math.sin(phi / 2)
  
      return math.array(
          [
              [1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
              [0, c, 0, 0, -s, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
              [0, 0, c, 0, 0, 0, 0, 0, -s, 0, 0, 0, 0, 0, 0, 0],
              [
                  0,
                  0,
                  0,
                  c**2,
                  0,
                  0,
                  c * s,
                  0,
                  0,
                  -c * s,
                  0,
                  0,
                  s**2,
                  0,
                  0,
                  0,
              ],
              [0, s, 0, 0, c, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
              [0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
              [
                  0,
                  0,
                  0,
                  -c * s,
                  0,
                  0,
                  c**2,
                  0,
                  0,
                  s**2,
                  0,
                  0,
                  c * s,
                  0,
                  0,
                  0,
              ],
              [0, 0, 0, 0, 0, 0, 0, c, 0, 0, 0, 0, 0, s, 0, 0],
              [0, 0, s, 0, 0, 0, 0, 0, c, 0, 0, 0, 0, 0, 0, 0],
              [
                  0,
                  0,
                  0,
                  c * s,
                  0,
                  0,
                  s**2,
                  0,
                  0,
                  c**2,
                  0,
                  0,
                  -c * s,
                  0,
                  0,
                  0,
              ],
              [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0],
              [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, c, 0, 0, s, 0],
              [
                  0,
                  0,
                  0,
                  s**2,
                  0,
                  0,
                  -c * s,
                  0,
                  0,
                  c * s,
                  0,
                  0,
                  c**2,
                  0,
                  0,
                  0,
              ],
              [0, 0, 0, 0, 0, 0, 0, -s, 0, 0, 0, 0, 0, c, 0, 0],
              [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -s, 0, 0, c, 0],
              [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1],
          ]
      )
  
imports: import pennylane
  from pennylane import numpy
  from pennylane import math
label_desc: Use image node_id 29 for calling a global function with example usage: OrbitalRotation(phi) and returns: math