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