Blame view

fvn_fnlib/d9admp.f 28.8 KB
38581db0c   daniau   git-svn-id: https...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
        subroutine d9admp (x, ampl, phi)
  c july 1980 edition.  w. fullerton, bell labs. revised nov 1983.
  c
  c evaluate the derivative of the airy function modulus and phase
  c for x .le. -1.0 .
  c
        double precision x, ampl, phi, an22cs(74), an21cs(60), an20cs(57),
       1  aph2cs(72), aph1cs(58), aph0cs(53), xsml, z, sqrtx, pi34,
       2  d1mach, dcsevl, dsqrt
        external d1mach, dcsevl, initds
  c
  c series for an22 on the interval -1.00000e+00 to -1.25000e-01
  c                                        with weighted error   4.05e-32
  c                                         log weighted error  31.39
  c                               significant figures required  29.87
  c                                    decimal places required  32.33
  c
        data an22cs(  1) /  0.0537418629 6297943290 9110336091 7783d0/
        data an22cs(  2) / -0.0126661435 8598831934 6631208503 6450d0/
        data an22cs(  3) / -0.0011924334 1065930068 4084891691 3681d0/
        data an22cs(  4) / -0.0002032327 6272756545 5268715517 6363d0/
        data an22cs(  5) / -0.0000446468 9630751639 7951616490 5945d0/
        data an22cs(  6) / -0.0000113359 0360531234 9041699789 3086d0/
        data an22cs(  7) / -0.0000031641 3523785461 0735667135 5827d0/
        data an22cs(  8) / -0.0000009446 7088861489 3912088853 2442d0/
        data an22cs(  9) / -0.0000002966 5622364717 6552790090 5456d0/
        data an22cs( 10) / -0.0000000969 1188920243 6779990866 1433d0/
        data an22cs( 11) / -0.0000000326 8225386532 7409153307 2559d0/
        data an22cs( 12) / -0.0000000113 1446189635 8386590044 7294d0/
        data an22cs( 13) / -0.0000000040 0426910017 4150173827 8050d0/
        data an22cs( 14) / -0.0000000014 4403336839 0742377852 2199d0/
        data an22cs( 15) / -0.0000000005 2928537461 5261158566 3541d0/
        data an22cs( 16) / -0.0000000001 9677633737 0788952824 5726d0/
        data an22cs( 17) / -0.0000000000 7408000957 5584985881 6731d0/
        data an22cs( 18) / -0.0000000000 2820163142 9466198284 2740d0/
        data an22cs( 19) / -0.0000000000 1084400664 6312833133 7590d0/
        data an22cs( 20) / -0.0000000000 0420748006 8264423692 0617d0/
        data an22cs( 21) / -0.0000000000 0164591496 7063481972 4739d0/
        data an22cs( 22) / -0.0000000000 0064868267 0512101889 6077d0/
        data an22cs( 23) / -0.0000000000 0025740950 0335410583 2300d0/
        data an22cs( 24) / -0.0000000000 0010278890 2940782213 2143d0/
        data an22cs( 25) / -0.0000000000 0004128458 2719522272 0128d0/
        data an22cs( 26) / -0.0000000000 0001667110 2933286250 9726d0/
        data an22cs( 27) / -0.0000000000 0000676566 9616560802 3403d0/
        data an22cs( 28) / -0.0000000000 0000275854 4823269357 6823d0/
        data an22cs( 29) / -0.0000000000 0000112963 9791529716 8938d0/
        data an22cs( 30) / -0.0000000000 0000046448 4822545731 4333d0/
        data an22cs( 31) / -0.0000000000 0000019171 9803503391 2928d0/
        data an22cs( 32) / -0.0000000000 0000007941 9757011189 3530d0/
        data an22cs( 33) / -0.0000000000 0000003301 1649230036 8930d0/
        data an22cs( 34) / -0.0000000000 0000001376 5805772654 9714d0/
        data an22cs( 35) / -0.0000000000 0000000575 7809372001 2791d0/
        data an22cs( 36) / -0.0000000000 0000000241 5270085863 2017d0/
        data an22cs( 37) / -0.0000000000 0000000101 5930170093 3666d0/
        data an22cs( 38) / -0.0000000000 0000000042 8443495533 0055d0/
        data an22cs( 39) / -0.0000000000 0000000018 1134405216 8016d0/
        data an22cs( 40) / -0.0000000000 0000000007 6760204561 9422d0/
        data an22cs( 41) / -0.0000000000 0000000003 2602634675 8614d0/
        data an22cs( 42) / -0.0000000000 0000000001 3877380668 2627d0/
        data an22cs( 43) / -0.0000000000 0000000000 5919162710 3729d0/
        data an22cs( 44) / -0.0000000000 0000000000 2529725643 1944d0/
        data an22cs( 45) / -0.0000000000 0000000000 1083207729 3819d0/
        data an22cs( 46) / -0.0000000000 0000000000 0464667488 0404d0/
        data an22cs( 47) / -0.0000000000 0000000000 0199679778 3865d0/
        data an22cs( 48) / -0.0000000000 0000000000 0085952410 8705d0/
        data an22cs( 49) / -0.0000000000 0000000000 0037058415 2073d0/
        data an22cs( 50) / -0.0000000000 0000000000 0016002750 3479d0/
        data an22cs( 51) / -0.0000000000 0000000000 0006920812 4999d0/
        data an22cs( 52) / -0.0000000000 0000000000 0002997444 8994d0/
        data an22cs( 53) / -0.0000000000 0000000000 0001300035 6362d0/
        data an22cs( 54) / -0.0000000000 0000000000 0000564610 0942d0/
        data an22cs( 55) / -0.0000000000 0000000000 0000245534 1103d0/
        data an22cs( 56) / -0.0000000000 0000000000 0000106911 9686d0/
        data an22cs( 57) / -0.0000000000 0000000000 0000046609 5090d0/
        data an22cs( 58) / -0.0000000000 0000000000 0000020344 1579d0/
        data an22cs( 59) / -0.0000000000 0000000000 0000008890 0866d0/
        data an22cs( 60) / -0.0000000000 0000000000 0000003889 1813d0/
        data an22cs( 61) / -0.0000000000 0000000000 0000001703 2637d0/
        data an22cs( 62) / -0.0000000000 0000000000 0000000746 7295d0/
        data an22cs( 63) / -0.0000000000 0000000000 0000000327 7097d0/
        data an22cs( 64) / -0.0000000000 0000000000 0000000143 9618d0/
        data an22cs( 65) / -0.0000000000 0000000000 0000000063 3031d0/
        data an22cs( 66) / -0.0000000000 0000000000 0000000027 8620d0/
        data an22cs( 67) / -0.0000000000 0000000000 0000000012 2743d0/
        data an22cs( 68) / -0.0000000000 0000000000 0000000005 4121d0/
        data an22cs( 69) / -0.0000000000 0000000000 0000000002 3884d0/
        data an22cs( 70) / -0.0000000000 0000000000 0000000001 0549d0/
        data an22cs( 71) / -0.0000000000 0000000000 0000000000 4663d0/
        data an22cs( 72) / -0.0000000000 0000000000 0000000000 2063d0/
        data an22cs( 73) / -0.0000000000 0000000000 0000000000 0913d0/
        data an22cs( 74) / -0.0000000000 0000000000 0000000000 0405d0/
  c
  c series for an21 on the interval -1.25000e-01 to -1.56250e-02
  c                                        with weighted error   3.55e-32
  c                                         log weighted error  31.45
  c                               significant figures required  29.46
  c                                    decimal places required  32.34
  c
        data an21cs(  1) /  0.0198313155 2631693944 2034248316 5643d0/
        data an21cs(  2) / -0.0029376249 0670875334 6059374559 4484d0/
        data an21cs(  3) / -0.0001136260 6959581955 4987261113 7182d0/
        data an21cs(  4) / -0.0000100554 4510871560 0975098164 5918d0/
        data an21cs(  5) / -0.0000013048 7871165632 5042178559 8252d0/
        data an21cs(  6) / -0.0000002123 8819931506 6483066607 9609d0/
        data an21cs(  7) / -0.0000000402 2708333842 6904034785 0109d0/
        data an21cs(  8) / -0.0000000084 9967459531 6179914220 1792d0/
        data an21cs(  9) / -0.0000000019 5148394261 7861409953 2934d0/
        data an21cs( 10) / -0.0000000004 7838653438 4038428299 2480d0/
        data an21cs( 11) / -0.0000000001 2367339920 9945050113 7105d0/
        data an21cs( 12) / -0.0000000000 3341374863 9875423221 9789d0/
        data an21cs( 13) / -0.0000000000 0937028235 4076632989 7780d0/
        data an21cs( 14) / -0.0000000000 0271301281 5613956468 7240d0/
        data an21cs( 15) / -0.0000000000 0080759538 0058347953 5949d0/
        data an21cs( 16) / -0.0000000000 0024632143 0470012525 2160d0/
        data an21cs( 17) / -0.0000000000 0007676556 8910932156 4410d0/
        data an21cs( 18) / -0.0000000000 0002438825 9880735491 9791d0/
        data an21cs( 19) / -0.0000000000 0000788314 6635876030 8462d0/
        data an21cs( 20) / -0.0000000000 0000258824 0099558586 4077d0/
        data an21cs( 21) / -0.0000000000 0000086194 5786294569 0828d0/
        data an21cs( 22) / -0.0000000000 0000029079 9473966312 8534d0/
        data an21cs( 23) / -0.0000000000 0000009928 4679612289 0484d0/
        data an21cs( 24) / -0.0000000000 0000003427 2022918777 4480d0/
        data an21cs( 25) / -0.0000000000 0000001195 1104820551 5026d0/
        data an21cs( 26) / -0.0000000000 0000000420 6972904367 8359d0/
        data an21cs( 27) / -0.0000000000 0000000149 3969776281 8400d0/
        data an21cs( 28) / -0.0000000000 0000000053 4898116158 9517d0/
        data an21cs( 29) / -0.0000000000 0000000019 2987757782 6238d0/
        data an21cs( 30) / -0.0000000000 0000000007 0131370101 8203d0/
        data an21cs( 31) / -0.0000000000 0000000002 5658573850 9682d0/
        data an21cs( 32) / -0.0000000000 0000000000 9447589456 2734d0/
        data an21cs( 33) / -0.0000000000 0000000000 3499640194 1465d0/
        data an21cs( 34) / -0.0000000000 0000000000 1303762246 6397d0/
        data an21cs( 35) / -0.0000000000 0000000000 0488333416 3346d0/
        data an21cs( 36) / -0.0000000000 0000000000 0183847758 6152d0/
        data an21cs( 37) / -0.0000000000 0000000000 0069552732 4058d0/
        data an21cs( 38) / -0.0000000000 0000000000 0026435191 0209d0/
        data an21cs( 39) / -0.0000000000 0000000000 0010091809 4655d0/
        data an21cs( 40) / -0.0000000000 0000000000 0003868892 4289d0/
        data an21cs( 41) / -0.0000000000 0000000000 0001489203 6525d0/
        data an21cs( 42) / -0.0000000000 0000000000 0000575434 2426d0/
        data an21cs( 43) / -0.0000000000 0000000000 0000223172 5971d0/
        data an21cs( 44) / -0.0000000000 0000000000 0000086860 7480d0/
        data an21cs( 45) / -0.0000000000 0000000000 0000033922 0403d0/
        data an21cs( 46) / -0.0000000000 0000000000 0000013291 0128d0/
        data an21cs( 47) / -0.0000000000 0000000000 0000005223 9309d0/
        data an21cs( 48) / -0.0000000000 0000000000 0000002059 4383d0/
        data an21cs( 49) / -0.0000000000 0000000000 0000000814 2614d0/
        data an21cs( 50) / -0.0000000000 0000000000 0000000322 8473d0/
        data an21cs( 51) / -0.0000000000 0000000000 0000000128 3529d0/
        data an21cs( 52) / -0.0000000000 0000000000 0000000051 1622d0/
        data an21cs( 53) / -0.0000000000 0000000000 0000000020 4451d0/
        data an21cs( 54) / -0.0000000000 0000000000 0000000008 1901d0/
        data an21cs( 55) / -0.0000000000 0000000000 0000000003 2886d0/
        data an21cs( 56) / -0.0000000000 0000000000 0000000001 3235d0/
        data an21cs( 57) / -0.0000000000 0000000000 0000000000 5338d0/
        data an21cs( 58) / -0.0000000000 0000000000 0000000000 2158d0/
        data an21cs( 59) / -0.0000000000 0000000000 0000000000 0874d0/
        data an21cs( 60) / -0.0000000000 0000000000 0000000000 0355d0/
  c
  c series for an20 on the interval -1.56250e-02 to  0.00000e+00
  c                                        with weighted error   3.18e-32
  c                                         log weighted error  31.50
  c                               significant figures required  29.30
  c                                    decimal places required  32.38
  c
        data an20cs(  1) /  0.0126732217 1457380271 5461075103 4240d0/
        data an20cs(  2) / -0.0005212847 0726156211 8478094230 9478d0/
        data an20cs(  3) / -0.0000052672 1111403704 2980907405 2969d0/
        data an20cs(  4) / -0.0000001628 2021850264 8375263246 0680d0/
        data an20cs(  5) / -0.0000000090 9914426873 7138632597 3075d0/
        data an20cs(  6) / -0.0000000007 4386471262 4219289068 5403d0/
        data an20cs(  7) / -0.0000000000 7954947515 9146948612 2822d0/
        data an20cs(  8) / -0.0000000000 1040509442 8830374280 3960d0/
        data an20cs(  9) / -0.0000000000 0159324255 9841455152 3990d0/
        data an20cs( 10) / -0.0000000000 0027706482 7234191394 6674d0/
        data an20cs( 11) / -0.0000000000 0005353426 2923760629 5104d0/
        data an20cs( 12) / -0.0000000000 0001130615 4178172831 4051d0/
        data an20cs( 13) / -0.0000000000 0000257721 9007894316 7788d0/
        data an20cs( 14) / -0.0000000000 0000062780 3311603248 5076d0/
        data an20cs( 15) / -0.0000000000 0000016212 9540018993 9757d0/
        data an20cs( 16) / -0.0000000000 0000004409 9298524067 5353d0/
        data an20cs( 17) / -0.0000000000 0000001256 5551655325 8972d0/
        data an20cs( 18) / -0.0000000000 0000000373 3690698801 5204d0/
        data an20cs( 19) / -0.0000000000 0000000115 2462692672 4671d0/
        data an20cs( 20) / -0.0000000000 0000000036 8308149909 9144d0/
        data an20cs( 21) / -0.0000000000 0000000012 1520696533 1797d0/
        data an20cs( 22) / -0.0000000000 0000000004 1291617772 4016d0/
        data an20cs( 23) / -0.0000000000 0000000001 4417736423 9347d0/
        data an20cs( 24) / -0.0000000000 0000000000 5163184287 5864d0/
        data an20cs( 25) / -0.0000000000 0000000000 1893124266 8250d0/
        data an20cs( 26) / -0.0000000000 0000000000 0709605466 8569d0/
        data an20cs( 27) / -0.0000000000 0000000000 0271540664 6904d0/
        data an20cs( 28) / -0.0000000000 0000000000 0105948697 9400d0/
        data an20cs( 29) / -0.0000000000 0000000000 0042103003 5685d0/
        data an20cs( 30) / -0.0000000000 0000000000 0017023378 1664d0/
        data an20cs( 31) / -0.0000000000 0000000000 0006996667 7028d0/
        data an20cs( 32) / -0.0000000000 0000000000 0002920664 3813d0/
        data an20cs( 33) / -0.0000000000 0000000000 0001237312 8203d0/
        data an20cs( 34) / -0.0000000000 0000000000 0000531587 1095d0/
        data an20cs( 35) / -0.0000000000 0000000000 0000231462 2618d0/
        data an20cs( 36) / -0.0000000000 0000000000 0000102077 9922d0/
        data an20cs( 37) / -0.0000000000 0000000000 0000045570 6227d0/
        data an20cs( 38) / -0.0000000000 0000000000 0000020583 1071d0/
        data an20cs( 39) / -0.0000000000 0000000000 0000009401 5189d0/
        data an20cs( 40) / -0.0000000000 0000000000 0000004340 5874d0/
        data an20cs( 41) / -0.0000000000 0000000000 0000002024 7792d0/
        data an20cs( 42) / -0.0000000000 0000000000 0000000953 9214d0/
        data an20cs( 43) / -0.0000000000 0000000000 0000000453 7234d0/
        data an20cs( 44) / -0.0000000000 0000000000 0000000217 8016d0/
        data an20cs( 45) / -0.0000000000 0000000000 0000000105 4823d0/
        data an20cs( 46) / -0.0000000000 0000000000 0000000051 5242d0/
        data an20cs( 47) / -0.0000000000 0000000000 0000000025 3763d0/
        data an20cs( 48) / -0.0000000000 0000000000 0000000012 5983d0/
        data an20cs( 49) / -0.0000000000 0000000000 0000000006 3030d0/
        data an20cs( 50) / -0.0000000000 0000000000 0000000003 1771d0/
        data an20cs( 51) / -0.0000000000 0000000000 0000000001 6131d0/
        data an20cs( 52) / -0.0000000000 0000000000 0000000000 8248d0/
        data an20cs( 53) / -0.0000000000 0000000000 0000000000 4246d0/
        data an20cs( 54) / -0.0000000000 0000000000 0000000000 2200d0/
        data an20cs( 55) / -0.0000000000 0000000000 0000000000 1147d0/
        data an20cs( 56) / -0.0000000000 0000000000 0000000000 0602d0/
        data an20cs( 57) / -0.0000000000 0000000000 0000000000 0318d0/
  c
  c series for aph2 on the interval -1.00000e+00 to -1.25000e-01
  c                                        with weighted error   5.44e-32
  c                                         log weighted error  31.26
  c                               significant figures required  30.31
  c                                    decimal places required  32.19
  c
        data aph2cs(  1) / -0.2057088719 7814651069 7364866560 2125d0/
        data aph2cs(  2) /  0.0422196961 3577719216 7311498036 9460d0/
        data aph2cs(  3) /  0.0020482560 5112072750 4266057781 3334d0/
        data aph2cs(  4) /  0.0002607800 7351650056 3118787992 2652d0/
        data aph2cs(  5) /  0.0000474824 2680047288 7538175051 9293d0/
        data aph2cs(  6) /  0.0000105102 7564316117 4347363002 6955d0/
        data aph2cs(  7) /  0.0000026353 5340146679 4510931404 1983d0/
        data aph2cs(  8) /  0.0000007208 8248634991 4729979078 3731d0/
        data aph2cs(  9) /  0.0000002103 2366644733 5285974947 7082d0/
        data aph2cs( 10) /  0.0000000644 9756345552 9559843736 2273d0/
        data aph2cs( 11) /  0.0000000205 8023772643 6850797811 6888d0/
        data aph2cs( 12) /  0.0000000067 8362739209 0642896351 3918d0/
        data aph2cs( 13) /  0.0000000022 9740152840 0940016834 3792d0/
        data aph2cs( 14) /  0.0000000007 9613067654 9118753488 3226d0/
        data aph2cs( 15) /  0.0000000002 8138606097 4159171900 3632d0/
        data aph2cs( 16) /  0.0000000001 0117490569 3197392284 1793d0/
        data aph2cs( 17) /  0.0000000000 3693067379 5247655909 7060d0/
        data aph2cs( 18) /  0.0000000000 1366150661 2709803177 8842d0/
        data aph2cs( 19) /  0.0000000000 0511427514 1604504511 9388d0/
        data aph2cs( 20) /  0.0000000000 0193516889 3170651624 7975d0/
        data aph2cs( 21) /  0.0000000000 0073936069 1649322421 7271d0/
        data aph2cs( 22) /  0.0000000000 0028497922 1922274359 7555d0/
        data aph2cs( 23) /  0.0000000000 0011072807 8245964833 5733d0/
        data aph2cs( 24) /  0.0000000000 0004334121 9937013463 3169d0/
        data aph2cs( 25) /  0.0000000000 0001708008 2526567036 7471d0/
        data aph2cs( 26) /  0.0000000000 0000677330 8019563111 4673d0/
        data aph2cs( 27) /  0.0000000000 0000270169 0478926241 4108d0/
        data aph2cs( 28) /  0.0000000000 0000108347 2075181078 2141d0/
        data aph2cs( 29) /  0.0000000000 0000043670 6031297028 6167d0/
        data aph2cs( 30) /  0.0000000000 0000017685 1173805336 6608d0/
        data aph2cs( 31) /  0.0000000000 0000007193 5921309364 5717d0/
        data aph2cs( 32) /  0.0000000000 0000002938 2361000293 3154d0/
        data aph2cs( 33) /  0.0000000000 0000001204 8281152584 8357d0/
        data aph2cs( 34) /  0.0000000000 0000000495 8665949109 1389d0/
        data aph2cs( 35) /  0.0000000000 0000000204 7943831584 7217d0/
        data aph2cs( 36) /  0.0000000000 0000000084 8601994441 0629d0/
        data aph2cs( 37) /  0.0000000000 0000000035 2735176538 4506d0/
        data aph2cs( 38) /  0.0000000000 0000000014 7056399680 4903d0/
        data aph2cs( 39) /  0.0000000000 0000000006 1481782690 2188d0/
        data aph2cs( 40) /  0.0000000000 0000000002 5773770656 5077d0/
        data aph2cs( 41) /  0.0000000000 0000000001 0832390359 0042d0/
        data aph2cs( 42) /  0.0000000000 0000000000 4563889802 4998d0/
        data aph2cs( 43) /  0.0000000000 0000000000 1927363540 3662d0/
        data aph2cs( 44) /  0.0000000000 0000000000 0815766856 9775d0/
        data aph2cs( 45) /  0.0000000000 0000000000 0346020282 8346d0/
        data aph2cs( 46) /  0.0000000000 0000000000 0147072648 2427d0/
        data aph2cs( 47) /  0.0000000000 0000000000 0062635607 4088d0/
        data aph2cs( 48) /  0.0000000000 0000000000 0026726129 2780d0/
        data aph2cs( 49) /  0.0000000000 0000000000 0011424694 8763d0/
        data aph2cs( 50) /  0.0000000000 0000000000 0004892346 0516d0/
        data aph2cs( 51) /  0.0000000000 0000000000 0002098580 7810d0/
        data aph2cs( 52) /  0.0000000000 0000000000 0000901661 8807d0/
        data aph2cs( 53) /  0.0000000000 0000000000 0000388012 9464d0/
        data aph2cs( 54) /  0.0000000000 0000000000 0000167228 2170d0/
        data aph2cs( 55) /  0.0000000000 0000000000 0000072179 0800d0/
        data aph2cs( 56) /  0.0000000000 0000000000 0000031198 2573d0/
        data aph2cs( 57) /  0.0000000000 0000000000 0000013503 5015d0/
        data aph2cs( 58) /  0.0000000000 0000000000 0000005852 4861d0/
        data aph2cs( 59) /  0.0000000000 0000000000 0000002539 7686d0/
        data aph2cs( 60) /  0.0000000000 0000000000 0000001103 5457d0/
        data aph2cs( 61) /  0.0000000000 0000000000 0000000480 0788d0/
        data aph2cs( 62) /  0.0000000000 0000000000 0000000209 0956d0/
        data aph2cs( 63) /  0.0000000000 0000000000 0000000091 1743d0/
        data aph2cs( 64) /  0.0000000000 0000000000 0000000039 7998d0/
        data aph2cs( 65) /  0.0000000000 0000000000 0000000017 3923d0/
        data aph2cs( 66) /  0.0000000000 0000000000 0000000007 6083d0/
        data aph2cs( 67) /  0.0000000000 0000000000 0000000003 3316d0/
        data aph2cs( 68) /  0.0000000000 0000000000 0000000001 4604d0/
        data aph2cs( 69) /  0.0000000000 0000000000 0000000000 6407d0/
        data aph2cs( 70) /  0.0000000000 0000000000 0000000000 2814d0/
        data aph2cs( 71) /  0.0000000000 0000000000 0000000000 1237d0/
        data aph2cs( 72) /  0.0000000000 0000000000 0000000000 0544d0/
  c
  c series for aph1 on the interval -1.25000e-01 to -1.56250e-02
  c                                        with weighted error   3.08e-32
  c                                         log weighted error  31.51
  c                               significant figures required  30.23
  c                                    decimal places required  32.39
  c
        data aph1cs(  1) / -0.1024172908 0775716940 2112332181 3917d0/
        data aph1cs(  2) /  0.0071697275 1465912480 4721164914 4704d0/
        data aph1cs(  3) /  0.0001209959 3631223285 8981385649 1397d0/
        data aph1cs(  4) /  0.0000073361 5128412199 1208029784 5684d0/
        data aph1cs(  5) /  0.0000007535 3829542716 0706998290 3869d0/
        data aph1cs(  6) /  0.0000001041 4781717413 0192688510 9155d0/
        data aph1cs(  7) /  0.0000000174 3587285185 4569185890 7606d0/
        data aph1cs(  8) /  0.0000000033 3997950333 4645166018 4961d0/
        data aph1cs(  9) /  0.0000000007 0730751743 6352708339 9508d0/
        data aph1cs( 10) /  0.0000000001 6191875151 8977326679 2272d0/
        data aph1cs( 11) /  0.0000000000 3945399818 8195488987 9668d0/
        data aph1cs( 12) /  0.0000000000 1011922817 3422713329 2631d0/
        data aph1cs( 13) /  0.0000000000 0270927782 5952033219 8030d0/
        data aph1cs( 14) /  0.0000000000 0075238064 1842254888 5854d0/
        data aph1cs( 15) /  0.0000000000 0021563687 3300896635 7328d0/
        data aph1cs( 16) /  0.0000000000 0006352827 7712606841 0174d0/
        data aph1cs( 17) /  0.0000000000 0001917569 7264150172 9345d0/
        data aph1cs( 18) /  0.0000000000 0000591430 7244646489 1558d0/
        data aph1cs( 19) /  0.0000000000 0000185971 2851727502 8357d0/
        data aph1cs( 20) /  0.0000000000 0000059504 4492394610 3668d0/
        data aph1cs( 21) /  0.0000000000 0000019342 2995643018 0252d0/
        data aph1cs( 22) /  0.0000000000 0000006378 4302148950 4324d0/
        data aph1cs( 23) /  0.0000000000 0000002131 2729008731 2393d0/
        data aph1cs( 24) /  0.0000000000 0000000720 8138065672 8500d0/
        data aph1cs( 25) /  0.0000000000 0000000246 5249414476 9247d0/
        data aph1cs( 26) /  0.0000000000 0000000085 1911057026 6154d0/
        data aph1cs( 27) /  0.0000000000 0000000029 7238446849 1170d0/
        data aph1cs( 28) /  0.0000000000 0000000010 4642664881 1446d0/
        data aph1cs( 29) /  0.0000000000 0000000003 7149303634 7327d0/
        data aph1cs( 30) /  0.0000000000 0000000001 3292324779 3472d0/
        data aph1cs( 31) /  0.0000000000 0000000000 4791283792 5909d0/
        data aph1cs( 32) /  0.0000000000 0000000000 1739061985 9336d0/
        data aph1cs( 33) /  0.0000000000 0000000000 0635358517 3501d0/
        data aph1cs( 34) /  0.0000000000 0000000000 0233564361 4263d0/
        data aph1cs( 35) /  0.0000000000 0000000000 0086364388 1606d0/
        data aph1cs( 36) /  0.0000000000 0000000000 0032112300 6944d0/
        data aph1cs( 37) /  0.0000000000 0000000000 0012003154 0983d0/
        data aph1cs( 38) /  0.0000000000 0000000000 0004509148 8699d0/
        data aph1cs( 39) /  0.0000000000 0000000000 0001702022 8580d0/
        data aph1cs( 40) /  0.0000000000 0000000000 0000645374 4630d0/
        data aph1cs( 41) /  0.0000000000 0000000000 0000245778 8564d0/
        data aph1cs( 42) /  0.0000000000 0000000000 0000093989 7684d0/
        data aph1cs( 43) /  0.0000000000 0000000000 0000036086 3150d0/
        data aph1cs( 44) /  0.0000000000 0000000000 0000013907 7884d0/
        data aph1cs( 45) /  0.0000000000 0000000000 0000005379 7184d0/
        data aph1cs( 46) /  0.0000000000 0000000000 0000002088 2551d0/
        data aph1cs( 47) /  0.0000000000 0000000000 0000000813 3371d0/
        data aph1cs( 48) /  0.0000000000 0000000000 0000000317 8080d0/
        data aph1cs( 49) /  0.0000000000 0000000000 0000000124 5700d0/
        data aph1cs( 50) /  0.0000000000 0000000000 0000000048 9742d0/
        data aph1cs( 51) /  0.0000000000 0000000000 0000000019 3099d0/
        data aph1cs( 52) /  0.0000000000 0000000000 0000000007 6349d0/
        data aph1cs( 53) /  0.0000000000 0000000000 0000000003 0269d0/
        data aph1cs( 54) /  0.0000000000 0000000000 0000000001 2032d0/
        data aph1cs( 55) /  0.0000000000 0000000000 0000000000 4795d0/
        data aph1cs( 56) /  0.0000000000 0000000000 0000000000 1915d0/
        data aph1cs( 57) /  0.0000000000 0000000000 0000000000 0767d0/
        data aph1cs( 58) /  0.0000000000 0000000000 0000000000 0308d0/
  c
  c series for aph0 on the interval -1.56250e-02 to  0.00000e+00
  c                                        with weighted error   2.81e-32
  c                                         log weighted error  31.55
  c                               significant figures required  30.18
  c                                    decimal places required  32.41
  c
        data aph0cs(  1) / -0.0855849241 1309332569 2012426017 9491d0/
        data aph0cs(  2) /  0.0011214378 8670652607 3578672247 1124d0/
        data aph0cs(  3) /  0.0000042721 0293536641 1395157374 2015d0/
        data aph0cs(  4) /  0.0000000817 6073814832 4364401806 2323d0/
        data aph0cs(  5) /  0.0000000033 9076450004 9272420781 6418d0/
        data aph0cs(  6) /  0.0000000002 2532644226 1911393984 5276d0/
        data aph0cs(  7) /  0.0000000000 2062842092 2901525125 6990d0/
        data aph0cs(  8) /  0.0000000000 0238587628 2813088762 7258d0/
        data aph0cs(  9) /  0.0000000000 0033016181 0588670548 0628d0/
        data aph0cs( 10) /  0.0000000000 0005270096 4850832858 1123d0/
        data aph0cs( 11) /  0.0000000000 0000945554 8220381349 2868d0/
        data aph0cs( 12) /  0.0000000000 0000187094 2695134483 6908d0/
        data aph0cs( 13) /  0.0000000000 0000040239 8004182539 2741d0/
        data aph0cs( 14) /  0.0000000000 0000009301 9287925898 3167d0/
        data aph0cs( 15) /  0.0000000000 0000002290 3863540237 9945d0/
        data aph0cs( 16) /  0.0000000000 0000000596 3435982208 3386d0/
        data aph0cs( 17) /  0.0000000000 0000000163 2027965940 3399d0/
        data aph0cs( 18) /  0.0000000000 0000000046 7114565886 1339d0/
        data aph0cs( 19) /  0.0000000000 0000000013 9233441536 3502d0/
        data aph0cs( 20) /  0.0000000000 0000000004 3064267028 5155d0/
        data aph0cs( 21) /  0.0000000000 0000000001 3778141631 8755d0/
        data aph0cs( 22) /  0.0000000000 0000000000 4547671048 0396d0/
        data aph0cs( 23) /  0.0000000000 0000000000 1544842020 3026d0/
        data aph0cs( 24) /  0.0000000000 0000000000 0538977055 1212d0/
        data aph0cs( 25) /  0.0000000000 0000000000 0192772673 7155d0/
        data aph0cs( 26) /  0.0000000000 0000000000 0070565932 0166d0/
        data aph0cs( 27) /  0.0000000000 0000000000 0026398508 4827d0/
        data aph0cs( 28) /  0.0000000000 0000000000 0010079130 1805d0/
        data aph0cs( 29) /  0.0000000000 0000000000 0003922892 8481d0/
        data aph0cs( 30) /  0.0000000000 0000000000 0001554742 2955d0/
        data aph0cs( 31) /  0.0000000000 0000000000 0000626830 6372d0/
        data aph0cs( 32) /  0.0000000000 0000000000 0000256856 3962d0/
        data aph0cs( 33) /  0.0000000000 0000000000 0000106885 8883d0/
        data aph0cs( 34) /  0.0000000000 0000000000 0000045134 7253d0/
        data aph0cs( 35) /  0.0000000000 0000000000 0000019326 7262d0/
        data aph0cs( 36) /  0.0000000000 0000000000 0000008386 5369d0/
        data aph0cs( 37) /  0.0000000000 0000000000 0000003685 7386d0/
        data aph0cs( 38) /  0.0000000000 0000000000 0000001639 6202d0/
        data aph0cs( 39) /  0.0000000000 0000000000 0000000737 9298d0/
        data aph0cs( 40) /  0.0000000000 0000000000 0000000335 8392d0/
        data aph0cs( 41) /  0.0000000000 0000000000 0000000154 4891d0/
        data aph0cs( 42) /  0.0000000000 0000000000 0000000071 8013d0/
        data aph0cs( 43) /  0.0000000000 0000000000 0000000033 7026d0/
        data aph0cs( 44) /  0.0000000000 0000000000 0000000015 9710d0/
        data aph0cs( 45) /  0.0000000000 0000000000 0000000007 6382d0/
        data aph0cs( 46) /  0.0000000000 0000000000 0000000003 6855d0/
        data aph0cs( 47) /  0.0000000000 0000000000 0000000001 7935d0/
        data aph0cs( 48) /  0.0000000000 0000000000 0000000000 8800d0/
        data aph0cs( 49) /  0.0000000000 0000000000 0000000000 4353d0/
        data aph0cs( 50) /  0.0000000000 0000000000 0000000000 2170d0/
        data aph0cs( 51) /  0.0000000000 0000000000 0000000000 1090d0/
        data aph0cs( 52) /  0.0000000000 0000000000 0000000000 0551d0/
        data aph0cs( 53) /  0.0000000000 0000000000 0000000000 0281d0/
  c
        data pi34 / 2.3561944 9019234492 8846982537 4596271631 3d0  /
        data nan20, nan21, nan22, naph0, naph1, naph2 / 6*0 /
        data xsml / 0.0d0 /
  c
        if (nan20.ne.0) go to 10
        eta = 0.1d0*sngl(d1mach(3))
        nan20 = initds (an20cs, 57, eta)
        nan21 = initds (an21cs, 60, eta)
        nan22 = initds (an22cs, 74, eta)
        naph0 = initds (aph0cs, 53, eta)
        naph1 = initds (aph1cs, 58, eta)
        naph2 = initds (aph2cs, 72, eta)
  c
        xsml = -(128.d0/d1mach(3))**0.3333d0
  c
   10   if (x.ge.(-4.0d0)) go to 20
        z = 1.0d0
        if (x.gt.xsml) z = 128.d0/x**3 + 1.0d0
        ampl = 0.3125d0 + dcsevl (z, an20cs, nan20)
        phi = -0.625d0 + dcsevl (z, aph0cs, naph0)
        go to 40
  c
   20   if (x.ge.(-2.d0)) go to 30
        z = (128.d0/x**3 + 9.0d0) / 7.0d0
        ampl = 0.3125d0 + dcsevl (z, an21cs, nan21)
        phi = -0.625 + dcsevl (z, aph1cs, naph1)
  c
        go to 40
  c
   30   if (x.gt.(-1.0d0)) call seteru (
       1  25hd9admp  x must be le -1.0, 25, 1, 2)
  c
        z = (16.d0/x**3 + 9.0d0) / 7.0d0
        ampl = 0.3125d0 + dcsevl (z, an22cs, nan22)
        phi = -0.625d0 + dcsevl (z, aph2cs, naph2)
  c
   40   sqrtx = dsqrt(-x)
        ampl = dsqrt (ampl*sqrtx)
        phi = pi34 - x*sqrtx*phi
  c
        return
        end