code: def mathieu_odd_coef(m, q):
      r"""Fourier coefficients for even Mathieu and modified Mathieu functions.
  
      The Fourier series of the odd solutions of the Mathieu differential
      equation are of the form
  
      .. math:: \mathrm{se}_{2n+1}(z, q) = \sum_{k=0}^{\infty} B_{(2n+1)}^{(2k+1)} \sin (2k+1)z
  
      .. math:: \mathrm{se}_{2n+2}(z, q) = \sum_{k=0}^{\infty} B_{(2n+2)}^{(2k+2)} \sin (2k+2)z
  
      This function returns the coefficients :math:`B_{(2n+2)}^{(2k+2)}` for even
      input m=2n+2, and the coefficients :math:`B_{(2n+1)}^{(2k+1)}` for odd
      input m=2n+1.
  
      Parameters
      ----------
      m : int
          Order of Mathieu functions.  Must be non-negative.
      q : float (>=0)
          Parameter of Mathieu functions.  Must be non-negative.
  
      Returns
      -------
      Bk : ndarray
          Even or odd Fourier coefficients, corresponding to even or odd m.
  
      References
      ----------
      .. [1] Zhang, Shanjie and Jin, Jianming. "Computation of Special
             Functions", John Wiley and Sons, 1996.
             https://people.sc.fsu.edu/~jburkardt/f_src/special_functions/special_functions.html
  
      """
      if not (isscalar(m) and isscalar(q)):
          raise ValueError("m and q must be scalars.")
      if q < 0:
          raise ValueError("q >=0")
      if (m != floor(m)) or (m <= 0):
          raise ValueError("m must be an integer > 0")
  
      if q <= 1:
          qm = 7.5 + 56.1 * sqrt(q) - 134.7 * q + 90.7 * sqrt(q) * q
      else:
          qm = 17.0 + 3.1 * sqrt(q) - 0.126 * q + 0.0037 * sqrt(q) * q
      km = int(qm + 0.5 * m)
      if km > 251:
          print("Warning, too many predicted coefficients.")
      kd = 4
      m = int(floor(m))
      if m % 2:
          kd = 3
  
      b = mathieu_b(m, q)
      fc = specfun.fcoef(kd, m, q, b)
      return fc[:km]
  
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 28 for calling a global function with example usage: mathieu_odd_coef(m, q) and returns: fc