double precision function de1 (x)
c july 1977 edition.  w. fullerton, c3, los alamos scientific lab.
      double precision x, ae10cs(50), ae11cs(60), ae12cs(41), e11cs(29),
     1  e12cs(25), ae13cs(50), ae14cs(64), xmax, d1mach, dcsevl,
     2  dexp, dlog
      external d1mach, dcsevl, initds
c
c series for ae10       on the interval -3.12500e-02 to  0.
c                                        with weighted error   4.62e-32
c                                         log weighted error  31.34
c                               significant figures required  29.70
c                                    decimal places required  32.18
c
      data ae10cs(  1) / +.3284394579 6166990878 7384420188 1 d-1      /
      data ae10cs(  2) / -.1669920452 0313628514 7618434338 7 d-1      /
      data ae10cs(  3) / +.2845284724 3613468074 2489985325 2 d-3      /
      data ae10cs(  4) / -.7563944358 5162064894 8786693853 3 d-5      /
      data ae10cs(  5) / +.2798971289 4508591575 0484318087 9 d-6      /
      data ae10cs(  6) / -.1357901828 5345310695 2556392625 5 d-7      /
      data ae10cs(  7) / +.8343596202 0404692558 5610290490 6 d-9      /
      data ae10cs(  8) / -.6370971727 6402484382 7524298853 2 d-10     /
      data ae10cs(  9) / +.6007247608 8118612357 6083156158 4 d-11     /
      data ae10cs( 10) / -.7022876174 6797735907 5062615008 8 d-12     /
      data ae10cs( 11) / +.1018302673 7036876930 9665234688 3 d-12     /
      data ae10cs( 12) / -.1761812903 4308800404 0630996642 2 d-13     /
      data ae10cs( 13) / +.3250828614 2353606942 4403035387 7 d-14     /
      data ae10cs( 14) / -.5071770025 5058186788 2487225904 4 d-15     /
      data ae10cs( 15) / +.1665177387 0432942981 7248608415 6 d-16     /
      data ae10cs( 16) / +.3166753890 7975144006 7700353655 5 d-16     /
      data ae10cs( 17) / -.1588403763 6641415151 3311834353 8 d-16     /
      data ae10cs( 18) / +.4175513256 1380188330 0303461848 4 d-17     /
      data ae10cs( 19) / -.2892347749 7071419067 1071447885 2 d-18     /
      data ae10cs( 20) / -.2800625903 3966081035 0634058966 9 d-18     /
      data ae10cs( 21) / +.1322938639 5392709037 0758002378 1 d-18     /
      data ae10cs( 22) / -.1804447444 1773016272 8388783355 7 d-19     /
      data ae10cs( 23) / -.7905384086 5226160762 9164481760 4 d-20     /
      data ae10cs( 24) / +.4435711366 3695701039 4623583802 7 d-20     /
      data ae10cs( 25) / -.4264103994 9781208688 6530920655 5 d-21     /
      data ae10cs( 26) / -.3920101766 9371175415 5371316204 8 d-21     /
      data ae10cs( 27) / +.1527378051 3439942663 4375232697 1 d-21     /
      data ae10cs( 28) / +.1024849527 0493723393 1030878311 7 d-22     /
      data ae10cs( 29) / -.2134907874 7714335762 6271140588 2 d-22     /
      data ae10cs( 30) / +.3239139475 1600282670 6169470036 6 d-23     /
      data ae10cs( 31) / +.2142183762 2998899547 6264316829 6 d-23     /
      data ae10cs( 32) / -.8234609419 6010184147 0034808231 2 d-24     /
      data ae10cs( 33) / -.1524652829 6458094796 1369440114 0 d-24     /
      data ae10cs( 34) / +.1378208282 4606391346 6848036432 5 d-24     /
      data ae10cs( 35) / +.2131311202 8339478795 2322499925 3 d-26     /
      data ae10cs( 36) / -.2012649651 5264841218 1746676312 7 d-25     /
      data ae10cs( 37) / +.1995535662 2633580161 0631178267 3 d-26     /
      data ae10cs( 38) / +.2798995808 9840034649 4868652031 9 d-26     /
      data ae10cs( 39) / -.5534511845 3896266376 4081927782 3 d-27     /
      data ae10cs( 40) / -.3884995396 1599688616 8254402614 6 d-27     /
      data ae10cs( 41) / +.1121304434 5073593828 5068035467 9 d-27     /
      data ae10cs( 42) / +.5566568152 4237409482 5656383351 4 d-28     /
      data ae10cs( 43) / -.2045482929 8104997004 4853393817 6 d-28     /
      data ae10cs( 44) / -.8453813992 7123362334 1145749367 4 d-29     /
      data ae10cs( 45) / +.3565758433 4312915628 1611111628 7 d-29     /
      data ae10cs( 46) / +.1383653872 1256347055 3994909887 1 d-29     /
      data ae10cs( 47) / -.6062167864 4513724365 8453376477 8 d-30     /
      data ae10cs( 48) / -.2447198043 9893132674 3765511918 9 d-30     /
      data ae10cs( 49) / +.1006850640 9339983480 1154818048 0 d-30     /
      data ae10cs( 50) / +.4623685555 0148690156 6434146167 4 d-31     /
c
c series for ae11       on the interval -1.25000e-01 to -3.12500e-02
c                                        with weighted error   2.22e-32
c                                         log weighted error  31.65
c                               significant figures required  30.75
c                                    decimal places required  32.54
c
      data ae11cs(  1) / +.2026315064 7078889499 4012365173 81 d+0     /
      data ae11cs(  2) / -.7365514099 1203130439 5368987280 34 d-1     /
      data ae11cs(  3) / +.6390934911 8361915862 7532838400 20 d-2     /
      data ae11cs(  4) / -.6079725270 5247911780 6531533639 99 d-3     /
      data ae11cs(  5) / -.7370649862 0176629330 6814114934 84 d-4     /
      data ae11cs(  6) / +.4873285744 9450183453 4649924880 76 d-4     /
      data ae11cs(  7) / -.2383706484 0448290766 5884894602 35 d-5     /
      data ae11cs(  8) / -.3051861262 8561521027 0273322461 21 d-5     /
      data ae11cs(  9) / +.1705033157 2564559009 6880329929 07 d-6     /
      data ae11cs( 10) / +.2383420452 7487747258 6015981364 03 d-6     /
      data ae11cs( 11) / +.1078177255 6163166562 5968723640 20 d-7     /
      data ae11cs( 12) / -.1795569284 7399102653 6426914465 99 d-7     /
      data ae11cs( 13) / -.4128407234 1950457727 9123946404 36 d-8     /
      data ae11cs( 14) / +.6862214858 8631968618 3468445266 64 d-9     /
      data ae11cs( 15) / +.5313018312 0506356147 6020096759 61 d-9     /
      data ae11cs( 16) / +.7879688026 1490694831 3050228935 15 d-10    /
      data ae11cs( 17) / -.2626176232 9356522290 3416752712 32 d-10    /
      data ae11cs( 18) / -.1548368763 6308261963 1257562941 00 d-10    /
      data ae11cs( 19) / -.2581896237 7261390492 8024051225 91 d-11    /
      data ae11cs( 20) / +.5954287919 1591072658 9035299593 52 d-12    /
      data ae11cs( 21) / +.4645140038 7681525833 7849193214 05 d-12    /
      data ae11cs( 22) / +.1155785502 3255861496 2880062037 31 d-12    /
      data ae11cs( 23) / -.1047523687 0835799012 3175471896 70 d-14    /
      data ae11cs( 24) / -.1189665350 2709004368 1044892609 29 d-13    /
      data ae11cs( 25) / -.4774907749 0261778752 6430193499 50 d-14    /
      data ae11cs( 26) / -.8107764961 5772777976 2497347541 35 d-15    /
      data ae11cs( 27) / +.1343556925 0031554199 3769879981 78 d-15    /
      data ae11cs( 28) / +.1413453002 2913106260 2488738812 87 d-15    /
      data ae11cs( 29) / +.4945159257 3953173115 5206632328 83 d-16    /
      data ae11cs( 30) / +.7988404848 0080665648 8585873993 67 d-17    /
      data ae11cs( 31) / -.1400863218 8089809829 2487119353 93 d-17    /
      data ae11cs( 32) / -.1481424695 8417372107 7228040016 80 d-17    /
      data ae11cs( 33) / -.5582617364 6025601904 0106939371 13 d-18    /
      data ae11cs( 34) / -.1144207454 2191647264 7830725445 98 d-18    /
      data ae11cs( 35) / +.2537182387 9566853500 5240184799 23 d-20    /
      data ae11cs( 36) / +.1320532815 4805359813 2788633890 97 d-19    /
      data ae11cs( 37) / +.6293026108 1586809166 2874267894 85 d-20    /
      data ae11cs( 38) / +.1768827042 4882713734 9992613325 48 d-20    /
      data ae11cs( 39) / +.2326618798 5146045209 6742968874 32 d-21    /
      data ae11cs( 40) / -.6780306081 1125233043 7738318441 13 d-22    /
      data ae11cs( 41) / -.5944087695 9676373802 8741505318 91 d-22    /
      data ae11cs( 42) / -.2361821453 1184415968 5325925034 66 d-22    /
      data ae11cs( 43) / -.6021449972 4601478214 1684787445 76 d-23    /
      data ae11cs( 44) / -.6551790647 4348299071 3704441446 39 d-24    /
      data ae11cs( 45) / +.2938875529 7497724587 0420386993 49 d-24    /
      data ae11cs( 46) / +.2260160620 0642115173 2157287585 10 d-24    /
      data ae11cs( 47) / +.8953436924 5958628745 0912068730 87 d-25    /
      data ae11cs( 48) / +.2401592347 1098457555 7720674577 06 d-25    /
      data ae11cs( 49) / +.3411837688 8907172955 6664230434 13 d-26    /
      data ae11cs( 50) / -.7161707169 4630342052 3550133452 79 d-27    /
      data ae11cs( 51) / -.7562039065 9281725157 9286519807 99 d-27    /
      data ae11cs( 52) / -.3377461215 7467324637 9529207808 00 d-27    /
      data ae11cs( 53) / -.1047932570 3300941711 5264303322 45 d-27    /
      data ae11cs( 54) / -.2165455025 2170342240 8548802013 86 d-28    /
      data ae11cs( 55) / -.7529712574 5288269994 6892984320 00 d-30    /
      data ae11cs( 56) / +.1910317939 2798935768 6380840004 26 d-29    /
      data ae11cs( 57) / +.1149210496 6530338547 7907288337 06 d-29    /
      data ae11cs( 58) / +.4389697058 2661751514 4103591936 00 d-30    /
      data ae11cs( 59) / +.1232088323 9205686471 6471577258 66 d-30    /
      data ae11cs( 60) / +.2222017445 7553175317 5385811626 66 d-31    /
c
c series for ae12       on the interval -2.50000e-01 to -1.25000e-01
c                                        with weighted error   5.19e-32
c                                         log weighted error  31.28
c                               significant figures required  30.82
c                                    decimal places required  32.09
c
      data ae12cs(  1) / +.6362958979 6747038767 1298878068 03 d+0     /
      data ae12cs(  2) / -.1308116867 5067634385 8126711211 35 d+0     /
      data ae12cs(  3) / -.8436741021 3053930014 4876621297 52 d-2     /
      data ae12cs(  4) / +.2656849153 1006685413 0294280689 06 d-2     /
      data ae12cs(  5) / +.3282272178 1658133778 7921701425 17 d-3     /
      data ae12cs(  6) / -.2378344777 1430248269 5798078510 50 d-4     /
      data ae12cs(  7) / -.1143980430 8100055514 4470767970 47 d-4     /
      data ae12cs(  8) / -.1440594343 3238338455 2397176993 23 d-5     /
      data ae12cs(  9) / +.5241595665 1148829963 7728180616 64 d-8     /
      data ae12cs( 10) / +.3840730640 7844323480 9792030597 16 d-7     /
      data ae12cs( 11) / +.8588024486 0267195879 6605157593 44 d-8     /
      data ae12cs( 12) / +.1021922662 5855003286 3399695539 11 d-8     /
      data ae12cs( 13) / +.2174913232 3289724542 8213398059 92 d-10    /
      data ae12cs( 14) / -.2209023814 2623144809 5235038117 41 d-10    /
      data ae12cs( 15) / -.6345753354 4928753294 3836222088 01 d-11    /
      data ae12cs( 16) / -.1083774656 6857661115 3405397329 19 d-11    /
      data ae12cs( 17) / -.1190982287 2222586730 2622004402 77 d-12    /
      data ae12cs( 18) / -.2843868238 9265590299 5087660086 61 d-14    /
      data ae12cs( 19) / +.2508032702 6686769668 5871954875 46 d-14    /
      data ae12cs( 20) / +.7872964152 8559842431 5977264212 65 d-15    /
      data ae12cs( 21) / +.1547506634 7785217148 4843346373 29 d-15    /
      data ae12cs( 22) / +.2257532283 1665075055 2726081972 90 d-16    /
      data ae12cs( 23) / +.2223335286 7266608760 2813808366 93 d-17    /
      data ae12cs( 24) / +.1696781956 3544153513 4641946623 99 d-19    /
      data ae12cs( 25) / -.5760831625 5947682105 3100873045 33 d-19    /
      data ae12cs( 26) / -.1759123577 4646878055 6253694088 53 d-19    /
      data ae12cs( 27) / -.3628605637 5103174394 7553286826 66 d-20    /
      data ae12cs( 28) / -.5923556979 7328991652 5581434880 00 d-21    /
      data ae12cs( 29) / -.7603038092 6310191114 4291368959 99 d-22    /
      data ae12cs( 30) / -.6254784352 1711763842 6414284799 99 d-23    /
      data ae12cs( 31) / +.2548336075 9307648606 0376064000 00 d-24    /
      data ae12cs( 32) / +.2559861573 1739857020 1688746666 66 d-24    /
      data ae12cs( 33) / +.7137623935 7899318800 2070528000 00 d-25    /
      data ae12cs( 34) / +.1470375993 9567568181 5789568000 00 d-25    /
      data ae12cs( 35) / +.2510552476 5386733555 1986346666 66 d-26    /
      data ae12cs( 36) / +.3588666638 7790890886 5836373333 33 d-27    /
      data ae12cs( 37) / +.3988603515 6771301763 3177599999 99 d-28    /
      data ae12cs( 38) / +.2176367694 7356220478 8053333333 33 d-29    /
      data ae12cs( 39) / -.4614699848 7618942367 6074666666 66 d-30    /
      data ae12cs( 40) / -.2071351787 7481987707 1530666666 66 d-30    /
      data ae12cs( 41) / -.5189037856 3534371596 9706666666 66 d-31    /
c
c series for e11        on the interval -4.00000e+00 to -1.00000e+00
c                                        with weighted error   8.49e-34
c                                         log weighted error  33.07
c                               significant figures required  34.13
c                                    decimal places required  33.80
c
      data e11 cs(  1) / -.1611346165 5571494025 7206639275 66180 d+2  /
      data e11 cs(  2) / +.7794072778 7426802769 2722458917 41497 d+1  /
      data e11 cs(  3) / -.1955405818 8631419507 1272838128 14491 d+1  /
      data e11 cs(  4) / +.3733729386 6277945611 5171908656 90209 d+0  /
      data e11 cs(  5) / -.5692503191 0929019385 2638922200 51166 d-1  /
      data e11 cs(  6) / +.7211077769 6600918537 8477248126 35813 d-2  /
      data e11 cs(  7) / -.7810490144 9841593997 7151840890 64148 d-3  /
      data e11 cs(  8) / +.7388093356 2621681878 9748813661 77858 d-4  /
      data e11 cs(  9) / -.6202861875 8082045134 3581336079 09712 d-5  /
      data e11 cs( 10) / +.4681600230 3176735524 4058238683 62657 d-6  /
      data e11 cs( 11) / -.3209288853 3298649524 0725530272 28719 d-7  /
      data e11 cs( 12) / +.2015199748 7404533394 8262622130 19548 d-8  /
      data e11 cs( 13) / -.1167368681 6697793105 3562716950 15419 d-9  /
      data e11 cs( 14) / +.6276270667 2039943397 7887483796 15573 d-11 /
      data e11 cs( 15) / -.3148154167 2275441045 2467818023 93600 d-12 /
      data e11 cs( 16) / +.1479904174 4493474210 8944722517 33333 d-13 /
      data e11 cs( 17) / -.6545709158 3979673774 2634015880 53333 d-15 /
      data e11 cs( 18) / +.2733687222 3137291142 5080127487 99999 d-16 /
      data e11 cs( 19) / -.1081352434 9754406876 7217276245 33333 d-17 /
      data e11 cs( 20) / +.4062832804 0434303295 3003485866 66666 d-19 /
      data e11 cs( 21) / -.1453553935 8960455858 9143722666 66666 d-20 /
      data e11 cs( 22) / +.4963274618 1648636830 1984426666 66666 d-22 /
      data e11 cs( 23) / -.1620861269 6636044604 8665600000 00000 d-23 /
      data e11 cs( 24) / +.5072144803 8607422226 4319999999 99999 d-25 /
      data e11 cs( 25) / -.1523581113 3372207813 9733333333 33333 d-26 /
      data e11 cs( 26) / +.4400151125 6103618696 5333333333 33333 d-28 /
      data e11 cs( 27) / -.1223614194 5416231594 6666666666 66666 d-29 /
      data e11 cs( 28) / +.3280921666 1066001066 6666666666 66666 d-31 /
      data e11 cs( 29) / -.8493345226 8306432000 0000000000 00000 d-33 /
c
c series for e12        on the interval -1.00000e+00 to  1.00000e+00
c                                        with weighted error   8.08e-33
c                                         log weighted error  32.09
c                        approx significant figures required  30.4
c                                    decimal places required  32.79
c
      data e12 cs(  1) / -.3739021479 22027951166 869820482 7 d-1      /
      data e12 cs(  2) / +.4272398606 2209577260 4917917652 8 d-1      /
      data e12 cs(  3) / -.1303182079 8497005441 5392055219 726 d+0    /
      data e12 cs(  4) / +.1441912402 4698890734 1095893982 137 d-1    /
      data e12 cs(  5) / -.1346170780 5106802211 6121527983 553 d-2    /
      data e12 cs(  6) / +.1073102925 3063779997 6115850970 073 d-3    /
      data e12 cs(  7) / -.7429999516 1194364961 0283062223 163 d-5    /
      data e12 cs(  8) / +.4537732569 0753713938 6383211511 827 d-6    /
      data e12 cs(  9) / -.2476417211 3906013184 6547423802 912 d-7    /
      data e12 cs( 10) / +.1220765813 7459095370 0228167846 102 d-8    /
      data e12 cs( 11) / -.5485141480 6409239382 1357398028 261 d-10   /
      data e12 cs( 12) / +.2263621421 3007879929 3688162377 002 d-11   /
      data e12 cs( 13) / -.8635897271 6980097940 4172916282 240 d-13   /
      data e12 cs( 14) / +.3062915536 6933299758 1032894881 279 d-14   /
      data e12 cs( 15) / -.1014857188 5594414755 7128906734 933 d-15   /
      data e12 cs( 16) / +.3154821740 3406987754 6855328426 666 d-17   /
      data e12 cs( 17) / -.9236042407 6924095448 4015923200 000 d-19   /
      data e12 cs( 18) / +.2555042679 7081400244 0435029333 333 d-20   /
      data e12 cs( 19) / -.6699128056 8456684721 7882453333 333 d-22   /
      data e12 cs( 20) / +.1669254054 3538731943 1987199999 999 d-23   /
      data e12 cs( 21) / -.3962549251 8437964185 6000000000 000 d-25   /
      data e12 cs( 22) / +.8981358965 9851133201 0666666666 666 d-27   /
      data e12 cs( 23) / -.1947633669 9301643332 2666666666 666 d-28   /
      data e12 cs( 24) / +.4048360190 2463003306 6666666666 666 d-30   /
      data e12 cs( 25) / -.8079815676 9984512000 0000000000 000 d-32   /
c
c series for ae13       on the interval  2.50000e-01 to  1.00000e+00
c                                        with weighted error   6.65e-32
c                                         log weighted error  31.18
c                               significant figures required  30.69
c                                    decimal places required  32.03
c
      data ae13cs(  1) / -.6057732466 4060345999 3193827377 47 d+0     /
      data ae13cs(  2) / -.1125352434 8366090030 6497688527 18 d+0     /
      data ae13cs(  3) / +.1343226624 7902779492 4878593294 14 d-1     /
      data ae13cs(  4) / -.1926845187 3811457249 2468389913 03 d-2     /
      data ae13cs(  5) / +.3091183377 2060318335 5867374753 68 d-3     /
      data ae13cs(  6) / -.5356413212 9618418776 3935597951 47 d-4     /
      data ae13cs(  7) / +.9827812880 2474923952 4918827172 37 d-5     /
      data ae13cs(  8) / -.1885368984 9165182826 9028919389 10 d-5     /
      data ae13cs(  9) / +.3749431935 6894735406 9640421905 31 d-6     /
      data ae13cs( 10) / -.7682345587 0552639273 7334656805 56 d-7     /
      data ae13cs( 11) / +.1614327056 7198777552 9563000608 68 d-7     /
      data ae13cs( 12) / -.3466802211 4907354566 3090602260 27 d-8     /
      data ae13cs( 13) / +.7587542091 9036277572 8897470541 14 d-9     /
      data ae13cs( 14) / -.1688643332 9881412573 5145266367 03 d-9     /
      data ae13cs( 15) / +.3814570674 9552265682 8042509272 72 d-10    /
      data ae13cs( 16) / -.8733026632 4446292706 8517182723 34 d-11    /
      data ae13cs( 17) / +.2023672864 5867960961 7943110643 30 d-11    /
      data ae13cs( 18) / -.4741328303 9555834655 2103408201 60 d-12    /
      data ae13cs( 19) / +.1122117204 8389864324 7317999289 20 d-12    /
      data ae13cs( 20) / -.2680422543 4840309912 8268090933 95 d-13    /
      data ae13cs( 21) / +.6457851441 7716530343 5803690672 12 d-14    /
      data ae13cs( 22) / -.1568276050 1666478830 3057028491 94 d-14    /
      data ae13cs( 23) / +.3836786539 9315404861 8215164414 08 d-15    /
      data ae13cs( 24) / -.9451717302 7579130478 8710489325 56 d-16    /
      data ae13cs( 25) / +.2343481228 8949573293 8966664391 33 d-16    /
      data ae13cs( 26) / -.5845866158 0214714576 1231944198 82 d-17    /
      data ae13cs( 27) / +.1466622986 7947778605 8736174191 95 d-17    /
      data ae13cs( 28) / -.3699392347 6444472706 5925382744 74 d-18    /
      data ae13cs( 29) / +.9379015993 6721242136 0142918178 13 d-19    /
      data ae13cs( 30) / -.2389367322 1937873136 3082240873 81 d-19    /
      data ae13cs( 31) / +.6115062462 9497608051 9342238378 66 d-20    /
      data ae13cs( 32) / -.1571858532 7554025507 7198532881 06 d-20    /
      data ae13cs( 33) / +.4057238728 5585397769 5192944913 06 d-21    /
      data ae13cs( 34) / -.1051402655 4738034990 5663671227 73 d-21    /
      data ae13cs( 35) / +.2734966493 0638667785 8060031317 33 d-22    /
      data ae13cs( 36) / -.7140160408 0205796099 3555742719 99 d-23    /
      data ae13cs( 37) / +.1870555243 2235079986 7569242111 99 d-23    /
      data ae13cs( 38) / -.4916746816 6870480520 4780209493 33 d-24    /
      data ae13cs( 39) / +.1296498811 9684031730 9160871253 33 d-24    /
      data ae13cs( 40) / -.3429251568 8362864461 6239404373 33 d-25    /
      data ae13cs( 41) / +.9097224164 3887034329 1048209066 66 d-26    /
      data ae13cs( 42) / -.2420211231 4316856489 9348479999 99 d-26    /
      data ae13cs( 43) / +.6456361293 4639510757 6704750933 33 d-27    /
      data ae13cs( 44) / -.1726913273 5340541122 3159876266 66 d-27    /
      data ae13cs( 45) / +.4630861165 9151500715 1942314666 66 d-28    /
      data ae13cs( 46) / -.1244870363 7214131241 7551701333 33 d-28    /
      data ae13cs( 47) / +.3354457409 0520678532 9070079999 99 d-29    /
      data ae13cs( 48) / -.9059886852 1070774437 5439359999 99 d-30    /
      data ae13cs( 49) / +.2452414705 1474238587 2732160000 00 d-30    /
      data ae13cs( 50) / -.6652817873 3552062817 1079679999 99 d-31    /
c
c series for ae14       on the interval  0.          to  2.50000e-01
c                                        with weighted error   5.07e-32
c                                         log weighted error  31.30
c                               significant figures required  30.40
c                                    decimal places required  32.20
c
      data ae14cs(  1) / -.1892918000 7530168254 9567994282 0 d+0      /
      data ae14cs(  2) / -.8648117855 2598714899 6881705682 4 d-1      /
      data ae14cs(  3) / +.7224101543 7465947470 2151483918 4 d-2      /
      data ae14cs(  4) / -.8097559457 5573861971 5965561018 1 d-3      /
      data ae14cs(  5) / +.1099913443 2661388671 7925115700 2 d-3      /
      data ae14cs(  6) / -.1717332998 9377673714 9535881448 7 d-4      /
      data ae14cs(  7) / +.2985627514 4792833228 2534249500 3 d-5      /
      data ae14cs(  8) / -.5659649145 7719300565 6016726715 5 d-6      /
      data ae14cs(  9) / +.1152680839 7141400192 2658350166 3 d-6      /
      data ae14cs( 10) / -.2495030440 2693382288 4212876506 5 d-7      /
      data ae14cs( 11) / +.5692324201 8337543670 3937036814 0 d-8      /
      data ae14cs( 12) / -.1359957664 8056003384 9003093917 6 d-8      /
      data ae14cs( 13) / +.3384662888 7608845901 8451292585 9 d-9      /
      data ae14cs( 14) / -.8737853904 4746819523 5084931658 0 d-10     /
      data ae14cs( 15) / +.2331588663 2226597186 1261340047 0 d-10     /
      data ae14cs( 16) / -.6411481049 2137859697 5316519632 6 d-11     /
      data ae14cs( 17) / +.1812246980 2048164333 8435948468 2 d-11     /
      data ae14cs( 18) / -.5253831761 5584606888 1940384046 6 d-12     /
      data ae14cs( 19) / +.1559218272 5919256988 5502860982 5 d-12     /
      data ae14cs( 20) / -.4729168297 0803987184 7642936946 6 d-13     /
      data ae14cs( 21) / +.1463761864 3932435020 7619949380 8 d-13     /
      data ae14cs( 22) / -.4617388988 7129241022 3217362360 4 d-14     /
      data ae14cs( 23) / +.1482710348 2893693237 8923966037 1 d-14     /
      data ae14cs( 24) / -.4841672496 2392291469 7316573441 7 d-15     /
      data ae14cs( 25) / +.1606215575 7002904081 1657196618 8 d-15     /
      data ae14cs( 26) / -.5408917538 9571709478 9502378425 2 d-16     /
      data ae14cs( 27) / +.1847470159 3468978813 7023140231 0 d-16     /
      data ae14cs( 28) / -.6395830792 7590944705 0061042505 0 d-17     /
      data ae14cs( 29) / +.2242780721 6997594572 5023327617 0 d-17     /
      data ae14cs( 30) / -.7961369173 9839475527 4455530864 6 d-18     /
      data ae14cs( 31) / +.2859308111 5401974598 0861992927 2 d-18     /
      data ae14cs( 32) / -.1038450244 7011371459 0069713744 6 d-18     /
      data ae14cs( 33) / +.3812040607 0979757808 6684100831 9 d-19     /
      data ae14cs( 34) / -.1413795417 7172007687 1756272369 6 d-19     /
      data ae14cs( 35) / +.5295367865 1827409583 0544259481 5 d-20     /
      data ae14cs( 36) / -.2002264245 0268259021 3721113143 9 d-20     /
      data ae14cs( 37) / +.7640262751 2751960147 3684861091 8 d-21     /
      data ae14cs( 38) / -.2941119006 8687878833 1126352336 2 d-21     /
      data ae14cs( 39) / +.1141823539 0789271930 3769148358 6 d-21     /
      data ae14cs( 40) / -.4469308475 9552984252 4702071848 9 d-22     /
      data ae14cs( 41) / +.1763262410 5717507706 3049140852 0 d-22     /
      data ae14cs( 42) / -.7009968187 9259023563 5151826234 0 d-23     /
      data ae14cs( 43) / +.2807573556 5583789222 8775750751 5 d-23     /
      data ae14cs( 44) / -.1132560944 9810864321 4188889156 2 d-23     /
      data ae14cs( 45) / +.4600574684 3750179461 5676423372 7 d-24     /
      data ae14cs( 46) / -.1881448598 9761334598 6460914810 8 d-24     /
      data ae14cs( 47) / +.7744916111 5077308454 4432847803 7 d-25     /
      data ae14cs( 48) / -.3208512760 5853689267 0270382626 1 d-25     /
      data ae14cs( 49) / +.1337445542 9108397606 1993042138 4 d-25     /
      data ae14cs( 50) / -.5608671881 8022170488 9477173521 0 d-26     /
      data ae14cs( 51) / +.2365839716 5285374837 1006947327 9 d-26     /
      data ae14cs( 52) / -.1003656195 0253053340 6583452685 6 d-26     /
      data ae14cs( 53) / +.4281490878 0941611312 8664255692 7 d-27     /
      data ae14cs( 54) / -.1836345261 8153181996 9132695825 0 d-27     /
      data ae14cs( 55) / +.7917798231 3495400000 9746867814 4 d-28     /
      data ae14cs( 56) / -.3431542358 7422203610 2501577523 1 d-28     /
      data ae14cs( 57) / +.1494705493 8971032374 7506600891 7 d-28     /
      data ae14cs( 58) / -.6542620279 8657054397 3904242005 3 d-29     /
      data ae14cs( 59) / +.2877581395 1991711143 4048735368 5 d-29     /
      data ae14cs( 60) / -.1271557211 7960247110 2798120004 2 d-29     /
      data ae14cs( 61) / +.5644615555 6487225223 8804462250 6 d-30     /
      data ae14cs( 62) / -.2516994994 2840951060 8061683029 3 d-30     /
      data ae14cs( 63) / +.1127259818 9275102063 7036880418 1 d-30     /
      data ae14cs( 64) / -.5069814875 8004608555 6258471936 0 d-31     /
c
      data ntae10, ntae11, ntae12, nte11, nte12, ntae13, ntae14 / 7*0 /
      data xmax / 0.d0 /
c
      if (ntae10.ne.0) go to 10
      eta = 0.1*sngl(d1mach(3))
      ntae10 = initds (ae10cs, 50, eta)
      ntae11 = initds (ae11cs, 60, eta)
      ntae12 = initds (ae12cs, 41, eta)
      nte11 = initds (e11cs, 29, eta)
      nte12 = initds (e12cs, 25, eta)
      ntae13 = initds (ae13cs, 50, eta)
      ntae14 = initds (ae14cs, 64, eta)
c
      xmax = -dlog (d1mach(1))
      xmax = xmax - dlog(xmax)
c
 10   if (x.gt.(-1.d0)) go to 50
      if (x.gt.(-32.d0)) go to 20
      de1 = dexp(-x)/x * (1.d0 + dcsevl (64.d0/x+1.d0, ae10cs, ntae10))
      return
c
 20   if (x.gt.(-8.d0)) go to 30
      de1 = dexp(-x)/x * (1.d0 + dcsevl ((64.d0/x+5.d0)/3.d0, ae11cs,
     1  ntae11))
      return
c
 30   if (x.gt.(-4.d0)) go to 40
      de1 = dexp(-x)/x * (1.d0 + dcsevl (16.d0/x+3.d0, ae12cs, ntae12))
      return
c
 40   de1 = -dlog(-x) + dcsevl ((2.d0*x+5.d0)/3.d0, e11cs, nte11)
      return
c
 50   if (x.gt.1.0d0) go to 60
      if (x.eq.0.d0) call seteru (14hde1     x is 0, 14, 2, 2)
      de1 = (-dlog(dabs(x)) - 0.6875d0 + x)  + dcsevl (x, e12cs, nte12)
      return
c
 60   if (x.gt.4.0d0) go to 70
      de1 = dexp(-x)/x * (1.d0 + dcsevl ((8.d0/x-5.d0)/3.d0, ae13cs,
     1  ntae13))
      return
c
 70   if (x.gt.xmax) go to 80
      de1 = dexp(-x)/x * (1.d0 + dcsevl (8.d0/x-1.d0, ae14cs, ntae14))
      return
c
 80   call seteru (30hde1     x so big e1 underflows, 30, 1, 0)
      de1 = 0.d0
      return
c
      end