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, dsqrt, 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