# Numbers plotted by ListPlot don't match the data

My mathematica is doing weird things. Consider the function

f[n_]:=Re[(1/2)*(I*Pi+(-1)^n*n!*Gamma[-n,-n])]

Let's make a list out of this

points:=Table[{n,f[n]},{n,60,63}]

when evaluated I get that the list points has these four entries

N@points

(* {{60., 0.0539289}, {61., 0.0534851}, {62., 0.00402142}, {63.,0.00395776}} *)

That is, the first two are above 0.05 and the last two are significantly smaller. Nonetheless, when I perform a ListPlot of points I get a plot where the four dots are way above 0.05. What is going on? Is this a bug?

ListPlot[points, PlotRange -> {0, 0.06}, PlotStyle -> Red]

• Your Table is missing an iterator (n) and Gamma[-60,-60] doesn't evaluate for me. Did you transcribe your code correctly? Commented Apr 20, 2017 at 13:52
• edited and Gamma[-60,-60] evaluates without any problems to me Commented Apr 20, 2017 at 14:04
• In order to see the decimal output you show, you had to use N, otherwise you would see something with very large integers multiplied by Gama[-60,-60], so I added that line to the post. So if you use ListPlot[N@points] then you see what you expect Commented Apr 20, 2017 at 14:04
• @JasonB. ok. But why does not my code work? It should work shouldnt it? Commented Apr 20, 2017 at 14:06
• @AnarchistBirdsWorshipFungus - you are multiplying very huge numbers (10^83) by very small numbers (10^-82), so you have to give some consideration to precision. I think that ListPlot is being smarter than N in this case, check the output from this: Block[{\$MaxExtraPrecision = 100}, N[points, {∞, 5}] ] Commented Apr 20, 2017 at 14:09