~ubuntu-branches/debian/squeeze/maxima/squeeze

« back to all changes in this revision

Viewing changes to src/numerical/slatec/d9b0mp.lisp

  • Committer: Bazaar Package Importer
  • Author(s): Camm Maguire
  • Date: 2006-10-18 14:52:42 UTC
  • mto: (1.1.5 upstream)
  • mto: This revision was merged to the branch mainline in revision 4.
  • Revision ID: james.westby@ubuntu.com-20061018145242-vzyrm5hmxr8kiosf
ImportĀ upstreamĀ versionĀ 5.10.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
;;; Compiled by f2cl version 2.0 beta 2002-05-06
 
1
;;; Compiled by f2cl version 2.0 beta Date: 2006/01/31 15:11:05 
 
2
;;; Using Lisp CMU Common Lisp Snapshot 2006-01 (19C)
2
3
;;; 
3
4
;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
4
5
;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
5
6
;;;           (:array-slicing nil) (:declare-common nil)
6
7
;;;           (:float-format double-float))
7
8
 
8
 
(in-package "SLATEC")
 
9
(in-package :slatec)
9
10
 
10
11
 
11
12
(let ((nbm0 0)
13
14
      (nbm02 0)
14
15
      (nbth0 0)
15
16
      (xmax 0.0)
16
 
      (bm0cs (make-array 37 :element-type 'double-float))
17
 
      (bth0cs (make-array 44 :element-type 'double-float))
18
 
      (bm02cs (make-array 40 :element-type 'double-float))
19
 
      (bt02cs (make-array 39 :element-type 'double-float))
 
17
      (bm0cs
 
18
       (make-array 37
 
19
                   :element-type 'double-float
 
20
                   :initial-contents '(0.09211656246827743
 
21
                                       -0.0010505909972719051
 
22
                                       1.4701598407687597e-5
 
23
                                       -5.058557606038554e-7
 
24
                                       2.7872545386324443e-8
 
25
                                       -2.062363611780915e-9
 
26
                                       1.8702143131388797e-10
 
27
                                       -1.9693309711356362e-11
 
28
                                       2.3259737939992754e-12
 
29
                                       -3.0095203449382503e-13
 
30
                                       4.1945213338506693e-14
 
31
                                       -6.219449312188446e-15
 
32
                                       9.71826041133607e-16
 
33
                                       -1.5884785857010752e-16
 
34
                                       2.7000721936713088e-17
 
35
                                       -4.7500923652340086e-18
 
36
                                       8.615128162604371e-19
 
37
                                       -1.6056086869561448e-19
 
38
                                       3.066513987314483e-20
 
39
                                       -5.987764223193956e-21
 
40
                                       1.1929712537482484e-21
 
41
                                       -2.4209691420448057e-22
 
42
                                       4.9967517605106164e-23
 
43
                                       -1.0474936393511585e-23
 
44
                                       2.2277868437974682e-24
 
45
                                       -4.801813239398163e-25
 
46
                                       1.04796272347096e-25
 
47
                                       -2.3138581656786152e-26
 
48
                                       5.164823088462674e-27
 
49
                                       -1.1646911918500655e-27
 
50
                                       2.651788486043319e-28
 
51
                                       -6.092559503825728e-29
 
52
                                       1.4118046861442593e-29
 
53
                                       -3.298094961231737e-30
 
54
                                       7.763931143074065e-31
 
55
                                       -1.8410313436614585e-31
 
56
                                       4.395880138594311e-32)))
 
57
      (bth0cs
 
58
       (make-array 44
 
59
                   :element-type 'double-float
 
60
                   :initial-contents '(-0.24901780862128936
 
61
                                       4.855029960962375e-4
 
62
                                       -5.4511837345017206e-6
 
63
                                       1.3558673059405963e-7
 
64
                                       -5.5691398902227624e-9
 
65
                                       3.2609031824994337e-10
 
66
                                       -2.491880786246134e-11
 
67
                                       2.344937742088252e-12
 
68
                                       -2.609653444431039e-13
 
69
                                       3.3353140420097393e-14
 
70
                                       -4.789000044057268e-15
 
71
                                       7.595617843619222e-16
 
72
                                       -1.3131556016891442e-16
 
73
                                       2.4483618345240857e-17
 
74
                                       -4.880572981061878e-18
 
75
                                       1.0327285029786316e-18
 
76
                                       -2.3057633815057217e-19
 
77
                                       5.404444300189269e-20
 
78
                                       -1.3240695194366572e-20
 
79
                                       3.378079562137197e-21
 
80
                                       -8.945762915711178e-22
 
81
                                       2.4519906889219317e-22
 
82
                                       -6.938842287686632e-23
 
83
                                       2.022827871489014e-23
 
84
                                       -6.062850000233549e-24
 
85
                                       1.8649748964037634e-24
 
86
                                       -5.878373238484989e-25
 
87
                                       1.8958591447999562e-25
 
88
                                       -6.248197937225886e-26
 
89
                                       2.1017901684551025e-26
 
90
                                       -7.208430093520926e-27
 
91
                                       2.5181363892474242e-27
 
92
                                       -8.951804225878578e-28
 
93
                                       3.23572374797623e-28
 
94
                                       -1.1883010519855353e-28
 
95
                                       4.4306286907358106e-29
 
96
                                       -1.676100964883483e-29
 
97
                                       6.429294692120746e-30
 
98
                                       -2.499226116697865e-30
 
99
                                       9.839979429952196e-31
 
100
                                       -3.9220375242408017e-31
 
101
                                       1.5818107030056521e-31
 
102
                                       -6.452550614489072e-32
 
103
                                       2.6611111369199356e-32)))
 
104
      (bm02cs
 
105
       (make-array 40
 
106
                   :element-type 'double-float
 
107
                   :initial-contents '(0.09500415145228382
 
108
                                       -3.801864682365671e-4
 
109
                                       2.258339301031481e-6
 
110
                                       -3.895725802372229e-8
 
111
                                       1.2468864165120817e-9
 
112
                                       -6.065949022102504e-11
 
113
                                       4.008461651421747e-12
 
114
                                       -3.3509981833980945e-13
 
115
                                       3.3771197165174173e-14
 
116
                                       -3.964585901635013e-15
 
117
                                       5.286111503883857e-16
 
118
                                       -7.852519083450852e-17
 
119
                                       1.2803005733866823e-17
 
120
                                       -2.26399629639143e-18
 
121
                                       4.3004969296567905e-19
 
122
                                       -8.705749805132587e-20
 
123
                                       1.8658627139620952e-20
 
124
                                       -4.210482486093065e-21
 
125
                                       9.9566769642284e-22
 
126
                                       -2.457357442805313e-22
 
127
                                       6.307692160762032e-23
 
128
                                       -1.6787736914407402e-23
 
129
                                       4.6202590646739044e-24
 
130
                                       -1.3117822668603088e-24
 
131
                                       3.834087564116303e-25
 
132
                                       -1.1514593240777412e-25
 
133
                                       3.547210007523339e-26
 
134
                                       -1.1192183858150046e-26
 
135
                                       3.611879427629838e-27
 
136
                                       -1.1906877659133332e-27
 
137
                                       4.005094059403968e-28
 
138
                                       -1.3731694224522124e-28
 
139
                                       4.7941990887425316e-29
 
140
                                       -1.7029656276241095e-29
 
141
                                       6.14951242893633e-30
 
142
                                       -2.2557668965818283e-30
 
143
                                       8.3997075092943e-31
 
144
                                       -3.1729975955626022e-31
 
145
                                       1.2152052988812985e-31
 
146
                                       -4.7158527497544386e-32)))
 
147
      (bt02cs
 
148
       (make-array 39
 
149
                   :element-type 'double-float
 
150
                   :initial-contents '(-0.24548295213424598
 
151
                                       0.0012544121039084616
 
152
                                       -3.1253950414871526e-5
 
153
                                       1.4709778249940832e-6
 
154
                                       -9.954348893795004e-8
 
155
                                       8.549316673320304e-9
 
156
                                       -8.698975952655434e-10
 
157
                                       1.0052099533559791e-10
 
158
                                       -1.2828230601708893e-11
 
159
                                       1.7731700781805131e-12
 
160
                                       -2.617457456948558e-13
 
161
                                       4.082835138997206e-14
 
162
                                       -6.675166823974272e-15
 
163
                                       1.136576139307163e-15
 
164
                                       -2.005118962064716e-16
 
165
                                       3.649797879476627e-17
 
166
                                       -6.8309637564582306e-18
 
167
                                       1.3107583145670756e-18
 
168
                                       -2.5723363101850606e-19
 
169
                                       5.152165744186396e-20
 
170
                                       -1.0513017563758802e-20
 
171
                                       2.1820381991194814e-21
 
172
                                       -4.600470121036216e-22
 
173
                                       9.840700692546682e-23
 
174
                                       -2.1334038035728374e-23
 
175
                                       4.683103642397336e-24
 
176
                                       -1.0400213691985747e-24
 
177
                                       2.334910567730151e-25
 
178
                                       -5.295682532331862e-26
 
179
                                       1.2126341952959756e-26
 
180
                                       -2.801889708228943e-27
 
181
                                       6.529267898701287e-28
 
182
                                       -1.5337980061873346e-28
 
183
                                       3.6305884306364537e-29
 
184
                                       -8.656075571362912e-30
 
185
                                       2.0779909972536286e-30
 
186
                                       -5.021117022141722e-31
 
187
                                       1.2208360279441715e-31
 
188
                                       -2.986005626703991e-32)))
20
189
      (pi4 0.7853981633974483)
21
 
      (first nil))
22
 
  (declare (type f2cl-lib:logical first)
 
190
      (first$ nil))
 
191
  (declare (type f2cl-lib:logical first$)
23
192
           (type (simple-array double-float (39)) bt02cs)
24
193
           (type (simple-array double-float (40)) bm02cs)
25
194
           (type (simple-array double-float (44)) bth0cs)
26
195
           (type (simple-array double-float (37)) bm0cs)
27
196
           (type double-float pi4 xmax)
28
197
           (type f2cl-lib:integer4 nbth0 nbm02 nbt02 nbm0))
29
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (1) ((1 37))) 0.09211656246827743)
30
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (2) ((1 37))) -0.0010505909972719051)
31
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (3) ((1 37))) 1.4701598407687597e-5)
32
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (4) ((1 37))) -5.058557606038554e-7)
33
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (5) ((1 37))) 2.7872545386324443e-8)
34
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (6) ((1 37))) -2.062363611780915e-9)
35
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (7) ((1 37))) 1.87021431313888e-10)
36
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (8) ((1 37))) -1.9693309711356366e-11)
37
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (9) ((1 37))) 2.3259737939992753e-12)
38
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (10) ((1 37))) -3.0095203449382507e-13)
39
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (11) ((1 37))) 4.194521333850669e-14)
40
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (12) ((1 37))) -6.219449312188445e-15)
41
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (13) ((1 37))) 9.718260411336069e-16)
42
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (14) ((1 37))) -1.5884785857010755e-16)
43
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (15) ((1 37))) 2.700072193671309e-17)
44
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (16) ((1 37))) -4.7500923652340094e-18)
45
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (17) ((1 37))) 8.615128162604372e-19)
46
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (18) ((1 37))) -1.6056086869561448e-19)
47
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (19) ((1 37))) 3.066513987314483e-20)
48
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (20) ((1 37))) -5.987764223193957e-21)
49
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (21) ((1 37))) 1.1929712537482481e-21)
50
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (22) ((1 37))) -2.4209691420448054e-22)
51
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (23) ((1 37))) 4.996751760510617e-23)
52
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (24) ((1 37))) -1.0474936393511587e-23)
53
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (25) ((1 37))) 2.2277868437974682e-24)
54
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (26) ((1 37))) -4.8018132393981616e-25)
55
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (27) ((1 37))) 1.0479627234709601e-25)
56
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (28) ((1 37))) -2.3138581656786159e-26)
57
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (29) ((1 37))) 5.164823088462674e-27)
58
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (30) ((1 37))) -1.1646911918500655e-27)
59
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (31) ((1 37))) 2.65178848604332e-28)
60
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (32) ((1 37))) -6.092559503825729e-29)
61
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (33) ((1 37))) 1.4118046861442593e-29)
62
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (34) ((1 37))) -3.298094961231737e-30)
63
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (35) ((1 37))) 7.763931143074065e-31)
64
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (36) ((1 37))) -1.8410313436614584e-31)
65
 
  (f2cl-lib:fset (f2cl-lib:fref bm0cs (37) ((1 37))) 4.3958801385943114e-32)
66
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (1) ((1 44))) -0.24901780862128936)
67
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (2) ((1 44))) 4.855029960962375e-4)
68
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (3) ((1 44))) -5.45118373450172e-6)
69
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (4) ((1 44))) 1.3558673059405965e-7)
70
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (5) ((1 44))) -5.569139890222763e-9)
71
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (6) ((1 44))) 3.260903182499434e-10)
72
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (7) ((1 44))) -2.4918807862461345e-11)
73
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (8) ((1 44))) 2.344937742088252e-12)
74
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (9) ((1 44))) -2.609653444431039e-13)
75
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (10) ((1 44))) 3.3353140420097394e-14)
76
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (11) ((1 44))) -4.789000044057268e-15)
77
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (12) ((1 44))) 7.595617843619222e-16)
78
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (13) ((1 44))) -1.313155601689144e-16)
79
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (14) ((1 44))) 2.4483618345240857e-17)
80
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (15) ((1 44))) -4.880572981061878e-18)
81
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (16) ((1 44))) 1.0327285029786318e-18)
82
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (17) ((1 44))) -2.305763381505722e-19)
83
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (18) ((1 44))) 5.404444300189268e-20)
84
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (19) ((1 44))) -1.3240695194366572e-20)
85
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (20) ((1 44))) 3.3780795621371973e-21)
86
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (21) ((1 44))) -8.94576291571118e-22)
87
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (22) ((1 44))) 2.4519906889219317e-22)
88
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (23) ((1 44))) -6.938842287686633e-23)
89
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (24) ((1 44))) 2.0228278714890138e-23)
90
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (25) ((1 44))) -6.062850000233548e-24)
91
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (26) ((1 44))) 1.8649748964037632e-24)
92
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (27) ((1 44))) -5.878373238484988e-25)
93
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (28) ((1 44))) 1.895859144799956e-25)
94
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (29) ((1 44))) -6.248197937225886e-26)
95
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (30) ((1 44))) 2.101790168455103e-26)
96
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (31) ((1 44))) -7.208430093520926e-27)
97
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (32) ((1 44))) 2.5181363892474246e-27)
98
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (33) ((1 44))) -8.951804225878578e-28)
99
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (34) ((1 44))) 3.2357237479762296e-28)
100
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (35) ((1 44))) -1.1883010519855355e-28)
101
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (36) ((1 44))) 4.43062869073581e-29)
102
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (37) ((1 44))) -1.676100964883483e-29)
103
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (38) ((1 44))) 6.429294692120746e-30)
104
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (39) ((1 44))) -2.499226116697865e-30)
105
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (40) ((1 44))) 9.839979429952197e-31)
106
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (41) ((1 44))) -3.922037524240802e-31)
107
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (42) ((1 44))) 1.5818107030056522e-31)
108
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (43) ((1 44))) -6.452550614489073e-32)
109
 
  (f2cl-lib:fset (f2cl-lib:fref bth0cs (44) ((1 44))) 2.6611111369199364e-32)
110
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (1) ((1 40))) 0.09500415145228382)
111
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (2) ((1 40))) -3.801864682365671e-4)
112
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (3) ((1 40))) 2.2583393010314814e-6)
113
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (4) ((1 40))) -3.895725802372228e-8)
114
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (5) ((1 40))) 1.2468864165120816e-9)
115
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (6) ((1 40))) -6.065949022102504e-11)
116
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (7) ((1 40))) 4.0084616514217464e-12)
117
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (8) ((1 40))) -3.350998183398094e-13)
118
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (9) ((1 40))) 3.377119716517417e-14)
119
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (10) ((1 40))) -3.9645859016350127e-15)
120
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (11) ((1 40))) 5.286111503883857e-16)
121
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (12) ((1 40))) -7.852519083450852e-17)
122
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (13) ((1 40))) 1.280300573386682e-17)
123
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (14) ((1 40))) -2.2639962963914298e-18)
124
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (15) ((1 40))) 4.30049692965679e-19)
125
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (16) ((1 40))) -8.705749805132588e-20)
126
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (17) ((1 40))) 1.8658627139620948e-20)
127
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (18) ((1 40))) -4.2104824860930656e-21)
128
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (19) ((1 40))) 9.956676964228401e-22)
129
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (20) ((1 40))) -2.4573574428053133e-22)
130
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (21) ((1 40))) 6.307692160762031e-23)
131
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (22) ((1 40))) -1.6787736914407403e-23)
132
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (23) ((1 40))) 4.6202590646739045e-24)
133
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (24) ((1 40))) -1.3117822668603088e-24)
134
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (25) ((1 40))) 3.8340875641163025e-25)
135
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (26) ((1 40))) -1.1514593240777413e-25)
136
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (27) ((1 40))) 3.5472100075233387e-26)
137
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (28) ((1 40))) -1.1192183858150047e-26)
138
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (29) ((1 40))) 3.611879427629838e-27)
139
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (30) ((1 40))) -1.1906877659133333e-27)
140
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (31) ((1 40))) 4.0050940594039686e-28)
141
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (32) ((1 40))) -1.3731694224522123e-28)
142
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (33) ((1 40))) 4.794199088742532e-29)
143
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (34) ((1 40))) -1.7029656276241095e-29)
144
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (35) ((1 40))) 6.14951242893633e-30)
145
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (36) ((1 40))) -2.2557668965818284e-30)
146
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (37) ((1 40))) 8.3997075092943e-31)
147
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (38) ((1 40))) -3.1729975955626016e-31)
148
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (39) ((1 40))) 1.2152052988812986e-31)
149
 
  (f2cl-lib:fset (f2cl-lib:fref bm02cs (40) ((1 40))) -4.7158527497544395e-32)
150
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (1) ((1 39))) -0.24548295213424598)
151
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (2) ((1 39))) 0.0012544121039084616)
152
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (3) ((1 39))) -3.1253950414871523e-5)
153
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (4) ((1 39))) 1.470977824994083e-6)
154
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (5) ((1 39))) -9.954348893795003e-8)
155
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (6) ((1 39))) 8.549316673320305e-9)
156
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (7) ((1 39))) -8.698975952655433e-10)
157
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (8) ((1 39))) 1.0052099533559791e-10)
158
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (9) ((1 39))) -1.2828230601708893e-11)
159
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (10) ((1 39))) 1.7731700781805132e-12)
160
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (11) ((1 39))) -2.6174574569485576e-13)
161
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (12) ((1 39))) 4.082835138997206e-14)
162
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (13) ((1 39))) -6.675166823974272e-15)
163
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (14) ((1 39))) 1.136576139307163e-15)
164
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (15) ((1 39))) -2.0051189620647164e-16)
165
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (16) ((1 39))) 3.6497978794766267e-17)
166
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (17) ((1 39))) -6.83096375645823e-18)
167
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (18) ((1 39))) 1.310758314567076e-18)
168
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (19) ((1 39))) -2.572336310185061e-19)
169
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (20) ((1 39))) 5.152165744186394e-20)
170
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (21) ((1 39))) -1.0513017563758804e-20)
171
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (22) ((1 39))) 2.1820381991194815e-21)
172
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (23) ((1 39))) -4.6004701210362153e-22)
173
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (24) ((1 39))) 9.840700692546682e-23)
174
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (25) ((1 39))) -2.1334038035728375e-23)
175
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (26) ((1 39))) 4.6831036423973366e-24)
176
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (27) ((1 39))) -1.0400213691985746e-24)
177
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (28) ((1 39))) 2.334910567730151e-25)
178
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (29) ((1 39))) -5.295682532331861e-26)
179
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (30) ((1 39))) 1.2126341952959759e-26)
180
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (31) ((1 39))) -2.801889708228943e-27)
181
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (32) ((1 39))) 6.529267898701289e-28)
182
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (33) ((1 39))) -1.5337980061873346e-28)
183
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (34) ((1 39))) 3.6305884306364533e-29)
184
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (35) ((1 39))) -8.656075571362912e-30)
185
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (36) ((1 39))) 2.0779909972536278e-30)
186
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (37) ((1 39))) -5.021117022141722e-31)
187
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (38) ((1 39))) 1.2208360279441714e-31)
188
 
  (f2cl-lib:fset (f2cl-lib:fref bt02cs (39) ((1 39))) -2.986005626703992e-32)
189
 
  (setq first f2cl-lib:%true%)
 
198
  (setq first$ f2cl-lib:%true%)
190
199
  (defun d9b0mp (x ampl theta)
191
200
    (declare (type double-float theta ampl x))
192
201
    (prog ((z 0.0) (eta 0.0f0))
193
202
      (declare (type single-float eta) (type double-float z))
194
203
      (cond
195
 
       (first (setf eta (* 0.1f0 (f2cl-lib:freal (f2cl-lib:d1mach 3))))
196
 
              (setf nbm0 (initds bm0cs 37 eta))
197
 
              (setf nbt02 (initds bt02cs 39 eta))
198
 
              (setf nbm02 (initds bm02cs 40 eta))
199
 
              (setf nbth0 (initds bth0cs 44 eta))
200
 
              (setf xmax (/ 1.0 (f2cl-lib:d1mach 4)))))
201
 
      (setf first f2cl-lib:%false%)
 
204
        (first$
 
205
         (setf eta (* 0.1f0 (f2cl-lib:freal (f2cl-lib:d1mach 3))))
 
206
         (setf nbm0 (initds bm0cs 37 eta))
 
207
         (setf nbt02 (initds bt02cs 39 eta))
 
208
         (setf nbm02 (initds bm02cs 40 eta))
 
209
         (setf nbth0 (initds bth0cs 44 eta))
 
210
         (setf xmax (/ 1.0 (f2cl-lib:d1mach 4)))))
 
211
      (setf first$ f2cl-lib:%false%)
202
212
      (if (< x 4.0) (xermsg "SLATEC" "D9B0MP" "X MUST BE GE 4" 1 2))
203
213
      (if (> x 8.0) (go label20))
204
214
      (setf z (/ (- (/ 128.0 (* x x)) 5.0) 3.0))