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