[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Axiom-developer] 20080209.01.tpd.patch (add Exponential Integral to boo
From: |
daly |
Subject: |
[Axiom-developer] 20080209.01.tpd.patch (add Exponential Integral to book) |
Date: |
Sun, 10 Feb 2008 01:12:28 -0600 |
Add Exponential Integral family of functions to the Jenks book.
====================================================================
diff --git a/changelog b/changelog
index c9a8e1a..c15d753 100644
--- a/changelog
+++ b/changelog
@@ -1,3 +1,4 @@
+20080109 tpd src/doc/book add Ei,En,Ei1,Ei2,Ei3,Ei4,Ei5,Ei6
20080107 tpd Makefile fix GCLOPTS-CUSTRELOC for macosxppc
20080107 tpd Makefile add Makefile.macosxppc stanza and GCLOPTS-CUSTRELOC
20080107 tpd src/algebra/Makefile make cp of upper/lower files conditional
diff --git a/src/doc/book.pamphlet b/src/doc/book.pamphlet
index 377d5c6..4fdb58c 100644
--- a/src/doc/book.pamphlet
+++ b/src/doc/book.pamphlet
@@ -17252,7 +17252,7 @@ $z=3$:
%\epsffile[0 0 295 295]{ps/newmap.ps}
% I think this is good to say here: it shows a lot of depth. RSS
-{\sloppy
+%{\sloppy
The {\tt CoordinateSystems} package exports the following
\index{coordinate system}
operations:
@@ -18231,12 +18231,12 @@ $Gamma(z)$ is the Euler gamma function,
{\bf Beta}: $F -> F$\hfill\newline
$Beta(u, v)$ is the Euler Beta function,
\index{function!Euler Beta}
- $B(u,v)$, defined by
+ $Beta(u,v)$, defined by
\index{Euler!Beta function}
- $$B(u,v) = \int_{0}^{1} t^{u-1} (1-t)^{v-1} dt.$$
+ $$Beta(u,v) = \int_{0}^{1} t^{u-1} (1-t)^{v-1} dt.$$
This is related to $\Gamma(z)$ by
- $$B(u,v) = \frac{\Gamma(u) \Gamma(v)}{\Gamma(u + v)}.$$
+ $$Beta(u,v) = \frac{\Gamma(u) \Gamma(v)}{\Gamma(u + v)}.$$
\noindent
{\bf logGamma}: $F -> F$\hfill\newline
@@ -18261,6 +18261,63 @@ is the function $\psi(z)$,
$\psi(z)$, written $\psi^{(n)}(z)$.
\noindent
+{\bf E1}: $(DoubleFloat) -> OnePointCompletion DoubleFloat$\hfill\newline
+ E1(x) is the Exponential Integral function
+ The current implementation is a piecewise approximation
+ involving one poly from $-4..4$ and a second poly for $x > 4$
+\index{function!E1}
+
+\noindent
+{\bf En}: $(PI, DFLOAT) -> OnePointCompletion DoubleFloat$\hfill\newline
+ En(PI,R) is the nth Exponential Integral
+\index{function!En}
+
+\noindent
+{\bf Ei}: $(OnePointCompletion DFLOAT) -> OnePointCompletion DFLOAT$
+\hfill\newline
+ Ei is the Exponential Integral function
+ This is computed using a 6 part piecewise approximation.
+ DoubleFloat can only preserve about 16 digits but the
+ Chebyshev approximation used can give 30 digits.
+\index{function!Ei}
+
+\noindent
+{\bf Ei1}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei1 is the first approximation of Ei where the result is
+ $x*e^-x*Ei(x)$ from -infinity to -10 (preserves digits)
+\index{function!Ei1}
+
+\noindent
+{\bf Ei2}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei2 is the first approximation of Ei where the result is
+ $x*e^-x*Ei(x)$ from -10 to -4 (preserves digits)
+\index{function!Ei2}
+
+\noindent
+{\bf Ei3}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei3 is the first approximation of Ei where the result is
+ $(Ei(x)-log |x| - gamma)/x$ from -4 to 4 (preserves digits)
+\index{function!Ei3}
+
+\noindent
+{\bf Ei4}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei4 is the first approximation of Ei where the result is
+ $x*e^-x*Ei(x)$ from 4 to 12 (preserves digits)
+\index{function!Ei4}
+
+\noindent
+{\bf Ei5}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei5 is the first approximation of Ei where the result is
+ $x*e^-x*Ei(x)$ from 12 to 32 (preserves digits)
+\index{function!Ei5}
+
+\noindent
+{\bf Ei6}: $(DoubleFloat) -> DoubleFloat$\hfill\newline
+ Ei6 is the first approximation of Ei where the result is
+ $x*e^-x*Ei(x)$ from 32 to infinity (preserves digits)
+\index{function!Ei6}
+
+\noindent
{\bf besselJ}: $(F,F) -> F$\hfill\newline
$besselJ(v,z)$ is the Bessel function of the first kind,
\index{function!Bessel}
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Axiom-developer] 20080209.01.tpd.patch (add Exponential Integral to book),
daly <=