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.