Blame view
fvn_fnlib/d9aimp.f
30.3 KB
38581db0c
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 |
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
|