python,numpy,scipy,vectorization,bessel-functions

You can write a cython function to speedup the calculation, the first thing you must to do is to get the address of the fortran function SPHJ, here is how to do this in Python: from scipy import special as sp sphj = sp.specfun.sphj import ctypes addr = ctypes.pythonapi.PyCObject_AsVoidPtr(ctypes.py_object(sphj._cpointer)) Then...

matlab,equation-solving,bessel-functions

You can use (http://nl.mathworks.com/help/matlab/ref/besselj.html, http://nl.mathworks.com/help/optim/ug/fsolve.html) x = fsolve(@(x)(besselj(0,a*x)-b),x0); or using a simpler solver as pointed in the comment below (http://nl.mathworks.com/help/optim/ug/fzero.html) x = fzero(@(x)(besselj(0,a*x)-b),x0); ...