code: def pbdv_seq(v, x):
      """Parabolic cylinder functions Dv(x) and derivatives.
  
      Parameters
      ----------
      v : float
          Order of the parabolic cylinder function
      x : float
          Value at which to evaluate the function and derivatives
  
      Returns
      -------
      dv : ndarray
          Values of D_vi(x), for vi=v-int(v), vi=1+v-int(v), ..., vi=v.
      dp : ndarray
          Derivatives D_vi'(x), for vi=v-int(v), vi=1+v-int(v), ..., vi=v.
  
      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(v) and isscalar(x)):
          raise ValueError("arguments must be scalars.")
      n = int(v)
      v0 = v - n
      if n < 1:
          n1 = 1
      else:
          n1 = n
      v1 = n1 + v0
      dv, dp, pdf, pdd = specfun.pbdv(v1, x)
      return dv[: n1 + 1], dp[: 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 39 for calling a global function with example usage: pbdv_seq(v, x) and returns: dv,  dp