SciPy

scipy.special.mathieu_b

scipy.special.mathieu_b(x1, x2[, out]) = <ufunc 'mathieu_b'>

lmbda=mathieu_b(m,q) returns the characteristic value for the odd solution, se_m(z,q), of Mathieu’s equation