The BesselJ function is an oscillating function that has an infinite number of zeros. There is no known closed form for the solutions for the equation, so solve() is not able to return a symbolic solution.
However I would also point out that you have used beselj(0,a) . That corresponds to
but the equation in the image is
which would correspond to besselj(1,a) rather than besselj(0,a) What can you do? Well, you cannot find a useful symbolic solution for
. Numeric solutions exist, and you can approximate how far apart they are, but you cannot give a formula for them except in terms of infinite summations or infinite integrals. vpasolve() can find numeric roots without complaining about it. Because there are an infinite number of roots, you should pass a starting point into vpasolve(), such as
vpasolve(besselj(1,a), 3)
vpasolve(besselj(1,a), 6)