Initial population.

1 files changed, 701 insertions, 0 deletions

diff --git a/src/input/nsfip.input.pamphlet b/src/input/nsfip.input.pamphlet
new file mode 100644
index 00000000..87c68a44
--- /dev/null
+++ b/src/input/nsfip.input.pamphlet
@@ -0,0 +1,701 @@
+\documentclass{article}
+\usepackage{axiom}
+\begin{document}
+\title{\$SPAD/src/input nsfip.input}
+\author{The Axiom Team}
+\maketitle
+\begin{abstract}
+\end{abstract}
+\eject
+\tableofcontents
+\eject
+\section{License}
+<<license>>=
+--Copyright The Numerical Algorithms Group Limited 1994.
+@
+<<*>>=
+<<license>>
+
+-- Note that the conversions of Results from DoubleFloat to Float
+-- will become unnecessary if outputGeneral is extended to apply to
+-- DoubleFloat quantities.
+
+
+outputGeneral 4
+
+-- DF here means DoubleFloat.
+-- Results converted to Float as outputGeneral not working on DF.
+
+-- nagExpInt : DF -> DF ;
+
+nagExpInt(2) :: Float
+
+--       0.0489
+
+nagExpInt(-1) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S13AAF: IFAIL =     1
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S13AAF. The error number (IFAIL value) is 1, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+-- nagSinInt : DF -> DF ;
+
+nagSinInt(0) :: Float
+
+--       0.0
+
+nagSinInt(0.2) :: Float
+
+--       0.1996
+
+nagSinInt(0.4) :: Float
+
+--       0.3965
+
+nagSinInt(0.6) :: Float
+
+--       0.5881
+
+nagSinInt(0.8) :: Float
+
+--       0.7721
+
+nagSinInt(1) :: Float
+
+--       0.9461
+
+-- nagCosInt : DF -> DF ;
+
+nagCosInt(0.2) :: Float
+
+--       - 1.042
+
+nagCosInt(0.4) :: Float
+
+--       - 0.3788
+
+nagCosInt(0.6) :: Float
+
+--       - 0.02227
+
+nagCosInt(0.8) :: Float
+
+--       0.1983
+
+nagCosInt(1) :: Float
+
+--       0.3374
+
+-- nagIncompleteGammaP : (DF, DF) -> DF ; (to machine precision)
+
+nagIncompleteGammaP(2,3) :: Float
+
+--       0.8009
+
+nagIncompleteGammaP(7,1) :: Float
+
+--       0.00008324
+
+nagIncompleteGammaP(0.5,99) :: Float
+
+--       1.0
+
+nagIncompleteGammaP(20,21) :: Float
+
+--       0.6157
+
+nagIncompleteGammaP(21,20) :: Float
+
+--       0.4409
+
+-- nagIncompleteGammaP : (DF, DF, DF) -> DF ; (to specified precision)
+
+nagIncompleteGammaP(7,1,0.1) :: Float
+
+--       0.00008313
+
+-- nagIncompleteGammaQ : (DF, DF) -> DF ; (to machine precision)
+
+nagIncompleteGammaQ(2,3) :: Float
+
+--       0.1991
+
+nagIncompleteGammaQ(7,1) :: Float
+
+--       0.9999
+
+nagIncompleteGammaQ(0.5,99) :: Float
+
+--       0.5705 E -44
+
+nagIncompleteGammaQ(20,21) :: Float
+
+--       0.3843
+
+nagIncompleteGammaQ(21,20) :: Float
+
+--       0.5591
+
+nagIncompleteGammaQ(25,14) :: Float
+
+--       0.995
+
+-- nagIncompleteGammaQ : (DF, DF, DF) -> DF ; (to specified precision)
+
+nagIncompleteGammaQ(25,14,0.1) :: Float
+
+--       0.9953
+
+-- nagErf : DF -> DF ;
+
+nagErf(-6) :: Float
+
+--       - 1.0
+
+nagErf(-4.5) :: Float
+
+--       - 1.0
+
+nagErf(-1) :: Float
+
+--       - 0.8427
+
+nagErf(1) :: Float
+
+--       0.8427
+
+nagErf(4.5) :: Float
+
+--       1.0
+
+nagErf(6) :: Float
+
+--       1.0
+
+-- nagErfC : DF -> DF ;
+
+nagErfC(-10) :: Float
+
+--       2.0
+
+nagErfC(-1) :: Float
+
+--       1.843
+
+nagErfC(0) :: Float
+
+--       1.0
+
+nagErfC(1) :: Float
+
+--       0.1573
+
+nagErfC(15) :: Float
+
+--       0.7213 E -99
+
+-- nagDAiryAi : DF -> DF ;
+
+nagDAiryAi(-10) :: Float
+
+--       0.9963
+
+nagDAiryAi(-1) :: Float
+
+--       - 0.01016
+
+nagDAiryAi(0) :: Float
+
+--       - 0.2588
+
+nagDAiryAi(1) :: Float
+
+--       - 0.1591
+
+nagDAiryAi(5) :: Float
+
+--       - 0.0002474
+
+nagDAiryAi(10) :: Float
+
+--       - 0.3521 E -9
+
+nagDAiryAi(20) :: Float
+
+--       - 0.7586 E -26
+
+-- nagDAiryAi : CDF -> CDF ;
+
+nagDAiryAi(0.3+0.4*%i) :: Complex Float
+
+--       - 0.2612 + 0.03848 %i
+
+-- nagDAiryBi : DF -> DF ;
+
+nagDAiryBi(-10) :: Float
+
+--       0.1194
+
+nagDAiryBi(-1) :: Float
+
+--       0.5924
+
+nagDAiryBi(0) :: Float
+
+--       0.4483
+
+nagDAiryBi(1) :: Float
+
+--       0.9324
+
+nagDAiryBi(5) :: Float
+
+--       1436.0
+
+nagDAiryBi(10) :: Float
+
+--       0.1429 E 10
+
+nagDAiryBi(20) :: Float
+
+--       0.9382 E 26
+
+-- nagDAiryBi : CDF -> CDF ;
+
+nagDAiryBi(0.3+0.4*%i) :: Complex Float
+
+--       0.4093 + 0.07966 %i
+
+-- nagScaledDAiryAi : CDF -> CDF ;
+
+nagScaledDAiryAi(0.3+0.4*%i) :: Complex Float
+
+--       - 0.2744 - 0.02356 %i
+
+-- nagScaledDAiryBi : CDF -> CDF ;
+
+nagScaledDAiryBi(0.3+0.4*%i) :: Complex Float
+
+--       0.3924 + 0.07638 %i
+
+-- nagHankelH1 : (DF, CDF, Int) -> List CDF ;
+
+nagHankelH1(0,0.3+0.4*%i,2) :: Matrix Complex Float
+
+--       [0.3466 - 0.5588 %i  - 0.7912 - 0.8178 %i]
+
+nagHankelH1(2.3,2,2) :: Matrix Complex Float
+
+--       [0.2721 - 0.7398 %i  0.08902 - 1.412 %i]
+
+nagHankelH1(2.12,-1,2) :: Matrix Complex Float
+
+--       [- 0.7722 - 1.693 %i  2.601 + 6.527 %i]
+
+-- nagHankelH2 : (DF, CDF, Int) -> List CDF ;
+
+nagHankelH2(6,3.1-1.6*%i,2) :: Matrix Complex Float
+
+--       [- 1.371 - 1.28 %i  - 1.491 - 5.993 %i]
+
+-- nagScaledHankelH1 : (DF, CDF, Int) -> List CDF ;
+
+nagScaledHankelH1(0,0.3+0.4*%i,2) :: Matrix Complex Float
+
+--       [0.2477 - 0.9492 %i  - 1.488 - 0.8166 %i]
+
+-- nagScaledHankelH2 : (DF, CDF, Int) -> List CDF ;
+
+nagScaledHankelH2(6,3.1-1.6*%i,2) :: Matrix Complex Float
+
+--       [7.05 + 6.052 %i  8.614 + 29.35 %i]
+
+-- nagKelvinBer : DF -> DF ;
+
+nagKelvinBer(0.1) :: Float
+
+--       1.0
+
+nagKelvinBer(1) :: Float
+
+--       0.9844
+
+nagKelvinBer(2.5) :: Float
+
+--       0.4
+
+nagKelvinBer(5) :: Float
+
+--       - 6.23
+
+nagKelvinBer(10) :: Float
+
+--       138.8
+
+nagKelvinBer(15) :: Float
+
+--       - 2967.0
+
+nagKelvinBer(60) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19AAF: IFAIL =     1
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19AAF. The error number (IFAIL value) is 1, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+nagKelvinBer(-1) :: Float
+
+--       0.9844
+
+-- nagKelvinBei : DF -> DF ;
+
+nagKelvinBei(0.1) :: Float
+
+--       0.0025
+
+nagKelvinBei(1) :: Float
+
+--       0.2496
+
+nagKelvinBei(2.5) :: Float
+
+--       1.457
+
+nagKelvinBei(5) :: Float
+
+--       0.116
+
+nagKelvinBei(10) :: Float
+
+--       56.37
+
+nagKelvinBei(15) :: Float
+
+--       - 2953.0
+
+nagKelvinBei(60) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ABF: IFAIL =     1
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ABF. The error number (IFAIL value) is 1, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+nagKelvinBei(-1) :: Float
+
+--       0.2496
+
+-- nagKelvinKer : DF -> DF ;
+
+nagKelvinKer(0) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ACF: IFAIL =     2
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ACF. The error number (IFAIL value) is 2, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+nagKelvinKer(0.1) :: Float
+
+--       2.42
+
+nagKelvinKer(1) :: Float
+
+--       0.2867
+
+nagKelvinKer(2.5) :: Float
+
+--       - 0.06969
+
+nagKelvinKer(5) :: Float
+
+--       - 0.01151
+
+nagKelvinKer(10) :: Float
+
+--       0.0001295
+
+nagKelvinKer(15) :: Float
+
+--       - 0.1514 E -7
+
+nagKelvinKer(1100) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ACF: IFAIL =     1
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ACF. The error number (IFAIL value) is 1, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+nagKelvinKer(-1) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ACF: IFAIL =     2
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ACF. The error number (IFAIL value) is 2, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+-- nagKelvinKei : DF -> DF ;
+
+nagKelvinKei(0) :: Float
+
+--       - 0.7854
+
+nagKelvinKei(0.1) :: Float
+
+--       - 0.7769
+
+nagKelvinKei(1) :: Float
+
+--       - 0.495
+
+nagKelvinKei(2.5) :: Float
+
+--       - 0.1107
+
+nagKelvinKei(5) :: Float
+
+--       0.01119
+
+nagKelvinKei(10) :: Float
+
+--       - 0.0003075
+
+nagKelvinKei(15) :: Float
+
+--       0.000007963
+
+nagKelvinKei(1100) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ADF: IFAIL =     1
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ADF. The error number (IFAIL value) is 1, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+nagKelvinKei(-1) :: Float
+
+-- ** ABNORMAL EXIT from NAG Library routine S19ADF: IFAIL =     2
+-- ** NAG soft failure - control returned
+-- 
+--   Error signalled from user code:
+--      An error was detected when calling the NAG Library routine 
+--      S19ADF. The error number (IFAIL value) is 2, please consult the 
+--      NAG manual via the Browser for diagnostic information.
+
+
+-- nagFresnelS : DF -> DF ;
+
+nagFresnelS(0) :: Float
+
+--       0.0
+
+nagFresnelS(0.5) :: Float
+
+--       0.06473
+
+nagFresnelS(1) :: Float
+
+--       0.4383
+
+nagFresnelS(2) :: Float
+
+--       0.3434
+
+nagFresnelS(4) :: Float
+
+--       0.4205
+
+nagFresnelS(5) :: Float
+
+--       0.4992
+
+nagFresnelS(6) :: Float
+
+--       0.447
+
+nagFresnelS(8) :: Float
+
+--       0.4602
+
+nagFresnelS(10) :: Float
+
+--       0.4682
+
+nagFresnelS(-1) :: Float
+
+--       - 0.4383
+
+nagFresnelS(1000) :: Float
+
+--       0.4997
+
+-- nagFresnelC : DF -> DF ;
+
+nagFresnelC(0) :: Float
+
+--       0.0
+
+nagFresnelC(0.5) :: Float
+
+--       0.4923
+
+nagFresnelC(1) :: Float
+
+--       0.7799
+
+nagFresnelC(2) :: Float
+
+--       0.4883
+
+nagFresnelC(4) :: Float
+
+--       0.4984
+
+nagFresnelC(5) :: Float
+
+--       0.5636
+
+nagFresnelC(6) :: Float
+
+--       0.4995
+
+nagFresnelC(8) :: Float
+
+--       0.4998
+
+nagFresnelC(10) :: Float
+
+--       0.4999
+
+nagFresnelC(-1) :: Float
+
+--       - 0.7799
+
+nagFresnelC(1000) :: Float
+
+--       0.5
+
+-- nagEllipticIntegralRC : (DF, DF) -> DF ;
+
+nagEllipticIntegralRC(0.5,1) :: Float
+
+--       1.111
+
+nagEllipticIntegralRC(1,1) :: Float
+
+--       1.0
+
+nagEllipticIntegralRC(1.5,1) :: Float
+
+--       0.9312
+
+-- nagEllipticIntegralRD : (DF, DF, DF) -> DF ;
+
+nagEllipticIntegralRD(0.5,0.5,1) :: Float
+
+--       1.479
+
+nagEllipticIntegralRD(0.5,1,1) :: Float
+
+--       1.211
+
+nagEllipticIntegralRD(0.5,1.5,1) :: Float
+
+--       1.061
+
+nagEllipticIntegralRD(1,1,1) :: Float
+
+--       1.0
+
+nagEllipticIntegralRD(1,1.5,1) :: Float
+
+--       0.8805
+
+nagEllipticIntegralRD(1.5,1.5,1) :: Float
+
+--       0.7775
+
+-- nagEllipticIntegralRF : (DF, DF, DF) -> DF ;
+
+nagEllipticIntegralRF(0.5,1,1.5) :: Float
+
+--       1.028
+
+nagEllipticIntegralRF(1,1.5,2) :: Float
+
+--       0.826
+
+nagEllipticIntegralRF(1.5,2,2.5) :: Float
+
+--       0.7116
+
+-- nagEllipticIntegralRJ : (DF, DF, DF, DF) -> DF ;
+
+nagEllipticIntegralRJ(0.5,0.5,0.5,2) :: Float
+
+--       1.118
+
+nagEllipticIntegralRJ(0.5,0.5,1,2) :: Float
+
+--       0.9221
+
+nagEllipticIntegralRJ(0.5,0.5,1.5,2) :: Float
+
+--       0.8115
+
+nagEllipticIntegralRJ(0.5,1,1,2) :: Float
+
+--       0.7671
+
+nagEllipticIntegralRJ(0.5,1,1.5,2) :: Float
+
+--       0.6784
+
+nagEllipticIntegralRJ(0.5,1.5,1.5,2) :: Float
+
+--       0.6017
+
+nagEllipticIntegralRJ(1,1,1,2) :: Float
+
+--       0.6438
+
+nagEllipticIntegralRJ(1,1,1.5,2) :: Float
+
+--       0.5722
+
+nagEllipticIntegralRJ(1,1.5,1.5,2) :: Float
+
+--       0.5101
+
+nagEllipticIntegralRJ(1.5,1.5,1.5,2) :: Float
+
+--       0.4561
+
+outputGeneral()
+
+output "End of tests"
+
+@
+\eject
+\begin{thebibliography}{99}
+\bibitem{1} nothing
+\end{thebibliography}
+\end{document}