I'm looking for a function (preferrably in Python) which takes a complex number $a$ as input, and outputs the zeros of the function $f:x \rightarrow J_x(a)$ where $J_x(a)$ is the Bessel function of the first kind. This problem is relevant to solving a Schroedinger equation in a 1D condensed matter system (N.B. I am not asking to find the zeros of the bessel function $f:x \rightarrow J_n(x)$).
Does anybody know of a programming language such as Python, Wolfram language, etc. which has this built-in; and if so what is that function? It would save me the effort of implementing some sort of Newton's method-esque function myself.
Please let me know if this belongs on stackoverflow.