Fabrice23 said:
What I don't get, why is the max value of BesselJ(n, eta)^2 rougly around eta?
It is an intrinsic property of the Jn(x) function.
In case of constant and large x, the behaviour of Jn(x) as a function of n is known for a long time. Roughly, three domains are distinguished :
First domain: relatively small n (compared to large x).
Jn(x) is oscillating with amplitude roughly = sqrt(2/(pi*x))
Second domain: The transition.
As n increases, the amplitude of oscillation increases. For n not far from x, the oscillations cease. Then, with n increasing, Jn(x) begins to decrease.
In the transition domain, the asymptotic formulas are very complicated. A good compilation can be found in "Handbook of Mathematical Functions", M.Abramowitz, I.A.Stegun, Dover Publications, N.-Y., 1972, “Bessel functions of integer order", Ch.9, pp.355-389, especially “Asymptotic expansions in the transition region” pp.367-368.
It should be too complicated to deal with it on a forum. For even more explanation about those formulas, it should be necessary to consult some specialized documents listed p.368
Third domain : relatively large n (compared to x)
As n increases, Jn(x) decreases very rapidly. An asymptotic formula is :
Jn(x) roughly = (1/sqrt(2*pi*n))*(e*x/(2*n))^n
This formula is not valid on the transition domain. More accurate formulas can be found in the book referenced above.
The behaviour of the function Jn(n) is different (because representing for each n a particular point in the transition domain). This is a continuously decreasing function. An asymptotic formula is :
Jn(n) roughly=(1/Gamma(2/3))* (2/(9*n))^(1/3)
All above is for positive n.
For negative n, see the book referenced above (especially about oscillations in case of large negative n)