[pstricks] postscript language to launch a rocket ?

[Herbert Voss]

Am 12.07.2010 09:01, schrieb mathias legrand:
> Hi all,
> 
> it looks like when above several digits, calculations within pstricks or
> postscript (?)
> become difficult. Is there a way to correct this ?

> \psplot[linewidth=1.25pt,algebraic=true,linecolor=red,plotpoints=500,linearc=.001]{0}{84.3}{(1/6402373705728000)*\Exp^(-(1/2)*x)*(x^18+6402373705728000-324*x^17+46818*x^16-3995136*x^15+224726400*x^14-8809274880*x^13+248127909120*x^12-5104345559040*x^11+77203226580480*x^10-857813628672000*x^9+6948290392243200*x^8-40426416827596800*x^7+165074535379353600*x^6-457129482588979200*x^5+816302647480320000*x^4-870722823979008000*x^3+489781588488192000*x^2-115242726703104000*x)}

use the exponential format of numbers:

\documentclass{book}
\usepackage{pstricks-add}

\begin{document}
\psset{xunit=0.08,yunit=1.3}
\begin{pspicture}(0,-1.3)(84.3,1.3)
\psframe[fillstyle=solid,fillcolor=gray!10,linestyle=none](0,-1.3)(84.3,1.3)
\psline[linecolor=gray!50](0,0)(84.3,0)
\psline[linecolor=gray!50](0,-1.29)(0,1.29)
\psplot[linewidth=1.25pt,algebraic,linecolor=red,plotpoints=1000,
  linearc=.001]{0}{84.3}{Euler^(-(1/2)*x)/6.402373705728e15*(
  x^18
  +6.402373705728000e15
  -324*x^17
  +4.6818e4*x^16
  -3.995136e6*x^15
  +2.24726400e8*x^14
  -8.809274880e9*x^13
  +2.48127909120e11*x^12
  -5.104345559040e12*x^11
  +7.7203226580480e13*x^10
  -8.57813628672000e14*x^9
  +6.948290392243200e15*x^8
  -4.0426416827596800e16*x^7
  +1.65074535379353600e17*x^6
  -4.57129482588979200e17*x^5
  +8.16302647480320000e17*x^4
  -8.70722823979008000e17*x^3
  +4.89781588488192000e17*x^2
  -1.15242726703104000e17*x)}
\end{pspicture}
\end{document}

Herbert