One image is worth one thousand words, so:

Question: At what options should I be looking, so that what I see in my notebook is what I get in my printed pdf ?
Here is the code I'm using:
deq = x''[t] == x[t] - x[t]^3 + x[t]^2 + a Cos[t]; genPoincareMapPoints[deq_] := Module[{}, Reap[ NDSolve[{deq, WhenEvent[ Mod[t, 2 \[Pi]] == 0, Sow[{x[t], x'[t]}]] }, {}, {t, 0, 5 10^4}, AccuracyGoal -> 10, PrecisionGoal -> 10, MaxSteps -> ∞]]][[-1, 1]]; Style[ Grid@{ Parallelize[ ListPlot[ genPoincareMapPoints[ {deq /. a -> #, x[0] == 0, x'[0] == 0.0}], AxesLabel -> {"x[t]", "x'[t]"}, PlotRange -> All, PlotLabel -> "Poincare Map for a=" <> ToString@#, PlotStyle -> AbsolutePointSize[0.5]] & /@ {0.03, 0.3, 3}] }, ImageSizeMultipliers -> 1] EDIT: Mac OSX 10.9, Mathematica 9.0.1. I am printing with 'File -> Print... -> PDF -> Save as PDF...'. I'm viewing the PDF with 'Preview' app of Mac OSX.

