#t = 0.4  > dsolve.log

cat >tmpfile <<\EOF
AppendTo[$Echo, "stdout"]
Needs["Graphics`Legend`"]
Plot[{4 Exp[-x^2],3(1-x^2)},{x,-1,1},
     PlotRange -> {Automatic,{-4,4}},
     PlotStyle -> {AbsoluteThickness[1],
	              {AbsoluteThickness[0.5], AbsoluteDashing[{4,4}]}},
     PlotLegend -> {"Theory","Expt"}]

Display["plotExample1.eps",%,"EPS"]

dataone = Table[N[{2 Pi i/20, Sin[2 Pi i/20]}], {i, 0, 40}];
datatwo = Table[N[{2 Pi i/20, Cos[2 Pi i/20]}], {i, 0, 40}];
fnone = Interpolation[dataone];
fntwo = Interpolation[datatwo];
Plot[{fnone[x], fntwo[x]}, {x, 0, 4 Pi},
	PlotStyle -> {
    {AbsoluteThickness[1], RGBColor[0, 0, 1]},
    {AbsoluteThickness[1], AbsoluteDashing[{4,4}]}},
	PlotLegend -> {" Sindata"," Cosdata"},
	LegendPosition -> {1.2, 0},
	LegendShadow -> {0, 0},
	LegendSize -> {1, 0.5}]

Display["plotExample2.eps",%,"EPS"]

KSubsets[l_List,0] := { {} }
KSubsets[l_List,1] := Partition[l,1]
KSubsets[l_List,k_Integer?Positive] := {l} /; (k == Length[l])
KSubsets[l_List,k_Integer?Positive] := {}  /; (k > Length[l])
KSubsets[l_List,k_Integer?Positive] :=
	Join[
		Map[(Prepend[#,First[l]])&, KSubsets[Rest[l],k-1]],
		KSubsets[Rest[l],k]
		]
KSubsets[{1, 2, 4, 3}, 1]
Vertices = Table[{Cos[2 Pi i/20], Sin[2 Pi i/20]}, {i, 1, 20}];
Show[Graphics[{
{Thickness[0.002], RGBColor[0, 0, 1], Map[Line, KSubsets[Vertices, 2]]}, 
{PointSize[0.03], RGBColor[1, 0, 0], Map[Point, Vertices]}}, 
AspectRatio -> 1]];
Display["plotExample3.eps",%,"EPS"]

randata = NestList[(# + Random[] - 0.5)&, 0, 1999];˜`
rstream = OpenWrite["rseries.dat"];
Do[
Write[rstream, randata[[i]]],{i,1,2000}];
Close[rstream];
ListPlot[randata]
Display["plotExample4.eps",%,"EPS"]


osccidata = Table[N[Sin[x] + Cos[Pi x] + 0.1*Random[]],
{x, 0, 10, 0.01}];
Needs["Utilities`BinaryFiles`"]
ostream = OpenWriteBinary["oseries.dat"];
Do[
WriteBinary[ostream, osccidata[[j]]],
{j, 1, 1001}];
Close[ostream]
ListPlot[osccidata];
Display["plotExample5.eps",%,"EPS"]
LocationReport[osccidata]
DispersionReport[osccidata]
{Min[osccidata], Max[osccidata]}
bindata = BinCounts[osccidata, {-4, 4, .1}]
BarChart[bindata, PlotRange -> All, BarValues -> True,
BarLabels -> Map[ToString, bins], 
DefaultFont -> {"Courier", 6}];
Display["plotExample6.eps",%,"EPS"]

Exit
\EOF

echo "starting mathematica in batch mode"
math <tmpfile 

ghostview plotExample3.eps