diff options
Diffstat (limited to 'Master/texmf-dist/dvips/pst-ode/pst-ode.pro')
| -rw-r--r-- | Master/texmf-dist/dvips/pst-ode/pst-ode.pro | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/Master/texmf-dist/dvips/pst-ode/pst-ode.pro b/Master/texmf-dist/dvips/pst-ode/pst-ode.pro index 848ac7d728c..df404a2de86 100644 --- a/Master/texmf-dist/dvips/pst-ode/pst-ode.pro +++ b/Master/texmf-dist/dvips/pst-ode/pst-ode.pro @@ -1,7 +1,7 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % PostScript prologue for pst-ode.tex. -% Version 0.8, 2017/06/12 +% Version 0.9, 2017/06/14 % % Alexander Grahn (C) 2012--today % @@ -117,10 +117,11 @@ end %scaling vector for step size adjustment (Numerical Recipies) eabsvect k1 eabsvect addvect {1e-30 add} forall xlength array astore % => [x(t)] [x(t+ddt) by RKF4] [err] [xscale] - %maximum rel. error - /errmax 0 def + /errmax 0 def %maximum rel. error edivvect {abs dup errmax gt {/errmax exch def}{pop} ifelse} forall errmax % => [x(t)] [x(t+ddt) by RKF4] errmax + %Sometimes, the rel. error is zero (Why is that???). So we add 1e-30 + 1e-30 add end } bind def /ODEINT { % performs integration over output step [t,t+dt] @@ -137,7 +138,7 @@ end exch pop pshrink exp 0.1 max sfty mul ddt mul /ddt exch def ode@dict tcur ddt add tcur end eq { % error: step size underflow in ODEINT - (!) odeprint + (! t=) odeprint tcur 256 string cvs odeprint true }{ (-) odeprint |