6 files changed, 45 insertions, 32 deletions
diff --git a/Master/texmf-dist/doc/generic/pst-func/Changes b/Master/texmf-dist/doc/generic/pst-func/Changes
index 6631b7f9601..e8af464d75a 100644
--- a/Master/texmf-dist/doc/generic/pst-func/Changes
+++ b/Master/texmf-dist/doc/generic/pst-func/Changes
@@ -1,4 +1,6 @@
 ..... pst-func.tex
+0.91  2018-11-15  - new algorithm for sin/cos integral to prevent problems
+                    with large values 
 0.90  2018-11-14  - fix bug with latest Ghostscript 
 0.89  2018-04-08  - added labelangle and labeldistance for \psZero, 
 0.88  2017-08-24  - use always URW fonts instead of PostScript, 
@@ -109,6 +111,7 @@
 
 
 ..... pst-func.pro
+0.16  2018-11-15  - new algorithm for Si/si/Ci/ci
 0.15  2013-06-10  - fix bug with wrong arrow in psBezier#
 0.14  2013-06-03  - fix bug with arrow in psBezier#
 0.13  2010-06-21  - run FindZeros in a local dictionary
diff --git a/Master/texmf-dist/doc/generic/pst-func/README.md b/Master/texmf-dist/doc/generic/pst-func/README.md
index e561341de8e..3297389d3e7 100644
--- a/Master/texmf-dist/doc/generic/pst-func/README.md
+++ b/Master/texmf-dist/doc/generic/pst-func/README.md
@@ -7,7 +7,7 @@
 - psFourier
 - psBessel
 - psSi and pssi (integral sin)
-- psCi (integral cosin)
+- psCi and \psci (integral cosin)
 - psIntegral
 - psCumIntegral
 - psConv
@@ -21,7 +21,7 @@
 - psTDist for Student's t-distribution
 - psBetaDist for Beta distribution
 - psWeibull for Weibull distribution 
-- psLame (Lam\`e\ Curve - a superellipse)
+- psLame (Lamé Curve - a superellipse)
 - psThomae (popcorn function)
 - psWeierstrass (original and modified)
 - psplotImp (plotting implicit defined functions)
@@ -44,7 +44,4 @@ PSTricks is PostScript Tricks, the documentation cannot be run
 with pdftex, use the sequence latex->dvips->ps2pdf or
 pdflatex with package auto-pst-pdf or xelatex.
 
-hvoss@tug.org
-
 %% $Id: README 897 2014-03-21 08:06:41Z herbert $
-
diff --git a/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.pdf b/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.pdf
index 94a83645b63..ea5c1d52bf5 100644
--- a/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.pdf
+++ b/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.pdf
diff --git a/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.tex b/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.tex
index bd34476ee19..3003c3c0856 100644
--- a/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.tex
+++ b/Master/texmf-dist/doc/generic/pst-func/pst-func-doc.tex
@@ -1184,11 +1184,11 @@ The integral sin  and cosin are defined as
 
 \begin{LTXexample}[pos=t]
 \psset{xunit=0.5}
-\begin{pspicture}(-14,-4.5)(14,2)
-  \psaxes[dx=1cm,Dx=2]{->}(0,0)(-13.1,-4)(13,2)
+\begin{pspicture}(-15,-4.5)(15,2)
+  \psaxes[dx=1cm,Dx=2]{->}(0,0)(-14.1,-4)(14,2)
   \psplot[plotpoints=1000]{-12.5}{12.5}{ x RadtoDeg sin x div }
-  \psSi[plotpoints=1500,linecolor=red,linewidth=1pt]{-12.5}{12.5}
-  \pssi[plotpoints=1500,linecolor=blue,linewidth=1pt]{-12.5}{12.5}
+  \psSi[plotpoints=1500,linecolor=red,linewidth=1pt]{-13.5}{13.5}
+  \pssi[plotpoints=1500,linecolor=blue,linewidth=1pt]{-13.5}{13.5}
   \rput(-5,1.5){\color{red}$Si(x)=\int\limits_{0}^x \frac{\sin(t)}{t}\dt$}  
   \rput(8,-1.5){\color{blue}$si(x)=-\int\limits_{x}^{\infty} \frac{\sin(t)}{t}\dt=Si(x)-\frac{\pi}{2}$}  
   \rput(8,.5){$f(x)= \frac{\sin(t)}{t}$}
diff --git a/Master/texmf-dist/dvips/pst-func/pst-func.pro b/Master/texmf-dist/dvips/pst-func/pst-func.pro
index 8be99845cd6..d5c66b3829e 100644
--- a/Master/texmf-dist/dvips/pst-func/pst-func.pro
+++ b/Master/texmf-dist/dvips/pst-func/pst-func.pro
@@ -16,7 +16,7 @@
 %%   `pst-func' is a PSTricks package to plot special math functions
 %%
 %%
-%% version 0.16 / 2014-05-07  Herbert Voss
+%% version 0.16 / 2018-11-15  Herbert Voss
 %
 /tx@FuncDict 100 dict def
 tx@FuncDict begin
@@ -127,19 +127,31 @@ tx@FuncDict begin
 %%
 /Si { % integral sin from 0 to x (arg on stack)
   10 dict begin				% hold all local
-  /arg exch def  % x
-  /Sum arg def   %
-  /sign -1 def   
-  /I 3 def
-  { %  a sequence of x - x^3/(3*3!) + x^5/(5*5!) -...+...
-    arg I Power dup abs 1e30 gt { pop exit } if
-    I factorial div I div sign mul 
-    dup abs eps8 lt { pop exit } if 
-    Sum add /Sum exch def
-    /sign sign neg def
-    /I I 2 add def
-  } loop
-  Sum
+  dup 0 eq 
+  { pop 0 } 
+  { 
+    /arg exch def  % x
+    /arg2 arg dup mul def
+    /Sum arg def   %
+    /sign -1 def   
+    /I 3 def
+    /Frac arg2 arg mul 6 div def
+    { %  a sequence of x - x^3/(3*3!) + x^5/(5*5!) -...+...
+      Frac I div sign mul
+      dup abs eps5 lt { pop exit } if
+      Sum add /Sum exch def
+      /sign sign neg def
+      /I I 2 add def
+      Frac arg2 mul I 1 sub I mul div /Frac ED
+%    arg I Power dup abs 1e30 gt { pop exit } if
+%    I factorial div I div sign mul 
+%    dup abs eps8 lt { pop exit } if 
+%    Sum add /Sum exch def
+%    /sign sign neg def
+%    /I I 2 add def
+   } loop
+   Sum
+  } ifelse
   end
 } def
 %
@@ -151,16 +163,17 @@ tx@FuncDict begin
   10 dict begin				% hold all local
   abs /arg exch def
   arg 0 eq { 0 } { 
+    /arg2 arg dup mul def 
     /Sum CEuler arg ln add def
     /sign -1 def
     /I 2 def
-    { 
-    arg I Power dup abs 1e30 gt { pop exit } if
-    I factorial div I div sign mul 
-    dup abs eps8 lt { pop exit } if 
-    Sum add /Sum exch def
-    /sign sign neg def
-    /I I 2 add def
+    /Frac arg2 2 div def % first fraction
+    { Frac I div sign mul
+      dup abs eps5 lt { pop exit } if
+      Sum add /Sum exch def
+     /sign sign neg def
+     /I I 2 add def
+     Frac arg2 mul I 1 sub I mul div /Frac ED
     } loop
     Sum
   } ifelse
diff --git a/Master/texmf-dist/tex/generic/pst-func/pst-func.tex b/Master/texmf-dist/tex/generic/pst-func/pst-func.tex
index f65c20d6c7b..a1490b48993 100644
--- a/Master/texmf-dist/tex/generic/pst-func/pst-func.tex
+++ b/Master/texmf-dist/tex/generic/pst-func/pst-func.tex
@@ -31,8 +31,8 @@
 % interface to the `xkeyval' package
 \pst@addfams{pst-func}
 %
-\def\fileversion{0.90}
-\def\filedate{2018/11/14}
+\def\fileversion{0.91}
+\def\filedate{2018/11/15}
 \message{`PST-func' v\fileversion, \filedate\space (hv)}
 %
 \pstheader{pst-func.pro}