Since MaxRecursion
is set to 0
, this should work:
{int, {evals}} =
Reap[NIntegrate[
0.00159155/(x^2 + 0.0000001 x + 0.000025^2), {x, -10, 10, 0.2},
Method -> {"TrapezoidalRule", "Points" -> 51,
"RombergQuadrature" -> False},
EvaluationMonitor :> Sow[x, "EVAL"],
MaxRecursion -> 0] // Quiet,
{"EVAL"}];
evals
(*
{{-10., -9.8, -9.6, -9.4, -9.2, -9., -8.8, -8.6, -8.4, -8.2, -8.,
-7.8, -7.6, -7.4, -7.2, -7., -6.8, -6.6, -6.4, -6.2, -6., -5.8, -5.6,
-5.4, -5.2, -5., -4.8, -4.6, -4.4, -4.2, -4., -3.8, -3.6, -3.4, -3.2,
-3., -2.8, -2.6, -2.4, -2.2, -2., -1.8, -1.6, -1.4, -1.2, -1., -0.8,
-0.6, -0.4, -0.2, 0., 0.2, 0.4, 0.6, 0.8, 1., 1.2, 1.4, 1.6, 1.8, 2.,
2.2, 2.4, 2.6, 2.8, 3., 3.2, 3.4, 3.6, 3.8, 4., 4.2, 4.4, 4.6, 4.8,
5., 5.2, 5.4, 5.6, 5.8, 6., 6.2, 6.4, 6.6, 6.8, 7., 7.2, 7.4, 7.6,
7.8, 8., 8.2, 8.4, 8.6, 8.8, 9., 9.2, 9.4, 9.6, 9.8, 10., 10.,
9.902, 9.804, 9.706, 9.608, 9.51, 9.412, 9.314, 9.216, 9.118, 9.02,
8.922, 8.824, 8.726, 8.628, 8.53, 8.432, 8.334, 8.236, 8.138, 8.04,
7.942, 7.844, 7.746, 7.648, 7.55, 7.452, 7.354, 7.256, 7.158, 7.06,
6.962, 6.864, 6.766, 6.668, 6.57, 6.472, 6.374, 6.276, 6.178, 6.08,
5.982, 5.884, 5.786, 5.688, 5.59, 5.492, 5.394, 5.296, 5.198, 5.1,
5.002, 4.904, 4.806, 4.708, 4.61, 4.512, 4.414, 4.316, 4.218, 4.12,
4.022, 3.924, 3.826, 3.728, 3.63, 3.532, 3.434, 3.336, 3.238, 3.14,
3.042, 2.944, 2.846, 2.748, 2.65, 2.552, 2.454, 2.356, 2.258, 2.16,
2.062, 1.964, 1.866, 1.768, 1.67, 1.572, 1.474, 1.376, 1.278, 1.18,
1.082, 0.984, 0.886, 0.788, 0.69, 0.592, 0.494, 0.396, 0.298, 0.2}}
*)
If you are surprised that the values go from -10.
to 10.
and then from 10.
to 0.2
, it's because your code {x, -10, 10, 0.2}
tells NIntegrate
to do that.