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