summaryrefslogtreecommitdiff
path: root/Master/texmf-dist/dvips/pst-ode/pst-ode.pro
diff options
context:
space:
mode:
Diffstat (limited to 'Master/texmf-dist/dvips/pst-ode/pst-ode.pro')
-rw-r--r--Master/texmf-dist/dvips/pst-ode/pst-ode.pro9
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