I am trying to plot this function
f[X_,Y_]=-(1/Sqrt[2])(√((1/2 - 0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2/((1/2 -
0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2 +
10000.3 BesselJ[1, 0.01 X]^2 Sin[Y]^2) + (
10000.3 BesselJ[1, 0.01 X]^2 Sin[
Y]^2)/((1/2 - 0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2 +
10000.3 BesselJ[1, 0.01 X]^2 Sin[
Y]^2) + (√((1/2 - 0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2 ((1/2 -
0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2 +
10000.3 BesselJ[1, 0.01 X]^2 Sin[Y]^2)))/((1/2 -
0.500013 BesselJ[0, 0.01 X] +
100.001 BesselJ[1, 0.01 X] Cos[Y])^2 +
10000.3 BesselJ[1, 0.01 X]^2 Sin[Y]^2)))
But I get this plot.
Plot[f[X, π/2], {X, 1, 10^4}]
I was expecting more regular oscillations, whereas I see in the lower part certain irregularities, and this make me think of some plotting error. Is there anything I can do to make the bessel function be plotted right?
Thanks!
Plot
probably you could try increasing the number of points the function is evaluated on. beware this will be slow:Plot[f[X, Pi/2], {X, 1, 10^4}, PlotRange -> {-1.1, -0.7}, PlotPoints -> 100, MaxRecursion -> 15]
$\endgroup$ – rhermans Oct 21 '15 at 14:53