subroutine d9aimp (x, ampl, theta)
c july 1977 edition.  w. fullerton, c3, los alamos scientific lab.
c
c evaluate the airy modulus and phase for x .le. -1.0
c
      double precision x, ampl, theta, am20cs(57), ath0cs(53),
     1  am21cs(60), ath1cs(58), am22cs(74), ath2cs(72), pi4, sqrtx,
     2  xsml, z,  d1mach, dcsevl, dsqrt
      external d1mach, dcsevl, initds
c
c series for am20       on the interval -1.56250e-02 to  0.
c                                        with weighted error   3.12e-32
c                                         log weighted error  31.51
c                               significant figures required  29.24
c                                    decimal places required  32.38
c
      data am20cs(  1) / +.1087167490 8656185661 5730588125 d-1        /
      data am20cs(  2) / +.3694892289 8266355509 1728665146 d-3        /
      data am20cs(  3) / +.4406801004 8468956366 7507001327 d-5        /
      data am20cs(  4) / +.1436867623 6191115392 9183952833 d-6        /
      data am20cs(  5) / +.8242755523 9007830867 0628855353 d-8        /
      data am20cs(  6) / +.6844267588 9366160617 3927278180 d-9        /
      data am20cs(  7) / +.7395666972 8273928773 1004740213 d-10       /
      data am20cs(  8) / +.9745956336 9682501763 8702600847 d-11       /
      data am20cs(  9) / +.1500768858 2940577565 0973119497 d-11       /
      data am20cs( 10) / +.2621479102 2152763420 6252854802 d-12       /
      data am20cs( 11) / +.5083541113 7648718035 7278966914 d-13       /
      data am20cs( 12) / +.1076847533 5881144049 2985997070 d-13       /
      data am20cs( 13) / +.2460912866 1843342933 5914062617 d-14       /
      data am20cs( 14) / +.6007863803 5865641843 6110373550 d-15       /
      data am20cs( 15) / +.1554491561 0238807115 0651388384 d-15       /
      data am20cs( 16) / +.4235351250 3557660442 6382780182 d-16       /
      data am20cs( 17) / +.1208621662 8929984015 4401109189 d-16       /
      data am20cs( 18) / +.3596096512 1465824086 1499706423 d-17       /
      data am20cs( 19) / +.1111342183 8639563826 1774604677 d-17       /
      data am20cs( 20) / +.3555595324 3236660989 3680289225 d-18       /
      data am20cs( 21) / +.1174330216 0013930999 8766947387 d-18       /
      data am20cs( 22) / +.3993974546 6107756138 9162200966 d-19       /
      data am20cs( 23) / +.1395766715 2891631042 5606325640 d-19       /
      data am20cs( 24) / +.5002400553 0923604139 3459280716 d-20       /
      data am20cs( 25) / +.1835527609 5813267918 4834866457 d-20       /
      data am20cs( 26) / +.6884909981 7920274319 7790112404 d-21       /
      data am20cs( 27) / +.2636310356 1141701235 9996885105 d-21       /
      data am20cs( 28) / +.1029248902 3733836028 7153563785 d-21       /
      data am20cs( 29) / +.4092469666 7159488548 9762960571 d-22       /
      data am20cs( 30) / +.1655585734 0673465103 9727903828 d-22       /
      data am20cs( 31) / +.6807974670 6303335611 6599685727 d-23       /
      data am20cs( 32) / +.2843265599 3407983241 9751134476 d-23       /
      data am20cs( 33) / +.1205073983 4896525509 7287818819 d-23       /
      data am20cs( 34) / +.5179612432 8750521797 6613610424 d-24       /
      data am20cs( 35) / +.2256226134 2756281630 3268640887 d-24       /
      data am20cs( 36) / +.9954188011 4774516883 2117078246 d-25       /
      data am20cs( 37) / +.4445516963 9734242430 8280582053 d-25       /
      data am20cs( 38) / +.2008651954 6150110142 5916097338 d-25       /
      data am20cs( 39) / +.9177863441 5177516597 3885645402 d-26       /
      data am20cs( 40) / +.4238729581 0558924066 1672197948 d-26       /
      data am20cs( 41) / +.1977892720 0784609237 0846251490 d-26       /
      data am20cs( 42) / +.9321163512 8462066568 0435253373 d-27       /
      data am20cs( 43) / +.4434821332 4991809995 5611379722 d-27       /
      data am20cs( 44) / +.2129456723 6557389559 4589552837 d-27       /
      data am20cs( 45) / +.1031585696 5107597755 2209344907 d-27       /
      data am20cs( 46) / +.5040237730 2259119915 7904590029 d-28       /
      data am20cs( 47) / +.2483013045 7015594530 4046541005 d-28       /
      data am20cs( 48) / +.1233017831 2856219605 4198238560 d-28       /
      data am20cs( 49) / +.6170334499 2052174612 1976730507 d-29       /
      data am20cs( 50) / +.3110926174 1591889723 3869792213 d-29       /
      data am20cs( 51) / +.1579830852 0170617301 5269071503 d-29       /
      data am20cs( 52) / +.8079319875 3828360767 8121339092 d-30       /
      data am20cs( 53) / +.4159973941 3866756272 2951360052 d-30       /
      data am20cs( 54) / +.2156109340 9771690047 1935862504 d-30       /
      data am20cs( 55) / +.1124688572 6586917829 6752823613 d-30       /
      data am20cs( 56) / +.5903315606 3283809112 3040811797 d-31       /
      data am20cs( 57) / +.3117356676 9292856204 6280505333 d-31       /
c
c series for ath0       on the interval -1.56250e-02 to  0.
c                                        with weighted error   2.75e-32
c                                         log weighted error  31.56
c                               significant figures required  30.17
c                                    decimal places required  32.42
c
      data ath0cs(  1) / -.8172601764 1616344998 4020870054 3 d-1      /
      data ath0cs(  2) / -.8004012824 7882732875 9648111306 8 d-3      /
      data ath0cs(  3) / -.3186525268 7821132037 9555362824 2 d-5      /
      data ath0cs(  4) / -.6688388266 4775093307 4169886503 3 d-7      /
      data ath0cs(  5) / -.2931759284 9945645165 0682246318 4 d-8      /
      data ath0cs(  6) / -.2011263760 8836216690 4903030718 6 d-9      /
      data ath0cs(  7) / -.1877522678 0559734260 7400816665 2 d-10     /
      data ath0cs(  8) / -.2199637137 7046012518 9900219984 8 d-11     /
      data ath0cs(  9) / -.3071616682 5922724490 2574660558 6 d-12     /
      data ath0cs( 10) / -.4936140553 6734183610 2560098538 9 d-13     /
      data ath0cs( 11) / -.8902833722 5836604169 3523696986 6 d-14     /
      data ath0cs( 12) / -.1768987764 6152726136 5681419946 7 d-14     /
      data ath0cs( 13) / -.3817868689 0322770146 7819960960 0 d-15     /
      data ath0cs( 14) / -.8851159014 8199475941 5628650998 4 d-16     /
      data ath0cs( 15) / -.2184818181 4143659531 4967767956 8 d-16     /
      data ath0cs( 16) / -.5700849046 9864523805 9944229511 9 d-17     /
      data ath0cs( 17) / -.1563121122 1778753925 1603179549 5 d-17     /
      data ath0cs( 18) / -.4481437996 7689950679 0668877635 3 d-18     /
      data ath0cs( 19) / -.1337794883 7361880220 4456604409 8 d-18     /
      data ath0cs( 20) / -.4143340036 8741144537 7685244544 2 d-19     /
      data ath0cs( 21) / -.1327263385 7188050250 8048116465 2 d-19     /
      data ath0cs( 22) / -.4385728589 1284405222 1575683595 5 d-20     /
      data ath0cs( 23) / -.1491360695 9528180676 8620174395 6 d-20     /
      data ath0cs( 24) / -.5208104738 6307113771 5423818877 3 d-21     /
      data ath0cs( 25) / -.1864382222 3904989238 7252660497 9 d-21     /
      data ath0cs( 26) / -.6830263751 1679690129 7543538188 1 d-22     /
      data ath0cs( 27) / -.2557117058 0293296292 9620759134 7 d-22     /
      data ath0cs( 28) / -.9770158640 2543002182 4690725404 6 d-23     /
      data ath0cs( 29) / -.3805161433 4166790840 6842825488 6 d-23     /
      data ath0cs( 30) / -.1509022750 7370540634 9392648299 5 d-23     /
      data ath0cs( 31) / -.6087551341 2424249290 0556801452 5 d-24     /
      data ath0cs( 32) / -.2495879513 8097114954 2598212405 8 d-24     /
      data ath0cs( 33) / -.1039157654 5819209489 0958808427 4 d-24     /
      data ath0cs( 34) / -.4390235913 9768465369 7459496905 1 d-25     /
      data ath0cs( 35) / -.1880790678 4479902116 7582682058 2 d-25     /
      data ath0cs( 36) / -.8165070764 1994629488 6302220575 3 d-26     /
      data ath0cs( 37) / -.3589944503 7497505142 6643558504 1 d-26     /
      data ath0cs( 38) / -.1597658126 6321328729 8129160870 8 d-26     /
      data ath0cs( 39) / -.7193250175 7038239691 1380283530 5 d-27     /
      data ath0cs( 40) / -.3274943012 7278565062 0935113272 1 d-27     /
      data ath0cs( 41) / -.1507042445 7836906658 1697504727 2 d-27     /
      data ath0cs( 42) / -.7006624198 3199047178 4396794914 0 d-28     /
      data ath0cs( 43) / -.3289907402 9837182265 2881567835 6 d-28     /
      data ath0cs( 44) / -.1559518084 3651465264 4532271149 6 d-28     /
      data ath0cs( 45) / -.7460690508 2082545828 3385111972 1 d-29     /
      data ath0cs( 46) / -.3600877034 8246620205 6327724943 1 d-29     /
      data ath0cs( 47) / -.1752851437 4737722573 5040221919 7 d-29     /
      data ath0cs( 48) / -.8603275775 1885129096 2377862872 4 d-30     /
      data ath0cs( 49) / -.4256432603 2269465346 6803948010 5 d-30     /
      data ath0cs( 50) / -.2122161865 0442629277 2365069820 6 d-30     /
      data ath0cs( 51) / -.1065996156 7048790524 7206079856 1 d-30     /
      data ath0cs( 52) / -.5393568608 8169491164 1068808689 2 d-31     /
      data ath0cs( 53) / -.2748174851 0439548222 7849651787 0 d-31     /
c
c series for am21       on the interval -1.25000e-01 to -1.56250e-02
c                                        with weighted error   3.40e-32
c                                         log weighted error  31.47
c                               significant figures required  29.02
c                                    decimal places required  32.36
c
      data am21cs(  1) / +.5927902667 2130958837 5717482814 d-2        /
      data am21cs(  2) / +.2005694053 9316518642 8695217690 d-2        /
      data am21cs(  3) / +.9110818502 6227589355 3072526291 d-4        /
      data am21cs(  4) / +.8498943063 7204715563 3172107475 d-5        /
      data am21cs(  5) / +.1132979089 7691307663 7929215494 d-5        /
      data am21cs(  6) / +.1875179461 0066649618 0950627804 d-6        /
      data am21cs(  7) / +.3593065190 1824583269 9035211192 d-7        /
      data am21cs(  8) / +.7657577140 7168386403 9093517470 d-8        /
      data am21cs(  9) / +.1769999671 6803917392 5953460744 d-8        /
      data am21cs( 10) / +.4362595556 5459893272 0546585535 d-9        /
      data am21cs( 11) / +.1132916413 3785323003 5520085219 d-9        /
      data am21cs( 12) / +.3072576909 8241924413 7868398126 d-10       /
      data am21cs( 13) / +.8644824164 8220107554 1200465766 d-11       /
      data am21cs( 14) / +.2510152500 6092440211 5104562212 d-11       /
      data am21cs( 15) / +.7491024967 6444037160 1802227751 d-12       /
      data am21cs( 16) / +.2289969284 8799407308 9565214432 d-12       /
      data am21cs( 17) / +.7151136589 2798769494 9327491175 d-13       /
      data am21cs( 18) / +.2276079249 5956684194 6395165061 d-13       /
      data am21cs( 19) / +.7369421427 6088651396 9953227782 d-14       /
      data am21cs( 20) / +.2423286752 6782749046 3991742006 d-14       /
      data am21cs( 21) / +.8081537745 4823986928 3406558403 d-15       /
      data am21cs( 22) / +.2730080798 0435608665 9174563386 d-15       /
      data am21cs( 23) / +.9332360708 9138531847 3519474326 d-16       /
      data am21cs( 24) / +.3225080996 8108462221 3867546973 d-16       /
      data am21cs( 25) / +.1125819323 4644454121 7757573416 d-16       /
      data am21cs( 26) / +.3966994639 8693882166 0259459530 d-17       /
      data am21cs( 27) / +.1410065679 4431950466 0865034527 d-17       /
      data am21cs( 28) / +.5053020865 3785121337 5537393032 d-18       /
      data am21cs( 29) / +.1824615232 1594514119 7999102789 d-18       /
      data am21cs( 30) / +.6635845682 6213046692 8029121642 d-19       /
      data am21cs( 31) / +.2429637316 3127617974 1747455826 d-19       /
      data am21cs( 32) / +.8952389151 2368780201 3669922963 d-20       /
      data am21cs( 33) / +.3318452893 5005079126 0229250755 d-20       /
      data am21cs( 34) / +.1237061961 8865831538 4437905922 d-20       /
      data am21cs( 35) / +.4636366770 1239084030 6767734243 d-21       /
      data am21cs( 36) / +.1746531359 4776447546 9758765989 d-21       /
      data am21cs( 37) / +.6611168102 3499117630 7910643111 d-22       /
      data am21cs( 38) / +.2514099189 9407248617 6125666459 d-22       /
      data am21cs( 39) / +.9602749955 7173256869 4034386998 d-23       /
      data am21cs( 40) / +.3683249522 8929639568 6436898078 d-23       /
      data am21cs( 41) / +.1418431382 6915913614 5535939553 d-23       /
      data am21cs( 42) / +.5483426742 7693583010 6345800990 d-24       /
      data am21cs( 43) / +.2127610546 2311880665 0372562616 d-24       /
      data am21cs( 44) / +.8284437008 4941859148 7734760953 d-25       /
      data am21cs( 45) / +.3236705639 2612700142 1028600927 d-25       /
      data am21cs( 46) / +.1268688829 6328605735 5055062493 d-25       /
      data am21cs( 47) / +.4988438189 9212162693 5068934362 d-26       /
      data am21cs( 48) / +.1967345844 6764939096 7119381790 d-26       /
      data am21cs( 49) / +.7781359710 2032695771 3212064836 d-27       /
      data am21cs( 50) / +.3086339414 9891115291 9192968451 d-27       /
      data am21cs( 51) / +.1227446470 4545311978 9338037234 d-27       /
      data am21cs( 52) / +.4894312791 3429220588 5241216204 d-28       /
      data am21cs( 53) / +.1956468798 0290982117 5925099724 d-28       /
      data am21cs( 54) / +.7839889529 2242617116 6311492266 d-29       /
      data am21cs( 55) / +.3148969140 0248422374 8298978099 d-29       /
      data am21cs( 56) / +.1267697631 3725068130 7067842559 d-29       /
      data am21cs( 57) / +.5114706919 0690014164 1632107724 d-30       /
      data am21cs( 58) / +.2068017097 9553877025 0900316706 d-30       /
      data am21cs( 59) / +.8378913447 6851900132 5996867583 d-31       /
      data am21cs( 60) / +.3401689919 7148980205 2339079577 d-31       /
c
c series for ath1       on the interval -1.25000e-01 to -1.56250e-02
c                                        with weighted error   2.94e-32
c                                         log weighted error  31.53
c                               significant figures required  30.08
c                                    decimal places required  32.41
c
      data ath1cs(  1) / -.6972849916 2088838458 8814841503 7 d-1      /
      data ath1cs(  2) / -.5108722790 6500449870 7344807796 1 d-2      /
      data ath1cs(  3) / -.8644335996 9897550945 2533474951 2 d-4      /
      data ath1cs(  4) / -.5604720044 2352635421 8869891612 5 d-5      /
      data ath1cs(  5) / -.6045735125 6238974091 5637664007 7 d-6      /
      data ath1cs(  6) / -.8639802632 4883343932 1972113849 9 d-7      /
      data ath1cs(  7) / -.1480809484 3099271571 4778248078 0 d-7      /
      data ath1cs(  8) / -.2885809334 5772360399 9944990871 2 d-8      /
      data ath1cs(  9) / -.6191631975 6656996093 0919123180 0 d-9      /
      data ath1cs( 10) / -.1431992808 8609578309 3136525987 9 d-9      /
      data ath1cs( 11) / -.3518141102 1372147215 0461687432 1 d-10     /
      data ath1cs( 12) / -.9084761919 9550782900 7033980805 1 d-11     /
      data ath1cs( 13) / -.2446171672 6885984493 4328366476 7 d-11     /
      data ath1cs( 14) / -.6826083203 2134462408 2899671026 4 d-12     /
      data ath1cs( 15) / -.1964579931 1949401712 7854625780 2 d-12     /
      data ath1cs( 16) / -.5808933227 1396931640 0919126585 6 d-13     /
      data ath1cs( 17) / -.1759042249 5274419927 9540095902 4 d-13     /
      data ath1cs( 18) / -.5440902932 7148966136 3253894531 9 d-14     /
      data ath1cs( 19) / -.1715247407 4868068026 2235851945 1 d-14     /
      data ath1cs( 20) / -.5500929233 5769915468 7110184716 1 d-15     /
      data ath1cs( 21) / -.1791878287 7393172594 9515263875 4 d-15     /
      data ath1cs( 22) / -.5920372520 0866941977 7841106223 1 d-16     /
      data ath1cs( 23) / -.1981713027 8764839624 7097220659 0 d-16     /
      data ath1cs( 24) / -.6713232347 0163522620 4998434379 0 d-17     /
      data ath1cs( 25) / -.2299450243 6582811161 2235861983 2 d-17     /
      data ath1cs( 26) / -.7957300928 2363765953 0463714563 4 d-18     /
      data ath1cs( 27) / -.2779994027 2917841571 7229023373 9 d-18     /
      data ath1cs( 28) / -.9798924361 3269852244 0679548081 4 d-19     /
      data ath1cs( 29) / -.3482717006 0615743867 0264556584 9 d-19     /
      data ath1cs( 30) / -.1247489122 5585990571 7330005808 4 d-19     /
      data ath1cs( 31) / -.4501210041 4782281134 8775182445 2 d-20     /
      data ath1cs( 32) / -.1635346244 0133521355 9611416466 7 d-20     /
      data ath1cs( 33) / -.5980102897 7803362680 9876226594 1 d-21     /
      data ath1cs( 34) / -.2200246286 2861234540 2819629547 5 d-21     /
      data ath1cs( 35) / -.8142463073 5150858974 0820529151 9 d-22     /
      data ath1cs( 36) / -.3029924773 6600425374 3233070967 4 d-22     /
      data ath1cs( 37) / -.1133390098 5746235377 2294396968 9 d-22     /
      data ath1cs( 38) / -.4260766024 7492957192 8304988979 1 d-23     /
      data ath1cs( 39) / -.1609363396 2781897187 9750063445 3 d-23     /
      data ath1cs( 40) / -.6106377190 8250262930 4533044428 7 d-24     /
      data ath1cs( 41) / -.2326954318 0216940618 3657788757 3 d-24     /
      data ath1cs( 42) / -.8903987877 4722526044 7412955818 6 d-25     /
      data ath1cs( 43) / -.3420558530 0056750241 1791475234 1 d-25     /
      data ath1cs( 44) / -.1319026715 2572726590 1721210060 7 d-25     /
      data ath1cs( 45) / -.5104899493 6120430913 1619117738 6 d-26     /
      data ath1cs( 46) / -.1982599478 4745474512 4244466346 6 d-26     /
      data ath1cs( 47) / -.7725702356 8808305356 3611185151 9 d-27     /
      data ath1cs( 48) / -.3020234733 6646801008 1577686357 3 d-27     /
      data ath1cs( 49) / -.1184379739 0741699937 1294638080 0 d-27     /
      data ath1cs( 50) / -.4658430227 9223085205 7325284010 6 d-28     /
      data ath1cs( 51) / -.1837554188 1003846471 5750200661 3 d-28     /
      data ath1cs( 52) / -.7268566894 4279909533 2187668480 0 d-29     /
      data ath1cs( 53) / -.2882863120 3914681355 2708987562 6 d-29     /
      data ath1cs( 54) / -.1146374629 4599063504 1759166463 9 d-29     /
      data ath1cs( 55) / -.4570031437 7485330581 7999168853 3 d-30     /
      data ath1cs( 56) / -.1826276602 0453461048 0993402879 9 d-30     /
      data ath1cs( 57) / -.7315349993 3852504691 1106635093 3 d-31     /
      data ath1cs( 58) / -.2936925599 9714297816 3781577386 6 d-31     /
c
c series for am22       on the interval -1.00000e+00 to -1.25000e-01
c                                        with weighted error   3.76e-32
c                                         log weighted error  31.42
c                               significant figures required  29.47
c                                    decimal places required  32.36
c
      data am22cs(  1) / -.1562844480 6253411275 3545828583 d-1        /
      data am22cs(  2) / +.7783364452 3968130701 8943100334 d-2        /
      data am22cs(  3) / +.8670577704 7718952840 6072812110 d-3        /
      data am22cs(  4) / +.1569662731 5611371946 9953482266 d-3        /
      data am22cs(  5) / +.3563962571 4328651132 4100666302 d-4        /
      data am22cs(  6) / +.9245983354 2504315449 5080090994 d-5        /
      data am22cs(  7) / +.2621101618 5042238952 3194982066 d-5        /
      data am22cs(  8) / +.7918822165 1601256148 9469982263 d-6        /
      data am22cs(  9) / +.2510415279 2101184780 3162690862 d-6        /
      data am22cs( 10) / +.8265223206 6540773447 2997712940 d-7        /
      data am22cs( 11) / +.2805711662 8130526439 6384290014 d-7        /
      data am22cs( 12) / +.9768210904 8468078667 4631273890 d-8        /
      data am22cs( 13) / +.3474079232 2771034328 7279035573 d-8        /
      data am22cs( 14) / +.1258281321 6983691421 9092738164 d-8        /
      data am22cs( 15) / +.4629882606 4189526449 7330784625 d-9        /
      data am22cs( 16) / +.1727282588 1360407246 8143128696 d-9        /
      data am22cs( 17) / +.6523192001 3115413514 8574124970 d-10       /
      data am22cs( 18) / +.2490471685 2098205601 9881087112 d-10       /
      data am22cs( 19) / +.9601568205 5376594807 8189890126 d-11       /
      data am22cs( 20) / +.3734480020 6772685697 4776596757 d-11       /
      data am22cs( 21) / +.1464175650 3205339172 2216189678 d-11       /
      data am22cs( 22) / +.5782654711 6851282547 5827881553 d-12       /
      data am22cs( 23) / +.2299154072 4470611856 0254184494 d-12       /
      data am22cs( 24) / +.9197807112 3199725715 0883662365 d-13       /
      data am22cs( 25) / +.3700600688 1309006580 7504045556 d-13       /
      data am22cs( 26) / +.1496757616 9867298782 3326345205 d-13       /
      data am22cs( 27) / +.6083611949 3846114872 0451399443 d-14       /
      data am22cs( 28) / +.2484040871 1512139763 5425326873 d-14       /
      data am22cs( 29) / +.1018624765 2676908072 7914465339 d-14       /
      data am22cs( 30) / +.4193838563 5275398942 9640310957 d-15       /
      data am22cs( 31) / +.1733189017 6293075614 9702493501 d-15       /
      data am22cs( 32) / +.7188219023 8850851782 0445406811 d-16       /
      data am22cs( 33) / +.2991236335 9840360771 2470896113 d-16       /
      data am22cs( 34) / +.1248689904 3323862785 5713110880 d-16       /
      data am22cs( 35) / +.5228293446 0948366192 8651193632 d-17       /
      data am22cs( 36) / +.2195329617 2471339659 5998454359 d-17       /
      data am22cs( 37) / +.9242983252 2977728115 4410024332 d-18       /
      data am22cs( 38) / +.3901577082 3609140782 5543197309 d-18       /
      data am22cs( 39) / +.1650938926 9386370721 3759030367 d-18       /
      data am22cs( 40) / +.7002218157 1599436756 5716554487 d-19       /
      data am22cs( 41) / +.2976518336 1678691557 3214963506 d-19       /
      data am22cs( 42) / +.1267965390 8690207257 1134261229 d-19       /
      data am22cs( 43) / +.5412434006 9707762868 7581725061 d-20       /
      data am22cs( 44) / +.2314873502 1815525229 6382133283 d-20       /
      data am22cs( 45) / +.9919202883 8656656346 2623851167 d-21       /
      data am22cs( 46) / +.4258030153 2373235715 8897608174 d-21       /
      data am22cs( 47) / +.1831018429 7302450167 8402003088 d-21       /
      data am22cs( 48) / +.7886787123 1107537556 4526811022 d-22       /
      data am22cs( 49) / +.3402546073 8622987495 6582997235 d-22       /
      data am22cs( 50) / +.1470208814 0571253079 1860892535 d-22       /
      data am22cs( 51) / +.6362110183 2491695773 3348071767 d-23       /
      data am22cs( 52) / +.2757070506 8098072191 9395987768 d-23       /
      data am22cs( 53) / +.1196458580 9010407135 6261780457 d-23       /
      data am22cs( 54) / +.5199125457 2924214798 1768210567 d-24       /
      data am22cs( 55) / +.2262176748 4710447526 0575286850 d-24       /
      data am22cs( 56) / +.9855261137 5443181944 8565068283 d-25       /
      data am22cs( 57) / +.4298706303 3250871722 3681286187 d-25       /
      data am22cs( 58) / +.1877236416 6158063982 9657670189 d-25       /
      data am22cs( 59) / +.8207219417 7284213726 8801052115 d-26       /
      data am22cs( 60) / +.3592146656 0461550781 2767944463 d-26       /
      data am22cs( 61) / +.1573905946 1277331561 1458940587 d-26       /
      data am22cs( 62) / +.6903297810 3933383496 5319153586 d-27       /
      data am22cs( 63) / +.3030920790 7896853460 7859331415 d-27       /
      data am22cs( 64) / +.1332049341 6048121918 5689121944 d-27       /
      data am22cs( 65) / +.5859788368 5152349011 7937981442 d-28       /
      data am22cs( 66) / +.2580168684 8948780633 8425080457 d-28       /
      data am22cs( 67) / +.1137124336 3728366722 3632182863 d-28       /
      data am22cs( 68) / +.5015925572 2606850923 6430548549 d-29       /
      data am22cs( 69) / +.2214458293 9550937332 2569708484 d-29       /
      data am22cs( 70) / +.9784702838 8650728998 4691416411 d-30       /
      data am22cs( 71) / +.4326954149 3418017011 2000952983 d-30       /
      data am22cs( 72) / +.1914972881 9399457061 2929860440 d-30       /
      data am22cs( 73) / +.8481646224 0239235417 1298331562 d-31       /
      data am22cs( 74) / +.3759470651 7395591994 7455052934 d-31       /
c
c series for ath2       on the interval -1.00000e+00 to -1.25000e-01
c                                        with weighted error   4.97e-32
c                                         log weighted error  31.30
c                               significant figures required  29.79
c                                    decimal places required  32.23
c
      data ath2cs(  1) / +.4405273458 7187789970 6112705777 5 d-2      /
      data ath2cs(  2) / -.3042919452 3184546084 8384423987 3 d-1      /
      data ath2cs(  3) / -.1385653283 7717937916 0269284265 3 d-2      /
      data ath2cs(  4) / -.1804443908 9549523026 7048691095 2 d-3      /
      data ath2cs(  5) / -.3380847108 3273086710 5746532361 8 d-4      /
      data ath2cs(  6) / -.7678183535 2290230552 5767681776 5 d-5      /
      data ath2cs(  7) / -.1967839443 7160353246 9093541707 7 d-5      /
      data ath2cs(  8) / -.5483727115 8777003615 8614365928 1 d-6      /
      data ath2cs(  9) / -.1625461550 5326124527 1269621225 8 d-6      /
      data ath2cs( 10) / -.5053049981 2688950152 7763784207 8 d-7      /
      data ath2cs( 11) / -.1631580701 1240668811 8385171561 7 d-7      /
      data ath2cs( 12) / -.5434204112 3485175079 6343669481 7 d-8      /
      data ath2cs( 13) / -.1857398556 4099003257 6385010963 0 d-8      /
      data ath2cs( 14) / -.6489512033 3261088162 1351364067 6 d-9      /
      data ath2cs( 15) / -.2310594885 8009447204 8299598707 9 d-9      /
      data ath2cs( 16) / -.8363282183 2044116828 1932954674 5 d-10     /
      data ath2cs( 17) / -.3071196844 8901914626 6066130389 1 d-10     /
      data ath2cs( 18) / -.1142367142 4327168194 0951457989 2 d-10     /
      data ath2cs( 19) / -.4298116066 3458030658 2247010897 1 d-11     /
      data ath2cs( 20) / -.1633898699 5967154406 0164608663 2 d-11     /
      data ath2cs( 21) / -.6269328620 0166194321 2344375407 6 d-12     /
      data ath2cs( 22) / -.2426052694 8162573573 5615920399 1 d-12     /
      data ath2cs( 23) / -.9461198321 6240390907 4252776505 2 d-13     /
      data ath2cs( 24) / -.3716060313 4115048068 4779828126 9 d-13     /
      data ath2cs( 25) / -.1469155684 0975267631 7013881030 9 d-13     /
      data ath2cs( 26) / -.5843694726 1409119445 5640136309 4 d-14     /
      data ath2cs( 27) / -.2337502595 5919512988 3267503493 4 d-14     /
      data ath2cs( 28) / -.9399231371 1714354011 6016735841 1 d-15     /
      data ath2cs( 29) / -.3798014669 3728945000 7633526371 5 d-15     /
      data ath2cs( 30) / -.1541731043 9849725248 8344368177 5 d-15     /
      data ath2cs( 31) / -.6285287079 5353071629 2566236520 2 d-16     /
      data ath2cs( 32) / -.2572731812 8114554247 5538399277 4 d-16     /
      data ath2cs( 33) / -.1057098119 3540178093 4097486655 5 d-16     /
      data ath2cs( 34) / -.4359080267 4026969666 9599269996 4 d-17     /
      data ath2cs( 35) / -.1803634315 9599780139 5317694554 0 d-17     /
      data ath2cs( 36) / -.7486838064 3805368217 1943167691 4 d-18     /
      data ath2cs( 37) / -.3117261367 3476046567 9959720998 5 d-18     /
      data ath2cs( 38) / -.1301687980 9277007347 9287162069 6 d-18     /
      data ath2cs( 39) / -.5450527587 5195224689 7388390990 9 d-19     /
      data ath2cs( 40) / -.2288293490 1142318722 6863593190 3 d-19     /
      data ath2cs( 41) / -.9631059503 8295386556 5506044008 8 d-20     /
      data ath2cs( 42) / -.4063281001 5246140890 9219541643 4 d-20     /
      data ath2cs( 43) / -.1718203980 9080267639 0041385851 0 d-20     /
      data ath2cs( 44) / -.7281574619 8925363674 1532247332 8 d-21     /
      data ath2cs( 45) / -.3092352652 6806431279 6068034579 0 d-21     /
      data ath2cs( 46) / -.1315917855 9654404903 8341702325 4 d-21     /
      data ath2cs( 47) / -.5610606786 0870555126 6490741266 8 d-22     /
      data ath2cs( 48) / -.2396621894 0863552060 2030433789 5 d-22     /
      data ath2cs( 49) / -.1025574332 3905812008 3295442392 4 d-22     /
      data ath2cs( 50) / -.4396264138 1436564764 0360732366 3 d-23     /
      data ath2cs( 51) / -.1887652998 3725773733 4250871945 0 d-23     /
      data ath2cs( 52) / -.8118140359 5768076035 7943323044 5 d-24     /
      data ath2cs( 53) / -.3496734274 3662868563 7595208921 4 d-24     /
      data ath2cs( 54) / -.1508402925 1568732151 7175147586 7 d-24     /
      data ath2cs( 55) / -.6516268284 7786710597 8777383434 1 d-25     /
      data ath2cs( 56) / -.2818945797 5292074245 0594211458 3 d-25     /
      data ath2cs( 57) / -.1221127596 5122627445 9809446450 5 d-25     /
      data ath2cs( 58) / -.5296674341 1698671686 2001170507 3 d-26     /
      data ath2cs( 59) / -.2300359270 7736734313 5887097174 4 d-26     /
      data ath2cs( 60) / -.1000279482 3553674947 8122034893 0 d-26     /
      data ath2cs( 61) / -.4354760404 1808793948 0689316217 9 d-27     /
      data ath2cs( 62) / -.1898056134 7414775225 1548282703 0 d-27     /
      data ath2cs( 63) / -.8282111868 7129746975 5400930931 5 d-28     /
      data ath2cs( 64) / -.3617815493 0665690065 8621348437 4 d-28     /
      data ath2cs( 65) / -.1582018896 1780036548 5894184363 6 d-28     /
      data ath2cs( 66) / -.6925068597 8022700117 7282038324 7 d-29     /
      data ath2cs( 67) / -.3034390239 7786291289 0862972733 5 d-29     /
      data ath2cs( 68) / -.1330889568 1667252247 6197744650 9 d-29     /
      data ath2cs( 69) / -.5842848522 1730901204 8760697170 6 d-30     /
      data ath2cs( 70) / -.2567488423 2383026311 2127435767 8 d-30     /
      data ath2cs( 71) / -.1129232322 2688821857 9150581915 1 d-30     /
      data ath2cs( 72) / -.4970947029 7533369165 5057010502 3 d-31     /
c
      data pi4 / 0.7853981633 9744830961 5660845819 88d0 /
      data nam20, nath0, nam21, nath1, nam22, nath2 / 6*0 /
      data xsml / 0.d0 /
c
      if (nam20.ne.0) go to 10
      eta = 0.1*sngl(d1mach(3))
      nam20 = initds (am20cs, 57, eta)
      nath0 = initds (ath0cs, 53, eta)
      nam21 = initds (am21cs, 60, eta)
      nath1 = initds (ath1cs, 58, eta)
      nam22 = initds (am22cs, 74, eta)
      nath2 = initds (ath2cs, 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, am20cs, nam20)
      theta = -0.625d0 + dcsevl (z, ath0cs, nath0)
      go to 40
c
 20   if (x.ge.(-2.0d0)) go to 30
      z = (128.d0/x**3 + 9.0d0)/7.0d0
      ampl = 0.3125d0 + dcsevl (z, am21cs, nam21)
      theta = -0.625d0 + dcsevl (z, ath1cs, nath1)
      go to 40
c
 30   if (x.gt.(-1.0d0)) call seteru (
     1  25hd9aimp  x must be le -1.0, 25, 1, 2)
c
      z = (16.d0/x**3 + 9.0d0)/7.0d0
      ampl = 0.3125d0 + dcsevl (z, am22cs, nam22)
      theta = -0.625d0 + dcsevl (z, ath2cs, nath2)
c
 40   sqrtx = dsqrt(-x)
      ampl = dsqrt (ampl/sqrtx)
      theta = pi4 - x*sqrtx*theta
c
      return
      end