← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/numerical/slatec/djairy.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:
src/numerical/slatec/djairy.lisp
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 ((n1 14)
20
21
      (con2 5.03154716196777)
21
22
      (con3 0.380004589867293)
22
23
      (con4 0.833333333333333)
23
 
      (con5 0.8660254037844392)
24
 
      (ak1 (make-array 14 :element-type 'double-float))
25
 
      (ak2 (make-array 23 :element-type 'double-float))
26
 
      (ak3 (make-array 14 :element-type 'double-float))
27
 
      (ajp (make-array 19 :element-type 'double-float))
28
 
      (ajn (make-array 19 :element-type 'double-float))
29
 
      (a (make-array 15 :element-type 'double-float))
30
 
      (b (make-array 15 :element-type 'double-float))
 
24
      (con5 0.866025403784439)
 
25
      (ak1
 
26
       (make-array 14
 
27
                   :element-type 'double-float
 
28
                   :initial-contents '(0.220423090987793 -0.1252902427877
 
29
                                       0.0103881163359194 8.22844152006343e-4
 
30
                                       -2.34614345891226e-4 1.63824280172116e-5
 
31
                                       3.06902589573189e-7 -1.29621999359332e-7
 
32
                                       8.22908158823668e-9 1.53963968623298e-11
 
33
                                       -3.39165465615682e-11
 
34
                                       2.03253257423626e-12
 
35
                                       -1.10679546097884e-14
 
36
                                       -5.1616949778508e-15)))
 
37
      (ak2
 
38
       (make-array 23
 
39
                   :element-type 'double-float
 
40
                   :initial-contents '(0.274366150869598 0.00539790969736903
 
41
                                       -0.0015733922062119 4.2742752824875e-4
 
42
                                       -1.12124917399925e-4 2.88763171318904e-5
 
43
                                       -7.36804225370554e-6 1.87290209741024e-6
 
44
                                       -4.75892793962291e-7 1.21130416955909e-7
 
45
                                       -3.09245374270614e-8 7.92454705282654e-9
 
46
                                       -2.03902447167914e-9
 
47
                                       5.26863056595742e-10
 
48
                                       -1.36704767639569e-10
 
49
                                       3.56141039013708e-11
 
50
                                       -9.3138829654843e-12
 
51
                                       2.44464450473635e-12
 
52
                                       -6.43840261990955e-13
 
53
                                       1.70106030559349e-13
 
54
                                       -4.50760104503281e-14
 
55
                                       1.19774799164811e-14
 
56
                                       -3.19077040865066e-15)))
 
57
      (ak3
 
58
       (make-array 14
 
59
                   :element-type 'double-float
 
60
                   :initial-contents '(0.280271447340791 -0.00178127042844379
 
61
                                       4.03422579628999e-5 -1.63249965269003e-6
 
62
                                       9.21181482476768e-8 -6.52294330229155e-9
 
63
                                       5.47138404576546e-10
 
64
                                       -5.2440825180026e-11
 
65
                                       5.60477904117209e-12
 
66
                                       -6.56375244639313e-13
 
67
                                       8.31285761966247e-14
 
68
                                       -1.12705134691063e-14
 
69
                                       1.62267976598129e-15
 
70
                                       -2.46480324312426e-16)))
 
71
      (ajp
 
72
       (make-array 19
 
73
                   :element-type 'double-float
 
74
                   :initial-contents '(0.0778952966437581 -0.184356363456801
 
75
                                       0.0301412605216174 0.0305342724277608
 
76
                                       -0.00495424702513079
 
77
                                       -0.00172749552563952 2.4313763783919e-4
 
78
                                       5.04564777517082e-5 -6.16316582695208e-6
 
79
                                       -9.03986745510768e-7 9.70243778355884e-8
 
80
                                       1.09639453305205e-8 -1.04716330588766e-9
 
81
                                       -9.60359441344646e-11
 
82
                                       8.25358789454134e-12
 
83
                                       6.36123439018768e-13
 
84
                                       -4.96629614116015e-14
 
85
                                       -3.29810288929615e-15
 
86
                                       2.35798252031104e-16)))
 
87
      (ajn
 
88
       (make-array 19
 
89
                   :element-type 'double-float
 
90
                   :initial-contents '(0.0380497887617242 -0.245319541845546
 
91
                                       0.165820623702696 0.0749330045818789
 
92
                                       -0.0263476288106641 -0.00592535597304981
 
93
                                       0.00144744409589804 2.18311831322215e-4
 
94
                                       -4.10662077680304e-5
 
95
                                       -4.66874994171766e-6 7.1521880727716e-7
 
96
                                       6.52964770854633e-8 -8.44284027565946e-9
 
97
                                       -6.44186158976978e-10
 
98
                                       7.20802286505285e-11
 
99
                                       4.72465431717846e-12
 
100
                                       -4.66022632547045e-13
 
101
                                       -2.67762710389189e-14
 
102
                                       2.36161316570019e-15)))
 
103
      (a
 
104
       (make-array 15
 
105
                   :element-type 'double-float
 
106
                   :initial-contents '(0.490275424742791 0.00157647277946204
 
107
                                       -9.66195963140306e-5 1.35916080268815e-7
 
108
                                       2.98157342654859e-7 -1.86824767559979e-8
 
109
                                       -1.03685737667141e-9
 
110
                                       3.28660818434328e-10
 
111
                                       -2.5709141063278e-11
 
112
                                       -2.32357655300677e-12
 
113
                                       9.57523279048255e-13
 
114
                                       -1.20340828049719e-13
 
115
                                       -2.90907716770715e-15
 
116
                                       4.55656454580149e-15
 
117
                                       -9.99003874810259e-16)))
 
118
      (b
 
119
       (make-array 15
 
120
                   :element-type 'double-float
 
121
                   :initial-contents '(0.278593552803079 -0.00352915691882584
 
122
                                       -2.31149677384994e-5 4.7131784226356e-6
 
123
                                       -1.12415907931333e-7
 
124
                                       -2.00100301184339e-8 2.60948075302193e-9
 
125
                                       -3.55098136101216e-11
 
126
                                       -3.50849978423875e-11
 
127
                                       5.83007187954202e-12
 
128
                                       -2.04644828753326e-13
 
129
                                       -1.10529179476742e-13
 
130
                                       2.87724778038775e-14
 
131
                                       -2.88205111009939e-15
 
132
                                       -3.32656311696166e-16)))
31
133
      (n1d 14)
32
134
      (n2d 24)
33
135
      (n3d 19)
36
138
      (m2d 22)
37
139
      (m3d 17)
38
140
      (m4d 13)
39
 
      (dak1 (make-array 14 :element-type 'double-float))
40
 
      (dak2 (make-array 24 :element-type 'double-float))
41
 
      (dak3 (make-array 14 :element-type 'double-float))
42
 
      (dajp (make-array 19 :element-type 'double-float))
43
 
      (dajn (make-array 19 :element-type 'double-float))
44
 
      (da (make-array 15 :element-type 'double-float))
45
 
      (db (make-array 15 :element-type 'double-float)))
 
141
      (dak1
 
142
       (make-array 14
 
143
                   :element-type 'double-float
 
144
                   :initial-contents '(0.204567842307887 -0.0661322739905664
 
145
                                       -0.00849845800989287 0.00312183491556289
 
146
                                       -2.70016489829432e-4
 
147
                                       -6.35636298679387e-6 3.02397712409509e-6
 
148
                                       -2.18311195330088e-7
 
149
                                       -5.36194289332826e-10 1.1309803562231e-9
 
150
                                       -7.43023834629073e-11
 
151
                                       4.28804170826891e-13
 
152
                                       2.23810925754539e-13
 
153
                                       -1.39140135641182e-14)))
 
154
      (dak2
 
155
       (make-array 24
 
156
                   :element-type 'double-float
 
157
                   :initial-contents '(0.29333234388323 -0.00806196784743112
 
158
                                       0.0024254017233314 -6.82297548850235e-4
 
159
                                       1.85786427751181e-4 -4.97457447684059e-5
 
160
                                       1.32090681239497e-5 -3.49528240444943e-6
 
161
                                       9.24362451078835e-7 -2.44732671521867e-7
 
162
                                       6.4930783764891e-8 -1.72717621501538e-8
 
163
                                       4.60725763604656e-9 -1.2324905529155e-9
 
164
                                       3.30620409488102e-10
 
165
                                       -8.89252099772401e-11
 
166
                                       2.39773319878298e-11
 
167
                                       -6.4801392115345e-12
 
168
                                       1.75510132023731e-12
 
169
                                       -4.76303829833637e-13
 
170
                                       1.2949824110081e-13 -3.5267962221043e-14
 
171
                                       9.62005151585923e-15
 
172
                                       -2.62786914342292e-15)))
 
173
      (dak3
 
174
       (make-array 14
 
175
                   :element-type 'double-float
 
176
                   :initial-contents '(0.284675828811349 0.0025307307261908
 
177
                                       -4.83481130337976e-5 1.84907283946343e-6
 
178
                                       -1.01418491178576e-7 7.05925634457153e-9
 
179
                                       -5.85325291400382e-10
 
180
                                       5.56357688831339e-11 -5.908890947795e-12
 
181
                                       6.88574353784436e-13
 
182
                                       -8.68588256452194e-14
 
183
                                       1.17374762617213e-14
 
184
                                       -1.68523146510923e-15
 
185
                                       2.55374773097056e-16)))
 
186
      (dajp
 
187
       (make-array 19
 
188
                   :element-type 'double-float
 
189
                   :initial-contents '(0.0653219131311457 -0.120262933688823
 
190
                                       0.00978010236263823 0.0167948429230505
 
191
                                       -0.00197146140182132
 
192
                                       -8.45560295098867e-4 9.42889620701976e-5
 
193
                                       2.25827860945475e-5 -2.29067870915987e-6
 
194
                                       -3.76343991136919e-7 3.45663933559565e-8
 
195
                                       4.29611332003007e-9
 
196
                                       -3.58673691214989e-10
 
197
                                       -3.57245881361895e-11
 
198
                                       2.72696091066336e-12
 
199
                                       2.26120653095771e-13
 
200
                                       -1.58763205238303e-14
 
201
                                       -1.12604374485125e-15
 
202
                                       7.31327529515367e-17)))
 
203
      (dajn
 
204
       (make-array 19
 
205
                   :element-type 'double-float
 
206
                   :initial-contents '(0.0108594539632967 0.0853313194857091
 
207
                                       -0.315277068113058 -0.0878420725294257
 
208
                                       0.0553251906976048 0.00941674060503241
 
209
                                       -0.00332187026018996
 
210
                                       -4.11157343156826e-4 1.01297326891346e-4
 
211
                                       9.87633682208396e-6 -1.87312969812393e-6
 
212
                                       -1.50798500131468e-7 2.32687669525394e-8
 
213
                                       1.59599917419225e-9
 
214
                                       -2.07665922668385e-10
 
215
                                       -1.24103350500302e-11
 
216
                                       1.39631765331043e-12 7.3940097115574e-14
 
217
                                       -7.328874756275e-15)))
 
218
      (da
 
219
       (make-array 15
 
220
                   :element-type 'double-float
 
221
                   :initial-contents '(0.491627321104601 0.00311164930427489
 
222
                                       8.23140762854081e-5 -4.61769776172142e-6
 
223
                                       -6.13158880534626e-8 2.8729580465652e-8
 
224
                                       -1.81959715372117e-9
 
225
                                       -1.44752826642035e-10
 
226
                                       4.53724043420422e-11
 
227
                                       -3.99655065847223e-12
 
228
                                       -3.24089119830323e-13
 
229
                                       1.62098952568741e-13
 
230
                                       -2.40765247974057e-14
 
231
                                       1.69384811284491e-16
 
232
                                       8.17900786477396e-16)))
 
233
      (db
 
234
       (make-array 15
 
235
                   :element-type 'double-float
 
236
                   :initial-contents '(-0.277571356944231 0.0044421283341992
 
237
                                       -8.42328522190089e-5 -2.5804031841871e-6
 
238
                                       3.42389720217621e-7 -6.24286894709776e-9
 
239
                                       -2.36377836844577e-9
 
240
                                       3.16991042656673e-10
 
241
                                       -4.40995691658191e-12
 
242
                                       -5.18674221093575e-12
 
243
                                       9.64874015137022e-13
 
244
                                       -4.9019057660871e-14
 
245
                                       -1.77253430678112e-14
 
246
                                       5.55950610442662e-15
 
247
                                       -7.1179333757953e-16))))
46
248
  (declare (type (simple-array double-float (24)) dak2)
47
249
           (type (simple-array double-float (15)) db da b a)
48
250
           (type (simple-array double-float (19)) dajn dajp ajn ajp)
50
252
           (type (simple-array double-float (14)) dak3 dak1 ak3 ak1)
51
253
           (type double-float con5 con4 con3 con2 fpi12)
52
254
           (type f2cl-lib:integer4 m4d m3d m2d m1d n4d n3d n2d n1d m4 m3 m2 m1
53
 
            n4 n3 n2 n1))
54
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (1) ((1 14))) 0.220423090987793)
55
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (2) ((1 14))) -0.1252902427877)
56
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (3) ((1 14))) 0.0103881163359194)
57
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (4) ((1 14))) 8.228441520063431e-4)
58
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (5) ((1 14))) -2.3461434589122598e-4)
59
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (6) ((1 14))) 1.6382428017211603e-5)
60
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (7) ((1 14))) 3.06902589573189e-7)
61
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (8) ((1 14))) -1.29621999359332e-7)
62
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (9) ((1 14))) 8.229081588236679e-9)
63
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (10) ((1 14))) 1.5396396862329798e-11)
64
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (11) ((1 14))) -3.3916546561568195e-11)
65
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (12) ((1 14))) 2.0325325742362604e-12)
66
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (13) ((1 14))) -1.10679546097884e-14)
67
 
  (f2cl-lib:fset (f2cl-lib:fref ak1 (14) ((1 14))) -5.1616949778508e-15)
68
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (1) ((1 23))) 0.274366150869598)
69
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (2) ((1 23))) 0.005397909697369031)
70
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (3) ((1 23))) -0.0015733922062118998)
71
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (4) ((1 23))) 4.274275282487501e-4)
72
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (5) ((1 23))) -1.12124917399925e-4)
73
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (6) ((1 23))) 2.8876317131890405e-5)
74
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (7) ((1 23))) -7.368042253705539e-6)
75
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (8) ((1 23))) 1.87290209741024e-6)
76
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (9) ((1 23))) -4.75892793962291e-7)
77
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (10) ((1 23))) 1.2113041695590901e-7)
78
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (11) ((1 23))) -3.09245374270614e-8)
79
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (12) ((1 23))) 7.924547052826542e-9)
80
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (13) ((1 23))) -2.0390244716791406e-9)
81
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (14) ((1 23))) 5.268630565957421e-10)
82
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (15) ((1 23))) -1.36704767639569e-10)
83
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (16) ((1 23))) 3.561410390137079e-11)
84
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (17) ((1 23))) -9.3138829654843e-12)
85
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (18) ((1 23))) 2.44464450473635e-12)
86
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (19) ((1 23))) -6.43840261990955e-13)
87
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (20) ((1 23))) 1.7010603055934903e-13)
88
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (21) ((1 23))) -4.5076010450328097e-14)
89
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (22) ((1 23))) 1.19774799164811e-14)
90
 
  (f2cl-lib:fset (f2cl-lib:fref ak2 (23) ((1 23))) -3.1907704086506605e-15)
91
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (1) ((1 14))) 0.280271447340791)
92
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (2) ((1 14))) -0.00178127042844379)
93
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (3) ((1 14))) 4.034225796289991e-5)
94
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (4) ((1 14))) -1.6324996526900298e-6)
95
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (5) ((1 14))) 9.21181482476768e-8)
96
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (6) ((1 14))) -6.522943302291551e-9)
97
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (7) ((1 14))) 5.471384045765461e-10)
98
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (8) ((1 14))) -5.244082518002599e-11)
99
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (9) ((1 14))) 5.60477904117209e-12)
100
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (10) ((1 14))) -6.56375244639313e-13)
101
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (11) ((1 14))) 8.31285761966247e-14)
102
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (12) ((1 14))) -1.12705134691063e-14)
103
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (13) ((1 14))) 1.6226797659812903e-15)
104
 
  (f2cl-lib:fset (f2cl-lib:fref ak3 (14) ((1 14))) -2.4648032431242603e-16)
105
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (1) ((1 19))) 0.0778952966437581)
106
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (2) ((1 19))) -0.184356363456801)
107
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (3) ((1 19))) 0.030141260521617398)
108
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (4) ((1 19))) 0.0305342724277608)
109
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (5) ((1 19))) -0.004954247025130791)
110
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (6) ((1 19))) -0.0017274955256395199)
111
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (7) ((1 19))) 2.4313763783918999e-4)
112
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (8) ((1 19))) 5.045647775170821e-5)
113
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (9) ((1 19))) -6.163165826952081e-6)
114
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (10) ((1 19))) -9.03986745510768e-7)
115
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (11) ((1 19))) 9.70243778355884e-8)
116
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (12) ((1 19))) 1.09639453305205e-8)
117
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (13) ((1 19))) -1.0471633058876599e-9)
118
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (14) ((1 19))) -9.603594413446459e-11)
119
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (15) ((1 19))) 8.253587894541341e-12)
120
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (16) ((1 19))) 6.36123439018768e-13)
121
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (17) ((1 19))) -4.9662961411601503e-14)
122
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (18) ((1 19))) -3.29810288929615e-15)
123
 
  (f2cl-lib:fset (f2cl-lib:fref ajp (19) ((1 19))) 2.3579825203110397e-16)
124
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (1) ((1 19))) 0.038049788761724206)
125
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (2) ((1 19))) -0.24531954184554602)
126
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (3) ((1 19))) 0.165820623702696)
127
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (4) ((1 19))) 0.0749330045818789)
128
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (5) ((1 19))) -0.0263476288106641)
129
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (6) ((1 19))) -0.00592535597304981)
130
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (7) ((1 19))) 0.00144744409589804)
131
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (8) ((1 19))) 2.18311831322215e-4)
132
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (9) ((1 19))) -4.1066207768030405e-5)
133
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (10) ((1 19))) -4.6687499417176603e-6)
134
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (11) ((1 19))) 7.152188072771599e-7)
135
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (12) ((1 19))) 6.52964770854633e-8)
136
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (13) ((1 19))) -8.442840275659461e-9)
137
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (14) ((1 19))) -6.441861589769781e-10)
138
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (15) ((1 19))) 7.208022865052851e-11)
139
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (16) ((1 19))) 4.72465431717846e-12)
140
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (17) ((1 19))) -4.6602263254704496e-13)
141
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (18) ((1 19))) -2.6776271038918903e-14)
142
 
  (f2cl-lib:fset (f2cl-lib:fref ajn (19) ((1 19))) 2.3616131657001907e-15)
143
 
  (f2cl-lib:fset (f2cl-lib:fref a (1) ((1 15))) 0.490275424742791)
144
 
  (f2cl-lib:fset (f2cl-lib:fref a (2) ((1 15))) 0.00157647277946204)
145
 
  (f2cl-lib:fset (f2cl-lib:fref a (3) ((1 15))) -9.66195963140306e-5)
146
 
  (f2cl-lib:fset (f2cl-lib:fref a (4) ((1 15))) 1.35916080268815e-7)
147
 
  (f2cl-lib:fset (f2cl-lib:fref a (5) ((1 15))) 2.98157342654859e-7)
148
 
  (f2cl-lib:fset (f2cl-lib:fref a (6) ((1 15))) -1.8682476755997898e-8)
149
 
  (f2cl-lib:fset (f2cl-lib:fref a (7) ((1 15))) -1.0368573766714102e-9)
150
 
  (f2cl-lib:fset (f2cl-lib:fref a (8) ((1 15))) 3.28660818434328e-10)
151
 
  (f2cl-lib:fset (f2cl-lib:fref a (9) ((1 15))) -2.5709141063278e-11)
152
 
  (f2cl-lib:fset (f2cl-lib:fref a (10) ((1 15))) -2.3235765530067695e-12)
153
 
  (f2cl-lib:fset (f2cl-lib:fref a (11) ((1 15))) 9.575232790482551e-13)
154
 
  (f2cl-lib:fset (f2cl-lib:fref a (12) ((1 15))) -1.2034082804971902e-13)
155
 
  (f2cl-lib:fset (f2cl-lib:fref a (13) ((1 15))) -2.9090771677071503e-15)
156
 
  (f2cl-lib:fset (f2cl-lib:fref a (14) ((1 15))) 4.5565645458014903e-15)
157
 
  (f2cl-lib:fset (f2cl-lib:fref a (15) ((1 15))) -9.99003874810259e-16)
158
 
  (f2cl-lib:fset (f2cl-lib:fref b (1) ((1 15))) 0.278593552803079)
159
 
  (f2cl-lib:fset (f2cl-lib:fref b (2) ((1 15))) -0.00352915691882584)
160
 
  (f2cl-lib:fset (f2cl-lib:fref b (3) ((1 15))) -2.31149677384994e-5)
161
 
  (f2cl-lib:fset (f2cl-lib:fref b (4) ((1 15))) 4.7131784226356e-6)
162
 
  (f2cl-lib:fset (f2cl-lib:fref b (5) ((1 15))) -1.12415907931333e-7)
163
 
  (f2cl-lib:fset (f2cl-lib:fref b (6) ((1 15))) -2.0010030118433902e-8)
164
 
  (f2cl-lib:fset (f2cl-lib:fref b (7) ((1 15))) 2.6094807530219305e-9)
165
 
  (f2cl-lib:fset (f2cl-lib:fref b (8) ((1 15))) -3.55098136101216e-11)
166
 
  (f2cl-lib:fset (f2cl-lib:fref b (9) ((1 15))) -3.50849978423875e-11)
167
 
  (f2cl-lib:fset (f2cl-lib:fref b (10) ((1 15))) 5.830071879542019e-12)
168
 
  (f2cl-lib:fset (f2cl-lib:fref b (11) ((1 15))) -2.04644828753326e-13)
169
 
  (f2cl-lib:fset (f2cl-lib:fref b (12) ((1 15))) -1.1052917947674198e-13)
170
 
  (f2cl-lib:fset (f2cl-lib:fref b (13) ((1 15))) 2.87724778038775e-14)
171
 
  (f2cl-lib:fset (f2cl-lib:fref b (14) ((1 15))) -2.8820511100993906e-15)
172
 
  (f2cl-lib:fset (f2cl-lib:fref b (15) ((1 15))) -3.32656311696166e-16)
173
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (1) ((1 14))) 0.204567842307887)
174
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (2) ((1 14))) -0.0661322739905664)
175
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (3) ((1 14))) -0.00849845800989287)
176
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (4) ((1 14))) 0.0031218349155628902)
177
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (5) ((1 14))) -2.70016489829432e-4)
178
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (6) ((1 14))) -6.356362986793869e-6)
179
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (7) ((1 14))) 3.0239771240950897e-6)
180
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (8) ((1 14))) -2.1831119533008797e-7)
181
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (9) ((1 14))) -5.361942893328261e-10)
182
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (10) ((1 14))) 1.1309803562231e-9)
183
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (11) ((1 14))) -7.43023834629073e-11)
184
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (12) ((1 14))) 4.2880417082689104e-13)
185
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (13) ((1 14))) 2.23810925754539e-13)
186
 
  (f2cl-lib:fset (f2cl-lib:fref dak1 (14) ((1 14))) -1.39140135641182e-14)
187
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (1) ((1 24))) 0.29333234388323)
188
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (2) ((1 24))) -0.00806196784743112)
189
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (3) ((1 24))) 0.0024254017233314)
190
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (4) ((1 24))) -6.82297548850235e-4)
191
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (5) ((1 24))) 1.85786427751181e-4)
192
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (6) ((1 24))) -4.974574476840591e-5)
193
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (7) ((1 24))) 1.3209068123949702e-5)
194
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (8) ((1 24))) -3.4952824044494296e-6)
195
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (9) ((1 24))) 9.24362451078835e-7)
196
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (10) ((1 24))) -2.4473267152186698e-7)
197
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (11) ((1 24))) 6.4930783764891e-8)
198
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (12) ((1 24))) -1.72717621501538e-8)
199
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (13) ((1 24))) 4.6072576360465606e-9)
200
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (14) ((1 24))) -1.2324905529154999e-9)
201
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (15) ((1 24))) 3.30620409488102e-10)
202
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (16) ((1 24))) -8.892520997724009e-11)
203
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (17) ((1 24))) 2.3977331987829795e-11)
204
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (18) ((1 24))) -6.480139211534499e-12)
205
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (19) ((1 24))) 1.75510132023731e-12)
206
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (20) ((1 24))) -4.7630382983363695e-13)
207
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (21) ((1 24))) 1.2949824110081e-13)
208
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (22) ((1 24))) -3.5267962221043e-14)
209
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (23) ((1 24))) 9.62005151585923e-15)
210
 
  (f2cl-lib:fset (f2cl-lib:fref dak2 (24) ((1 24))) -2.62786914342292e-15)
211
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (1) ((1 14))) 0.284675828811349)
212
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (2) ((1 14))) 0.0025307307261908003)
213
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (3) ((1 14))) -4.83481130337976e-5)
214
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (4) ((1 14))) 1.8490728394634298e-6)
215
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (5) ((1 14))) -1.01418491178576e-7)
216
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (6) ((1 14))) 7.059256344571531e-9)
217
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (7) ((1 14))) -5.85325291400382e-10)
218
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (8) ((1 14))) 5.563576888313391e-11)
219
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (9) ((1 14))) -5.908890947794999e-12)
220
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (10) ((1 14))) 6.88574353784436e-13)
221
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (11) ((1 14))) -8.68588256452194e-14)
222
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (12) ((1 14))) 1.17374762617213e-14)
223
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (13) ((1 14))) -1.68523146510923e-15)
224
 
  (f2cl-lib:fset (f2cl-lib:fref dak3 (14) ((1 14))) 2.5537477309705603e-16)
225
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (1) ((1 19))) 0.0653219131311457)
226
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (2) ((1 19))) -0.120262933688823)
227
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (3) ((1 19))) 0.00978010236263823)
228
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (4) ((1 19))) 0.016794842923050498)
229
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (5) ((1 19))) -0.00197146140182132)
230
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (6) ((1 19))) -8.45560295098867e-4)
231
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (7) ((1 19))) 9.428896207019759e-5)
232
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (8) ((1 19))) 2.25827860945475e-5)
233
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (9) ((1 19))) -2.29067870915987e-6)
234
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (10) ((1 19))) -3.7634399113691896e-7)
235
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (11) ((1 19))) 3.4566393355956504e-8)
236
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (12) ((1 19))) 4.2961133200300705e-9)
237
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (13) ((1 19))) -3.58673691214989e-10)
238
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (14) ((1 19))) -3.5724588136189495e-11)
239
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (15) ((1 19))) 2.7269609106633597e-12)
240
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (16) ((1 19))) 2.2612065309577098e-13)
241
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (17) ((1 19))) -1.5876320523830303e-14)
242
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (18) ((1 19))) -1.1260437448512499e-15)
243
 
  (f2cl-lib:fset (f2cl-lib:fref dajp (19) ((1 19))) 7.313275295153671e-17)
244
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (1) ((1 19))) 0.0108594539632967)
245
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (2) ((1 19))) 0.0853313194857091)
246
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (3) ((1 19))) -0.315277068113058)
247
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (4) ((1 19))) -0.0878420725294257)
248
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (5) ((1 19))) 0.0553251906976048)
249
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (6) ((1 19))) 0.009416740605032411)
250
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (7) ((1 19))) -0.00332187026018996)
251
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (8) ((1 19))) -4.1115734315682606e-4)
252
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (9) ((1 19))) 1.0129732689134603e-4)
253
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (10) ((1 19))) 9.876336822083961e-6)
254
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (11) ((1 19))) -1.8731296981239298e-6)
255
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (12) ((1 19))) -1.5079850013146798e-7)
256
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (13) ((1 19))) 2.3268766952539402e-8)
257
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (14) ((1 19))) 1.5959991741922502e-9)
258
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (15) ((1 19))) -2.0766592266838502e-10)
259
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (16) ((1 19))) -1.24103350500302e-11)
260
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (17) ((1 19))) 1.39631765331043e-12)
261
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (18) ((1 19))) 7.3940097115574e-14)
262
 
  (f2cl-lib:fset (f2cl-lib:fref dajn (19) ((1 19))) -7.328874756275e-15)
263
 
  (f2cl-lib:fset (f2cl-lib:fref da (1) ((1 15))) 0.49162732110460106)
264
 
  (f2cl-lib:fset (f2cl-lib:fref da (2) ((1 15))) 0.00311164930427489)
265
 
  (f2cl-lib:fset (f2cl-lib:fref da (3) ((1 15))) 8.23140762854081e-5)
266
 
  (f2cl-lib:fset (f2cl-lib:fref da (4) ((1 15))) -4.6176977617214193e-6)
267
 
  (f2cl-lib:fset (f2cl-lib:fref da (5) ((1 15))) -6.13158880534626e-8)
268
 
  (f2cl-lib:fset (f2cl-lib:fref da (6) ((1 15))) 2.8729580465652e-8)
269
 
  (f2cl-lib:fset (f2cl-lib:fref da (7) ((1 15))) -1.8195971537211703e-9)
270
 
  (f2cl-lib:fset (f2cl-lib:fref da (8) ((1 15))) -1.44752826642035e-10)
271
 
  (f2cl-lib:fset (f2cl-lib:fref da (9) ((1 15))) 4.5372404342042194e-11)
272
 
  (f2cl-lib:fset (f2cl-lib:fref da (10) ((1 15))) -3.99655065847223e-12)
273
 
  (f2cl-lib:fset (f2cl-lib:fref da (11) ((1 15))) -3.24089119830323e-13)
274
 
  (f2cl-lib:fset (f2cl-lib:fref da (12) ((1 15))) 1.6209895256874102e-13)
275
 
  (f2cl-lib:fset (f2cl-lib:fref da (13) ((1 15))) -2.4076524797405702e-14)
276
 
  (f2cl-lib:fset (f2cl-lib:fref da (14) ((1 15))) 1.6938481128449098e-16)
277
 
  (f2cl-lib:fset (f2cl-lib:fref da (15) ((1 15))) 8.17900786477396e-16)
278
 
  (f2cl-lib:fset (f2cl-lib:fref db (1) ((1 15))) -0.277571356944231)
279
 
  (f2cl-lib:fset (f2cl-lib:fref db (2) ((1 15))) 0.0044421283341992005)
280
 
  (f2cl-lib:fset (f2cl-lib:fref db (3) ((1 15))) -8.42328522190089e-5)
281
 
  (f2cl-lib:fset (f2cl-lib:fref db (4) ((1 15))) -2.5804031841870995e-6)
282
 
  (f2cl-lib:fset (f2cl-lib:fref db (5) ((1 15))) 3.4238972021762104e-7)
283
 
  (f2cl-lib:fset (f2cl-lib:fref db (6) ((1 15))) -6.24286894709776e-9)
284
 
  (f2cl-lib:fset (f2cl-lib:fref db (7) ((1 15))) -2.3637783684457703e-9)
285
 
  (f2cl-lib:fset (f2cl-lib:fref db (8) ((1 15))) 3.16991042656673e-10)
286
 
  (f2cl-lib:fset (f2cl-lib:fref db (9) ((1 15))) -4.4099569165819097e-12)
287
 
  (f2cl-lib:fset (f2cl-lib:fref db (10) ((1 15))) -5.186742210935751e-12)
288
 
  (f2cl-lib:fset (f2cl-lib:fref db (11) ((1 15))) 9.64874015137022e-13)
289
 
  (f2cl-lib:fset (f2cl-lib:fref db (12) ((1 15))) -4.9019057660870996e-14)
290
 
  (f2cl-lib:fset (f2cl-lib:fref db (13) ((1 15))) -1.77253430678112e-14)
291
 
  (f2cl-lib:fset (f2cl-lib:fref db (14) ((1 15))) 5.55950610442662e-15)
292
 
  (f2cl-lib:fset (f2cl-lib:fref db (15) ((1 15))) -7.117933375795299e-16)
 
255
                                   n4 n3 n2 n1))
293
256
  (defun djairy (x rx c ai dai)
294
257
    (declare (type double-float dai ai c rx x))
295
258
    (prog ((ccv 0.0) (cv 0.0) (ec 0.0) (e1 0.0) (e2 0.0) (f1 0.0) (f2 0.0)
296
 
           (rtrx 0.0) (scv 0.0) (t_ 0.0) (temp1 0.0) (temp2 0.0) (tt 0.0) (i 0)
 
259
           (rtrx 0.0) (scv 0.0) (t$ 0.0) (temp1 0.0) (temp2 0.0) (tt 0.0) (i 0)
297
260
           (j 0))
298
261
      (declare (type f2cl-lib:integer4 j i)
299
 
               (type double-float tt temp2 temp1 t_ scv rtrx f2 f1 e2 e1 ec cv
300
 
                ccv))
 
262
               (type double-float tt temp2 temp1 t$ scv rtrx f2 f1 e2 e1 ec cv
 
263
                                  ccv))
301
264
      (if (< x 0.0) (go label90))
302
265
      (if (> c 5.0) (go label60))
303
266
      (if (> x 1.2) (go label30))
304
 
      (setf t_ (* (- (+ x x) 1.2) con4))
305
 
      (setf tt (+ t_ t_))
 
267
      (setf t$ (* (- (+ x x) 1.2) con4))
 
268
      (setf tt (+ t$ t$))
306
269
      (setf j n1)
307
270
      (setf f1 (f2cl-lib:fref ak1 (j) ((1 14))))
308
271
      (setf f2 0.0)
314
277
          (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref ak1 (j) ((1 14)))))
315
278
          (setf f2 temp1)
316
279
         label10))
317
 
      (setf ai (+ (- (* t_ f1) f2) (f2cl-lib:fref ak1 (1) ((1 14)))))
 
280
      (setf ai (+ (- (* t$ f1) f2) (f2cl-lib:fref ak1 (1) ((1 14)))))
318
281
      (setf j n1d)
319
282
      (setf f1 (f2cl-lib:fref dak1 (j) ((1 14))))
320
283
      (setf f2 0.0)
326
289
          (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dak1 (j) ((1 14)))))
327
290
          (setf f2 temp1)
328
291
         label20))
329
 
      (setf dai (- (+ (- (* t_ f1) f2) (f2cl-lib:fref dak1 (1) ((1 14))))))
 
292
      (setf dai (- (+ (- (* t$ f1) f2) (f2cl-lib:fref dak1 (1) ((1 14))))))
330
293
      (go end_label)
331
294
     label30
332
 
      (setf t_ (* (- (+ x x) con2) con3))
333
 
      (setf tt (+ t_ t_))
 
295
      (setf t$ (* (- (+ x x) con2) con3))
 
296
      (setf tt (+ t$ t$))
334
297
      (setf j n2)
335
298
      (setf f1 (f2cl-lib:fref ak2 (j) ((1 23))))
336
299
      (setf f2 0.0)
345
308
      (setf rtrx (f2cl-lib:fsqrt rx))
346
309
      (setf ec (exp (- c)))
347
310
      (setf ai
348
 
              (/ (* ec (+ (- (* t_ f1) f2) (f2cl-lib:fref ak2 (1) ((1 23)))))
 
311
              (/ (* ec (+ (- (* t$ f1) f2) (f2cl-lib:fref ak2 (1) ((1 23)))))
349
312
                 rtrx))
350
313
      (setf j n2d)
351
314
      (setf f1 (f2cl-lib:fref dak2 (j) ((1 24))))
360
323
         label50))
361
324
      (setf dai
362
325
              (* (- ec)
363
 
                 (+ (- (* t_ f1) f2) (f2cl-lib:fref dak2 (1) ((1 24))))
 
326
                 (+ (- (* t$ f1) f2) (f2cl-lib:fref dak2 (1) ((1 24))))
364
327
                 rtrx))
365
328
      (go end_label)
366
329
     label60
367
 
      (setf t_ (- (/ 10.0 c) 1.0))
368
 
      (setf tt (+ t_ t_))
 
330
      (setf t$ (- (/ 10.0 c) 1.0))
 
331
      (setf tt (+ t$ t$))
369
332
      (setf j n1)
370
333
      (setf f1 (f2cl-lib:fref ak3 (j) ((1 14))))
371
334
      (setf f2 0.0)
380
343
      (setf rtrx (f2cl-lib:fsqrt rx))
381
344
      (setf ec (exp (- c)))
382
345
      (setf ai
383
 
              (/ (* ec (+ (- (* t_ f1) f2) (f2cl-lib:fref ak3 (1) ((1 14)))))
 
346
              (/ (* ec (+ (- (* t$ f1) f2) (f2cl-lib:fref ak3 (1) ((1 14)))))
384
347
                 rtrx))
385
348
      (setf j n1d)
386
349
      (setf f1 (f2cl-lib:fref dak3 (j) ((1 14))))
396
359
      (setf dai
397
360
              (* (- rtrx)
398
361
                 ec
399
 
                 (+ (- (* t_ f1) f2) (f2cl-lib:fref dak3 (1) ((1 14))))))
 
362
                 (+ (- (* t$ f1) f2) (f2cl-lib:fref dak3 (1) ((1 14))))))
400
363
      (go end_label)
401
364
     label90
402
365
      (if (> c 5.0) (go label120))
403
 
      (setf t_ (- (* 0.4 c) 1.0))
404
 
      (setf tt (+ t_ t_))
 
366
      (setf t$ (- (* 0.4 c) 1.0))
 
367
      (setf tt (+ t$ t$))
405
368
      (setf j n3)
406
369
      (setf f1 (f2cl-lib:fref ajp (j) ((1 19))))
407
370
      (setf e1 (f2cl-lib:fref ajn (j) ((1 19))))
419
382
          (setf e2 temp2)
420
383
         label100))
421
384
      (setf ai
422
 
              (- (+ (- (* t_ e1) e2) (f2cl-lib:fref ajn (1) ((1 19))))
423
 
                 (* x (+ (- (* t_ f1) f2) (f2cl-lib:fref ajp (1) ((1 19)))))))
 
385
              (- (+ (- (* t$ e1) e2) (f2cl-lib:fref ajn (1) ((1 19))))
 
386
                 (* x (+ (- (* t$ f1) f2) (f2cl-lib:fref ajp (1) ((1 19)))))))
424
387
      (setf j n3d)
425
388
      (setf f1 (f2cl-lib:fref dajp (j) ((1 19))))
426
389
      (setf e1 (f2cl-lib:fref dajn (j) ((1 19))))
438
401
          (setf e2 temp2)
439
402
         label110))
440
403
      (setf dai
441
 
              (+ (* x x (+ (- (* t_ f1) f2) (f2cl-lib:fref dajp (1) ((1 19)))))
442
 
                 (+ (- (* t_ e1) e2) (f2cl-lib:fref dajn (1) ((1 19))))))
 
404
              (+ (* x x (+ (- (* t$ f1) f2) (f2cl-lib:fref dajp (1) ((1 19)))))
 
405
                 (+ (- (* t$ e1) e2) (f2cl-lib:fref dajn (1) ((1 19))))))
443
406
      (go end_label)
444
407
     label120
445
 
      (setf t_ (- (/ 10.0 c) 1.0))
446
 
      (setf tt (+ t_ t_))
 
408
      (setf t$ (- (/ 10.0 c) 1.0))
 
409
      (setf tt (+ t$ t$))
447
410
      (setf j n4)
448
411
      (setf f1 (f2cl-lib:fref a (j) ((1 15))))
449
412
      (setf e1 (f2cl-lib:fref b (j) ((1 15))))
460
423
          (setf f2 temp1)
461
424
          (setf e2 temp2)
462
425
         label130))
463
 
      (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref a (1) ((1 15)))))
464
 
      (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref b (1) ((1 15)))))
 
426
      (setf temp1 (+ (- (* t$ f1) f2) (f2cl-lib:fref a (1) ((1 15)))))
 
427
      (setf temp2 (+ (- (* t$ e1) e2) (f2cl-lib:fref b (1) ((1 15)))))
465
428
      (setf rtrx (f2cl-lib:fsqrt rx))
466
429
      (setf cv (- c fpi12))
467
430
      (setf ccv (cos cv))
483
446
          (setf f2 temp1)
484
447
          (setf e2 temp2)
485
448
         label140))
486
 
      (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref da (1) ((1 15)))))
487
 
      (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref db (1) ((1 15)))))
 
449
      (setf temp1 (+ (- (* t$ f1) f2) (f2cl-lib:fref da (1) ((1 15)))))
 
450
      (setf temp2 (+ (- (* t$ e1) e2) (f2cl-lib:fref db (1) ((1 15)))))
488
451
      (setf e1 (+ (* ccv con5) (* 0.5 scv)))
489
452
      (setf e2 (- (* scv con5) (* 0.5 ccv)))
490
453
      (setf dai (* (- (* temp1 e1) (* temp2 e2)) rtrx))