1 DOUBLE PRECISION FUNCTION DE1(X) 2C***BEGIN PROLOGUE DE1 3C***DATE WRITTEN 770701 (YYMMDD) 4C***REVISION DATE 820801 (YYMMDD) 5C***CATEGORY NO. C5 6C***KEYWORDS DOUBLE PRECISION,EXPONENTIAL INTEGRAL,SPECIAL FUNCTION 7C***AUTHOR FULLERTON, W., (LANL) 8C***PURPOSE Computes the d.p. exponential integral E1(X). 9C***DESCRIPTION 10C 11C DE1(X) calculates the double precision exponential integral 12C E1(X) for positive double precision argument X and the 13C Cauchy principal value for negative X. If principal values are 14C used everywhere, then for all X 15C 16C Ei(X) = -E1(-X) or 17C E1(X) = -Ei(-X). 18C 19C 20C Series for AE10 on the interval -3.12500E-02 to 0. 21C with weighted error 4.62E-32 22C log weighted error 31.34 23C significant figures required 29.70 24C decimal places required 32.18 25C 26C 27C Series for AE11 on the interval -1.25000E-01 to -3.12500E-02 28C with weighted error 2.22E-32 29C log weighted error 31.65 30C significant figures required 30.75 31C decimal places required 32.54 32C 33C 34C Series for AE12 on the interval -2.50000E-01 to -1.25000E-01 35C with weighted error 5.19E-32 36C log weighted error 31.28 37C significant figures required 30.82 38C decimal places required 32.09 39C 40C 41C Series for E11 on the interval -4.00000E+00 to -1.00000E+00 42C with weighted error 8.49E-34 43C log weighted error 33.07 44C significant figures required 34.13 45C decimal places required 33.80 46C 47C 48C Series for E12 on the interval -1.00000E+00 to 1.00000E+00 49C with weighted error 8.08E-33 50C log weighted error 32.09 51C approx significant figures required 30.4 52C decimal places required 32.79 53C 54C 55C Series for AE13 on the interval 2.50000E-01 to 1.00000E+00 56C with weighted error 6.65E-32 57C log weighted error 31.18 58C significant figures required 30.69 59C decimal places required 32.03 60C 61C 62C Series for AE14 on the interval 0. to 2.50000E-01 63C with weighted error 5.07E-32 64C log weighted error 31.30 65C significant figures required 30.40 66C decimal places required 32.20 67C***REFERENCES (NONE) 68C***ROUTINES CALLED D1MACH,DCSEVL,INITDS,XERROR 69C***END PROLOGUE DE1 70 DOUBLE PRECISION X, AE10CS(50), AE11CS(60), AE12CS(41), E11CS(29), 71 1 E12CS(25), AE13CS(50), AE14CS(64), XMAX, D1MACH, DCSEVL 72 DATA AE10CS( 1) / +.3284394579 6166990878 7384420188 1 D-1 / 73 DATA AE10CS( 2) / -.1669920452 0313628514 7618434338 7 D-1 / 74 DATA AE10CS( 3) / +.2845284724 3613468074 2489985325 2 D-3 / 75 DATA AE10CS( 4) / -.7563944358 5162064894 8786693853 3 D-5 / 76 DATA AE10CS( 5) / +.2798971289 4508591575 0484318087 9 D-6 / 77 DATA AE10CS( 6) / -.1357901828 5345310695 2556392625 5 D-7 / 78 DATA AE10CS( 7) / +.8343596202 0404692558 5610290490 6 D-9 / 79 DATA AE10CS( 8) / -.6370971727 6402484382 7524298853 2 D-10 / 80 DATA AE10CS( 9) / +.6007247608 8118612357 6083156158 4 D-11 / 81 DATA AE10CS( 10) / -.7022876174 6797735907 5062615008 8 D-12 / 82 DATA AE10CS( 11) / +.1018302673 7036876930 9665234688 3 D-12 / 83 DATA AE10CS( 12) / -.1761812903 4308800404 0630996642 2 D-13 / 84 DATA AE10CS( 13) / +.3250828614 2353606942 4403035387 7 D-14 / 85 DATA AE10CS( 14) / -.5071770025 5058186788 2487225904 4 D-15 / 86 DATA AE10CS( 15) / +.1665177387 0432942981 7248608415 6 D-16 / 87 DATA AE10CS( 16) / +.3166753890 7975144006 7700353655 5 D-16 / 88 DATA AE10CS( 17) / -.1588403763 6641415151 3311834353 8 D-16 / 89 DATA AE10CS( 18) / +.4175513256 1380188330 0303461848 4 D-17 / 90 DATA AE10CS( 19) / -.2892347749 7071419067 1071447885 2 D-18 / 91 DATA AE10CS( 20) / -.2800625903 3966081035 0634058966 9 D-18 / 92 DATA AE10CS( 21) / +.1322938639 5392709037 0758002378 1 D-18 / 93 DATA AE10CS( 22) / -.1804447444 1773016272 8388783355 7 D-19 / 94 DATA AE10CS( 23) / -.7905384086 5226160762 9164481760 4 D-20 / 95 DATA AE10CS( 24) / +.4435711366 3695701039 4623583802 7 D-20 / 96 DATA AE10CS( 25) / -.4264103994 9781208688 6530920655 5 D-21 / 97 DATA AE10CS( 26) / -.3920101766 9371175415 5371316204 8 D-21 / 98 DATA AE10CS( 27) / +.1527378051 3439942663 4375232697 1 D-21 / 99 DATA AE10CS( 28) / +.1024849527 0493723393 1030878311 7 D-22 / 100 DATA AE10CS( 29) / -.2134907874 7714335762 6271140588 2 D-22 / 101 DATA AE10CS( 30) / +.3239139475 1600282670 6169470036 6 D-23 / 102 DATA AE10CS( 31) / +.2142183762 2998899547 6264316829 6 D-23 / 103 DATA AE10CS( 32) / -.8234609419 6010184147 0034808231 2 D-24 / 104 DATA AE10CS( 33) / -.1524652829 6458094796 1369440114 0 D-24 / 105 DATA AE10CS( 34) / +.1378208282 4606391346 6848036432 5 D-24 / 106 DATA AE10CS( 35) / +.2131311202 8339478795 2322499925 3 D-26 / 107 DATA AE10CS( 36) / -.2012649651 5264841218 1746676312 7 D-25 / 108 DATA AE10CS( 37) / +.1995535662 2633580161 0631178267 3 D-26 / 109 DATA AE10CS( 38) / +.2798995808 9840034649 4868652031 9 D-26 / 110 DATA AE10CS( 39) / -.5534511845 3896266376 4081927782 3 D-27 / 111 DATA AE10CS( 40) / -.3884995396 1599688616 8254402614 6 D-27 / 112 DATA AE10CS( 41) / +.1121304434 5073593828 5068035467 9 D-27 / 113 DATA AE10CS( 42) / +.5566568152 4237409482 5656383351 4 D-28 / 114 DATA AE10CS( 43) / -.2045482929 8104997004 4853393817 6 D-28 / 115 DATA AE10CS( 44) / -.8453813992 7123362334 1145749367 4 D-29 / 116 DATA AE10CS( 45) / +.3565758433 4312915628 1611111628 7 D-29 / 117 DATA AE10CS( 46) / +.1383653872 1256347055 3994909887 1 D-29 / 118 DATA AE10CS( 47) / -.6062167864 4513724365 8453376477 8 D-30 / 119 DATA AE10CS( 48) / -.2447198043 9893132674 3765511918 9 D-30 / 120 DATA AE10CS( 49) / +.1006850640 9339983480 1154818048 0 D-30 / 121 DATA AE10CS( 50) / +.4623685555 0148690156 6434146167 4 D-31 / 122 DATA AE11CS( 1) / +.2026315064 7078889499 4012365173 81 D+0 / 123 DATA AE11CS( 2) / -.7365514099 1203130439 5368987280 34 D-1 / 124 DATA AE11CS( 3) / +.6390934911 8361915862 7532838400 20 D-2 / 125 DATA AE11CS( 4) / -.6079725270 5247911780 6531533639 99 D-3 / 126 DATA AE11CS( 5) / -.7370649862 0176629330 6814114934 84 D-4 / 127 DATA AE11CS( 6) / +.4873285744 9450183453 4649924880 76 D-4 / 128 DATA AE11CS( 7) / -.2383706484 0448290766 5884894602 35 D-5 / 129 DATA AE11CS( 8) / -.3051861262 8561521027 0273322461 21 D-5 / 130 DATA AE11CS( 9) / +.1705033157 2564559009 6880329929 07 D-6 / 131 DATA AE11CS( 10) / +.2383420452 7487747258 6015981364 03 D-6 / 132 DATA AE11CS( 11) / +.1078177255 6163166562 5968723640 20 D-7 / 133 DATA AE11CS( 12) / -.1795569284 7399102653 6426914465 99 D-7 / 134 DATA AE11CS( 13) / -.4128407234 1950457727 9123946404 36 D-8 / 135 DATA AE11CS( 14) / +.6862214858 8631968618 3468445266 64 D-9 / 136 DATA AE11CS( 15) / +.5313018312 0506356147 6020096759 61 D-9 / 137 DATA AE11CS( 16) / +.7879688026 1490694831 3050228935 15 D-10 / 138 DATA AE11CS( 17) / -.2626176232 9356522290 3416752712 32 D-10 / 139 DATA AE11CS( 18) / -.1548368763 6308261963 1257562941 00 D-10 / 140 DATA AE11CS( 19) / -.2581896237 7261390492 8024051225 91 D-11 / 141 DATA AE11CS( 20) / +.5954287919 1591072658 9035299593 52 D-12 / 142 DATA AE11CS( 21) / +.4645140038 7681525833 7849193214 05 D-12 / 143 DATA AE11CS( 22) / +.1155785502 3255861496 2880062037 31 D-12 / 144 DATA AE11CS( 23) / -.1047523687 0835799012 3175471896 70 D-14 / 145 DATA AE11CS( 24) / -.1189665350 2709004368 1044892609 29 D-13 / 146 DATA AE11CS( 25) / -.4774907749 0261778752 6430193499 50 D-14 / 147 DATA AE11CS( 26) / -.8107764961 5772777976 2497347541 35 D-15 / 148 DATA AE11CS( 27) / +.1343556925 0031554199 3769879981 78 D-15 / 149 DATA AE11CS( 28) / +.1413453002 2913106260 2488738812 87 D-15 / 150 DATA AE11CS( 29) / +.4945159257 3953173115 5206632328 83 D-16 / 151 DATA AE11CS( 30) / +.7988404848 0080665648 8585873993 67 D-17 / 152 DATA AE11CS( 31) / -.1400863218 8089809829 2487119353 93 D-17 / 153 DATA AE11CS( 32) / -.1481424695 8417372107 7228040016 80 D-17 / 154 DATA AE11CS( 33) / -.5582617364 6025601904 0106939371 13 D-18 / 155 DATA AE11CS( 34) / -.1144207454 2191647264 7830725445 98 D-18 / 156 DATA AE11CS( 35) / +.2537182387 9566853500 5240184799 23 D-20 / 157 DATA AE11CS( 36) / +.1320532815 4805359813 2788633890 97 D-19 / 158 DATA AE11CS( 37) / +.6293026108 1586809166 2874267894 85 D-20 / 159 DATA AE11CS( 38) / +.1768827042 4882713734 9992613325 48 D-20 / 160 DATA AE11CS( 39) / +.2326618798 5146045209 6742968874 32 D-21 / 161 DATA AE11CS( 40) / -.6780306081 1125233043 7738318441 13 D-22 / 162 DATA AE11CS( 41) / -.5944087695 9676373802 8741505318 91 D-22 / 163 DATA AE11CS( 42) / -.2361821453 1184415968 5325925034 66 D-22 / 164 DATA AE11CS( 43) / -.6021449972 4601478214 1684787445 76 D-23 / 165 DATA AE11CS( 44) / -.6551790647 4348299071 3704441446 39 D-24 / 166 DATA AE11CS( 45) / +.2938875529 7497724587 0420386993 49 D-24 / 167 DATA AE11CS( 46) / +.2260160620 0642115173 2157287585 10 D-24 / 168 DATA AE11CS( 47) / +.8953436924 5958628745 0912068730 87 D-25 / 169 DATA AE11CS( 48) / +.2401592347 1098457555 7720674577 06 D-25 / 170 DATA AE11CS( 49) / +.3411837688 8907172955 6664230434 13 D-26 / 171 DATA AE11CS( 50) / -.7161707169 4630342052 3550133452 79 D-27 / 172 DATA AE11CS( 51) / -.7562039065 9281725157 9286519807 99 D-27 / 173 DATA AE11CS( 52) / -.3377461215 7467324637 9529207808 00 D-27 / 174 DATA AE11CS( 53) / -.1047932570 3300941711 5264303322 45 D-27 / 175 DATA AE11CS( 54) / -.2165455025 2170342240 8548802013 86 D-28 / 176 DATA AE11CS( 55) / -.7529712574 5288269994 6892984320 00 D-30 / 177 DATA AE11CS( 56) / +.1910317939 2798935768 6380840004 26 D-29 / 178 DATA AE11CS( 57) / +.1149210496 6530338547 7907288337 06 D-29 / 179 DATA AE11CS( 58) / +.4389697058 2661751514 4103591936 00 D-30 / 180 DATA AE11CS( 59) / +.1232088323 9205686471 6471577258 66 D-30 / 181 DATA AE11CS( 60) / +.2222017445 7553175317 5385811626 66 D-31 / 182 DATA AE12CS( 1) / +.6362958979 6747038767 1298878068 03 D+0 / 183 DATA AE12CS( 2) / -.1308116867 5067634385 8126711211 35 D+0 / 184 DATA AE12CS( 3) / -.8436741021 3053930014 4876621297 52 D-2 / 185 DATA AE12CS( 4) / +.2656849153 1006685413 0294280689 06 D-2 / 186 DATA AE12CS( 5) / +.3282272178 1658133778 7921701425 17 D-3 / 187 DATA AE12CS( 6) / -.2378344777 1430248269 5798078510 50 D-4 / 188 DATA AE12CS( 7) / -.1143980430 8100055514 4470767970 47 D-4 / 189 DATA AE12CS( 8) / -.1440594343 3238338455 2397176993 23 D-5 / 190 DATA AE12CS( 9) / +.5241595665 1148829963 7728180616 64 D-8 / 191 DATA AE12CS( 10) / +.3840730640 7844323480 9792030597 16 D-7 / 192 DATA AE12CS( 11) / +.8588024486 0267195879 6605157593 44 D-8 / 193 DATA AE12CS( 12) / +.1021922662 5855003286 3399695539 11 D-8 / 194 DATA AE12CS( 13) / +.2174913232 3289724542 8213398059 92 D-10 / 195 DATA AE12CS( 14) / -.2209023814 2623144809 5235038117 41 D-10 / 196 DATA AE12CS( 15) / -.6345753354 4928753294 3836222088 01 D-11 / 197 DATA AE12CS( 16) / -.1083774656 6857661115 3405397329 19 D-11 / 198 DATA AE12CS( 17) / -.1190982287 2222586730 2622004402 77 D-12 / 199 DATA AE12CS( 18) / -.2843868238 9265590299 5087660086 61 D-14 / 200 DATA AE12CS( 19) / +.2508032702 6686769668 5871954875 46 D-14 / 201 DATA AE12CS( 20) / +.7872964152 8559842431 5977264212 65 D-15 / 202 DATA AE12CS( 21) / +.1547506634 7785217148 4843346373 29 D-15 / 203 DATA AE12CS( 22) / +.2257532283 1665075055 2726081972 90 D-16 / 204 DATA AE12CS( 23) / +.2223335286 7266608760 2813808366 93 D-17 / 205 DATA AE12CS( 24) / +.1696781956 3544153513 4641946623 99 D-19 / 206 DATA AE12CS( 25) / -.5760831625 5947682105 3100873045 33 D-19 / 207 DATA AE12CS( 26) / -.1759123577 4646878055 6253694088 53 D-19 / 208 DATA AE12CS( 27) / -.3628605637 5103174394 7553286826 66 D-20 / 209 DATA AE12CS( 28) / -.5923556979 7328991652 5581434880 00 D-21 / 210 DATA AE12CS( 29) / -.7603038092 6310191114 4291368959 99 D-22 / 211 DATA AE12CS( 30) / -.6254784352 1711763842 6414284799 99 D-23 / 212 DATA AE12CS( 31) / +.2548336075 9307648606 0376064000 00 D-24 / 213 DATA AE12CS( 32) / +.2559861573 1739857020 1688746666 66 D-24 / 214 DATA AE12CS( 33) / +.7137623935 7899318800 2070528000 00 D-25 / 215 DATA AE12CS( 34) / +.1470375993 9567568181 5789568000 00 D-25 / 216 DATA AE12CS( 35) / +.2510552476 5386733555 1986346666 66 D-26 / 217 DATA AE12CS( 36) / +.3588666638 7790890886 5836373333 33 D-27 / 218 DATA AE12CS( 37) / +.3988603515 6771301763 3177599999 99 D-28 / 219 DATA AE12CS( 38) / +.2176367694 7356220478 8053333333 33 D-29 / 220 DATA AE12CS( 39) / -.4614699848 7618942367 6074666666 66 D-30 / 221 DATA AE12CS( 40) / -.2071351787 7481987707 1530666666 66 D-30 / 222 DATA AE12CS( 41) / -.5189037856 3534371596 9706666666 66 D-31 / 223 DATA E11 CS( 1) / -.1611346165 5571494025 7206639275 66180 D+2 / 224 DATA E11 CS( 2) / +.7794072778 7426802769 2722458917 41497 D+1 / 225 DATA E11 CS( 3) / -.1955405818 8631419507 1272838128 14491 D+1 / 226 DATA E11 CS( 4) / +.3733729386 6277945611 5171908656 90209 D+0 / 227 DATA E11 CS( 5) / -.5692503191 0929019385 2638922200 51166 D-1 / 228 DATA E11 CS( 6) / +.7211077769 6600918537 8477248126 35813 D-2 / 229 DATA E11 CS( 7) / -.7810490144 9841593997 7151840890 64148 D-3 / 230 DATA E11 CS( 8) / +.7388093356 2621681878 9748813661 77858 D-4 / 231 DATA E11 CS( 9) / -.6202861875 8082045134 3581336079 09712 D-5 / 232 DATA E11 CS( 10) / +.4681600230 3176735524 4058238683 62657 D-6 / 233 DATA E11 CS( 11) / -.3209288853 3298649524 0725530272 28719 D-7 / 234 DATA E11 CS( 12) / +.2015199748 7404533394 8262622130 19548 D-8 / 235 DATA E11 CS( 13) / -.1167368681 6697793105 3562716950 15419 D-9 / 236 DATA E11 CS( 14) / +.6276270667 2039943397 7887483796 15573 D-11 / 237 DATA E11 CS( 15) / -.3148154167 2275441045 2467818023 93600 D-12 / 238 DATA E11 CS( 16) / +.1479904174 4493474210 8944722517 33333 D-13 / 239 DATA E11 CS( 17) / -.6545709158 3979673774 2634015880 53333 D-15 / 240 DATA E11 CS( 18) / +.2733687222 3137291142 5080127487 99999 D-16 / 241 DATA E11 CS( 19) / -.1081352434 9754406876 7217276245 33333 D-17 / 242 DATA E11 CS( 20) / +.4062832804 0434303295 3003485866 66666 D-19 / 243 DATA E11 CS( 21) / -.1453553935 8960455858 9143722666 66666 D-20 / 244 DATA E11 CS( 22) / +.4963274618 1648636830 1984426666 66666 D-22 / 245 DATA E11 CS( 23) / -.1620861269 6636044604 8665600000 00000 D-23 / 246 DATA E11 CS( 24) / +.5072144803 8607422226 4319999999 99999 D-25 / 247 DATA E11 CS( 25) / -.1523581113 3372207813 9733333333 33333 D-26 / 248 DATA E11 CS( 26) / +.4400151125 6103618696 5333333333 33333 D-28 / 249 DATA E11 CS( 27) / -.1223614194 5416231594 6666666666 66666 D-29 / 250 DATA E11 CS( 28) / +.3280921666 1066001066 6666666666 66666 D-31 / 251 DATA E11 CS( 29) / -.8493345226 8306432000 0000000000 00000 D-33 / 252 DATA E12 CS( 1) / -.3739021479 22027951166 869820482 7 D-1 / 253 DATA E12 CS( 2) / +.4272398606 2209577260 4917917652 8 D-1 / 254 DATA E12 CS( 3) / -.1303182079 8497005441 5392055219 726 D+0 / 255 DATA E12 CS( 4) / +.1441912402 4698890734 1095893982 137 D-1 / 256 DATA E12 CS( 5) / -.1346170780 5106802211 6121527983 553 D-2 / 257 DATA E12 CS( 6) / +.1073102925 3063779997 6115850970 073 D-3 / 258 DATA E12 CS( 7) / -.7429999516 1194364961 0283062223 163 D-5 / 259 DATA E12 CS( 8) / +.4537732569 0753713938 6383211511 827 D-6 / 260 DATA E12 CS( 9) / -.2476417211 3906013184 6547423802 912 D-7 / 261 DATA E12 CS( 10) / +.1220765813 7459095370 0228167846 102 D-8 / 262 DATA E12 CS( 11) / -.5485141480 6409239382 1357398028 261 D-10 / 263 DATA E12 CS( 12) / +.2263621421 3007879929 3688162377 002 D-11 / 264 DATA E12 CS( 13) / -.8635897271 6980097940 4172916282 240 D-13 / 265 DATA E12 CS( 14) / +.3062915536 6933299758 1032894881 279 D-14 / 266 DATA E12 CS( 15) / -.1014857188 5594414755 7128906734 933 D-15 / 267 DATA E12 CS( 16) / +.3154821740 3406987754 6855328426 666 D-17 / 268 DATA E12 CS( 17) / -.9236042407 6924095448 4015923200 000 D-19 / 269 DATA E12 CS( 18) / +.2555042679 7081400244 0435029333 333 D-20 / 270 DATA E12 CS( 19) / -.6699128056 8456684721 7882453333 333 D-22 / 271 DATA E12 CS( 20) / +.1669254054 3538731943 1987199999 999 D-23 / 272 DATA E12 CS( 21) / -.3962549251 8437964185 6000000000 000 D-25 / 273 DATA E12 CS( 22) / +.8981358965 9851133201 0666666666 666 D-27 / 274 DATA E12 CS( 23) / -.1947633669 9301643332 2666666666 666 D-28 / 275 DATA E12 CS( 24) / +.4048360190 2463003306 6666666666 666 D-30 / 276 DATA E12 CS( 25) / -.8079815676 9984512000 0000000000 000 D-32 / 277 DATA AE13CS( 1) / -.6057732466 4060345999 3193827377 47 D+0 / 278 DATA AE13CS( 2) / -.1125352434 8366090030 6497688527 18 D+0 / 279 DATA AE13CS( 3) / +.1343226624 7902779492 4878593294 14 D-1 / 280 DATA AE13CS( 4) / -.1926845187 3811457249 2468389913 03 D-2 / 281 DATA AE13CS( 5) / +.3091183377 2060318335 5867374753 68 D-3 / 282 DATA AE13CS( 6) / -.5356413212 9618418776 3935597951 47 D-4 / 283 DATA AE13CS( 7) / +.9827812880 2474923952 4918827172 37 D-5 / 284 DATA AE13CS( 8) / -.1885368984 9165182826 9028919389 10 D-5 / 285 DATA AE13CS( 9) / +.3749431935 6894735406 9640421905 31 D-6 / 286 DATA AE13CS( 10) / -.7682345587 0552639273 7334656805 56 D-7 / 287 DATA AE13CS( 11) / +.1614327056 7198777552 9563000608 68 D-7 / 288 DATA AE13CS( 12) / -.3466802211 4907354566 3090602260 27 D-8 / 289 DATA AE13CS( 13) / +.7587542091 9036277572 8897470541 14 D-9 / 290 DATA AE13CS( 14) / -.1688643332 9881412573 5145266367 03 D-9 / 291 DATA AE13CS( 15) / +.3814570674 9552265682 8042509272 72 D-10 / 292 DATA AE13CS( 16) / -.8733026632 4446292706 8517182723 34 D-11 / 293 DATA AE13CS( 17) / +.2023672864 5867960961 7943110643 30 D-11 / 294 DATA AE13CS( 18) / -.4741328303 9555834655 2103408201 60 D-12 / 295 DATA AE13CS( 19) / +.1122117204 8389864324 7317999289 20 D-12 / 296 DATA AE13CS( 20) / -.2680422543 4840309912 8268090933 95 D-13 / 297 DATA AE13CS( 21) / +.6457851441 7716530343 5803690672 12 D-14 / 298 DATA AE13CS( 22) / -.1568276050 1666478830 3057028491 94 D-14 / 299 DATA AE13CS( 23) / +.3836786539 9315404861 8215164414 08 D-15 / 300 DATA AE13CS( 24) / -.9451717302 7579130478 8710489325 56 D-16 / 301 DATA AE13CS( 25) / +.2343481228 8949573293 8966664391 33 D-16 / 302 DATA AE13CS( 26) / -.5845866158 0214714576 1231944198 82 D-17 / 303 DATA AE13CS( 27) / +.1466622986 7947778605 8736174191 95 D-17 / 304 DATA AE13CS( 28) / -.3699392347 6444472706 5925382744 74 D-18 / 305 DATA AE13CS( 29) / +.9379015993 6721242136 0142918178 13 D-19 / 306 DATA AE13CS( 30) / -.2389367322 1937873136 3082240873 81 D-19 / 307 DATA AE13CS( 31) / +.6115062462 9497608051 9342238378 66 D-20 / 308 DATA AE13CS( 32) / -.1571858532 7554025507 7198532881 06 D-20 / 309 DATA AE13CS( 33) / +.4057238728 5585397769 5192944913 06 D-21 / 310 DATA AE13CS( 34) / -.1051402655 4738034990 5663671227 73 D-21 / 311 DATA AE13CS( 35) / +.2734966493 0638667785 8060031317 33 D-22 / 312 DATA AE13CS( 36) / -.7140160408 0205796099 3555742719 99 D-23 / 313 DATA AE13CS( 37) / +.1870555243 2235079986 7569242111 99 D-23 / 314 DATA AE13CS( 38) / -.4916746816 6870480520 4780209493 33 D-24 / 315 DATA AE13CS( 39) / +.1296498811 9684031730 9160871253 33 D-24 / 316 DATA AE13CS( 40) / -.3429251568 8362864461 6239404373 33 D-25 / 317 DATA AE13CS( 41) / +.9097224164 3887034329 1048209066 66 D-26 / 318 DATA AE13CS( 42) / -.2420211231 4316856489 9348479999 99 D-26 / 319 DATA AE13CS( 43) / +.6456361293 4639510757 6704750933 33 D-27 / 320 DATA AE13CS( 44) / -.1726913273 5340541122 3159876266 66 D-27 / 321 DATA AE13CS( 45) / +.4630861165 9151500715 1942314666 66 D-28 / 322 DATA AE13CS( 46) / -.1244870363 7214131241 7551701333 33 D-28 / 323 DATA AE13CS( 47) / +.3354457409 0520678532 9070079999 99 D-29 / 324 DATA AE13CS( 48) / -.9059886852 1070774437 5439359999 99 D-30 / 325 DATA AE13CS( 49) / +.2452414705 1474238587 2732160000 00 D-30 / 326 DATA AE13CS( 50) / -.6652817873 3552062817 1079679999 99 D-31 / 327 DATA AE14CS( 1) / -.1892918000 7530168254 9567994282 0 D+0 / 328 DATA AE14CS( 2) / -.8648117855 2598714899 6881705682 4 D-1 / 329 DATA AE14CS( 3) / +.7224101543 7465947470 2151483918 4 D-2 / 330 DATA AE14CS( 4) / -.8097559457 5573861971 5965561018 1 D-3 / 331 DATA AE14CS( 5) / +.1099913443 2661388671 7925115700 2 D-3 / 332 DATA AE14CS( 6) / -.1717332998 9377673714 9535881448 7 D-4 / 333 DATA AE14CS( 7) / +.2985627514 4792833228 2534249500 3 D-5 / 334 DATA AE14CS( 8) / -.5659649145 7719300565 6016726715 5 D-6 / 335 DATA AE14CS( 9) / +.1152680839 7141400192 2658350166 3 D-6 / 336 DATA AE14CS( 10) / -.2495030440 2693382288 4212876506 5 D-7 / 337 DATA AE14CS( 11) / +.5692324201 8337543670 3937036814 0 D-8 / 338 DATA AE14CS( 12) / -.1359957664 8056003384 9003093917 6 D-8 / 339 DATA AE14CS( 13) / +.3384662888 7608845901 8451292585 9 D-9 / 340 DATA AE14CS( 14) / -.8737853904 4746819523 5084931658 0 D-10 / 341 DATA AE14CS( 15) / +.2331588663 2226597186 1261340047 0 D-10 / 342 DATA AE14CS( 16) / -.6411481049 2137859697 5316519632 6 D-11 / 343 DATA AE14CS( 17) / +.1812246980 2048164333 8435948468 2 D-11 / 344 DATA AE14CS( 18) / -.5253831761 5584606888 1940384046 6 D-12 / 345 DATA AE14CS( 19) / +.1559218272 5919256988 5502860982 5 D-12 / 346 DATA AE14CS( 20) / -.4729168297 0803987184 7642936946 6 D-13 / 347 DATA AE14CS( 21) / +.1463761864 3932435020 7619949380 8 D-13 / 348 DATA AE14CS( 22) / -.4617388988 7129241022 3217362360 4 D-14 / 349 DATA AE14CS( 23) / +.1482710348 2893693237 8923966037 1 D-14 / 350 DATA AE14CS( 24) / -.4841672496 2392291469 7316573441 7 D-15 / 351 DATA AE14CS( 25) / +.1606215575 7002904081 1657196618 8 D-15 / 352 DATA AE14CS( 26) / -.5408917538 9571709478 9502378425 2 D-16 / 353 DATA AE14CS( 27) / +.1847470159 3468978813 7023140231 0 D-16 / 354 DATA AE14CS( 28) / -.6395830792 7590944705 0061042505 0 D-17 / 355 DATA AE14CS( 29) / +.2242780721 6997594572 5023327617 0 D-17 / 356 DATA AE14CS( 30) / -.7961369173 9839475527 4455530864 6 D-18 / 357 DATA AE14CS( 31) / +.2859308111 5401974598 0861992927 2 D-18 / 358 DATA AE14CS( 32) / -.1038450244 7011371459 0069713744 6 D-18 / 359 DATA AE14CS( 33) / +.3812040607 0979757808 6684100831 9 D-19 / 360 DATA AE14CS( 34) / -.1413795417 7172007687 1756272369 6 D-19 / 361 DATA AE14CS( 35) / +.5295367865 1827409583 0544259481 5 D-20 / 362 DATA AE14CS( 36) / -.2002264245 0268259021 3721113143 9 D-20 / 363 DATA AE14CS( 37) / +.7640262751 2751960147 3684861091 8 D-21 / 364 DATA AE14CS( 38) / -.2941119006 8687878833 1126352336 2 D-21 / 365 DATA AE14CS( 39) / +.1141823539 0789271930 3769148358 6 D-21 / 366 DATA AE14CS( 40) / -.4469308475 9552984252 4702071848 9 D-22 / 367 DATA AE14CS( 41) / +.1763262410 5717507706 3049140852 0 D-22 / 368 DATA AE14CS( 42) / -.7009968187 9259023563 5151826234 0 D-23 / 369 DATA AE14CS( 43) / +.2807573556 5583789222 8775750751 5 D-23 / 370 DATA AE14CS( 44) / -.1132560944 9810864321 4188889156 2 D-23 / 371 DATA AE14CS( 45) / +.4600574684 3750179461 5676423372 7 D-24 / 372 DATA AE14CS( 46) / -.1881448598 9761334598 6460914810 8 D-24 / 373 DATA AE14CS( 47) / +.7744916111 5077308454 4432847803 7 D-25 / 374 DATA AE14CS( 48) / -.3208512760 5853689267 0270382626 1 D-25 / 375 DATA AE14CS( 49) / +.1337445542 9108397606 1993042138 4 D-25 / 376 DATA AE14CS( 50) / -.5608671881 8022170488 9477173521 0 D-26 / 377 DATA AE14CS( 51) / +.2365839716 5285374837 1006947327 9 D-26 / 378 DATA AE14CS( 52) / -.1003656195 0253053340 6583452685 6 D-26 / 379 DATA AE14CS( 53) / +.4281490878 0941611312 8664255692 7 D-27 / 380 DATA AE14CS( 54) / -.1836345261 8153181996 9132695825 0 D-27 / 381 DATA AE14CS( 55) / +.7917798231 3495400000 9746867814 4 D-28 / 382 DATA AE14CS( 56) / -.3431542358 7422203610 2501577523 1 D-28 / 383 DATA AE14CS( 57) / +.1494705493 8971032374 7506600891 7 D-28 / 384 DATA AE14CS( 58) / -.6542620279 8657054397 3904242005 3 D-29 / 385 DATA AE14CS( 59) / +.2877581395 1991711143 4048735368 5 D-29 / 386 DATA AE14CS( 60) / -.1271557211 7960247110 2798120004 2 D-29 / 387 DATA AE14CS( 61) / +.5644615555 6487225223 8804462250 6 D-30 / 388 DATA AE14CS( 62) / -.2516994994 2840951060 8061683029 3 D-30 / 389 DATA AE14CS( 63) / +.1127259818 9275102063 7036880418 1 D-30 / 390 DATA AE14CS( 64) / -.5069814875 8004608555 6258471936 0 D-31 / 391 DATA NTAE10, NTAE11, NTAE12, NTE11, NTE12, NTAE13, NTAE14 / 7*0 / 392 DATA XMAX / 0.D0 / 393C***FIRST EXECUTABLE STATEMENT DE1 394 IF (NTAE10.NE.0) GO TO 10 395 ETA = 0.1*SNGL(D1MACH(3)) 396 NTAE10 = INITDS (AE10CS, 50, ETA) 397 NTAE11 = INITDS (AE11CS, 60, ETA) 398 NTAE12 = INITDS (AE12CS, 41, ETA) 399 NTE11 = INITDS (E11CS, 29, ETA) 400 NTE12 = INITDS (E12CS, 25, ETA) 401 NTAE13 = INITDS (AE13CS, 50, ETA) 402 NTAE14 = INITDS (AE14CS, 64, ETA) 403C 404 XMAX = -DLOG(D1MACH(1)) 405 XMAX = XMAX - DLOG(XMAX) 406C 407 10 IF (X.GT.(-1.D0)) GO TO 50 408 IF (X.GT.(-32.D0)) GO TO 20 409 DE1 = DEXP(-X)/X * (1.D0 + DCSEVL (64.D0/X+1.D0, AE10CS, NTAE10)) 410 RETURN 411C 412 20 IF (X.GT.(-8.D0)) GO TO 30 413 DE1 = DEXP(-X)/X * (1.D0 + DCSEVL ((64.D0/X+5.D0)/3.D0, AE11CS, 414 1 NTAE11)) 415 RETURN 416C 417 30 IF (X.GT.(-4.D0)) GO TO 40 418 DE1 = DEXP(-X)/X * (1.D0 + DCSEVL (16.D0/X+3.D0, AE12CS, NTAE12)) 419 RETURN 420C 421 40 DE1 = -DLOG(-X) + DCSEVL ((2.D0*X+5.D0)/3.D0, E11CS, NTE11) 422 RETURN 423C 424 50 IF (X.GT.1.0D0) GO TO 60 425 IF (X.EQ.0.D0) CALL XERROR ( 'DE1 X IS 0', 14, 2, 2) 426 DE1 = (-DLOG(DABS(X)) - 0.6875D0 + X) + DCSEVL (X, E12CS, NTE12) 427 RETURN 428C 429 60 IF (X.GT.4.0D0) GO TO 70 430 DE1 = DEXP(-X)/X * (1.D0 + DCSEVL ((8.D0/X-5.D0)/3.D0, AE13CS, 431 1 NTAE13)) 432 RETURN 433C 434 70 IF (X.GT.XMAX) GO TO 80 435 DE1 = DEXP(-X)/X * (1.D0 + DCSEVL (8.D0/X-1.D0, AE14CS, NTAE14)) 436 RETURN 437C 438 80 CALL XERROR ( 'DE1 X SO BIG E1 UNDERFLOWS', 30, 1, 1) 439 DE1 = 0.D0 440 RETURN 441C 442 END 443