REAL             FUNCTION SEI (XIN)
c Copyright (c) 1996 California Institute of Technology, Pasadena, CA.
c ALL RIGHTS RESERVED.
c Based on Government Sponsored Research NAS7-03001.
c>> 1996-04-27 SEI Krogh  Changes to use .C. and C%%.
c>> 1996-03-30 SEI Krogh  Added external statements.
C>> 1995-11-28 SEI Krogh  GO TO's => blodk IF's, removed multiple entry.
C>> 1995-11-14 SEI Krogh  Changes to simplify conversion to C.
C>> 1995-11-03 SEI Krogh  Removed blanks in numbers for C conversion.
C>> 1994-10-20 SEI Krogh  Changes to use M77CON
C>> 1994-04-20 SEI CLL  Edited type stmts to make DP & SP files similar.
C>> 1992-03-13 SEI FTK  Removed implicit statements.
C>> 1991-01-14 CLL SE1 Changed to generic name ABS
C>> 1990-11-29 CLL
C>> 1985-08-02 SEI    Lawson  Initial code.
C
C     JULY 1977 EDITION. W. FULLERTON, C3,
C     LOS ALAMOS SCIENTIFIC LAB.
C
C     REORGANIZATION OF FULLERTON'S SEI & SE1,
C     C.L. LAWSON & S.CHAN, JUNE 1983, JPL.
C
C     ----------------------------------------------------------------
C     MACHINE DEPENDENT VALUES ARE SET ON THE FIRST ENTRY
C     TO THIS CODE. EXAMPLES OF THESE VALUES FOLLOW:
C
C     SYSTEM         NTAE10   NTAE11   NTAE12   NTAE13   NTAE14
C     ------         ------   ------   ------   ------   ------
C     UNIVAC S.P.      **       16       14       13       12
C     UNIVAC D.P.      21       34       24       29       32
C
C     SYSTEM         NTE11    NTE12    XMAX     XMIN     DELTA
C     ------         -----    -----    ----     ----     -----
C     UNIVAC S.P.     13       10      84.8    -92.41    4.48
C     UNIVAC D.P.     19       17     703.2    -714.9    6.56
C
C     SYSTEM         EARG1    EARG2
C     ------         -----    -----
C     UNIVAC S.P.    89.42    88.03
C     UNIVAC D.P.    710.5    709.1
C
C     ----------------------------------------------------------------
c--S replaces "?": ?E1, ?EI, ?INITS, ?CSEVL, ?ERM1, ?ERV1
c     Also uses ERMSG
c     ------------------------------------------------------------------
      EXTERNAL SE1
      REAL             XIN
      REAL             SE1
c
      SEI = -SE1(-XIN)
      return
      end
c
      REAL             FUNCTION SE1 (XIN)
      EXTERNAL R1MACH, SCSEVL
      INTEGER NTAE10,NTAE11,NTAE12,NTAE13,NTAE14,NTE11,NTE12
      REAL             R1MACH,SCSEVL
      REAL             AE10CS(50),AE11CS(60),AE12CS(41),AE13CS(50)
      REAL             AE14CS(64),BIGNUM,DELTA
      REAL             E11CS(29),E12CS(25),EARG1,EARG2
      REAL             ETA,FAC,X,XIN,XMAX,XMIN,ZERO
C
      CHARACTER * 33 MSG1
      CHARACTER * 38 MSG2
      CHARACTER * 34 MSG3
      CHARACTER * 10 NAME
C
      SAVE DELTA,EARG1,EARG2,NTAE11,NTAE12,NTE11,NTE12,
     *     NTAE13,NTAE14,XMAX,XMIN,NTAE10
C
      DATA NTAE10, NTAE11, NTAE12, NTE11, NTE12, NTAE13, NTAE14 / 7*0 /
      DATA XMAX,ZERO / 2*0.E0 /
      DATA FAC / 0.999E0 /
C
      DATA MSG1 / '|X| SO LARGE SE1 or SEI OVERFLOWS' /
C
      DATA MSG2 / 'X = ZERO, SE1(0) or SEI(0) not defined' /
C
      DATA MSG3 / '|X| SO LARGE SE1 or SEI UNDERFLOWS' /
C
      DATA NAME / 'SE1 or SEI' /
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
c++ Save data by elements if ~.C.
      DATA AE10CS(  1) / +.3284394579616699087873844201881E-1      /
      DATA AE10CS(  2) / -.1669920452031362851476184343387E-1      /
      DATA AE10CS(  3) / +.2845284724361346807424899853252E-3      /
      DATA AE10CS(  4) / -.7563944358516206489487866938533E-5      /
      DATA AE10CS(  5) / +.2798971289450859157504843180879E-6      /
      DATA AE10CS(  6) / -.1357901828534531069525563926255E-7      /
      DATA AE10CS(  7) / +.8343596202040469255856102904906E-9      /
      DATA AE10CS(  8) / -.6370971727640248438275242988532E-10     /
      DATA AE10CS(  9) / +.6007247608811861235760831561584E-11     /
      DATA AE10CS( 10) / -.7022876174679773590750626150088E-12     /
      DATA AE10CS( 11) / +.1018302673703687693096652346883E-12     /
      DATA AE10CS( 12) / -.1761812903430880040406309966422E-13     /
      DATA AE10CS( 13) / +.3250828614235360694244030353877E-14     /
      DATA AE10CS( 14) / -.5071770025505818678824872259044E-15     /
      DATA AE10CS( 15) / +.1665177387043294298172486084156E-16     /
      DATA AE10CS( 16) / +.3166753890797514400677003536555E-16     /
      DATA AE10CS( 17) / -.1588403763664141515133118343538E-16     /
      DATA AE10CS( 18) / +.4175513256138018833003034618484E-17     /
      DATA AE10CS( 19) / -.2892347749707141906710714478852E-18     /
      DATA AE10CS( 20) / -.2800625903396608103506340589669E-18     /
      DATA AE10CS( 21) / +.1322938639539270903707580023781E-18     /
      DATA AE10CS( 22) / -.1804447444177301627283887833557E-19     /
      DATA AE10CS( 23) / -.7905384086522616076291644817604E-20     /
      DATA AE10CS( 24) / +.4435711366369570103946235838027E-20     /
      DATA AE10CS( 25) / -.4264103994978120868865309206555E-21     /
      DATA AE10CS( 26) / -.3920101766937117541553713162048E-21     /
      DATA AE10CS( 27) / +.1527378051343994266343752326971E-21     /
      DATA AE10CS( 28) / +.1024849527049372339310308783117E-22     /
      DATA AE10CS( 29) / -.2134907874771433576262711405882E-22     /
      DATA AE10CS( 30) / +.3239139475160028267061694700366E-23     /
      DATA AE10CS( 31) / +.2142183762299889954762643168296E-23     /
      DATA AE10CS( 32) / -.8234609419601018414700348082312E-24     /
      DATA AE10CS( 33) / -.1524652829645809479613694401140E-24     /
      DATA AE10CS( 34) / +.1378208282460639134668480364325E-24     /
      DATA AE10CS( 35) / +.2131311202833947879523224999253E-26     /
      DATA AE10CS( 36) / -.2012649651526484121817466763127E-25     /
      DATA AE10CS( 37) / +.1995535662263358016106311782673E-26     /
      DATA AE10CS( 38) / +.2798995808984003464948686520319E-26     /
      DATA AE10CS( 39) / -.5534511845389626637640819277823E-27     /
      DATA AE10CS( 40) / -.3884995396159968861682544026146E-27     /
      DATA AE10CS( 41) / +.1121304434507359382850680354679E-27     /
      DATA AE10CS( 42) / +.5566568152423740948256563833514E-28     /
      DATA AE10CS( 43) / -.2045482929810499700448533938176E-28     /
      DATA AE10CS( 44) / -.8453813992712336233411457493674E-29     /
      DATA AE10CS( 45) / +.3565758433431291562816111116287E-29     /
      DATA AE10CS( 46) / +.1383653872125634705539949098871E-29     /
      DATA AE10CS( 47) / -.6062167864451372436584533764778E-30     /
      DATA AE10CS( 48) / -.2447198043989313267437655119189E-30     /
      DATA AE10CS( 49) / +.1006850640933998348011548180480E-30     /
      DATA AE10CS( 50) / +.4623685555014869015664341461674E-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
c++ Save data by elements if ~.C.
      DATA AE11CS(  1) / +.20263150647078889499401236517381E+0     /
      DATA AE11CS(  2) / -.73655140991203130439536898728034E-1     /
      DATA AE11CS(  3) / +.63909349118361915862753283840020E-2     /
      DATA AE11CS(  4) / -.60797252705247911780653153363999E-3     /
      DATA AE11CS(  5) / -.73706498620176629330681411493484E-4     /
      DATA AE11CS(  6) / +.48732857449450183453464992488076E-4     /
      DATA AE11CS(  7) / -.23837064840448290766588489460235E-5     /
      DATA AE11CS(  8) / -.30518612628561521027027332246121E-5     /
      DATA AE11CS(  9) / +.17050331572564559009688032992907E-6     /
      DATA AE11CS( 10) / +.23834204527487747258601598136403E-6     /
      DATA AE11CS( 11) / +.10781772556163166562596872364020E-7     /
      DATA AE11CS( 12) / -.17955692847399102653642691446599E-7     /
      DATA AE11CS( 13) / -.41284072341950457727912394640436E-8     /
      DATA AE11CS( 14) / +.68622148588631968618346844526664E-9     /
      DATA AE11CS( 15) / +.53130183120506356147602009675961E-9     /
      DATA AE11CS( 16) / +.78796880261490694831305022893515E-10    /
      DATA AE11CS( 17) / -.26261762329356522290341675271232E-10    /
      DATA AE11CS( 18) / -.15483687636308261963125756294100E-10    /
      DATA AE11CS( 19) / -.25818962377261390492802405122591E-11    /
      DATA AE11CS( 20) / +.59542879191591072658903529959352E-12    /
      DATA AE11CS( 21) / +.46451400387681525833784919321405E-12    /
      DATA AE11CS( 22) / +.11557855023255861496288006203731E-12    /
      DATA AE11CS( 23) / -.10475236870835799012317547189670E-14    /
      DATA AE11CS( 24) / -.11896653502709004368104489260929E-13    /
      DATA AE11CS( 25) / -.47749077490261778752643019349950E-14    /
      DATA AE11CS( 26) / -.81077649615772777976249734754135E-15    /
      DATA AE11CS( 27) / +.13435569250031554199376987998178E-15    /
      DATA AE11CS( 28) / +.14134530022913106260248873881287E-15    /
      DATA AE11CS( 29) / +.49451592573953173115520663232883E-16    /
      DATA AE11CS( 30) / +.79884048480080665648858587399367E-17    /
      DATA AE11CS( 31) / -.14008632188089809829248711935393E-17    /
      DATA AE11CS( 32) / -.14814246958417372107722804001680E-17    /
      DATA AE11CS( 33) / -.55826173646025601904010693937113E-18    /
      DATA AE11CS( 34) / -.11442074542191647264783072544598E-18    /
      DATA AE11CS( 35) / +.25371823879566853500524018479923E-20    /
      DATA AE11CS( 36) / +.13205328154805359813278863389097E-19    /
      DATA AE11CS( 37) / +.62930261081586809166287426789485E-20    /
      DATA AE11CS( 38) / +.17688270424882713734999261332548E-20    /
      DATA AE11CS( 39) / +.23266187985146045209674296887432E-21    /
      DATA AE11CS( 40) / -.67803060811125233043773831844113E-22    /
      DATA AE11CS( 41) / -.59440876959676373802874150531891E-22    /
      DATA AE11CS( 42) / -.23618214531184415968532592503466E-22    /
      DATA AE11CS( 43) / -.60214499724601478214168478744576E-23    /
      DATA AE11CS( 44) / -.65517906474348299071370444144639E-24    /
      DATA AE11CS( 45) / +.29388755297497724587042038699349E-24    /
      DATA AE11CS( 46) / +.22601606200642115173215728758510E-24    /
      DATA AE11CS( 47) / +.89534369245958628745091206873087E-25    /
      DATA AE11CS( 48) / +.24015923471098457555772067457706E-25    /
      DATA AE11CS( 49) / +.34118376888907172955666423043413E-26    /
      DATA AE11CS( 50) / -.71617071694630342052355013345279E-27    /
      DATA AE11CS( 51) / -.75620390659281725157928651980799E-27    /
      DATA AE11CS( 52) / -.33774612157467324637952920780800E-27    /
      DATA AE11CS( 53) / -.10479325703300941711526430332245E-27    /
      DATA AE11CS( 54) / -.21654550252170342240854880201386E-28    /
      DATA AE11CS( 55) / -.75297125745288269994689298432000E-30    /
      DATA AE11CS( 56) / +.19103179392798935768638084000426E-29    /
      DATA AE11CS( 57) / +.11492104966530338547790728833706E-29    /
      DATA AE11CS( 58) / +.43896970582661751514410359193600E-30    /
      DATA AE11CS( 59) / +.12320883239205686471647157725866E-30    /
      DATA AE11CS( 60) / +.22220174457553175317538581162666E-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
c++ Save data by elements if ~.C.
      DATA AE12CS(  1) / +.63629589796747038767129887806803E+0     /
      DATA AE12CS(  2) / -.13081168675067634385812671121135E+0     /
      DATA AE12CS(  3) / -.84367410213053930014487662129752E-2     /
      DATA AE12CS(  4) / +.26568491531006685413029428068906E-2     /
      DATA AE12CS(  5) / +.32822721781658133778792170142517E-3     /
      DATA AE12CS(  6) / -.23783447771430248269579807851050E-4     /
      DATA AE12CS(  7) / -.11439804308100055514447076797047E-4     /
      DATA AE12CS(  8) / -.14405943433238338455239717699323E-5     /
      DATA AE12CS(  9) / +.52415956651148829963772818061664E-8     /
      DATA AE12CS( 10) / +.38407306407844323480979203059716E-7     /
      DATA AE12CS( 11) / +.85880244860267195879660515759344E-8     /
      DATA AE12CS( 12) / +.10219226625855003286339969553911E-8     /
      DATA AE12CS( 13) / +.21749132323289724542821339805992E-10    /
      DATA AE12CS( 14) / -.22090238142623144809523503811741E-10    /
      DATA AE12CS( 15) / -.63457533544928753294383622208801E-11    /
      DATA AE12CS( 16) / -.10837746566857661115340539732919E-11    /
      DATA AE12CS( 17) / -.11909822872222586730262200440277E-12    /
      DATA AE12CS( 18) / -.28438682389265590299508766008661E-14    /
      DATA AE12CS( 19) / +.25080327026686769668587195487546E-14    /
      DATA AE12CS( 20) / +.78729641528559842431597726421265E-15    /
      DATA AE12CS( 21) / +.15475066347785217148484334637329E-15    /
      DATA AE12CS( 22) / +.22575322831665075055272608197290E-16    /
      DATA AE12CS( 23) / +.22233352867266608760281380836693E-17    /
      DATA AE12CS( 24) / +.16967819563544153513464194662399E-19    /
      DATA AE12CS( 25) / -.57608316255947682105310087304533E-19    /
      DATA AE12CS( 26) / -.17591235774646878055625369408853E-19    /
      DATA AE12CS( 27) / -.36286056375103174394755328682666E-20    /
      DATA AE12CS( 28) / -.59235569797328991652558143488000E-21    /
      DATA AE12CS( 29) / -.76030380926310191114429136895999E-22    /
      DATA AE12CS( 30) / -.62547843521711763842641428479999E-23    /
      DATA AE12CS( 31) / +.25483360759307648606037606400000E-24    /
      DATA AE12CS( 32) / +.25598615731739857020168874666666E-24    /
      DATA AE12CS( 33) / +.71376239357899318800207052800000E-25    /
      DATA AE12CS( 34) / +.14703759939567568181578956800000E-25    /
      DATA AE12CS( 35) / +.25105524765386733555198634666666E-26    /
      DATA AE12CS( 36) / +.35886666387790890886583637333333E-27    /
      DATA AE12CS( 37) / +.39886035156771301763317759999999E-28    /
      DATA AE12CS( 38) / +.21763676947356220478805333333333E-29    /
      DATA AE12CS( 39) / -.46146998487618942367607466666666E-30    /
      DATA AE12CS( 40) / -.20713517877481987707153066666666E-30    /
      DATA AE12CS( 41) / -.51890378563534371596970666666666E-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
c++ Save data by elements if ~.C.
      DATA E11CS(  1) / -.16113461655571494025720663927566180E+2  /
      DATA E11CS(  2) / +.77940727787426802769272245891741497E+1  /
      DATA E11CS(  3) / -.19554058188631419507127283812814491E+1  /
      DATA E11CS(  4) / +.37337293866277945611517190865690209E+0  /
      DATA E11CS(  5) / -.56925031910929019385263892220051166E-1  /
      DATA E11CS(  6) / +.72110777696600918537847724812635813E-2  /
      DATA E11CS(  7) / -.78104901449841593997715184089064148E-3  /
      DATA E11CS(  8) / +.73880933562621681878974881366177858E-4  /
      DATA E11CS(  9) / -.62028618758082045134358133607909712E-5  /
      DATA E11CS( 10) / +.46816002303176735524405823868362657E-6  /
      DATA E11CS( 11) / -.32092888533298649524072553027228719E-7  /
      DATA E11CS( 12) / +.20151997487404533394826262213019548E-8  /
      DATA E11CS( 13) / -.11673686816697793105356271695015419E-9  /
      DATA E11CS( 14) / +.62762706672039943397788748379615573E-11 /
      DATA E11CS( 15) / -.31481541672275441045246781802393600E-12 /
      DATA E11CS( 16) / +.14799041744493474210894472251733333E-13 /
      DATA E11CS( 17) / -.65457091583979673774263401588053333E-15 /
      DATA E11CS( 18) / +.27336872223137291142508012748799999E-16 /
      DATA E11CS( 19) / -.10813524349754406876721727624533333E-17 /
      DATA E11CS( 20) / +.40628328040434303295300348586666666E-19 /
      DATA E11CS( 21) / -.14535539358960455858914372266666666E-20 /
      DATA E11CS( 22) / +.49632746181648636830198442666666666E-22 /
      DATA E11CS( 23) / -.16208612696636044604866560000000000E-23 /
      DATA E11CS( 24) / +.50721448038607422226431999999999999E-25 /
      DATA E11CS( 25) / -.15235811133372207813973333333333333E-26 /
      DATA E11CS( 26) / +.44001511256103618696533333333333333E-28 /
      DATA E11CS( 27) / -.12236141945416231594666666666666666E-29 /
      DATA E11CS( 28) / +.32809216661066001066666666666666666E-31 /
      DATA E11CS( 29) / -.84933452268306432000000000000000000E-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
c++ Save data by elements if ~.C.
      DATA E12CS(  1) / -.3739021479220279511668698204827E-1      /
      DATA E12CS(  2) / +.4272398606220957726049179176528E-1      /
      DATA E12CS(  3) / -.130318207984970054415392055219726E+0    /
      DATA E12CS(  4) / +.144191240246988907341095893982137E-1    /
      DATA E12CS(  5) / -.134617078051068022116121527983553E-2    /
      DATA E12CS(  6) / +.107310292530637799976115850970073E-3    /
      DATA E12CS(  7) / -.742999951611943649610283062223163E-5    /
      DATA E12CS(  8) / +.453773256907537139386383211511827E-6    /
      DATA E12CS(  9) / -.247641721139060131846547423802912E-7    /
      DATA E12CS( 10) / +.122076581374590953700228167846102E-8    /
      DATA E12CS( 11) / -.548514148064092393821357398028261E-10   /
      DATA E12CS( 12) / +.226362142130078799293688162377002E-11   /
      DATA E12CS( 13) / -.863589727169800979404172916282240E-13   /
      DATA E12CS( 14) / +.306291553669332997581032894881279E-14   /
      DATA E12CS( 15) / -.101485718855944147557128906734933E-15   /
      DATA E12CS( 16) / +.315482174034069877546855328426666E-17   /
      DATA E12CS( 17) / -.923604240769240954484015923200000E-19   /
      DATA E12CS( 18) / +.255504267970814002440435029333333E-20   /
      DATA E12CS( 19) / -.669912805684566847217882453333333E-22   /
      DATA E12CS( 20) / +.166925405435387319431987199999999E-23   /
      DATA E12CS( 21) / -.396254925184379641856000000000000E-25   /
      DATA E12CS( 22) / +.898135896598511332010666666666666E-27   /
      DATA E12CS( 23) / -.194763366993016433322666666666666E-28   /
      DATA E12CS( 24) / +.404836019024630033066666666666666E-30   /
      DATA E12CS( 25) / -.807981567699845120000000000000000E-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
c++ Save data by elements if ~.C.
      DATA AE13CS(  1) / -.60577324664060345999319382737747E+0     /
      DATA AE13CS(  2) / -.11253524348366090030649768852718E+0     /
      DATA AE13CS(  3) / +.13432266247902779492487859329414E-1     /
      DATA AE13CS(  4) / -.19268451873811457249246838991303E-2     /
      DATA AE13CS(  5) / +.30911833772060318335586737475368E-3     /
      DATA AE13CS(  6) / -.53564132129618418776393559795147E-4     /
      DATA AE13CS(  7) / +.98278128802474923952491882717237E-5     /
      DATA AE13CS(  8) / -.18853689849165182826902891938910E-5     /
      DATA AE13CS(  9) / +.37494319356894735406964042190531E-6     /
      DATA AE13CS( 10) / -.76823455870552639273733465680556E-7     /
      DATA AE13CS( 11) / +.16143270567198777552956300060868E-7     /
      DATA AE13CS( 12) / -.34668022114907354566309060226027E-8     /
      DATA AE13CS( 13) / +.75875420919036277572889747054114E-9     /
      DATA AE13CS( 14) / -.16886433329881412573514526636703E-9     /
      DATA AE13CS( 15) / +.38145706749552265682804250927272E-10    /
      DATA AE13CS( 16) / -.87330266324446292706851718272334E-11    /
      DATA AE13CS( 17) / +.20236728645867960961794311064330E-11    /
      DATA AE13CS( 18) / -.47413283039555834655210340820160E-12    /
      DATA AE13CS( 19) / +.11221172048389864324731799928920E-12    /
      DATA AE13CS( 20) / -.26804225434840309912826809093395E-13    /
      DATA AE13CS( 21) / +.64578514417716530343580369067212E-14    /
      DATA AE13CS( 22) / -.15682760501666478830305702849194E-14    /
      DATA AE13CS( 23) / +.38367865399315404861821516441408E-15    /
      DATA AE13CS( 24) / -.94517173027579130478871048932556E-16    /
      DATA AE13CS( 25) / +.23434812288949573293896666439133E-16    /
      DATA AE13CS( 26) / -.58458661580214714576123194419882E-17    /
      DATA AE13CS( 27) / +.14666229867947778605873617419195E-17    /
      DATA AE13CS( 28) / -.36993923476444472706592538274474E-18    /
      DATA AE13CS( 29) / +.93790159936721242136014291817813E-19    /
      DATA AE13CS( 30) / -.23893673221937873136308224087381E-19    /
      DATA AE13CS( 31) / +.61150624629497608051934223837866E-20    /
      DATA AE13CS( 32) / -.15718585327554025507719853288106E-20    /
      DATA AE13CS( 33) / +.40572387285585397769519294491306E-21    /
      DATA AE13CS( 34) / -.10514026554738034990566367122773E-21    /
      DATA AE13CS( 35) / +.27349664930638667785806003131733E-22    /
      DATA AE13CS( 36) / -.71401604080205796099355574271999E-23    /
      DATA AE13CS( 37) / +.18705552432235079986756924211199E-23    /
      DATA AE13CS( 38) / -.49167468166870480520478020949333E-24    /
      DATA AE13CS( 39) / +.12964988119684031730916087125333E-24    /
      DATA AE13CS( 40) / -.34292515688362864461623940437333E-25    /
      DATA AE13CS( 41) / +.90972241643887034329104820906666E-26    /
      DATA AE13CS( 42) / -.24202112314316856489934847999999E-26    /
      DATA AE13CS( 43) / +.64563612934639510757670475093333E-27    /
      DATA AE13CS( 44) / -.17269132735340541122315987626666E-27    /
      DATA AE13CS( 45) / +.46308611659151500715194231466666E-28    /
      DATA AE13CS( 46) / -.12448703637214131241755170133333E-28    /
      DATA AE13CS( 47) / +.33544574090520678532907007999999E-29    /
      DATA AE13CS( 48) / -.90598868521070774437543935999999E-30    /
      DATA AE13CS( 49) / +.24524147051474238587273216000000E-30    /
      DATA AE13CS( 50) / -.66528178733552062817107967999999E-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
c++ Save data by elements if ~.C.
      DATA AE14CS(  1) / -.1892918000753016825495679942820E+0      /
      DATA AE14CS(  2) / -.8648117855259871489968817056824E-1      /
      DATA AE14CS(  3) / +.7224101543746594747021514839184E-2      /
      DATA AE14CS(  4) / -.8097559457557386197159655610181E-3      /
      DATA AE14CS(  5) / +.1099913443266138867179251157002E-3      /
      DATA AE14CS(  6) / -.1717332998937767371495358814487E-4      /
      DATA AE14CS(  7) / +.2985627514479283322825342495003E-5      /
      DATA AE14CS(  8) / -.5659649145771930056560167267155E-6      /
      DATA AE14CS(  9) / +.1152680839714140019226583501663E-6      /
      DATA AE14CS( 10) / -.2495030440269338228842128765065E-7      /
      DATA AE14CS( 11) / +.5692324201833754367039370368140E-8      /
      DATA AE14CS( 12) / -.1359957664805600338490030939176E-8      /
      DATA AE14CS( 13) / +.3384662888760884590184512925859E-9      /
      DATA AE14CS( 14) / -.8737853904474681952350849316580E-10     /
      DATA AE14CS( 15) / +.2331588663222659718612613400470E-10     /
      DATA AE14CS( 16) / -.6411481049213785969753165196326E-11     /
      DATA AE14CS( 17) / +.1812246980204816433384359484682E-11     /
      DATA AE14CS( 18) / -.5253831761558460688819403840466E-12     /
      DATA AE14CS( 19) / +.1559218272591925698855028609825E-12     /
      DATA AE14CS( 20) / -.4729168297080398718476429369466E-13     /
      DATA AE14CS( 21) / +.1463761864393243502076199493808E-13     /
      DATA AE14CS( 22) / -.4617388988712924102232173623604E-14     /
      DATA AE14CS( 23) / +.1482710348289369323789239660371E-14     /
      DATA AE14CS( 24) / -.4841672496239229146973165734417E-15     /
      DATA AE14CS( 25) / +.1606215575700290408116571966188E-15     /
      DATA AE14CS( 26) / -.5408917538957170947895023784252E-16     /
      DATA AE14CS( 27) / +.1847470159346897881370231402310E-16     /
      DATA AE14CS( 28) / -.6395830792759094470500610425050E-17     /
      DATA AE14CS( 29) / +.2242780721699759457250233276170E-17     /
      DATA AE14CS( 30) / -.7961369173983947552744555308646E-18     /
      DATA AE14CS( 31) / +.2859308111540197459808619929272E-18     /
      DATA AE14CS( 32) / -.1038450244701137145900697137446E-18     /
      DATA AE14CS( 33) / +.3812040607097975780866841008319E-19     /
      DATA AE14CS( 34) / -.1413795417717200768717562723696E-19     /
      DATA AE14CS( 35) / +.5295367865182740958305442594815E-20     /
      DATA AE14CS( 36) / -.2002264245026825902137211131439E-20     /
      DATA AE14CS( 37) / +.7640262751275196014736848610918E-21     /
      DATA AE14CS( 38) / -.2941119006868787883311263523362E-21     /
      DATA AE14CS( 39) / +.1141823539078927193037691483586E-21     /
      DATA AE14CS( 40) / -.4469308475955298425247020718489E-22     /
      DATA AE14CS( 41) / +.1763262410571750770630491408520E-22     /
      DATA AE14CS( 42) / -.7009968187925902356351518262340E-23     /
      DATA AE14CS( 43) / +.2807573556558378922287757507515E-23     /
      DATA AE14CS( 44) / -.1132560944981086432141888891562E-23     /
      DATA AE14CS( 45) / +.4600574684375017946156764233727E-24     /
      DATA AE14CS( 46) / -.1881448598976133459864609148108E-24     /
      DATA AE14CS( 47) / +.7744916111507730845444328478037E-25     /
      DATA AE14CS( 48) / -.3208512760585368926702703826261E-25     /
      DATA AE14CS( 49) / +.1337445542910839760619930421384E-25     /
      DATA AE14CS( 50) / -.5608671881802217048894771735210E-26     /
      DATA AE14CS( 51) / +.2365839716528537483710069473279E-26     /
      DATA AE14CS( 52) / -.1003656195025305334065834526856E-26     /
      DATA AE14CS( 53) / +.4281490878094161131286642556927E-27     /
      DATA AE14CS( 54) / -.1836345261815318199691326958250E-27     /
      DATA AE14CS( 55) / +.7917798231349540000097468678144E-28     /
      DATA AE14CS( 56) / -.3431542358742220361025015775231E-28     /
      DATA AE14CS( 57) / +.1494705493897103237475066008917E-28     /
      DATA AE14CS( 58) / -.6542620279865705439739042420053E-29     /
      DATA AE14CS( 59) / +.2877581395199171114340487353685E-29     /
      DATA AE14CS( 60) / -.1271557211796024711027981200042E-29     /
      DATA AE14CS( 61) / +.5644615555648722522388044622506E-30     /
      DATA AE14CS( 62) / -.2516994994284095106080616830293E-30     /
      DATA AE14CS( 63) / +.1127259818927510206370368804181E-30     /
      DATA AE14CS( 64) / -.5069814875800460855562584719360E-31     /
C
C     ------------------------------------------------------------------
C
      X = XIN
      IF (NTAE10 .eq. 0) THEN
         BIGNUM = R1MACH(2)
         ETA = 0.1E0 * R1MACH(3)
         call SINITS(AE10CS, 50, ETA, NTAE10)
         call SINITS(AE11CS, 60, ETA, NTAE11)
         call SINITS(AE12CS, 41, ETA, NTAE12)
         call SINITS( E11CS, 29, ETA, NTE11 )
         call SINITS( E12CS, 25, ETA, NTE12 )
         call SINITS(AE13CS, 50, ETA, NTAE13)
         call SINITS(AE14CS, 64, ETA, NTAE14)
C
C     SETTING XMAX TO AVOID UNDERFLOW.
C
         EARG1 = -log(R1MACH(1))
         XMAX = (EARG1 - log(EARG1)) * FAC
C
C     SETTING XMIN TO AVOID OVERFLOW.
C
         EARG2 = log(BIGNUM)
         DELTA = log(EARG2)
         XMIN = -(EARG2 + DELTA) * FAC
      END IF
C
      IF (X.LE.(-1.E0)) THEN
         IF (X.LE.(-32.E0)) THEN
            IF (X .LT. XMIN) THEN
              CALL SERM1( NAME,1,0,MSG1,'X',XIN,',' )
              CALL SERV1('Limit |X|', XMIN,'.')
              SE1 = -BIGNUM
            ELSE
              SE1 = exp(-X-DELTA) * (EARG2/X) *
     *           (1.E0+SCSEVL(64.E0/X+1.E0,AE10CS,NTAE10))
            END IF
         ELSE IF (X.LE.(-8.E0)) THEN
            SE1 = exp(-X)/X * (1.E0+SCSEVL((64.E0/X+5.E0)/3.E0, AE11CS,
     1         NTAE11))
         ELSE IF (X.LE.(-4.E0)) THEN
            SE1 = exp(-X)/X*(1.E0+SCSEVL(16.E0/X+3.E0, AE12CS, NTAE12))
         ELSE
            SE1 = -log(-X) + SCSEVL((2.E0*X+5.E0)/3.E0, E11CS, NTE11)
         END IF
      ELSE IF (X .LE. 1.0E0) THEN
         IF (X .EQ. ZERO) THEN
           CALL ERMSG( NAME,2,0,MSG2,'.' )
           SE1 = BIGNUM
         ELSE
           SE1 = (-log(abs(X)) - 0.6875E0 + X) + SCSEVL(X,E12CS,NTE12)
         END IF
      ELSE IF (X.LE.4.0E0) THEN
         SE1 = exp(-X)/X * (1.E0 + SCSEVL((8.E0/X-5.E0)/3.E0, AE13CS,
     1      NTAE13))
      ELSE IF (X.LE.XMAX) THEN
         SE1 = exp(-X)/X * (1.E0 + SCSEVL(8.E0/X-1.E0, AE14CS, NTAE14))
      ELSE
         CALL SERM1( NAME,3,0,MSG3,'X',XIN,',' )
         CALL SERV1('Limit |X|', XMAX,'.')
         SE1 = ZERO
      END IF
C
      RETURN
C
      END

.