code: def ControlledPhaseShift(phi):
      r"""Controlled phase shift.
  
      Args:
          phi (float): rotation angle
  
      Returns:
          array: the two-wire controlled-phase matrix
      """
      return math.diag([1, 1, 1, math.exp(1j * phi)])
  
imports: import pennylane
  from pennylane import numpy
  from pennylane import math
label_desc: Use image node_id 19 for calling a global function with example usage: ControlledPhaseShift(phi) and returns: math