cat >tmpfile <<\EOF
AppendTo[$Echo, "stdout"]
NDSolve[{y'[x]^2 == 1/y[x]-1, y[0]==0.01}, y, {x,0,3},MaxSteps->Infinity]
Plot[Evaluate[y[x] /. %], {x, 0, 3}]
Display["dsolve.eps",%,"EPS"]
\EOF
math <tmpfile 
ghostview dsolve.eps &