(* this file defines a bunch of plots, puts them all on
one page, then saves the result to a postscript file called
ptot.ps *)

(* an example using the NDSolve function *)
e = NDSolve[{y'[x] == y[x], y[3] == 1}, y, {x, 0, 2}]
p1 = Plot[Evaluate[y[x] /. e], {x, 0, 2}]

(* simple sine curve plot *)
p2 = Plot[Sin[x], {x, 0, 3 Pi}]

(* plot with frame, gridlines *)
p3 = Plot[BesselJ[3,x], {x, 0, 10}, Frame -> True,          
         GridLines -> Automatic]

(* contour plot *)
p4 = ContourPlot[Sin[x] Sin[y], {x, -2, 2}, {y, -2, 2}]

(* 3-dimensional plot *)
p5 = Plot3D[Sin[x y], {x, 0, 3}, {y, 0, 3}]

(* add axes labels to the previous plot *)
p6 = Show[p5, AxesLabel -> {"X-AXIS", "Y-AXIS", "Z_AXIS"},
         FaceGrids -> All]

(* a parametric plot *)
p7 = ParametricPlot[{Sin[t], Sin[2t]}, {t, 0, 2 Pi}]

(* a couple 3-dimensional parametric plots *)
p8 = ParametricPlot3D[{Sin[t], Cos[t], t/3}, {t, 0, 15}]

p9 = ParametricPlot3D[{Cos[t] Cos[u], Sin[t] Cos[u], Sin[u]},
         {t, 0, 2 Pi}, {u, -Pi/2, Pi/2}]

(* add figure labels to all the plots above *)
pn1 = Show[p1, PlotLabel -> "FIGURE 1"];
pn2 = Show[p2, PlotLabel -> "FIGURE 2"];
pn3 = Show[p3, PlotLabel -> "FIGURE 3"];
pn4 = Show[p4, PlotLabel -> "FIGURE 4"];
pn5 = Show[p5, PlotLabel -> "FIGURE 5"];
pn6 = Show[p6, PlotLabel -> "FIGURE 6"];
pn7 = Show[p7, PlotLabel -> "FIGURE 7"];
pn8 = Show[p8, PlotLabel -> "FIGURE 8"];
pn9 = Show[p9, PlotLabel -> "FIGURE 9"];

(* put all the plots above together on one page using the GraphicsArray
command *)

ptot = Show[GraphicsArray[{{pn1, pn2},{ pn3, pn4}, {pn5, pn6}, {pn7, pn8}, {pn9}}]]

(* write the result to a postscript file called ptot.ps *)
Display["!psfix > ptot.ps", ptot]