code: def pbdn_seq(n, z):
      """Parabolic cylinder functions Dn(z) and derivatives.
  
      Parameters
      ----------
      n : int
          Order of the parabolic cylinder function
      z : complex
          Value at which to evaluate the function and derivatives
  
      Returns
      -------
      dv : ndarray
          Values of D_i(z), for i=0, ..., i=n.
      dp : ndarray
          Derivatives D_i'(z), for i=0, ..., i=n.
  
      References
      ----------
      .. [1] Zhang, Shanjie and Jin, Jianming. "Computation of Special
             Functions", John Wiley and Sons, 1996, chapter 13.
             https://people.sc.fsu.edu/~jburkardt/f_src/special_functions/special_functions.html
  
      """
      if not (isscalar(n) and isscalar(z)):
          raise ValueError("arguments must be scalars.")
      if floor(n) != n:
          raise ValueError("n must be an integer.")
      if abs(n) <= 1:
          n1 = 1
      else:
          n1 = n
      cpb, cpd = specfun.cpbdn(n1, z)
      return cpb[: n1 + 1], cpd[: n1 + 1]
  
imports: import operator
  import numpy
  import math
  from numpy import pi, asarray, floor, isscalar, iscomplex, real, imag, sqrt, where, mgrid, sin, place, issubdtype, extract, inexact, nan, zeros, sinc
  from .None import _ufuncs
  from ._ufuncs import mathieu_a, mathieu_b, iv, jv, gamma, psi, hankel1, hankel2, yv, kv, ndtri, poch, binom, hyp0f1
  from .None import specfun
  from .None import orthogonal
  from ._comb import _comb_int
label_desc: Use image node_id 41 for calling a global function with example usage: pbdn_seq(n, z) and returns: cpb,  cpd