SUBROUTINE DUMKA3(N,T,TEND,H0,FUN,Y,Z, & Z1,Y_EN,Z_EN,OLDEIGENVECTOR,ATOL,RTOL,IWORK) c ------- ABSTRACT ----------------------------------- c DUMKA3 - integrates initial value problems for systems of first c order ordinary differential equations y'=f(y,t). c It is based on a family of explicit Runge-Kutta-Chebyshev formulas c of order three. c It uses optimal third order accuracy stability polynomials c with the largest stability region along the negative real axis. c c c The algorithm dumka3 and dumka4 is described in the article: c Medovikov A.A. High order explicit methods for parabolic c equations. BIT, V38,No2,pp.372-390,1998 c c c ------ INPUT PARAMETERS -------------------------- c c In ordert to use dumka3.f two external procedures are required c c double precision FUNCTION RHO (N,T,Y)- returns UPPER estimate of the c c spectral radius of the Jacobian matrix of f(t,y) at (t,y). c Must have the following form c double precision function rho(n,t,y) c double precision y(n),t c rho=... c return c end c for example for y'=-50*y, you can use any number c larger than 50 (for example rho = 51.), or if you have c working Euler method with the step size h_{euler}, you can c try to use rho = 2./h_{euler} c also N.b. Gerschgorin's theorem can be helpful. c c SUBROUTINE FUN(N,T,Y,F) - right hand side of the ODE: y'=f(y,t) c c N - number of variables, c T - current time, C TEND - time interval, c H0 - initial step size, c FUN - external procedure for right hand evaluation, c Y - vector of variables, c ATOL(*) : Absolute and relative error tolerances c RTOL(*) can be both scalar (double precision) c or vectors of length N. c IWORK(*): Integer array of length 2 that gives information c on how the problem is to be solved c c IWORK(1): =0 DUMKA3 attempts to compute the spectral radius c internally. But you should define a c dummy function c c double precision function rho(n,t,y) c double precision y(n),t c integer n c rho = 0 c return c end c c =1 RHO is used to obtain upper estimate on c spectral radius. c c IWORK(2): =0 The Jacobian is not constant. c =1 The Jacobian is constant, c the function rho is called only once c IWORK(3): =0 dumka3 integrates solution until t=tend c =1 the code returns after each step t_i chosen c automatically between [t,tend] (solution c at t_i is in y(*) ). c To continue call dumka3 again without changing c any arguments.. c IWORK(4): =0 Atol and rtol are scalar (use arrays rtol(1) and atol(1)). c =1 Atol and rtol are array of length n (use arrays rtol(n) and atol(n)). c c Z,Z1,Y_EN,Z_EN,OLDEIGENVECTOR - work arrays size n c c ------ OUTPUT PARAMETERS -------------------------- c c t - T-value for which the solution has been computed c c y(n) - Numerical solution at tend. c c c IWORK(5) - Number of function evaluations. c IWORK(6) - Number of steps. c IWORK(7) - Number of accepted steps. c IWORK(8) - Number of rejected steps. c IWORK(9) - Number of evaluations of f used c to estimate the spectral radius c (equal to zero if iwork(1)=1). c IWORK(10) - Maximum number of stages used. c IWORK(11) - Maximum value of the estimated c bound for the spectral radius c (rounded to the nearest integer). c IWORK(12) - Minimum value of the estimated c bound for the spectral radius c (rounded to the nearest integer). c c This program was written by Alexei A. Medovikov. c c E-mail: amedovik@math.tulane.edu c c Many thanks Prof. Lebedev V.I. for exellent algorithm of ordering c roots of stability polynomials, and Prof. Wanner G. and Hairer E. c for step size control algorithm. c c --------------------------------------------------- PARAMETER(N_P=16) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DOUBLE PRECISION STAB_REG(N_P),COEF(736,6),C_2,C_3,C_4, & A_21,A_31,A_41,A_32,A_42,A_43,T,T2,T3,T4,R,COU,H_NEW EXTERNAL FUN INTEGER N_DEG(N_P),INDEX_FIRST(N_P),INDEX_LAST(N_P) DOUBLE PRECISION Y(*),Z(*),Z1(*),Y_EN(*),Z_EN(*),OLDEIGENVECTOR(*) DOUBLE PRECISION ATOL(*),RTOL(*) INTEGER IWORK(12) DATA N_DEG/3,6,9,15,21,27,36,48,63,81,135,189,243,324,432,576/ DATA INDEX_FIRST/1,2,4,7,12,19,28,40,56,77,104,149,212,293, & 401,545/ DATA INDEX_LAST/1,3,6,11,18,27,39,55,76,103,148,211,292,400, & 544,736/ DATA STAB_REG/ 2.5005127005D0, 9.784756428464169D0, & 23.818760475282560d0, & 68.932817124349140D0,136.648186730571300D0, & 226.8897061613812D0,404.7232537379578d0,720.7401653373073D0, & 1337.1268312643200D0, &3266.0271134242240D0,9072.0255133253400D0,17784.544093341440D0, &27323.320790121520985D0,50666.52463738415826D0, &88860.62470326659968122D0,142602.586088039999594D0/ DATA (COEF(1,J),J=1,6)/ &0.924574112262462D0,0.085748762860431D0,0.28771294386877D0, &0.085748762860431D0,0.28771294386877D0,0.6265382932708D0/ DATA ((COEF(I,J),J=1,6),I=2,3)/ c &.1231817161250036D+00,.1231817161250036D+0,.1042566436020026D+00, c &.1231817161250036D+00,.1042566436020026D+0,.3892933924117777D+00, c &-.3166077728088352D+00,.9929505828357067D0,-.7107751811698341D00, c &.9929505828357067D+00,-.7107751811698341D0,.1010928461953435D+00/ &.10425664360200260D0, 0.10425664360200260D0, .38929339241177780D0, &.10425664360200260D0, 0.38929339241177780D0, .12318171612500360D0, &-0.43147899243339550D0,.49334989598763150D0,-.29116420359694450D0, &.49334989598763150D0,-0.29116420359694450D0, .18108255547052920D0/ DATA ((COEF(I,J),J=1,6),I=4,6)/ &0.5178985446356218D-1,0.5178985446356218D-1,0.1499671381392263D+0, &0.5178985446356218D-1,0.1499671381392263D+0,0.8987058729109126D-1, &0.4518984000985679D-1,0.4518984000985679D-1,0.3655365960800432D+0, &0.4518984000985679D-1,0.3655365960800432D+0,0.4232408050425635D-1, &-.5350188768628940D+0,0.3882284364912002D+0,-.1974370120534677D+0, &0.3882284364912002D+0,-.1974370120534677D+0,0.6453047907423139D-1/ DATA ((COEF(I,J),J=1,6),I=7,11)/ &.01454747477255407D0, 0.01454747477255407D0,.02079304726135118D0, &.01454747477255407D0, 0.02079304726135118D0,.01487782470642270D0, &.07700788949112707D0, 0.07700788949112707D0,.01836387448690710D0, &.07700788949112707D0, 0.01836387448690710D0,.35446354514643800D0, &.02432895596370560D0, 0.02432895596370560D0,.05145775978136601D0, &.02432895596370560D0, 0.05145775978136601D0,.01669108438830406D0, &.13626416224605870D0, 0.13626416224605870D0,.01556953144025836D0, &.13626416224605870D0, 0.01556953144025836D0,.03779049284450304D0, &-.74030990899340230D0,.31849582639491420D0,-.13024769688516420D0, &.31849582639491420D0,-0.13024769688516420D0,.02959622796125402D0/ DATA ((COEF(I,J),J=1,6),I=12,18)/ &.01477637829587406D0,.01477637829587406D0, 0.00784871042703017D0, &.01477637829587406D0,.00784871042703017D0, 0.13270367742571640D0, &.00872655659409207D0,.00872655659409207D0, 0.04842768909048031D0, &.00872655659409207D0,.04842768909048031D0, 0.01136228953669463D0, &.02112866077563068D0,.02112866077563068D0, 0.00741183976173058D0, &.02112866077563068D0,.00741183976173058D0, 0.35151303099059050D0, &.01282729908575624D0,.01282729908575624D0, 0.01742474299312111D0, &.01282729908575624D0,.01742474299312111D0, 0.00822345466591715D0, &.07384130062447007D0,.07384130062447007D0, 0.00758257251237930D0, &.07384130062447007D0,.00758257251237930D0, 0.02651410461353078D0, &.00938667991391738D0,.00938667991391738D0, 0.03477132786973448D0, &.00938667991391738D0,.03477132786973448D0, 0.00732840358129688D0, &-.71406485063527700D0,.31962173072924080D0,-.13166563685102140D0, &.31962173072924080D0,-.13166563685102140D0, .01024518736381746D0/ DATA ((COEF(I,J),J=1,6),I=19,27)/ &.5107503045891083D-2,.5107503045891083D-02,0.3362065448136852D-1, &.5107503045891083D-2,.3362065448136852D-01,0.7188013295116514D-2, &.1156154505769375D-1,.1156154505769375D-01,0.4502527186119195D-2, &.1156154505769375D-1,.4502527186119195D-02,0.4725209571486408D-1, &.8864532407434520D-2,.8864532407434520D-02,0.4725743401141243D-2, &.8864532407434520D-2,.4725743401141243D-02,0.7259492479863046D-1, &.5372265460448424D-2,.5372265460448424D-02,0.2536472039688646D-1, &.5372265460448424D-2,.2536472039688646D-01,0.6586314262087265D-2, &.1355407127646629D-1,.1355407127646629D-01,0.4441351277928004D-2, &.1355407127646629D-1,.4441351277928004D-02,0.3502521986497042D+0, &.6097135119221143D-2,.6097135119221143D-02,0.1623684839319424D-1, &.6097135119221143D-2,.1623684839319424D-01,0.4596423609977717D-2, &.1312876556408217D+0,.1312876556408217D+00,0.5697996302599818D-2, &.1312876556408217D+0,.5697996302599818D-02,0.1996561282544545D-1, &.7933063094502557D-2,.7933063094502557D-02,0.4411179932825468D-2, &.7933063094502557D-2,.4411179932825468D-02,0.1004394121139120D-1, &-.6551731045007839D0,.3298286840565381D+00,-.1419813471117554D0, &.3298286840565381D+0,-.1419813471117554D0,0.4894346213457858D-02/ DATA ((COEF(I,J),J=1,6),I=28,39)/ &.2472004197149452D-2,.2472004197149452D-02,0.2845804896378570D-2, &.2472004197149452D-2,.2845804896378570D-02,0.3384488891281092D-2, &.3578890608230160D-2,.3578890608230160D-02,0.4381540707760226D-2, &.3578890608230160D-2,.4381540707760226D-02,0.1548938684256209D-1, &.4748598631674122D-2,.4748598631674122D-02,0.7174737492322378D-1, &.4748598631674122D-2,.7174737492322378D-01,0.5184597369665364D-2, &.1279818555220185D-1,.1279818555220185D-01,0.2481471514472173D-2, &.1279818555220185D-1,.2481471514472173D-02,0.5706570136850913D-2, &.1303049568819145D+0,.1303049568819145D+00,0.2500551935542241D-2, &.1303049568819145D+0,.2500551935542241D-02,0.1922188870896997D-1, &.1079337211241642D-1,.1079337211241642D-01,0.4070498864709592D-2, &.1079337211241642D-1,.4070498864709592D-02,0.2529541836892036D-2, &.3495107918546922D+0,.3495107918546922D+00,0.2681439019200542D-2, &.3495107918546922D+0,.2681439019200542D-02,0.7107926900931334D-2, &.4646569953972695D-1,.4646569953972695D-01,0.2568898066542698D-2, &.4646569953972695D-1,.2568898066542698D-02,0.8061631024330287D-2, &.3286244688464888D-1,.3286244688464888D-01,0.2619253475244260D-2, &.3286244688464888D-1,.2619253475244260D-02,0.9260036321170762D-2, &.2461846538572608D-1,.2461846538572608D-01,0.2950959867147192D-2, &.2461846538572608D-1,.2950959867147192D-02,0.6337320049319617D-2, &.3074015316944754D-2,.3074015316944754D-02,0.3217487048615952D-2, &.3074015316944754D-2,.3217487048615952D-02,0.2756513790498370D-2, &-.6576608049720105D0,.3283593691684084D+00,-.1404975080448568D+0, &.3283593691684084D+0,-.1404975080448568D0,0.3805529719813893D-2/ DATA ((COEF(I,J),J=1,6),I=40,55)/ &.2263873446504803D-2,0.2263873446504803D-02,.1396806966273506D-2, &.2263873446504803D-2,0.1396806966273506D-02,.2874096214080187D-2, &.1405856786430858D-2,0.1405856786430858D-02,.3082967237408208D-2, &.1405856786430858D-2,0.3082967237408208D-02,.3922287445663546D-2, &.1433480445994300D-2,0.1433480445994300D-02,.1508985502511977D-1, &.1433480445994300D-2,0.1508985502511977D-1,.1452300237076606D-2, &.4299540037293374D-2,0.4299540037293374D-2,.1474676421856554D-2, &.4299540037293374D-2,0.1474676421856554D-2,.7128613707930451D-1, &.1500817396882159D-2,0.1500817396882159D-2,.4744800792404933D-2, &.1500817396882159D-2,0.4744800792404933D-2,.1239756641917538D-1, &.1698486544452855D-2,0.1698486544452855D-2,.5275134958587701D-2, &.1698486544452855D-2,0.5275134958587701D-2,.1297655017833546D+0, &.1754329508563495D-2,0.1754329508563495D-2,.1882152472414553D-1, &.1754329508563495D-2,0.1882152472414553D-1,.1039021150600575D-1, &.3600074660496816D-2,0.3600074660496816D-2,.1530972551686858D-2, &.3600074660496816D-2,0.1530972551686858D-2,.3490549853839776D+0, &.3322903844895304D-2,0.3322903844895304D-2,.6690572174716599D-2, &.3322903844895304D-2,0.6690572174716599D-2,.1565437511457520D-2, &.4603762606912074D-1,0.4603762606912074D-1,.1887080972842806D-2, &.4603762606912074D-1,0.1887080972842806D-2,.7649980453576852D-2, &.1604560571449509D-2,0.1604560571449509D-2,.3245045171225396D-1, &.1604560571449509D-2,0.3245045171225396D-1,.1965669974385304D-2, &.8853146753713483D-2,0.8853146753713483D-2,.2421450066631322D-1, &.8853146753713483D-2,0.2421450066631322D-1,.2389190339801504D-2, &.5913364519720735D-2,0.5913364519720735D-2,.2530806273930749D-2, &.5913364519720735D-2,0.2530806273930749D-2,.1418045398325408D-2, &.2691365777559749D-2,0.2691365777559749D-2,.1390816912287954D-2, &.2691365777559749D-2,0.1390816912287954D-2,.2053776804223216D-2, &.1387834738590928D-2,0.1387834738590928D-2,.1648750573442188D-2, &.1387834738590928D-2,0.1648750573442188D-2,.2152673169773036D-2, &-.6637857547745568D0,.3265004113646613D0,-.1386521871511077D+00, &.3265004113646613D+0,-.1386521871511077D0,.1816937001325805D-02/ DATA ((COEF(I,J),J=1,6),I=56,76)/ &.1497462268771500D-2,.1497462268771500D-2,0.8016352455861625D-03, &.1497462268771500D-2,.8016352455861625D-3,0.1134428570591172D-01, &.9107244117098982D-3,.9107244117098982D-3,0.4209062910423610D-02, &.9107244117098982D-3,.4209062910423610D-2,0.1115182220312366D-02, &.2278506030972934D-2,.2278506030972934D-2,0.7536034930991046D-03, &.2278506030972934D-2,.7536034930991046D-3,0.1220779831587775D+00, &.8718837367637042D-3,.8718837367637042D-3,0.5300278416578640D-02, &.8718837367637042D-3,.5300278416578640D-2,0.1200547418986776D-02, &.1989485460853075D-2,.1989485460853075D-2,0.7621317815369420D-03, &.1989485460853075D-2,.7621317815369420D-3,0.4279900105852680D-01, &.1361585867641604D-2,.1361585867641604D-2,0.1663462555411732D-02, &.1361585867641604D-2,.1663462555411732D-2,0.8255956252741248D-03, &.8041038269364196D-2,.8041038269364196D-2,0.7825383597586035D-03, &.8041038269364196D-2,.7825383597586035D-3,0.1733501096812087D-01, &.9570875386251864D-3,.9570875386251864D-3,0.3439141925532424D-02, &.9570875386251864D-3,.3439141925532424D-2,0.1043622477050726D-02, &.2649703125972344D-2,.2649703125972344D-2,0.7489205255105934D-03, &.2649703125972344D-2,.7489205255105934D-3,0.3663467233197790D+00, &.8549622570811375D-3,.8549622570811375D-3,0.4704548039083622D-02, &.8549622570811375D-3,.4704548039083622D-2,0.1155963691837909D-02, &.1868538294053988D-2,.1868538294053988D-2,0.7573780373303584D-03, &.1868538294053988D-2,.7573780373303584D-3,0.3006933585815754D-01, &.1302844396906396D-2,.1302844396906396D-2,0.1576206717853303D-02, &.1302844396906396D-2,.1576206717853303D-2,0.8129768269505178D-03, &.6919669611338770D-2,.6919669611338770D-2,0.7746793699320062D-03, &.6919669611338770D-2,.7746793699320062D-3,0.1385421077615084D-01, &.9328931714716148D-3,.9328931714716148D-3,0.3136728923895398D-02, &.9328931714716148D-3,.3136728923895398D-2,0.1012249102774941D-02, &.2452066725852897D-2,.2452066725852897D-2,0.7479885938962553D-03, &.2452066725852897D-2,.7479885938962553D-3,0.6654408106210975D-01, &.8904320252826596D-3,.8904320252826596D-3,0.3792559920683190D-02, &.8904320252826596D-3,.3792559920683190D-2,0.1249348534492854D-02, &.2125317026738749D-2,.2125317026738749D-2,0.7678889941591767D-03, &.2125317026738749D-2,.7678889941591767D-3,0.2236984568081352D-01, &.1426211069499304D-2,.1426211069499304D-2,0.1760431216809798D-02, &.1426211069499304D-2,.1760431216809798D-2,0.8395635973886352D-03, &.6025220837783022D-2,.6025220837783022D-2,0.7915075709375068D-03, &.6025220837783022D-2,.7915075709375068D-3,0.9473449290746352D-02, &.1077840803925363D-2,.1077840803925363D-2,0.2876007491770444D-02, &.1077840803925363D-2,.2876007491770444D-2,0.7507890407502578D-03, &-.7025868366046315D0,.3410573174186340D+0,-.1408407318719962D+00, &.3410573174186340D+0,-.1408407318719962D0,0.9834760180518672D-03/ DATA ((COEF(I,J),J=1,6),I=77,103)/ &.3602252890687476D-3,.3602252890687476D-3,.2053141684033504D-2, &.3602252890687476D-3,.2053141684033504D-2,.4831209951340756D-3, &.8380603820115802D-3,.8380603820115802D-3,.3111037886540830D-3, &.8380603820115802D-3,.3111037886540830D-3,.2057135989638816D-1, &.5687593269538123D-3,.5687593269538123D-3,.6644986450860436D-3, &.5687593269538123D-3,.6644986450860436D-3,.3356309293764894D-3, &.3663093727840670D-3,.3663093727840670D-3,.2783006991234028D-2, &.3663093727840670D-3,.2783006991234028D-2,.5137328896150039D-3, &.8829983161130827D-3,.8829983161130827D-3,.3146625430082048D-3, &.8829983161130827D-3,.3146625430082048D-3,.1176114948662877D-1, &.5899945657034291D-3,.5899945657034291D-3,.7251651912695405D-3, &.5899945657034291D-3,.7251651912695405D-3,.3443798208440408D-3, &.3525769387698612D-2,.3525769387698612D-2,.3219448100672548D-3, &.3525769387698612D-2,.3219448100672548D-3,.6371551060177030D-2, &.3872985816371484D-3,.3872985816371484D-3,.1468179847166060D-2, &.3872985816371484D-3,.1468179847166060D-2,.4332912305927524D-3, &.1046913570305507D-2,.1046913570305507D-2,.3069038328618102D-3, &.1046913570305507D-2,.3069038328618102D-3,.8676048193961264D-1, &.5393647957387747D-2,.5393647957387747D-2,.3249207644000284D-3, &.5393647957387747D-2,.3249207644000284D-3,.9370103401532334D-2, &.4038633059540496D-3,.4038633059540496D-3,.1875584086669360D-2, &.4038633059540496D-3,.1875584086669360D-2,.4445014352820028D-3, &.1270617970958242D-2,.1270617970958242D-2,.3075985493257752D-3, &.1270617970958242D-2,.3075985493257752D-3,.2947128668608194D-1, &.6129127967380504D-3,.6129127967380504D-3,.3281863544835643D-3, &.6129127967380504D-3,.3281863544835643D-3,.4628932190009235D-2, &.3728320056096484D-3,.3728320056096484D-3,.1721383865094658D-2, &.3728320056096484D-3,.1721383865094658D-2,.4565007875900184D-3, &.9323756157185432D-3,.9323756157185432D-3,.3085281054284440D-3, &.9323756157185432D-3,.3085281054284440D-3,.4630801535879760D-1, &.3545562738952262D-3,.3545562738952262D-3,.2259014678545103D-2, &.3545562738952262D-3,.2259014678545103D-2,.4978866383126903D-3, &.7970559006761336D-3,.7970559006761336D-3,.5490570182311138D-3, &.7970559006761336D-3,.5490570182311138D-3,.6935763377230479D-3, &.3127578179191516D-3,.3127578179191516D-3,.1523284967991621D-1, &.3127578179191516D-3,.1523284967991621D-1,.3398351883455558D-3, &.3120231281205192D-2,.3120231281205192D-2,.3192487615768242D-3, &.3120231281205192D-2,.3192487615768242D-3,.7650648974430584D-2, &.3953019564372914D-3,.3953019564372914D-3,.1363527979037656D-2, &.3953019564372914D-3,.1363527979037656D-2,.4228143427761484D-3, &.1113584996470221D-2,.1113584996470221D-2,.3064418510275818D-3, &.1113584996470221D-2,.3064418510275818D-3,.4441014793556532D+0, &.3492809420252234D-3,.3492809420252234D-3,.2499554821215150D-2, &.3492809420252234D-3,.2499554821215150D-2,.4693511141335305D-3, &.7595495666968666D-3,.7595495666968666D-3,.3168239173466168D-3, &.7595495666968666D-3,.3168239173466168D-3,.4019457444625864D-2, &.5307544056836153D-3,.5307544056836153D-3,.6376835005816956D-3, &.5307544056836153D-3,.6376835005816956D-3,.3317524080451016D-3, &.1586624743617795D-2,.1586624743617795D-2,.3798189696643980D-3, &.1586624743617795D-2,.3798189696643980D-3,.9867829884641076D-3, &.3096953289201503D-3,.3096953289201503D-3,.1187767235014202D-2, &.3096953289201503D-3,.1187767235014202D-2,.3062112079598173D-3, &-.4595292409547097D0,.5537618913875574D+0,-.2982973149029844D0, &.5537618913875574D0,-.2982973149029844D+0,.4130202781675714D-3/ DATA ((COEF(I,J),J=1,6),I=104,120)/ &.1300979717625048D-3,.1300979717625048D-3,.7233069985422061D-3, &.1300979717625048D-3,.7233069985422061D-3,.1479865979466228D-3, &.4325905695142040D-3,.4325905695142040D-3,.1118913719696372D-3, &.4325905695142040D-3,.1118913719696372D-3,.7585567260231847D-2, &.1161063343045094D-3,.1161063343045094D-3,.2191575521772395D-2, &.1161063343045094D-3,.2191575521772395D-2,.1567383447080588D-3, &.3718918170158877D-3,.3718918170158877D-3,.1205361090670254D-3, &.3718918170158877D-3,.1205361090670254D-3,.1293982340612670D-2, &.2370803631413284D-3,.2370803631413284D-3,.1388585518045632D-3, &.2370803631413284D-3,.1388585518045632D-3,.5354856462012968D-3, &.2061424702090348D-3,.2061424702090348D-3,.1106817938650610D-3, &.2061424702090348D-3,.1106817938650610D-3,.2939862772511351D-1, &.2205316668063107D-3,.2205316668063107D-3,.1181459936611318D-3, &.2205316668063107D-3,.1181459936611318D-3,.1652950036594520D-2, &.1342074954405650D-3,.1342074954405650D-3,.6180865016739393D-3, &.1342074954405650D-3,.6180865016739393D-3,.1643005259557435D-3, &.3352779765046980D-3,.3352779765046980D-3,.1110729852040027D-3, &.3352779765046980D-3,.1110729852040027D-3,.1516714664210126D-1, &.1276330323761358D-3,.1276330323761358D-3,.8103090623879464D-3, &.1276330323761358D-3,.8103090623879464D-3,.1791822232033586D-3, &.2866899751783916D-3,.2866899751783916D-3,.1125948734121003D-3, &.2866899751783916D-3,.1125948734121003D-3,.5328164964739980D-2, &.1521857067779652D-3,.1521857067779652D-3,.4003043817587651D-3, &.1521857067779652D-3,.4003043817587651D-3,.1149303257924375D-3, &.2716175522656955D-2,.2716175522656955D-2,.1422901581051835D-3, &.2716175522656955D-2,.1422901581051835D-3,.4899227193665029D-3, &.1223369806032192D-3,.1223369806032192D-3,.1117406689032896D-2, &.1223369806032192D-3,.1117406689032896D-2,.1975790859531569D-3, &.2495178271696324D-3,.2495178271696324D-3,.1103223234115986D-3, &.2495178271696324D-3,.1103223234115986D-3,.4621346400269236D-1, &.1288392090552591D-3,.1288392090552591D-3,.7648510145779159D-3, &.1288392090552591D-3,.7648510145779159D-3,.1759521226738778D-3, &.2953654408369954D-3,.2953654408369954D-3,.1122270952009704D-3, &.2953654408369954D-3,.1122270952009704D-3,.6306291044327727D-2/ DATA ((COEF(I,J),J=1,6),I=121,148)/ &.1544159445083552D-3,.1544159445083552D-3,.3856535572687610D-3, &.1544159445083552D-3,.3856535572687610D-3,.1155004409094379D-3, &.2432424549003238D-2,.2432424549003238D-2,.1405396089644810D-3, &.2432424549003238D-2,.1405396089644810D-3,.5118749472635047D-3, &.1214150913210891D-3,.1214150913210891D-3,.1200672646554655D-2, &.1214150913210891D-3,.1200672646554655D-2,.2017653111614368D-3, &.2431468972385023D-3,.2431468972385023D-3,.1104120448884408D-3, &.2431468972385023D-3,.1104120448884408D-3,.8667516963670382D-1, &.1327802415082384D-3,.1327802415082384D-3,.6502766499608423D-3, &.1327802415082384D-3,.6502766499608423D-3,.1670349237499396D-3, &.3244400365212383D-3,.3244400365212383D-3,.1113146353560780D-3, &.3244400365212383D-3,.1113146353560780D-3,.1169603413901868D-1, &.2257889771491832D-3,.2257889771491832D-3,.1189025598809207D-3, &.2257889771491832D-3,.1189025598809207D-3,.1517834231694204D-2, &.2155136066576732D-3,.2155136066576732D-3,.1174283304930101D-3, &.2155136066576732D-3,.1174283304930101D-3,.1807493902559267D-2, &.1356948411014018D-3,.1356948411014018D-3,.5883857080272527D-3, &.1356948411014018D-3,.5883857080272527D-3,.1616761741067229D-3, &.3467639370112126D-3,.3467639370112126D-3,.1102625631020175D-3, &.3467639370112126D-3,.1102625631020175D-3,.2050372318892889D-1, &.1264778138423625D-3,.1264778138423625D-3,.8601876168349558D-3, &.1264778138423625D-3,.8601876168349558D-3,.1825520166041481D-3, &.2784656473621648D-3,.2784656473621648D-3,.1129951170496302D-3, &.2784656473621648D-3,.1129951170496302D-3,.4563202656757644D-2, &.2562128709546526D-3,.2562128709546526D-3,.1233029345562954D-3, &.2562128709546526D-3,.1233029345562954D-3,.1042793125539889D-2, &.1143953290968708D-3,.1143953290968708D-3,.3053689202336192D-2, &.1143953290968708D-3,.3053689202336192D-2,.1441129024472512D-3, &.4694781542091936D-3,.4694781542091936D-3,.1935737268930150D-3, &.4694781542091936D-3,.1935737268930150D-3,.1500437960573724D-3, &.4159214629811822D-3,.4159214629811822D-3,.1108621149877430D-3, &.4159214629811822D-3,.1108621149877430D-3,.4439828224987684D+0, &.1314110422564394D-3,.1314110422564394D-3,.6852411356133314D-3, &.1314110422564394D-3,.6852411356133314D-3,.1860686704835823D-3, &.2706625914249481D-3,.2706625914249481D-3,.1134282750606420D-3, &.2706625914249481D-3,.1134282750606420D-3,.3953467060562924D-2, &.1167487125030652D-3,.1167487125030652D-3,.1985382810604559D-2, &.1167487125030652D-3,.1985382810604559D-2,.1591569794186695D-3, &.3589495713626626D-3,.3589495713626626D-3,.1243141796164223D-3, &.3589495713626626D-3,.1243141796164223D-3,.9756737148526784D-3, &.1460107047877242D-3,.1460107047877242D-3,.4504072408797582D-3, &.1460107047877242D-3,.4504072408797582D-3,.1105318249567443D-3, &.9305108217882648D-2,.9305108217882648D-2,.2107213670342112D-3, &.9305108217882648D-2,.2107213670342112D-3,.1253720144824314D-3, &.9150779394173306D-3,.9150779394173306D-3,.2313000070835802D-3, &.9150779394173306D-3,.2313000070835802D-3,.1698845695088106D-3, &.3142031705565816D-3,.3142031705565816D-3,.1102326991347965D-3, &.3142031705565816D-3,.1102326991347965D-3,.3459507712468638D-2, &.1372444333085824D-3,.1372444333085824D-3,.5609250137857607D-3, &.1372444333085824D-3,.5609250137857607D-3,.1897398437653978D-3, &.2632533968883294D-3,.2632533968883294D-3,.1196989417613488D-3, &.2632533968883294D-3,.1196989417613488D-3,.1399018943806244D-2, &.1138948359758609D-3,.1138948359758609D-3,.3045246339149248D-3, &.1138948359758609D-3,.3045246339149248D-3,.1115873307883615D-3, &-.4760906730816906D0,.5431696227154399D+0,-.2877420590063439D0, &.5431696227154399D0,-.2877420590063439D+0,.1728550020493715D-3/ DATA ((COEF(I,J),J=1,6),I=149,169)/ &0.8794603386687293D-4,.8794603386687293D-4,.1559841060693178D-3, &0.8794603386687293D-4,.1559841060693178D-3,.6645585617088351D-4, &0.3657659216717562D-3,.3657659216717562D-3,.5639304508059442D-4, &0.3657659216717562D-3,.5639304508059442D-4,.3934735221765134D-2, &0.1347849952295428D-3,.1347849952295428D-3,.6142308606226574D-4, &0.1347849952295428D-3,.6142308606226574D-4,.6661893561026761D-3, &0.1054668600601382D-3,.1054668600601382D-3,.5927180544859897D-4, &0.1054668600601382D-3,.5927180544859897D-4,.1098680337850833D-2, &0.7071001772514381D-4,.7071001772514381D-4,.2747790909914622D-3, &0.7071001772514381D-4,.2747790909914622D-3,.7461855624800195D-4, &0.1886652725419461D-3,.1886652725419461D-3,.5644767750830879D-4, &0.1886652725419461D-3,.5644767750830879D-4,.1514703871410358D-1, &0.8479417947822094D-4,.8479417947822094D-4,.1669934028989259D-3, &0.8479417947822094D-4,.1669934028989259D-3,.6793473813586763D-4, &0.3266313227220002D-3,.3266313227220002D-3,.5674548854559881D-4, &0.3266313227220002D-3,.5674548854559881D-4,.6287282677721673D-2, &0.1106368401950715D-3,.1106368401950715D-3,.6053995731181062D-4, &0.1106368401950715D-3,.6053995731181062D-4,.7913976625843851D-3, &0.1143781572478465D-3,.1143781572478465D-3,.6000295999830059D-4, &0.1143781572478465D-3,.6000295999830059D-4,.8962483030328786D-3, &0.6899633034686719D-4,.6899633034686719D-4,.3041328611869304D-3, &0.6899633034686719D-4,.3041328611869304D-3,.8284258584704493D-4, &0.1751178732564210D-3,.1751178732564210D-3,.5658065044245652D-4, &0.1751178732564210D-3,.5658065044245652D-4,.9285720823459838D-2, &0.6468137427324778D-4,.6468137427324778D-4,.4308155322950510D-3, &0.6468137427324778D-4,.4308155322950510D-3,.9261260660914518D-4, &0.1431875236770320D-3,.1431875236770320D-3,.5757936051594736D-4, &0.1431875236770320D-3,.5757936051594736D-4,.2413816096453975D-2, &0.9930547940897894D-4,.9930547940897894D-4,.1296752264585964D-3, &0.9930547940897894D-4,.1296752264585964D-3,.6275338091083291D-4, &0.5416548200909192D-3,.5416548200909192D-3,.5842966551264406D-4, &0.5416548200909192D-3,.5842966551264406D-4,.1499203227632296D-2, &0.7323983531898397D-4,.7323983531898397D-4,.2422612758651356D-3, &0.7323983531898397D-4,.2422612758651356D-3,.7684141711636654D-4, &0.2097437725855342D-3,.2097437725855342D-3,.5625288509274553D-4, &0.2097437725855342D-3,.5625288509274553D-4,.2048301765309298D-1, &0.9020914763643894D-4,.9020914763643894D-4,.1493391786229736D-3, &0.9020914763643894D-4,.1493391786229736D-3,.6554149429275607D-4, &0.3961861287168930D-3,.3961861287168930D-3,.5729959745091044D-4, &0.3961861287168930D-3,.5729959745091044D-4,.3035041907808402D-2, &0.1022883561329008D-3,.1022883561329008D-3,.1249183946724889D-3, &0.1022883561329008D-3,.1249183946724889D-3,.6206573381626393D-4/ DATA ((COEF(I,J),J=1,6),I=170,190)/ &0.5989298808691080D-3,.5989298808691080D-3,.5883221466827894D-4, &0.5989298808691080D-3,.5883221466827894D-4,.1275310400484650D-2, &0.7193836788723518D-4,.7193836788723518D-4,.2576813057629360D-3, &0.7193836788723518D-4,.2576813057629360D-3,.7843378028023659D-4, &0.1987311001382434D-3,.1987311001382434D-3,.5630733649577963D-4, &0.1987311001382434D-3,.5630733649577963D-4,.4619111642398583D-1, &0.1124763955756592D-3,.1124763955756592D-3,.6026639769285611D-4, &0.1124763955756592D-3,.6026639769285611D-4,.8413183246610148D-3, &0.6845771781688004D-4,.6845771781688004D-4,.3150592204659011D-3, &0.6845771781688004D-4,.3150592204659011D-3,.8380407140692748D-4, &0.1709714508940318D-3,.1709714508940318D-3,.5665906501426691D-4, &0.1709714508940318D-3,.5665906501426691D-4,.7566397972122756D-2, &0.6510480557593150D-4,.6510480557593150D-4,.4129281267273862D-3, &0.6510480557593150D-4,.4129281267273862D-3,.9139271207113828D-4, &0.1462047251430903D-3,.1462047251430903D-3,.5743525528477778D-4, &0.1462047251430903D-3,.5743525528477778D-4,.2697552273731426D-2, &0.7762611102681718D-4,.7762611102681718D-4,.2041119533398660D-3, &0.7762611102681718D-4,.2041119533398660D-3,.5862637134354665D-4, &0.1380369292929180D-2,.1380369292929180D-2,.7257971136921318D-4, &0.1380369292929180D-2,.7257971136921318D-4,.2497756000560482D-3, &0.6240382512391533D-4,.6240382512391533D-4,.5691733826928257D-3, &0.6240382512391533D-4,.5691733826928257D-3,.1007734014166686D-3, &0.1272546822394078D-3,.1272546822394078D-3,.5627621280667196D-4, &0.1272546822394078D-3,.5627621280667196D-4,.8665955106287616D-1, &0.6173889013943212D-4,.6173889013943212D-4,.6311735522739983D-3, &0.6173889013943212D-4,.6311735522739983D-3,.1038521938386863D-3, &0.1226626631810164D-3,.1226626631810164D-3,.5904731686728005D-4, &0.1226626631810164D-3,.5904731686728005D-4,.1181975061792600D-2, &0.6599174040940792D-4,.6599174040940792D-4,.3804939876683146D-3, &0.6599174040940792D-4,.3804939876683146D-3,.8906065524080525D-4, &0.1525968531033595D-3,.1525968531033595D-3,.5717230974591810D-4, &0.1525968531033595D-3,.5717230974591810D-4,.3440825120605058D-2, &0.7131529866360054D-4,.7131529866360054D-4,.2660058367917730D-3, &0.7131529866360054D-4,.2660058367917730D-3,.7926513114741066D-4, &0.1935866985995343D-3,.1935866985995343D-3,.5634627339372908D-4, &0.1935866985995343D-3,.5634627339372908D-4,.2937709118909205D-1, &0.6082380692150844D-4,.6082380692150844D-4,.7458951270439268D-3, &0.6082380692150844D-4,.7458951270439268D-3,.1088569720659201D-3, &0.5974948342390019D-4,.5974948342390019D-4,.9568811950877288D-3, &0.5974948342390019D-4,.9568811950877288D-3,.1163447819683455D-3, &0.6742700987040936D-4,.6742700987040936D-4,.3389008440021692D-3, &0.6742700987040936D-4,.3389008440021692D-3,.8581386747399124D-4/ DATA ((COEF(I,J),J=1,6),I=191,211)/ &.1631749068467616D-3,.1631749068467616D-3,.5683996941383108D-4, &.1631749068467616D-3,.5683996941383108D-4,.5309276516506640D-2, &.6955097277938596D-4,.6955097277938596D-4,.2938052740290110D-3, &.6955097277938596D-4,.2938052740290110D-3,.8190880691002491D-4, &.1794421148111645D-3,.1794421148111645D-3,.5623734044645920D-4, &.1794421148111645D-3,.5623734044645920D-4,.1167634665834735D-1, &.6427091219301584D-4,.6427091219301584D-4,.4499550311805648D-3, &.6427091219301584D-4,.4499550311805648D-3,.9387014957107914D-4, &.1402819656449207D-3,.1402819656449207D-3,.5773199552171135D-4, &.1402819656449207D-3,.5773199552171135D-4,.2172974126167264D-2, &.9788282959394230D-4,.9788282959394230D-4,.1321839464002566D-3, &.9788282959394230D-4,.1321839464002566D-3,.6311462715182011D-4, &.5161554327051546D-3,.5161554327051546D-3,.5824198176735158D-4, &.5161554327051546D-3,.5824198176735158D-4,.1634333682182328D-2, &.7391926710804335D-4,.7391926710804335D-4,.2351131056298787D-3, &.7391926710804335D-4,.2351131056298787D-3,.7607902158886092D-4, &.2156422078856028D-3,.2156422078856028D-3,.5651020103131087D-4, &.2156422078856028D-3,.5651020103131087D-4,.4439592935638661D+0, &.8012090066203704D-4,.8012090066203704D-4,.3519247427973640D-3, &.8012090066203704D-4,.3519247427973640D-3,.8686413653407661D-4, &.1595077141859470D-3,.1595077141859470D-3,.5950581421224924D-4, &.1595077141859470D-3,.5950581421224924D-4,.1024035077128965D-2, &.1071344014445986D-3,.1071344014445986D-3,.6111812207724783D-4, &.1071344014445986D-3,.6111812207724783D-4,.7043043018724323D-3, &.1183790732789484D-3,.1183790732789484D-3,.5789324788946512D-4, &.1183790732789484D-3,.5789324788946512D-4,.1966781886872593D-2, &.7012205845840767D-4,.7012205845840767D-4,.2840336883698481D-3, &.7012205845840767D-4,.2840336883698481D-3,.8100185815976544D-4, &.2283078667002492D-3,.2283078667002492D-3,.5694256063267605D-4, &.2283078667002492D-3,.5694256063267605D-4,.4544403922678898D-2, &.6387314249555420D-4,.6387314249555420D-4,.4704660392145826D-3, &.6387314249555420D-4,.4704660392145826D-3,.1204839912812346D-3, &.5622957027120577D-4,.5622957027120577D-4,.4924825703519851D-3, &.5622957027120577D-4,.4924825703519851D-3,.5806321052489571D-4, &.1788887759333359D-2,.1788887759333359D-2,.7533827542425808D-4, &.1788887759333359D-2,.7533827542425808D-4,.2218241531500721D-3, &.9650377610676052D-4,.9650377610676052D-4,.1374827771253620D-3, &.9650377610676052D-4,.1374827771253620D-3,.6348779927329075D-4, &.9516672295803058D-4,.9516672295803058D-4,.1839542961703630D-3, &.9516672295803058D-4,.1839542961703630D-3,.5705331991573311D-4, &-.4557755622772656D0,.5559513892316170D+0,-.3005291574163789D0, &.5559513892316170D+0,-.3005291574163789D0,.6693416657141503D-4/ DATA ((COEF(I,J),J=1,6),I=212,232)/ &.55044886435908315D-4,.55044886435908315D-4,.10924521648384124D-3, &.55044886435908315D-4,.10924521648384124D-3,.44292243472068958D-4, &.21081952808812550D-3,.21081952808812550D-3,.36922232342692445D-4, &.21081952808812550D-3,.36922232342692445D-4,.42241719042147375D-2, &.72270051539673256D-4,.72270051539673256D-4,.39364616900087130D-4, &.72270051539673256D-4,.39364616900087130D-4,.52159834391255833D-3, &.78260909122679251D-4,.78260909122679251D-4,.39092798248711554D-4, &.78260909122679251D-4,.39092798248711554D-4,.57453143481445710D-3, &.44829588031624682D-4,.44829588031624682D-4,.19944100669939391D-3, &.44829588031624682D-4,.19944100669939391D-3,.54057102347705466D-4, &.11335614381928313D-3,.11335614381928313D-3,.36838793161261356D-4, &.11335614381928313D-3,.36838793161261356D-4,.57019347219706357D-2, &.42160613396809558D-4,.42160613396809558D-4,.27763204223781116D-3, &.42160613396809558D-4,.27763204223781116D-3,.60099349890311774D-4, &.93618971762938059D-4,.93618971762938059D-4,.37456578948844298D-4, &.93618971762938059D-4,.37456578948844298D-4,.16048221513304045D-2, &.41350770892630515D-4,.41350770892630515D-4,.35620561354022184D-3, &.41350770892630515D-4,.35620561354022184D-3,.64843118821321667D-4, &.84041566329960294D-4,.84041566329960294D-4,.38059247291629783D-4, &.84041566329960294D-4,.38059247291629783D-4,.95613058032920062D-3, &.47573572154190413D-4,.47573572154190413D-4,.14481206784523667D-3, &.47573572154190413D-4,.14481206784523667D-3,.50126274738492760D-4, &.13566508717168722D-3,.13566508717168722D-3,.36617514598916815D-4, &.13566508717168722D-3,.36617514598916815D-4,.12542671127637614D-1, &.42595984002896317D-4,.42595984002896317D-4,.26012414579348620D-3, &.42595984002896317D-4,.26012414579348620D-3,.96731928333512905D-4, &.58882888812587741D-4,.58882888812587741D-4,.37314966422694551D-4, &.58882888812587741D-4,.37314966422694551D-4,.19164386842987064D-2, &.66352792392554855D-4,.66352792392554855D-4,.81634285847049018D-4, &.66352792392554855D-4,.81634285847049018D-4,.40446038923794491D-4, &.38514732079839335D-3,.38514732079839335D-3,.38263011764130959D-4, &.38514732079839335D-3,.38263011764130959D-4,.84330830946168194D-3, &.46914819471622353D-4,.46914819471622353D-4,.16651447420864009D-3, &.46914819471622353D-4,.16651447420864009D-3,.50932253901212012D-4, &.15871412208610971D-3,.15871412208610971D-3,.36645077318639100D-4, &.15871412208610971D-3,.36645077318639100D-4,.31525691503432547D-1, &.71359273261155408D-4,.71359273261155408D-4,.39506554892796836D-4, &.71359273261155408D-4,.39506554892796836D-4,.49789559876827936D-3, &.75147767889215610D-4,.75147767889215610D-4,.38962816908880226D-4, &.75147767889215610D-4,.38962816908880226D-4,.60415210686553525D-3, &.55555970582991249D-4,.55555970582991249D-4,.10728640856105868D-3, &.55555970582991249D-4,.10728640856105868D-3,.44032685304033269D-4/ DATA ((COEF(I,J),J=1,6),I=233,253)/ &.21690526098680022D-3,.21690526098680022D-3,.36968696383592138D-4, &.21690526098680022D-3,.36968696383592138D-4,.36930831868967075D-2, &.45107612735098642D-4,.45107612735098642D-4,.19411809374750768D-3, &.45107612735098642D-4,.19411809374750768D-3,.53579829051401325D-4, &.11551384946628285D-3,.11551384946628285D-3,.36801789799676630D-4, &.11551384946628285D-3,.36801789799676630D-4,.67531692564439796D-2, &.40278706209645575D-4,.40278706209645575D-4,.40101134190648153D-3, &.40278706209645575D-4,.40101134190648153D-3,.67135960845910716D-4, &.80479245192672697D-4,.80479245192672697D-4,.38370327137740935D-4, &.80479245192672697D-4,.38370327137740935D-4,.79434353108711545D-3, &.58295721844918572D-4,.58295721844918572D-4,.98359430938613988D-4, &.58295721844918572D-4,.98359430938613988D-4,.42821589586023634D-4, &.25201589495371673D-3,.25201589495371673D-3,.37249125017575036D-4, &.25201589495371673D-3,.37249125017575036D-4,.21077907206935179D-2, &.46596377303863417D-4,.46596377303863417D-4,.17065386629320977D-3, &.46596377303863417D-4,.17065386629320977D-3,.51349025179459740D-4, &.12745113731596933D-3,.12745113731596933D-3,.36663467835956841D-4, &.12745113731596933D-3,.36663467835956841D-4,.21993424768517619D-1, &.60729432490863395D-4,.60729432490863395D-4,.92129976450777384D-4, &.60729432490863395D-4,.92129976450777384D-4,.41950668318758499D-4, &.28709338719081200D-3,.28709338719081200D-3,.37532398944882181D-4, &.28709338719081200D-3,.37532398944882181D-4,.14769845904712357D-2, &.64115508737857090D-4,.64115508737857090D-4,.85296137295142612D-4, &.64115508737857090D-4,.85296137295142612D-4,.40974986626216359D-4, &.34298536597357092D-3,.34298536597357092D-3,.37962725699828389D-4, &.34298536597357092D-3,.37962725699828389D-4,.10213699517661091D-2, &.47914199043072588D-4,.47914199043072588D-4,.15503706720996344D-3, &.47914199043072588D-4,.15503706720996344D-3,.49736643464922743D-4, &.13860378641239400D-3,.13860378641239400D-3,.36608333169139703D-4, &.13860378641239400D-3,.36608333169139703D-4,.92616992733712711D-1, &.73204537631936135D-4,.73204537631936135D-4,.39226715087758024D-4, &.73204537631936135D-4,.39226715087758024D-4,.54708259399162754D-3, &.44557838618547991D-4,.44557838618547991D-4,.20500331687559948D-3, &.44557838618547991D-4,.20500331687559948D-3,.37127258749107615D-4, &.11126755519498782D-3,.11126755519498782D-3,.36878936154620595D-4, &.11126755519498782D-3,.36878936154620595D-4,.48796679992481346D-2, &.45682991386666161D-4,.45682991386666161D-4,.18413747082473981D-3, &.45682991386666161D-4,.18413747082473981D-3,.52657158378245856D-4, &.12004887953998150D-3,.12004887953998150D-3,.36737153103038579D-4, &.12004887953998150D-3,.36737153103038579D-4,.99750351623657738D-2, &.42375689081510392D-4,.42375689081510392D-4,.26865301832779213D-3, &.42375689081510392D-4,.26865301832779213D-3,.59483969134597476D-4/ DATA ((COEF(I,J),J=1,6),I=254,274)/ &.95152379786018104D-4,.95152379786018104D-4,.37384109789850171D-4, &.95152379786018104D-4,.37384109789850171D-4,.17501860760722667D-2, &.65588699440111823D-4,.65588699440111823D-4,.82821346366277902D-4, &.65588699440111823D-4,.82821346366277902D-4,.40617818401824765D-4, &.37023704028143822D-3,.37023704028143822D-3,.38159330801636719D-4, &.37023704028143822D-3,.38159330801636719D-4,.89702536770757465D-3, &.47240499200164707D-4,.47240499200164707D-4,.16253732969068874D-3, &.47240499200164707D-4,.16253732969068874D-3,.50524742147609331D-4, &.13282986471746172D-3,.13282986471746172D-3,.36629761282483102D-4, &.13282986471746172D-3,.36629761282483102D-4,.49523516319748058D-1, &.43052600165222811D-4,.43052600165222811D-4,.24430112030043063D-3, &.43052600165222811D-4,.24430112030043063D-3,.57722094894444483D-4, &.10003679199991142D-3,.10003679199991142D-3,.37186562943089909D-4, &.10003679199991142D-3,.37186562943089909D-4,.23295724479988284D-2, &.67938790329243319D-4,.67938790329243319D-4,.79355134124617122D-4, &.67938790329243319D-4,.79355134124617122D-4,.40115755691604336D-4, &.41791214535637412D-3,.41791214535637412D-3,.38481315326362216D-4, &.41791214535637412D-3,.38481315326362216D-4,.74958632710527928D-3, &.46285021878322113D-4,.46285021878322113D-4,.17496438636607363D-3, &.46285021878322113D-4,.17496438636607363D-3,.51775278414661420D-4, &.12489893252886369D-3,.12489893252886369D-3,.36684938997666118D-4, &.12489893252886369D-3,.36684938997666118D-4,.16268620504023119D-1, &.70471458526692076D-4,.70471458526692076D-4,.76158162954803431D-4, &.70471458526692076D-4,.76158162954803431D-4,.39652582043053260D-4, &.47581212484097262D-3,.47581212484097262D-3,.38836723278455275D-4, &.47581212484097262D-3,.38836723278455275D-4,.63617993499509150D-3, &.43779050421102076D-4,.43779050421102076D-4,.22327735787388441D-3, &.43779050421102076D-4,.22327735787388441D-3,.56078915250456852D-4, &.10538855123368351D-3,.10538855123368351D-3,.37018344031328365D-4, &.10538855123368351D-3,.37018344031328365D-4,.32567191661860306D-2, &.45392037577823572D-4,.45392037577823572D-4,.18902109657543427D-3, &.45392037577823572D-4,.18902109657543427D-3,.53113268763947815D-4, &.11774369699559304D-3,.11774369699559304D-3,.36767913584558802D-4, &.11774369699559304D-3,.36767913584558802D-4,.81277945739512466D-2, &.41745767947950939D-4,.41745767947950939D-4,.29707215126818212D-3, &.41745767947950939D-4,.29707215126818212D-3,.61374633322977217D-4, &.90683746950053756D-4,.90683746950053756D-4,.37611596059445999D-4, &.90683746950053756D-4,.37611596059445999D-4,.13639620993797774D-2, &.63405347793084039D-4,.63405347793084039D-4,.86586305652146307D-4, &.63405347793084039D-4,.86586305652146307D-4,.41160515344792060D-4, &.33051501724888990D-3,.33051501724888990D-3,.37869731882423331D-4, &.33051501724888990D-3,.37869731882423331D-4,.10936233473859405D-2/ DATA ((COEF(I,J),J=1,6),I=275,292)/ &.48262545802113032D-4,.48262545802113032D-4,.15149886782325615D-3, &.48262545802113032D-4,.15149886782325615D-3,.52211243556735441D-4, &.14165099780195356D-3,.14165099780195356D-3,.36602213921238660D-4, &.14165099780195356D-3,.36602213921238660D-4,.43108450894259814D+0, &.57161647814334589D-4,.57161647814334589D-4,.10176601227740090D-3, &.57161647814334589D-4,.10176601227740090D-3,.43289113041943893D-4, &.23695483314780820D-3,.23695483314780820D-3,.37071192156164656D-4, &.23695483314780820D-3,.37071192156164656D-4,.25886117209878556D-2, &.89278712781927145D-4,.89278712781927145D-4,.68761889807490860D-4, &.89278712781927145D-4,.68761889807490860D-4,.39957124836697671D-4, &.43594201756465864D-3,.43594201756465864D-3,.38596016222135642D-4, &.43594201756465864D-3,.38596016222135642D-4,.70856790445854181D-3, &.45980607095937280D-4,.45980607095937280D-4,.17945553213931007D-3, &.45980607095937280D-4,.17945553213931007D-3,.48983089863766725D-4, &.12243276984397130D-3,.12243276984397130D-3,.37694197830374995D-4, &.12243276984397130D-3,.37694197830374995D-4,.12635486428370778D-2, &.39802753143808619D-4,.39802753143808619D-4,.45520375316940385D-3, &.39802753143808619D-4,.45520375316940385D-3,.69605892788794196D-4, &.77195571056441645D-4,.77195571056441645D-4,.38714471225651343D-4, &.77195571056441645D-4,.38714471225651343D-4,.67088302517673629D-3, &.43531228608653492D-4,.43531228608653492D-4,.22995399918782116D-3, &.43531228608653492D-4,.22995399918782116D-3,.56614033067770802D-4, &.10354919481243404D-3,.10354919481243404D-3,.36709498005207836D-4, &.10354919481243404D-3,.36709498005207836D-4,.28937703030015149D-2, &.48618783763881142D-4,.48618783763881142D-4,.14809267754754818D-3, &.48618783763881142D-4,.14809267754754818D-3,.62712134095548161D-4, &.87913375680110403D-4,.87913375680110403D-4,.41545829033424093D-4, &.87913375680110403D-4,.41545829033424093D-4,.31873896959308106D-3, &.37780233067604510D-4,.37780233067604510D-4,.11739344185279312D-2, &.37780233067604510D-4,.11739344185279312D-2,.62035384410473408D-4, &.30760668684523079D-3,.30760668684523079D-3,.74163507005584869D-4, &.30760668684523079D-3,.74163507005584869D-4,.13009337335596148D-3, &.36599154808666998D-4,.36599154808666998D-4,.54545360898317646D-4, &.36599154808666998D-4,.54545360898317646D-4,.49355646840265301D-4, &-.45363872495324625D0,.53279315738121825D+0,-.28639465815156346D0, &.53279315738121825D0,-.28639465815156346D+0,.40794111347048040D-4/ DATA ((COEF(I,J),J=1,6),I=293,310)/ &.2854939231264389D-4,.2854939231264389D-4,.6395210936050434D-4, &.2854939231264389D-4,.6395210936050434D-4,.2051785507702231D-4, &.5192983657517473D-3,.5192983657517473D-3,.2453709186636645D-4, &.5192983657517473D-3,.2453709186636645D-4,.1009174410957348D-3, &.2201130479865847D-4,.2201130479865847D-4,.1911126917488530D-3, &.2201130479865847D-4,.1911126917488530D-3,.3491779138465434D-4, &.4540290894573798D-4,.4540290894573798D-4,.1982169078787312D-4, &.4540290894573798D-4,.1982169078787312D-4,.4675003044911219D-2, &.2567778799291618D-4,.2567778799291618D-4,.8532749792941040D-4, &.2567778799291618D-4,.8532749792941040D-4,.2141378101648924D-4, &.2522185292731416D-3,.2522185292731416D-3,.3245142945627677D-4, &.2522185292731416D-3,.3245142945627677D-4,.5038181206330221D-4, &.1997363804031099D-4,.1997363804031099D-4,.1673041031667479D-2, &.1997363804031099D-4,.1673041031667479D-2,.2700505422077120D-4, &.7334814066832329D-4,.7334814066832329D-4,.2091875954309112D-4, &.7334814066832329D-4,.2091875954309112D-4,.3496830547243688D-3, &.3034888539262456D-4,.3034888539262456D-4,.5645387202965999D-4, &.3034888539262456D-4,.5645387202965999D-4,.2020448595755236D-4, &.8548820219728689D-3,.8548820219728689D-3,.2355839650054794D-4, &.8548820219728689D-3,.2355839650054794D-4,.1217136553896909D-3, &.2272183313047774D-4,.2272183313047774D-4,.1503040207111383D-3, &.2272183313047774D-4,.1503040207111383D-3,.3782585931086922D-4, &.4127664462817298D-4,.4127664462817298D-4,.1974629500149606D-4, &.4127664462817298D-4,.1974629500149606D-4,.4762481045181510D-1, &.2749400836943466D-4,.2749400836943466D-4,.6996845413363466D-4, &.2749400836943466D-4,.6996845413363466D-4,.2077502684646142D-4, &.3954133908703511D-3,.3954133908703511D-3,.2971797853369174D-4, &.3954133908703511D-3,.2971797853369174D-4,.5877494732055006D-4, &.2029955580276558D-4,.2029955580276558D-4,.7134956235495635D-3, &.2029955580276558D-4,.7134956235495635D-3,.2527814796849340D-4, &.9005880719852324D-4,.9005880719852324D-4,.2160109487553139D-4, &.9005880719852324D-4,.2160109487553139D-4,.2288452406575260D-3, &.3322968059810997D-4,.3322968059810997D-4,.4861416174116007D-4, &.3322968059810997D-4,.4861416174116007D-4,.1991436345914009D-4, &.2228708701482198D-2,.2228708701482198D-2,.2386794010935285D-4, &.2228708701482198D-2,.2386794010935285D-4,.1140704806409944D-3/ DATA ((COEF(I,J),J=1,6),I=311,330)/ &.2247168724990004D-4,.2247168724990004D-4,.1622514229624652D-3, &.2247168724990004D-4,.1622514229624652D-3,.3680174637706332D-4, &.4256932689326857D-4,.4256932689326857D-4,.1976301652904962D-4, &.4256932689326857D-4,.1976301652904962D-4,.1561203063451092D-1, &.2654015737880450D-4,.2654015737880450D-4,.7701214865746648D-4, &.2654015737880450D-4,.7701214865746648D-4,.2107293492156357D-4, &.3115769717259390D-3,.3115769717259390D-3,.3101335989061020D-4, &.3115769717259390D-3,.3101335989061020D-4,.5429014967269275D-4, &.2011857371932694D-4,.2011857371932694D-4,.1043408875341346D-2, &.2011857371932694D-4,.1043408875341346D-2,.2419392786277494D-4, &.1071692960215918D-3,.1071692960215918D-3,.2223500461256568D-4, &.1071692960215918D-3,.2223500461256568D-4,.1757597481708787D-3, &.3583362589366815D-4,.3583362589366815D-4,.4394258643733910D-4, &.3583362589366815D-4,.4394258643733910D-4,.1978813424166864D-4, &.7792331902620007D-2,.7792331902620007D-2,.2489821758432978D-4, &.7792331902620007D-2,.2489821758432978D-4,.9523629672676908D-4, &.2180014037073642D-4,.2180014037073642D-4,.2086615007631790D-3, &.2180014037073642D-4,.2086615007631790D-3,.3405084772151086D-4, &.4695745804197193D-4,.4695745804197193D-4,.1986374327159400D-4, &.4695745804197193D-4,.1986374327159400D-4,.3118002553877718D-2, &.2609810963864186D-4,.2609810963864186D-4,.8099287864400658D-4, &.2609810963864186D-4,.8099287864400658D-4,.2123783898272738D-4, &.2794900573031592D-3,.2794900573031592D-3,.3171346342409830D-4, &.2794900573031592D-3,.3171346342409830D-4,.5227018014573681D-4, &.2004166894529176D-4,.2004166894529176D-4,.1302653253130018D-2, &.2004166894529176D-4,.1302653253130018D-2,.2800832058428110D-4, &.6684474030911918D-4,.6684474030911918D-4,.2064147295859748D-4, &.6684474030911918D-4,.2064147295859748D-4,.4509463765114682D-3, &.2911873469281630D-4,.2911873469281630D-4,.6126863629398420D-4, &.2911873469281630D-4,.6126863629398420D-4,.2040395050934914D-4, &.6047427249358356D-3,.6047427249358356D-3,.2326461452656810D-4, &.6047427249358356D-3,.2326461452656810D-4,.1302086423226916D-3, &.2298595784140845D-4,.2298595784140845D-4,.1396863855598545D-3, &.2298595784140845D-4,.1396863855598545D-3,.3891001700544618D-4, &.4005865665173498D-4,.4005865665173498D-4,.1974258166862778D-4, &.4005865665173498D-4,.1974258166862778D-4,.8891921904368920D-1, &.2540912478004451D-4,.2540912478004451D-4,.8843472319843332D-4, &.2540912478004451D-4,.8843472319843332D-4,.2153737172612036D-4, &.2362505573914286D-3,.2362505573914286D-3,.3296562953781371D-4, &.2362505573914286D-3,.3296562953781371D-4,.4919058993052970D-4/ DATA ((COEF(I,J),J=1,6),I=331,350)/ &.1993315533691752D-4,.1993315533691752D-4,.2016000558169837D-2, &.1993315533691752D-4,.2016000558169837D-2,.2732827231549935D-4, &.7106525136178148D-4,.7106525136178148D-4,.2082179417909269D-4, &.7106525136178148D-4,.2082179417909269D-4,.3792025374418422D-3, &.2992467172057884D-4,.2992467172057884D-4,.5798286660460991D-4, &.2992467172057884D-4,.5798286660460991D-4,.2026683837686716D-4, &.7564156416098478D-3,.7564156416098478D-3,.2430636973929888D-4, &.7564156416098478D-3,.2430636973929888D-4,.1050173417816384D-3, &.2215901822636955D-4,.2215901822636955D-4,.1806566513351597D-3, &.2215901822636955D-4,.1806566513351597D-3,.3552271192937976D-4, &.4441934791292662D-4,.4441934791292662D-4,.1979837958814266D-4, &.4441934791292662D-4,.1979837958814266D-4,.6472868964056924D-2, &.2669252635982589D-4,.2669252635982589D-4,.7575731199074937D-4, &.2669252635982589D-4,.7575731199074937D-4,.2102036564396927D-4, &.3235405491587992D-3,.3235405491587992D-3,.3078801479336957D-4, &.3235405491587992D-3,.3078801479336957D-4,.5499477801229046D-4, &.2014620323386758D-4,.2014620323386758D-4,.9741216171055970D-3, &.2014620323386758D-4,.9741216171055970D-3,.2376296646961465D-4, &.1165307259924030D-3,.1165307259924030D-3,.2255354857872450D-4, &.1165307259924030D-3,.2255354857872450D-4,.1581078941528752D-3, &.3713670332761178D-4,.3713670332761178D-4,.4212978202468662D-4, &.3713670332761178D-4,.4212978202468662D-4,.1975651147000684D-4, &.2111621489849421D-1,.2111621489849421D-1,.2595564154722690D-4, &.2111621489849421D-1,.2595564154722690D-4,.8239644225017478D-4, &.2129524391762934D-4,.2129524391762934D-4,.2699148191260962D-3, &.2129524391762934D-4,.2699148191260962D-3,.3195512958833712D-4, &.5162666369478089D-4,.5162666369478089D-4,.2001801212202455D-4, &.5162666369478089D-4,.2001801212202455D-4,.1411039637077915D-2, &.2477579579571362D-4,.2477579579571362D-4,.9707095281092766D-4, &.2477579579571362D-4,.9707095281092766D-4,.2186916201648878D-4, &.2025428832074265D-3,.2025428832074265D-3,.3433459433490981D-4, &.2025428832074265D-3,.3433459433490981D-4,.4642833276575716D-4, &.1984877765826999D-4,.1984877765826999D-4,.3536595766607164D-2, &.1984877765826999D-4,.3536595766607164D-2,.2818563227211816D-4, &.6585599355623209D-4,.6585599355623209D-4,.2059917439958627D-4, &.6585599355623209D-4,.2059917439958627D-4,.4721260125077478D-3, &.2892571462762092D-4,.2892571462762092D-4,.6214113188896829D-4, &.2892571462762092D-4,.6214113188896829D-4,.2044085238718792D-4, &.5740284523167527D-3,.5740284523167527D-3,.2336083075337178D-4, &.5740284523167527D-3,.2336083075337178D-4,.1272747360373466D-3/ DATA ((COEF(I,J),J=1,6),I=351,370)/ &.2289632972349478D-4,.2289632972349478D-4,.1430903119882924D-3, &.2289632972349478D-4,.1430903119882924D-3,.2581554721843980D-4, &.8384078838081793D-4,.8384078838081793D-4,.2135388714249066D-4, &.8384078838081793D-4,.2135388714249066D-4,.2608355278136665D-3, &.3220108838311426D-4,.3220108838311426D-4,.5099734182177761D-4, &.3220108838311426D-4,.5099734182177761D-4,.1999533660427757D-4, &.1533635221452478D-2,.1533635221452478D-2,.2502272827030977D-4, &.1533635221452478D-2,.2502272827030977D-4,.9345752384129992D-4, &.2173246503325741D-4,.2173246503325741D-4,.2150719485912417D-3, &.2173246503325741D-4,.2150719485912417D-3,.3377218530504890D-4, &.4749792673493065D-4,.4749792673493065D-4,.1987966072924208D-4, &.4749792673493065D-4,.1987966072924208D-4,.2769856585202598D-2, &.2783398127395672D-4,.2783398127395672D-4,.6785913862031401D-4, &.2783398127395672D-4,.6785913862031401D-4,.2068487535288822D-4, &.4311809975149026D-3,.4311809975149026D-3,.2931507606773188D-4, &.4311809975149026D-3,.2931507606773188D-4,.6041720974299877D-4, &.2036810513094914D-4,.2036810513094914D-4,.6380221452089648D-3, &.2036810513094914D-4,.6380221452089648D-3,.2408339401971605D-4, &.1093932831160238D-3,.1093932831160238D-3,.2231243326343716D-4, &.1093932831160238D-3,.2231243326343716D-4,.1710672521125244D-3, &.3615034685234670D-4,.3615034685234670D-4,.4347549305662360D-4, &.3615034685234670D-4,.4347549305662360D-4,.1977882642445252D-4, &.9565881508594270D-2,.9565881508594270D-2,.2684747706778030D-4, &.9565881508594270D-2,.2684747706778030D-4,.7453637502627767D-4, &.2096897738042850D-4,.2096897738042850D-4,.3362229901262604D-3, &.2096897738042850D-4,.3362229901262604D-3,.3854169394789898D-4, &.4045669366319275D-4,.4045669366319275D-4,.1973979727999000D-4, &.4045669366319275D-4,.1973979727999000D-4,.4386222444610512D+0, &.3056654872244681D-4,.3056654872244681D-4,.5571585118489064D-4, &.3056654872244681D-4,.5571585118489064D-4,.2017483879376742D-4, &.9115606298064600D-3,.9115606298064600D-3,.2345875022619936D-4, &.9115606298064600D-3,.2345875022619936D-4,.1244446834401361D-3, &.2280829497854017D-4,.2280829497854017D-4,.1466271911395111D-3, &.2280829497854017D-4,.1466271911395111D-3,.3818036276923079D-4, &.4086262373442708D-4,.4086262373442708D-4,.1975093797731458D-4, &.4086262373442708D-4,.1975093797731458D-4,.3029017035433408D-1, &.2624299109483879D-4,.2624299109483879D-4,.7962858990618148D-4, &.2624299109483879D-4,.7962858990618148D-4,.2118166028719211D-4, &.2895980169015215D-3,.2895980169015215D-3,.3147600264954836D-4, &.2895980169015215D-3,.3147600264954836D-4,.5292830857310305D-4/ DATA ((COEF(I,J),J=1,6),I=371,390)/ &.2006631162197865D-4,.2006631162197865D-4,.1206359857654832D-02, &.2006631162197865D-4,.1206359857654832D-2,.2465543095699396D-04, &.9896384244937268D-4,.9896384244937268D-4,.2193954502023498D-04, &.9896384244937268D-4,.2193954502023498D-4,.1966986633729502D-03, &.3462353714478440D-4,.3462353714478440D-4,.4591024747423779D-04, &.3462353714478440D-4,.4591024747423779D-4,.1983476104934234D-04, &.4046094940818183D-2,.4046094940818183D-2,.2766256115849852D-04, &.4046094940818183D-2,.2766256115849852D-4,.6890006995395185D-04, &.2072939031926428D-4,.2072939031926428D-4,.4127066876776266D-03, &.2072939031926428D-4,.4127066876776266D-3,.2951480228392970D-04, &.5958619166455701D-4,.5958619166455701D-4,.2033330917839852D-04, &.5958619166455701D-4,.2033330917839852D-4,.6741606675715499D-03, &.2365979368802508D-4,.2365979368802508D-4,.1190771004598618D-03, &.2365979368802508D-4,.1190771004598618D-3,.2263692417691780D-04, &.1541282652727628D-3,.1541282652727628D-3,.3647300796836789D-04, &.1541282652727628D-3,.3647300796836789D-4,.4301781984257157D-04, &.1977045437989326D-4,.1977045437989326D-4,.1203189847779800D-01, &.1977045437989326D-4,.1203189847779800D-1,.2716530364689954D-04, &.7219146436412886D-4,.7219146436412886D-4,.2086970182172562D-04, &.7219146436412886D-4,.2086970182172562D-4,.3639856346175968D-03, &.3013495051182238D-4,.3013495051182238D-4,.5720936290522985D-04, &.3013495051182238D-4,.5720936290522985D-4,.2023515049246316D-04, &.8033703081975564D-3,.8033703081975564D-3,.2397474046842758D-04, &.8033703081975564D-3,.2397474046842758D-4,.1116925382739447D-03, &.2239132154776664D-4,.2239132154776664D-4,.1665679591598375D-03, &.2239132154776664D-4,.1665679591598375D-3,.3521747563310610D-04, &.4490603376079471D-4,.4490603376079471D-4,.1980956439884101D-04, &.4490603376079471D-4,.1980956439884101D-4,.5464009038574440D-02, &.2639032648639676D-4,.2639032648639676D-4,.7830213758507075D-04, &.2639032648639676D-4,.7830213758507075D-4,.2112669608345274D-04, &.3002789467902159D-3,.3002789467902159D-3,.3124266235529018D-04, &.3002789467902159D-3,.3124266235529018D-4,.5360148197894120D-04, &.2009194489801654D-4,.2009194489801654D-4,.1120423413505869D-02, &.2009194489801654D-4,.1120423413505869D-2,.2442074804906302D-04, &.1029343573766687D-3,.1029343573766687D-3,.2208445714291506D-04, &.1029343573766687D-3,.2208445714291506D-4,.1857699901750056D-03, &.3270624494111179D-4,.3270624494111179D-4,.4977968625559511D-04, &.3270624494111179D-4,.4977968625559511D-4,.1832482820013657D-02, &.2514936056600602D-4,.2514936056600602D-4,.9173240138754972D-04, &.2514936056600602D-4,.9173240138754972D-4,.2166612131400504D-04/ DATA ((COEF(I,J),J=1,6),I=391,400)/ &.22179306659979075D-3,.22179306659979075D-3,.33498498213184707D-4, &.22179306659979075D-3,.33498498213184707D-4,.48050052915205633D-4, &.19896533055280618D-4,.19896533055280618D-4,.24771592091872084D-2, &.19896533055280618D-4,.24771592091872084D-2,.25542326129265426D-4, &.13325152973885215D-3,.13325152973885215D-3,.21474938213123001D-4, &.13325152973885215D-3,.21474938213123001D-4,.24403297737851022D-3, &.23077200292764951D-4,.23077200292764951D-4,.13640884641872292D-3, &.23077200292764951D-4,.13640884641872292D-3,.39668321070634110D-4, &.64892053354505223D-4,.64892053354505223D-4,.20557971187156988D-4, &.64892053354505223D-4,.20557971187156988D-4,.54562222154862965D-3, &.28735953952505441D-4,.28735953952505441D-4,.63035384203333504D-4, &.28735953952505441D-4,.63035384203333504D-4,.20478818065726965D-4, &.23170078533532433D-4,.23170078533532433D-4,.49485815205661701D-3, &.23170078533532433D-4,.49485815205661701D-3,.19739797279989982D-4, &.37478024330443648D-4,.37478024330443648D-4,.22391321547766630D-4, &.37478024330443648D-4,.22391321547766630D-4,.41698960476035002D-4, &.19737013414938452D-4,.19737013414938452D-4,.28365970872542253D-4, &.19737013414938452D-4,.28365970872542253D-4,.86858229336522362D-4, &-.48121500846540849D0,.53017893189276518D+0,-.27846926883426556D0, &.53017893189276518D+0,-.2784692688342656D0,.19952912272757273D-4/ DATA ((COEF(I,J),J=1,6),I=401,420)/ &.16007316152806260D-4,.16007316152806260D-4,.37898134658626442D-4, &.16007316152806260D-4,.37898134658626442D-4,.11760217815002902D-4, &.26140255127179914D-3,.26140255127179914D-3,.16561221964136950D-4, &.26140255127179914D-3,.16561221964136950D-4,.35117366401543006D-4, &.11641727695428183D-4,.11641727695428183D-4,.33783187149825353D-3, &.11641727695428183D-4,.33783187149825353D-3,.14789203495903856D-4, &.47078621823147551D-4,.47078621823147551D-4,.12146232270348437D-4, &.47078621823147551D-4,.12146232270348437D-4,.15318863871675953D-3, &.18237273904380468D-4,.18237273904380468D-4,.29389980864747500D-4, &.18237273904380468D-4,.29389980864747500D-4,.11412213439691687D-4, &.81133763457012190D-3,.81133763457012190D-3,.14170050502889764D-4, &.81133763457012190D-3,.14170050502889764D-4,.54684915429518837D-4, &.12444631715316333D-4,.12444631715316333D-4,.11761906221029793D-3, &.12444631715316333D-4,.11761906221029793D-3,.19475067539216037D-4, &.26659377274690299D-4,.26659377274690299D-4,.11322638062013902D-4, &.26659377274690299D-4,.11322638062013902D-4,.18542334934794633D-2, &.15857898364668542D-4,.15857898364668542D-4,.38762846995482153D-4, &.15857898364668542D-4,.38762846995482153D-4,.11797219507375151D-4, &.24436615992378031D-3,.24436615992378031D-3,.16728783227142380D-4, &.24436615992378031D-3,.16728783227142380D-4,.34387012618777227D-4, &.11610958137780429D-4,.11610958137780429D-4,.36597610429905524D-3, &.11610958137780429D-4,.36597610429905524D-3,.13411031841498931D-4, &.69966868087496798D-4,.69966868087496798D-4,.21642378117264469D-4, &.69966868087496798D-4,.21642378117264469D-4,.84890392664090298D-4, &.12973860356423407D-4,.12973860356423407D-4,.23445386697297003D-4, &.12973860356423407D-4,.23445386697297003D-4,.11258080444400271D-4, &.30587620489677170D-1,.30587620489677170D-1,.15296621771754272D-4, &.30587620489677170D-1,.15296621771754272D-4,.42582100866682079D-4, &.11959692184299121D-4,.11959692184299121D-4,.19070269751147130D-3, &.11959692184299121D-4,.19070269751147130D-3,.17443846819263622D-4, &.31714666837516854D-4,.31714666837516854D-4,.11501204764404463D-4, &.31714666837516854D-4,.11501204764404463D-4,.52341179633617668D-3, &.13771047638167434D-4,.13771047638167434D-4,.61569376409481674D-4, &.13771047638167434D-4,.61569376409481674D-4,.20498506642835584D-4, &.99178293579936950D-4,.99178293579936950D-4,.12694365731300556D-4, &.99178293579936950D-4,.12694365731300556D-4,.24953887848146697D-4, &.11280739060005732D-4,.11280739060005732D-4,.47335914832266311D-2, &.11280739060005732D-4,.47335914832266311D-2,.14729297060972596D-4, &.47696145622415329D-4,.47696145622415329D-4,.12171309316668505D-4, &.47696145622415329D-4,.12171309316668505D-4,.14930883616158901D-3/ DATA ((COEF(I,J),J=1,6),I=421,440)/ &.18342447724854856D-4,.18342447724854856D-4,.29120892803552926D-4, &.18342447724854856D-4,.29120892803552926D-4,.11402516886760947D-4, &.86354542635543670D-3,.86354542635543670D-3,.14222844936537749D-4, &.86354542635543670D-3,.14222844936537749D-4,.77537235278883832D-4, &.12415405331967423D-4,.12415405331967423D-4,.12029551324102535D-3, &.12415405331967423D-4,.12029551324102535D-3,.19354828959913159D-4, &.26888034316789703D-4,.26888034316789703D-4,.11329241050089539D-4, &.26888034316789703D-4,.11329241050089539D-4,.16926722395854549D-2, &.15932133864980196D-4,.15932133864980196D-4,.38326326213263181D-4, &.15932133864980196D-4,.38326326213263181D-4,.11778541133267408D-4, &.25266572741798143D-3,.25266572741798143D-3,.16644466242546458D-4, &.25266572741798143D-3,.16644466242546458D-4,.34748851932579366D-4, &.11626173929208896D-4,.11626173929208896D-4,.35147704939134382D-3, &.11626173929208896D-4,.35147704939134382D-3,.13368606870758973D-4, &.71144768633401469D-4,.71144768633401469D-4,.21794687281656701D-4, &.71144768633401469D-4,.21794687281656701D-4,.83333915009290921D-4, &.13011000403904132D-4,.13011000403904132D-4,.23269225705938609D-4, &.13011000403904132D-4,.23269225705938609D-4,.11256591874385288D-4, &.47956953252117185D-1,.47956953252117185D-1,.14911304940714996D-4, &.47956953252117185D-1,.14911304940714996D-4,.45882618712941353D-4, &.12097268132291097D-4,.12097268132291097D-4,.16142922288697284D-3, &.12097268132291097D-4,.16142922288697284D-3,.22589976495968650D-4, &.29941861997064765D-4,.29941861997064765D-4,.11432551251897504D-4, &.29941861997064765D-4,.11432551251897504D-4,.72024662700823430D-3, &.14066452135306724D-4,.14066452135306724D-4,.56284674769103852D-4, &.14066452135306724D-4,.56284674769103852D-4,.12504388836955019D-4, &.11253612378299275D-3,.11253612378299275D-3,.19720503163723066D-4, &.11253612378299275D-3,.19720503163723066D-4,.26212792485563862D-4, &.11310346181792830D-4,.11310346181792830D-4,.22557019150327750D-2, &.11310346181792830D-4,.22557019150327750D-2,.15712217119245395D-4, &.39661739950219036D-4,.39661739950219036D-4,.11835649934469707D-4, &.39661739950219036D-4,.11835649934469707D-4,.22896631865186663D-3, &.16900697128642001D-4,.16900697128642001D-4,.33682735102511740D-4, &.16900697128642001D-4,.33682735102511740D-4,.11581533768627200D-4, &.39783487136906817D-3,.39783487136906817D-3,.13497563473985797D-4, &.39783487136906817D-3,.13497563473985797D-4,.67702462409557527D-4, &.21344194299835096D-4,.21344194299835096D-4,.88144191242198688D-4, &.21344194299835096D-4,.88144191242198688D-4,.12901076747669186D-4, &.23805664102448591D-4,.23805664102448591D-4,.11261951955008705D-4, &.23805664102448591D-4,.11261951955008705D-4,.15794640167391974D-1/ DATA ((COEF(I,J),J=1,6),I=441,460)/ &.15100323172298703D-4,.15100323172298703D-4,.44180912121592551D-4, &.15100323172298703D-4,.44180912121592551D-4,.12026753880090559D-4, &.17513132618029667D-3,.17513132618029667D-3,.17731690947072178D-4, &.17513132618029667D-3,.17731690947072178D-4,.30805480097326241D-4, &.11465436502139331D-4,.11465436502139331D-4,.61001843336396552D-3, &.11465436502139331D-4,.61001843336396552D-3,.13915918251617916D-4, &.58831124719412376D-4,.58831124719412376D-4,.20101468547197988D-4, &.58831124719412376D-4,.20101468547197988D-4,.10552835978900794D-3, &.12597341032982432D-4,.12597341032982432D-4,.25568680553915424D-4, &.12597341032982432D-4,.25568680553915424D-4,.11294184014571873D-4, &.31566355704869074D-2,.31566355704869074D-2,.15500485957318728D-4, &.31566355704869074D-2,.15500485957318728D-4,.41078135266080606D-4, &.11896012896407498D-4,.11896012896407498D-4,.20849935404662715D-3, &.11896012896407498D-4,.20849935404662715D-3,.17167007650242880D-4, &.32672598443218864D-4,.32672598443218864D-4,.11539891138007469D-4, &.32672598443218864D-4,.11539891138007469D-4,.45412749093781419D-3, &.13631654724944208D-4,.13631654724944208D-4,.64519077795486060D-4, &.13631654724944208D-4,.64519077795486060D-4,.20912451573097704D-4, &.93406229264489778D-4,.93406229264489778D-4,.12795572461485830D-4, &.93406229264489778D-4,.12795572461485830D-4,.24366735284720006D-4, &.11269998458881092D-4,.11269998458881092D-4,.78892160406046560D-2, &.11269998458881092D-4,.78892160406046560D-2,.14554043818743806D-4, &.49631409827654899D-4,.49631409827654899D-4,.12248956208453310D-4, &.49631409827654899D-4,.12248956208453310D-4,.13853586548275427D-3, &.18666467952404610D-4,.18666467952404610D-4,.28339884858442229D-4, &.18666467952404610D-4,.28339884858442229D-4,.11375304184904433D-4, &.10546138155858938D-2,.10546138155858938D-2,.14385302964692950D-4, &.10546138155858938D-2,.14385302964692950D-4,.51699455564304448D-4, &.12330292166756799D-4,.12330292166756799D-4,.12891785026199358D-3, &.12330292166756799D-4,.12891785026199358D-3,.19003713111479593D-4, &.27596360831608721D-4,.27596360831608721D-4,.11350886768711248D-4, &.27596360831608721D-4,.11350886768711248D-4,.13173387505836570D-2, &.16160573146372408D-4,.16160573146372408D-4,.37065915726437290D-4, &.16160573146372408D-4,.37065915726437290D-4,.11724628104735532D-4, &.28031586813720758D-3,.28031586813720758D-3,.16397891476262916D-4, &.28031586813720758D-3,.16397891476262916D-4,.35875073960823866D-4, &.11673856108098277D-4,.11673856108098277D-4,.31284629386502040D-3, &.11673856108098277D-4,.31284629386502040D-3,.13326734742444352D-4, &.72354598359952615D-4,.72354598359952615D-4,.21949192443455349D-4, &.72354598359952615D-4,.21949192443455349D-4,.81822002901579859D-4/ DATA ((COEF(I,J),J=1,6),I=461,480)/ &.13048645704248677D-4,.13048645704248677D-4,.23095651195431883D-4, &.13048645704248677D-4,.23095651195431883D-4,.11255401207321695D-4, &.88450005450447108D-1,.88450005450447108D-1,.22426283289812671D-4, &.88450005450447108D-1,.22426283289812671D-4,.43636956715395695D-4, &.12004019913628395D-4,.12004019913628395D-4,.18009808526378438D-3, &.12004019913628395D-4,.18009808526378438D-3,.17634489547962153D-4, &.31103327838724223D-4,.31103327838724223D-4,.11477037250100687D-4, &.31103327838724223D-4,.11477037250100687D-4,.57888671921843212D-3, &.13867008160790467D-4,.13867008160790467D-4,.59721642785232498D-4, &.13867008160790467D-4,.59721642785232498D-4,.20231988652192797D-4, &.10334279540396721D-3,.10334279540396721D-3,.12629224810200599D-4, &.10334279540396721D-3,.12629224810200599D-4,.25360577472314292D-4, &.11289401139380008D-4,.11289401139380008D-4,.35806469605106314D-2, &.11289401139380008D-4,.35806469605106314D-2,.15431672347673133D-4, &.41569382742129848D-4,.41569382742129848D-4,.11916867663722096D-4, &.41569382742129848D-4,.11916867663722096D-4,.20229451206493489D-3, &.17258094581805305D-4,.17258094581805305D-4,.32347664386630801D-4, &.17258094581805305D-4,.32347664386630801D-4,.11526669154415143D-4, &.47559629588350498D-3,.47559629588350498D-3,.13677521452935276D-4, &.47559629588350498D-3,.13677521452935276D-4,.63511034636790895D-4, &.20772537613658592D-4,.20772537613658592D-4,.95270441539504057D-4, &.20772537613658592D-4,.95270441539504057D-4,.12761365209932244D-4, &.24559480258621557D-4,.24559480258621557D-4,.11273278861947211D-4, &.24559480258621557D-4,.11273278861947211D-4,.65538177522978689D-2, &.14849869633331065D-4,.14849869633331065D-4,.46474235100298159D-4, &.14849869633331065D-4,.46474235100298159D-4,.12121552847836862D-4, &.15722589047113722D-3,.15722589047113722D-3,.18133480237631760D-4, &.15722589047113722D-3,.18133480237631760D-4,.29663604111381711D-4, &.11422224534057410D-4,.11422224534057410D-4,.76374780781249111D-3, &.11422224534057410D-4,.76374780781249111D-3,.14117922115315409D-4, &.55475412387603307D-4,.55475412387603307D-4,.12474291669763599D-4, &.55475412387603307D-4,.12474291669763599D-4,.11503396795255235D-3, &.19596949898135698D-4,.19596949898135698D-4,.26434320778341245D-4, &.19596949898135698D-4,.26434320778341245D-4,.11316339996595517D-4, &.20401837230112729D-2,.20401837230112729D-2,.15784596940989115D-4, &.20401837230112729D-2,.15784596940989115D-4,.39207911192319167D-4, &.11816255053780308D-4,.11816255053780308D-4,.23647515072161961D-3, &.11816255053780308D-4,.23647515072161961D-3,.16814188296198815D-4, &.34031691214172833D-4,.34031691214172833D-4,.11596078633360840D-4, &.34031691214172833D-4,.11596078633360840D-4,.38140181157899950D-3/ DATA ((COEF(I,J),J=1,6),I=481,500)/ &.13454015420147789D-4,.13454015420147789D-4,.68819786724504369D-4, &.13454015420147789D-4,.68819786724504369D-4,.21492226363361663D-4, &.86493200455835064D-4,.86493200455835064D-4,.12937220722645432D-4, &.86493200455835064D-4,.12937220722645432D-4,.23624182870468202D-4, &.11259867074901479D-4,.11259867074901479D-4,.21349492877649794D-1, &.11259867074901479D-4,.21349492877649794D-1,.14611726272817009D-4, &.48972138468133110D-4,.48972138468133110D-4,.12222668227838749D-4, &.48972138468133110D-4,.12222668227838749D-4,.14198977912752294D-3, &.18557022513661480D-4,.18557022513661480D-4,.28595937728523175D-4, &.18557022513661480D-4,.28595937728523175D-4,.11384063253525696D-4, &.98439355699152763D-3,.98439355699152763D-3,.14330463181987305D-4, &.98439355699152763D-3,.14330463181987305D-4,.52420401678078233D-4, &.12358239037847129D-4,.12358239037847129D-4,.12594015297857897D-3, &.12358239037847129D-4,.12594015297857897D-3,.19119177473405337D-4, &.27356448911004784D-4,.27356448911004784D-4,.11343364648227503D-4, &.27356448911004784D-4,.11343364648227503D-4,.14271756314916167D-2, &.16083458180344656D-4,.16083458180344656D-4,.37478064842959685D-4, &.16083458180344656D-4,.37478064842959685D-4,.11742247481138901D-4, &.27060789505677980D-3,.27060789505677980D-3,.16479035307633176D-4, &.27060789505677980D-3,.16479035307633176D-4,.35492718070706667D-4, &.11657621165083480D-4,.11657621165083480D-4,.32497480394128410D-3, &.11657621165083480D-4,.32497480394128410D-3,.15036521803337257D-4, &.44736293155516890D-4,.44736293155516890D-4,.12049871429151540D-4, &.44736293155516890D-4,.12049871429151540D-4,.17037170560559931D-3, &.17830171035396909D-4,.17830171035396909D-4,.30512693430340267D-4, &.17830171035396909D-4,.30512693430340267D-4,.11454156064390183D-4, &.64374974867375406D-3,.64374974867375406D-3,.13965457511795541D-4, &.64374974867375406D-3,.13965457511795541D-4,.57961899637490740D-4, &.13164662808033991D-4,.13164662808033991D-4,.10778688525780870D-3, &.13164662808033991D-4,.10778688525780870D-3,.19972733866662441D-4, &.25780040245567086D-4,.25780040245567086D-4,.11299268780647165D-4, &.25780040245567086D-4,.11299268780647165D-4,.28039536434331432D-2, &.15570173533831119D-4,.15570173533831119D-4,.40596612182183847D-4, &.15570173533831119D-4,.40596612182183847D-4,.11875526577528947D-4, &.21499994440361605D-3,.21499994440361605D-3,.17077090343869544D-4, &.21499994440361605D-3,.17077090343869544D-4,.33003330473576275D-4, &.11553441568870537D-4,.11553441568870537D-4,.43409199560445663D-3, &.11553441568870537D-4,.43409199560445663D-3,.13586376996307400D-4, &.65553061112370861D-4,.65553061112370861D-4,.21054342559815009D-4, &.65553061112370861D-4,.21054342559815009D-4,.91598607063187845D-4/ DATA ((COEF(I,J),J=1,6),I=501,520)/ &.12830257167514376D-4,.12830257167514376D-4,.24176887336394940D-4, &.12830257167514376D-4,.24176887336394940D-4,.11267017378707181D-4, &.96835365956925070D-2,.96835365956925070D-2,.14670141030560282D-4, &.96835365956925070D-2,.14670141030560282D-4,.48327186533404357D-4, &.12196786959431372D-4,.12196786959431372D-4,.14557840502588769D-3, &.12196786959431372D-4,.14557840502588769D-3,.18449023259109498D-4, &.28856243152464556D-4,.28856243152464556D-4,.11393133828569922D-4, &.28856243152464556D-4,.11393133828569922D-4,.92098909263034392D-3, &.14276313191419222D-4,.14276313191419222D-4,.53157950321663077D-4, &.14276313191419222D-4,.53157950321663077D-4,.12386608943438503D-4, &.12306767231462230D-3,.12306767231462230D-3,.19236207604298414D-4, &.12306767231462230D-3,.19236207604298414D-4,.27120366135768555D-4, &.11336149666001600D-4,.11336149666001600D-4,.15514086287870884D-2, &.11336149666001600D-4,.15514086287870884D-2,.13285409791087234D-4, &.73597517739401692D-4,.73597517739401692D-4,.22105933080185416D-4, &.73597517739401692D-4,.22105933080185416D-4,.80352975020505813D-4, &.13086801182298633D-4,.13086801182298633D-4,.22924615607925218D-4, &.13086801182298633D-4,.22924615607925218D-4,.11254508317224070D-4, &.43739131907417234D+0,.43739131907417234D+0,.14973519058439783D-4, &.43739131907417234D+0,.14973519058439783D-4,.45303418660508086D-4, &.12073375261901773D-4,.12073375261901773D-4,.16580786728958996D-3, &.12073375261901773D-4,.16580786728958996D-3,.17929949294045827D-4, &.30224856742924364D-4,.30224856742924364D-4,.11443194710845914D-4, &.30224856742924364D-4,.11443194710845914D-4,.68037853892465115D-3, &.14015633058427847D-4,.14015633058427847D-4,.57113297580298551D-4, &.14015633058427847D-4,.57113297580298551D-4,.12534926924897215D-4, &.11012165128645170D-3,.11012165128645170D-3,.19845755048735988D-4, &.11012165128645170D-3,.19845755048735988D-4,.25994721885928634D-4, &.11304655978549092D-4,.11304655978549092D-4,.25074219696264405D-2, &.11304655978549092D-4,.25074219696264405D-2,.15640746655926177D-4, &.40124561647965478D-4,.40124561647965478D-4,.11855406357496826D-4, &.40124561647965478D-4,.11855406357496826D-4,.22181537852856433D-3, &.16988325710888718D-4,.16988325710888718D-4,.33339996134818099D-4, &.16988325710888718D-4,.33339996134818099D-4,.11567321936322081D-4, &.41536496754785151D-3,.41536496754785151D-3,.13541681973911795D-4, &.41536496754785151D-3,.13541681973911795D-4,.66613878415321134D-4, &.12865423759464592D-4,.12865423759464592D-4,.89845311500692961D-4, &.12865423759464592D-4,.89845311500692961D-4,.21198245050312666D-4, &.23989881480309603D-4,.23989881480309603D-4,.11264335305469155D-4, &.23989881480309603D-4,.11264335305469155D-4,.12177163761232982D-1/ DATA ((COEF(I,J),J=1,6),I=521,540)/ &.15230362692577967D-4,.15230362692577967D-4,.43104118871167988D-4, &.15230362692577967D-4,.43104118871167988D-4,.11981666880841925D-4, &.18528417483361254D-3,.18528417483361254D-3,.17538547744014483D-4, &.18528417483361254D-3,.17538547744014483D-4,.31406350910385002D-4, &.11488959570965605D-4,.11488959570965605D-4,.55009406795549171D-3, &.11488959570965605D-4,.55009406795549171D-3,.13818720246547420D-4, &.60634151262476702D-4,.60634151262476702D-4,.12661565024468307D-4, &.60634151262476702D-4,.12661565024468307D-4,.10122708715214158D-3, &.20364324392084799D-4,.20364324392084799D-4,.25155667284389784D-4, &.20364324392084799D-4,.25155667284389784D-4,.11284919646577782D-4, &.40966882697210921D-2,.40966882697210921D-2,.15363721338586612D-4, &.40966882697210921D-2,.15363721338586612D-4,.42070614909451986D-4, &.11938093277398216D-4,.11938093277398216D-4,.19636774943497851D-3, &.11938093277398216D-4,.19636774943497851D-3,.17350368427104985D-4, &.76187536575496499D-4,.76187536575496499D-4,.11513774167370870D-4, &.76187536575496499D-4,.11513774167370870D-4,.49863851220432727D-3, &.13723983584732002D-4,.13723983584732002D-4,.62528074757355629D-4, &.13723983584732002D-4,.62528074757355629D-4,.20634566963387381D-4, &.97193626506792791D-4,.97193626506792791D-4,.12727631059022767D-4, &.97193626506792791D-4,.12727631059022767D-4,.24755178549705670D-4, &.11276858935731043D-4,.11276858935731043D-4,.55325014130113443D-2, &.11276858935731043D-4,.55325014130113443D-2,.14497084858994598D-4, &.50305424146135986D-4,.50305424146135986D-4,.12275654046010378D-4, &.50305424146135986D-4,.12275654046010378D-4,.13521000922831801D-3, &.18777382507839015D-4,.18777382507839015D-4,.28087995298969676D-4, &.18777382507839015D-4,.28087995298969676D-4,.11253615521569498D-4, &.11326640136322700D-2,.11326640136322700D-2,.14440840740754619D-4, &.11326640136322700D-2,.14440840740754619D-4,.50994620768592278D-4, &.12302764944248121D-4,.12302764944248121D-4,.13200596154049282D-3, &.12302764944248121D-4,.13200596154049282D-3,.18889789592626019D-4, &.27840182160185107D-4,.27840182160185107D-4,.11358716833039650D-4, &.27840182160185107D-4,.11358716833039650D-4,.12197544065627674D-2, &.16238674501646657D-4,.16238674501646657D-4,.36661492479174181D-4, &.16238674501646657D-4,.36661492479174181D-4,.11707357702859559D-4, &.29056374235603041D-3,.29056374235603041D-3,.16317775954446859D-4, &.29056374235603041D-3,.16317775954446859D-4,.36264606252638589D-4, &.11690434336100448D-4,.11690434336100448D-4,.30139231062301542D-3, &.11690434336100448D-4,.30139231062301542D-3,.13244626449494689D-4, &.74874740497964123D-4,.74874740497964123D-4,.13125471850028978D-4, &.74874740497964123D-4,.13125471850028978D-4,.78925228643769764D-4/ DATA ((COEF(I,J),J=1,6),I=541,544)/ &.22264949584552506D-4,.22264949584552506D-4,.22756072496823693D-4, &.22264949584552506D-4,.22756072496823693D-4,.13204379247046970D-4, &.11253913109628138D-4,.11253913109628138D-4,.53912611607795680D-4, &.11253913109628138D-4,.53912611607795680D-4,.18031045600407457D-4, &.32028396527394753D-4,.32028396527394753D-4,.11366855680626436D-4, &.32028396527394753D-4,.11366855680626436D-4,.15164933345976474D-4, &-.48184039296639574D0,.52752694704749011D00,-.27658313468001291D0, &.52752694704749011D+0,-.27658313468001291D0,.12565909708973445D-4/ DATA ((COEF(I,J),J=1,6),I=545,564)/ &.94245302992258385D-5,.94245302992258385D-5,.27401646727681066D-4, &.94245302992258385D-5,.27401646727681066D-4,.74889408954686915D-5, &.11025345819531954D-3,.11025345819531954D-3,.11026340575782703D-4, &.11025345819531954D-3,.11026340575782703D-4,.19264734253559409D-4, &.71471972532000718D-5,.71471972532000718D-5,.37240434325882008D-3, &.71471972532000718D-5,.37240434325882008D-3,.86599896242439190D-5, &.36864101809909466D-4,.36864101809909466D-4,.78572452093085833D-5, &.36864101809909466D-4,.78572452093085833D-5,.65235099186273704D-4, &.12556156548832625D-4,.12556156548832625D-4,.15883596047313187D-4, &.12556156548832625D-4,.15883596047313187D-4,.70366583662371115D-5, &.20520433040591450D-2,.20520433040591450D-2,.96426986374640598D-5, &.20520433040591450D-2,.96426986374640598D-5,.25710356360477812D-4, &.74176592926136943D-5,.74176592926136943D-5,.12842208688423301D-3, &.74176592926136943D-5,.12842208688423301D-3,.10718480696933845D-4, &.20282562027137344D-4,.20282562027137344D-4,.71877909795840991D-5, &.20282562027137344D-4,.71877909795840991D-5,.28773361307177762D-3, &.85050076995070280D-5,.85050076995070280D-5,.39964109473645595D-4, &.85050076995070280D-5,.39964109473645595D-4,.79653228529724083D-5, &.58630217827105670D-4,.58630217827105670D-4,.12998286703993232D-4, &.58630217827105670D-4,.12998286703993232D-4,.15228343379816742D-4, &.70234720358766300D-5,.70234720358766300D-5,.45352858358460842D-2, &.70234720358766300D-5,.45352858358460842D-2,.91095621765011335D-5, &.30464138639694080D-4,.30464138639694080D-4,.76143305005617474D-5, &.30464138639694080D-4,.76143305005617474D-5,.88739567435258671D-4, &.11554962480698175D-4,.11554962480698175D-4,.17838877918461762D-4, &.11554962480698175D-4,.17838877918461762D-4,.70945000098789796D-5, &.60754176507620710D-3,.60754176507620710D-3,.89255350908901812D-5, &.60754176507620710D-3,.89255350908901812D-5,.32720221925833957D-4, &.77030330399604522D-5,.77030330399604522D-5,.78239630649879277D-4, &.77030330399604522D-5,.78239630649879277D-4,.11922772308954658D-4, &.17027905105850960D-4,.17027905105850960D-4,.70678713228235503D-5, &.17027905105850960D-4,.70678713228235503D-5,.89692453731226012D-3, &.99453082637208940D-5,.99453082637208940D-5,.23781034288673099D-4, &.99453082637208940D-5,.23781034288673099D-4,.73353059579299221D-5, &.15940637589098751D-3,.15940637589098751D-3,.10352148184514540D-4, &.15940637589098751D-3,.10352148184514540D-4,.21738217141065871D-4, &.72482793012295026D-5,.72482793012295026D-5,.21568173679355765D-3, &.72482793012295026D-5,.21568173679355765D-3,.83205975672288570D-5, &.44609868004199732D-4,.44609868004199732D-4,.81163486787489765D-5, &.44609868004199732D-4,.81163486787489765D-5,.51567315825099454D-4/ DATA ((COEF(I,J),J=1,6),I=565,584)/ &.13616795060014110D-4,.13616795060014110D-4,.14458907055078780D-4, &.13616795060014110D-4,.14458907055078780D-4,.70140733950457911D-5, &.33846660188522594D-1,.33846660188522594D-1,.93945250377549229D-5, &.33846660188522594D-1,.93945250377549229D-5,.27658489885916262D-4, &.74996549280403107D-5,.74996549280403107D-5,.10798235780023396D-3, &.74996549280403107D-5,.10798235780023396D-3,.11072063924395267D-4, &.19126733369429586D-4,.19126733369429586D-4,.71418506386781355D-5, &.19126733369429586D-4,.71418506386781355D-5,.38752052606405053D-3, &.86829938959941017D-5,.86829938959941017D-5,.36452991417792442D-4, &.86829938959941017D-5,.36452991417792442D-4,.78424505327552497D-5, &.66273110462759295D-4,.66273110462759295D-4,.12495579165627886D-4, &.66273110462759295D-4,.12495579165627886D-4,.15981605198605262D-4, &.70389641258598842D-5,.70389641258598842D-5,.18731075029499479D-2, &.70389641258598842D-5,.18731075029499479D-2,.96750614363991077D-5, &.25483080260037917D-4,.25483080260037917D-4,.74079991975094170D-5, &.25483080260037917D-4,.74079991975094170D-5,.13138835380406750D-3, &.10676186759525129D-4,.10676186759525129D-4,.20435756202142354D-4, &.10676186759525129D-4,.20435756202142354D-4,.71940468499965623D-5, &.27805440959107494D-3,.27805440959107494D-3,.84837103068016240D-5, &.27805440959107494D-3,.84837103068016240D-5,.40441154738071411D-4, &.79814206450391445D-5,.79814206450391445D-5,.57772537032910487D-4, &.79814206450391445D-5,.57772537032910487D-4,.13064132913252675D-4, &.15138948540196529D-4,.15138948540196529D-4,.70220089226168414D-5, &.15138948540196529D-4,.70220089226168414D-5,.52404590875400475D-2, &.91368617050579777D-5,.91368617050579777D-5,.30162754430056536D-4, &.91368617050579777D-5,.30162754430056536D-4,.76022350808836604D-5, &.90416096171038157D-4,.90416096171038157D-4,.11504490674281419D-4, &.90416096171038157D-4,.11504490674281419D-4,.17960524373824040D-4, &.70987390854364978D-5,.70987390854364978D-5,.57798468948443182D-3, &.70987390854364978D-5,.57798468948443182D-3,.89002286340518929D-5, &.33064872506312161D-4,.33064872506312161D-4,.77162915000931510D-5, &.33064872506312161D-4,.77162915000931510D-5,.76897600432806318D-4, &.11977464459378985D-4,.11977464459378985D-4,.16917577989426801D-4, &.11977464459378985D-4,.16917577989426801D-4,.70644994463284859D-5, &.95475427102895019D-3,.95475427102895019D-3,.99805553551910107D-5, &.95475427102895019D-3,.99805553551910107D-5,.23581893718276244D-4, &.73267840114268009D-5,.73267840114268009D-5,.16354005803943878D-3, &.73267840114268009D-5,.16354005803943878D-3,.10313369261739703D-4, &.21911220761294489D-4,.21911220761294489D-4,.72555878631595786D-5, &.21911220761294489D-4,.72555878631595786D-5,.20940512607146622D-3/ DATA ((COEF(I,J),J=1,6),I=585,604)/ &.83010931721345929D-5,.83010931721345929D-5,.45178996330148403D-4, &.83010931721345929D-5,.45178996330148403D-4,.81340012823988125D-5, &.50865948519564435D-4,.50865948519564435D-4,.13689259467368197D-4, &.50865948519564435D-4,.13689259467368197D-4,.14378089346564405D-4, &.70135517401683392D-5,.70135517401683392D-5,.53095487758849161D-1, &.70135517401683392D-5,.53095487758849161D-1,.93061914660707206D-5, &.28453632753053526D-4,.28453632753053526D-4,.75326086930357818D-5, &.28453632753053526D-4,.75326086930357818D-5,.10158361231153129D-3, &.11211958929549878D-4,.11211958929549878D-4,.18723169480270177D-4, &.11211958929549878D-4,.18723169480270177D-4,.71264821718346257D-5, &.43887530175994454D-3,.43887530175994454D-3,.87533429794514489D-5, &.43887530175994454D-3,.87533429794514489D-5,.35263202473010327D-4, &.77990108357659675D-5,.77990108357659675D-5,.69546582491103191D-4, &.77990108357659675D-5,.69546582491103191D-4,.12317547380782591D-4, &.16282601462695181D-4,.16282601462695181D-4,.70465172963973866D-5, &.16282601462695181D-4,.70465172963973866D-5,.14573954141527583D-2, &.97740091228072002D-5,.97740091228072002D-5,.24821237940107113D-4, &.97740091228072002D-5,.24821237940107113D-4,.73797894577225617D-5, &.14094394402508064D-3,.14094394402508064D-3,.10551733504650674D-4, &.14094394402508064D-3,.10551733504650674D-4,.20907781137971538D-4, &.72135079471586990D-5,.72135079471586990D-5,.25179830460726353D-3, &.72135079471586990D-5,.25179830460726353D-3,.84210452417342069D-5, &.41928481348814915D-4,.41928481348814915D-4,.80307262879160231D-5, &.41928481348814915D-4,.80307262879160231D-5,.55314317815738251D-4, &.13265892678035781D-4,.13265892678035781D-4,.14876755338549496D-4, &.13265892678035781D-4,.14876755338549496D-4,.70182484970768530D-5, &.87336489979141906D-2,.87336489979141906D-2,.92203263971508705D-5, &.87336489979141906D-2,.92203263971508705D-5,.29287542202657544D-4, &.75667952733664597D-5,.75667952733664597D-5,.95749706711456304D-4, &.75667952733664597D-5,.95749706711456304D-4,.11356051949487734D-4, &.18334673417044500D-4,.18334673417044500D-4,.71121142774758928D-5, &.18334673417044500D-4,.71121142774758928D-5,.50123471265910548D-3, &.88257367227356355D-5,.88257367227356355D-5,.34135224455316730D-4, &.88257367227356355D-5,.34135224455316730D-4,.77569665592606447D-5, &.73078759103808899D-4,.73078759103808899D-4,.12144907708551192D-4, &.73078759103808899D-4,.12144907708551192D-4,.16594424451616243D-4, &.70550275382997734D-5,.70550275382997734D-5,.11663938319774565D-2, &.70550275382997734D-5,.11663938319774565D-2,.98758113216642850D-5, &.24188045316340235D-4,.24188045316340235D-4,.73527229262935955D-5, &.24188045316340235D-4,.73527229262935955D-5,.15160236526448107D-3/ DATA ((COEF(I,J),J=1,6),I=605,624)/ &.10430831984330381D-4,.10430831984330381D-4,.21399249408204877D-4, &.10430831984330381D-4,.21399249408204877D-4,.72340177915616730D-5, &.22912277313547641D-3,.22912277313547641D-3,.83601869302891388D-5, &.22912277313547641D-3,.83601869302891388D-5,.43505328126211834D-4, &.80815753966592059D-5,.80815753966592059D-5,.53016674842238096D-4, &.80815753966592059D-5,.53016674842238096D-4,.13474180435492761D-4, &.14623255782452738D-4,.14623255782452738D-4,.70154299398635239D-5, &.14623255782452738D-4,.70154299398635239D-5,.17480968059050604D-1, &.94853980782479993D-5,.94853980782479993D-5,.26899769600761121D-4, &.94853980782479993D-5,.26899769600761121D-4,.74679144293547372D-5, &.11502124544576709D-3,.11502124544576709D-3,.10936225823121623D-4, &.11502124544576709D-3,.10936225823121623D-4,.19546132169864770D-4, &.71582281224251119D-5,.71582281224251119D-5,.34472497192009961D-3, &.71582281224251119D-5,.34472497192009961D-3,.86146380507171119D-5, &.37709163695009371D-4,.37709163695009371D-4,.78873129357672457D-5, &.37709163695009371D-4,.78873129357672457D-5,.63233711507255913D-4, &.12679208442049206D-4,.12679208442049206D-4,.15690959900448468D-4, &.12679208442049206D-4,.15690959900448468D-4,.70323639556081438D-5, &.24966702088980862D-2,.24966702088980862D-2,.95788845267725159D-5, &.24966702088980862D-2,.95788845267725159D-5,.26175302881899980D-4, &.74373684853268159D-5,.74373684853268159D-5,.12278856060901337D-3, &.74373684853268159D-5,.12278856060901337D-3,.10804309753069211D-4, &.19982182850850370D-4,.19982182850850370D-4,.71756236516959891D-5, &.19982182850850370D-4,.71756236516959891D-5,.30868674317730292D-3, &.85482260221103670D-5,.85482260221103670D-5,.39036723128523021D-4, &.85482260221103670D-5,.39036723128523021D-4,.79336264509158694D-5, &.60406618278828360D-4,.60406618278828360D-4,.12868650851841925D-4, &.60406618278828360D-4,.12868650851841925D-4,.15410216025711804D-4, &.70267132335712992D-5,.70267132335712992D-5,.34943747654845032D-2, &.70267132335712992D-5,.34943747654845032D-2,.90557319780320561D-5, &.31082017043330456D-4,.31082017043330456D-4,.76389494170612945D-5, &.31082017043330456D-4,.76389494170612945D-5,.85527200409329156D-4, &.11657429736260115D-4,.11657429736260115D-4,.17600044519983241D-4, &.11657429736260115D-4,.17600044519983241D-4,.70863492905706581D-5, &.67392192453364925D-3,.67392192453364925D-3,.89768743992750964D-5, &.67392192453364925D-3,.89768743992750964D-5,.32048310561331380D-4, &.76769605011983092D-5,.76769605011983092D-5,.81034949995067838D-4, &.76769605011983092D-5,.81034949995067838D-4,.11815031096909232D-4, &.17252595968224813D-4,.17252595968224813D-4,.70749381442147086D-5, &.17252595968224813D-4,.70749381442147086D-5,.79602335477132200D-3/ DATA ((COEF(I,J),J=1,6),I=625,644)/ &.10088332308522055D-4,.10088332308522055D-4,.23001286250929036D-4, &.10088332308522055D-4,.23001286250929036D-4,.82249764599011137D-5, &.48205883477429849D-4,.48205883477429849D-4,.73019578496566915D-5, &.48205883477429849D-4,.73019578496566915D-5,.19215239946299297D-3, &.82064153236885917D-5,.82064153236885917D-5,.48849902175383898D-4, &.82064153236885917D-5,.48849902175383898D-4,.72939272744513248D-5, &.18182190342841469D-3,.18182190342841469D-3,.10016139405254361D-4, &.18182190342841469D-3,.10016139405254361D-4,.23385590541189348D-4, &.73183856983155347D-5,.73183856983155347D-5,.16783919194119574D-3, &.73183856983155347D-5,.16783919194119574D-3,.10274960376545449D-4, &.22086627779660424D-4,.22086627779660424D-4,.72630155699443653D-5, &.22086627779660424D-4,.72630155699443653D-5,.20340159268741136D-3, &.82817803177761428D-5,.82817803177761428D-5,.45759771311496335D-4, &.82817803177761428D-5,.45759771311496335D-4,.81518329224876888D-5, &.50179535404707896D-4,.50179535404707896D-4,.13762509127210122D-4, &.50179535404707896D-4,.13762509127210122D-4,.14298159431690026D-4, &.70131344535112218D-5,.70131344535112218D-5,.98634587896770917D-1, &.70131344535112218D-5,.98634587896770917D-1,.93353582670693143D-5, &.28184397267047336D-4,.28184397267047336D-4,.75214881081098402D-5, &.28184397267047336D-4,.75214881081098402D-5,.10365028950798750D-3, &.11164867915526471D-4,.11164867915526471D-4,.18855979678537097D-4, &.11164867915526471D-4,.18855979678537097D-4,.71314934142952811D-5, &.42067100571978452D-3,.42067100571978452D-3,.87296686819650397D-5, &.42067100571978452D-3,.87296686819650397D-5,.35652713481509513D-4, &.78133343632378617D-5,.78133343632378617D-5,.68427960247206701D-4, &.78133343632378617D-5,.68427960247206701D-4,.12376282247981658D-4, &.16181090458597078D-4,.16181090458597078D-4,.70438934433909418D-5, &.16181090458597078D-4,.70438934433909418D-5,.15790507605649101D-2, &.97407135998734465D-5,.97407135998734465D-5,.25038585897253080D-4, &.97407135998734465D-5,.25038585897253080D-4,.73890648898516850D-5, &.13764401574836472D-3,.13764401574836472D-3,.10592817565160546D-4, &.13764401574836472D-3,.10592817565160546D-4,.20748329445045336D-4, &.72069049358893230D-5,.72069049358893230D-5,.26011726125703448D-3, &.72069049358893230D-5,.26011726125703448D-3,.84417307141811128D-5, &.41423099456057091D-4,.41423099456057091D-4,.80141211846759126D-5, &.41423099456057091D-4,.80141211846759126D-5,.56115163147971960D-4, &.13197926710413674D-4,.13197926710413674D-4,.14963168741271402D-4, &.13197926710413674D-4,.14963168741271402D-4,.70193972310481613D-5, &.72555322731565177D-2,.72555322731565177D-2,.91922414094686697D-5, &.72555322731565177D-2,.91922414094686697D-5,.29574559276065508D-4/ DATA ((COEF(I,J),J=1,6),I=645,664)/ &.75784682617187188D-5,.75784682617187188D-5,.93919162108882341D-4, &.75784682617187188D-5,.93919162108882341D-4,.11405041178614201D-4, &.18208397283132284D-4,.18208397283132284D-4,.71075459521606686D-5, &.18208397283132284D-4,.71075459521606686D-5,.52501700766228478D-3, &.88503315258487724D-5,.88503315258487724D-5,.33772233840768839D-4, &.88503315258487724D-5,.33772233840768839D-4,.77432573344545658D-5, &.74318367835533234D-4,.74318367835533234D-4,.12088525306870367D-4, &.74318367835533234D-4,.12088525306870367D-4,.16700857621420006D-4, &.70580777470986616D-5,.70580777470986616D-5,.10886140693684510D-2, &.70580777470986616D-5,.10886140693684510D-2,.99103947101303707D-5, &.23983066353754984D-4,.23983066353754984D-4,.73439520794113249D-5, &.23983066353754984D-4,.73439520794113249D-5,.15542976405522179D-3, &.10391301090882721D-4,.10391301090882721D-4,.21567574007585653D-4, &.10391301090882721D-4,.21567574007585653D-4,.72410894270395114D-5, &.22224824181431041D-3,.22224824181431041D-3,.83402949858280055D-5, &.22224824181431041D-3,.83402949858280055D-5,.44052077999326403D-4, &.80988738112618116D-5,.80988738112618116D-5,.52284073767360611D-4, &.80988738112618116D-5,.52284073767360611D-4,.13545105468989272D-4, &.14540624983851848D-4,.14540624983851848D-4,.70146994491895719D-5, &.14540624983851848D-4,.70146994491895719D-5,.23625790134455604D-1, &.94548204000864880D-5,.94548204000864880D-5,.27148766722798232D-4, &.94548204000864880D-5,.27148766722798232D-4,.74783609602150804D-5, &.11259866791773670D-3,.11259866791773670D-3,.10981062873346978D-4, &.11259866791773670D-3,.10981062873346978D-4,.19404523773056786D-4, &.71526563029992670D-5,.71526563029992670D-5,.35816114697944561D-3, &.71526563029992670D-5,.35816114697944561D-3,.86372049525145922D-5, &.37282764194806265D-4,.37282764194806265D-4,.78721989930082050D-5, &.37282764194806265D-4,.78721989930082050D-5,.64222234400661054D-4, &.12617363646236412D-4,.12617363646236412D-4,.15786719828827913D-4, &.12617363646236412D-4,.15786719828827913D-4,.70344583533336761D-5, &.22579862980373561D-2,.22579862980373561D-2,.96106406589243601D-5, &.22579862980373561D-2,.96106406589243601D-5,.25941074241710623D-4, &.74274488440372355D-5,.74274488440372355D-5,.12555700557001037D-3, &.74274488440372355D-5,.12555700557001037D-3,.10761186839147333D-4, &.20131382371663484D-4,.20131382371663484D-4,.71816500397613673D-5, &.20131382371663484D-4,.71816500397613673D-5,.29793179660272853D-3, &.85265123743031309D-5,.85265123743031309D-5,.39496039871398632D-4, &.85265123743031309D-5,.39496039871398632D-4,.79493918347442568D-5, &.59508033563101788D-4,.59508033563101788D-4,.12933128954801887D-4, &.59508033563101788D-4,.12933128954801887D-4,.15318760912728098D-4/ DATA ((COEF(I,J),J=1,6),I=665,684)/ &.70250401083509866D-5,.70250401083509866D-5,.39638924996044439D-2, &.70250401083509866D-5,.39638924996044439D-2,.90825197095203892D-5, &.30770522747231520D-4,.30770522747231520D-4,.76265683383493165D-5, &.30770522747231520D-4,.76265683383493165D-5,.87110527348713333D-4, &.11605940169065357D-4,.11605940169065357D-4,.17718725838517398D-4, &.11605940169065357D-4,.17718725838517398D-4,.70903701615630901D-5, &.63943641432683820D-3,.63943641432683820D-3,.89510829638859337D-5, &.63943641432683820D-3,.89510829638859337D-5,.32381411378787013D-4, &.76899230050432565D-5,.76899230050432565D-5,.79618298973567185D-4, &.76899230050432565D-5,.79618298973567185D-4,.11868630050263832D-4, &.17139570779396539D-4,.17139570779396539D-4,.70713508189351903D-5, &.17139570779396539D-4,.70713508189351903D-5,.84420964046508901D-3, &.98076217842823239D-5,.98076217842823239D-5,.24607072280116662D-4, &.98076217842823239D-5,.24607072280116662D-4,.73706410425901031D-5, &.14436616879689001D-3,.14436616879689001D-3,.10511044039624024D-4, &.14436616879689001D-3,.10511044039624024D-4,.21069392642262451D-4, &.72202274807020163D-5,.72202274807020163D-5,.24387578434883469D-3, &.72202274807020163D-5,.24387578434883469D-3,.84005605057323301D-5, &.42443804373619273D-4,.42443804373619273D-4,.80475028767157237D-5, &.42443804373619273D-4,.80475028767157237D-5,.54531300940871876D-4, &.13334583871302231D-4,.13334583871302231D-4,.14791307731010641D-4, &.13334583871302231D-4,.14791307731010641D-4,.70172044109755499D-5, &.10719461948765566D-1,.10719461948765566D-1,.93648019150224895D-5, &.10719461948765566D-1,.93648019150224895D-5,.27919378559385924D-4, &.75105037610891027D-5,.75105037610891027D-5,.10578227804152673D-3, &.75105037610891027D-5,.10578227804152673D-3,.11118237988784151D-4, &.18990491509672226D-4,.18990491509672226D-4,.71366161323380111D-5, &.18990491509672226D-4,.71366161323380111D-5,.40358254327943184D-3, &.87062196221959018D-5,.87062196221959018D-5,.36049253407009406D-4, &.87062196221959018D-5,.36049253407009406D-4,.78278139254281424D-5, &.67337106363353220D-4,.67337106363353220D-4,.12435623641506669D-4, &.67337106363353220D-4,.12435623641506669D-4,.16080764241600419D-4, &.70413757701322685D-5,.70413757701322685D-5,.17166490951628401D-2, &.70413757701322685D-5,.17166490951628401D-2,.97077320712827280D-5, &.25259178653803265D-4,.25259178653803265D-4,.73984679347512398D-5, &.25259178653803265D-4,.73984679347512398D-5,.13446062651788587D-3, &.10634300526976350D-4,.10634300526976350D-4,.20591000036522397D-4, &.10634300526976350D-4,.20591000036522397D-4,.72004180374146885D-5, &.26885955534806699D-3,.26885955534806699D-3,.84626185273662012D-5, &.26885955534806699D-3,.84626185273662012D-5,.40927405610447576D-4/ DATA ((COEF(I,J),J=1,6),I=685,704)/ &.79976863655926030D-5,.79976863655926030D-5,.56934381165342989D-4, &.79976863655926030D-5,.56934381165342989D-4,.13130676544442687D-4, &.15050561777770844D-4,.15050561777770844D-4,.70206506813289957D-5, &.15050561777770844D-4,.70206506813289957D-5,.61249354916809752D-2, &.91644206528268454D-5,.91644206528268454D-5,.29866262793537184D-4, &.91644206528268454D-5,.29866262793537184D-4,.75902812679486336D-5, &.92141976125271525D-4,.92141976125271525D-4,.11454518845034693D-4, &.92141976125271525D-4,.11454518845034693D-4,.18083689299517499D-4, &.71030876449894805D-5,.71030876449894805D-5,.55054212553741478D-3, &.71030876449894805D-5,.55054212553741478D-3,.88751614756059939D-5, &.33415497169647286D-4,.33415497169647286D-4,.77296992919324384D-5, &.33415497169647286D-4,.77296992919324384D-5,.75590922723302442D-4, &.12032713192709885D-4,.12032713192709885D-4,.16808568879072246D-4, &.12032713192709885D-4,.16808568879072246D-4,.70612349865270662D-5, &.10183831709339949D-2,.10183831709339949D-2,.10052063883366224D-4, &.10183831709339949D-2,.10052063883366224D-4,.23192071920993803D-4, &.73101104857711391D-5,.73101104857711391D-5,.17231272574459348D-3, &.73101104857711391D-5,.17231272574459348D-3,.10236917639439980D-4, &.22264482092676085D-4,.22264482092676085D-4,.72705628868469607D-5, &.22264482092676085D-4,.72705628868469607D-5,.19765552189361646D-3, &.82626575407165351D-5,.82626575407165351D-5,.46352511844574516D-4, &.82626575407165351D-5,.46352511844574516D-4,.81698449165062775D-5, &.49507655406011214D-4,.49507655406011214D-4,.13836554656717775D-4, &.49507655406011214D-4,.13836554656717775D-4,.14219105102117412D-4, &.70128215102423398D-5,.70128215102423398D-5,.41743498478793817D0, &.70128215102423398D-5,.41743498478793817D0,.92772989197892976D-5, &.28727174057397125D-4,.28727174057397125D-4,.75438662503433974D-5, &.28727174057397125D-4,.75438662503433974D-5,.99579607020824966D-4, &.11259516335534382D-4,.11259516335534382D-4,.18592033104304197D-4, &.11259516335534382D-4,.18592033104304197D-4,.14063575429701773D-4, &.45829477017769808D-3,.45829477017769808D-3,.87772444441947399D-5, &.45829477017769808D-3,.87772444441947399D-5,.34880556161053211D-4, &.77848423461079303D-5,.77848423461079303D-5,.70693922388411455D-4, &.77848423461079303D-5,.70693922388411455D-4,.12259411557498746D-4, &.16385315205394168D-4,.16385315205394168D-4,.70492474863571056D-5, &.16385315205394168D-4,.70492474863571056D-5,.13493094613431448D-2, &.92486780710376462D-5,.92486780710376462D-5,.29005112589185273D-4, &.92486780710376462D-5,.29005112589185273D-4,.75552615251588716D-5, &.97635766229904323D-4,.97635766229904323D-4,.11307545519343828D-4, &.97635766229904323D-4,.11307545519343828D-4,.18462543310874277D-4/ DATA ((COEF(I,J),J=1,6),I=705,724)/ &.71167928983991944D-5,.71167928983991944D-5,.47904006005742867D-3, &.71167928983991944D-5,.47904006005742867D-3,.88013750314859325D-5, &.34504615100383551D-4,.34504615100383551D-4,.77708279107186617D-5, &.34504615100383551D-4,.77708279107186617D-5,.71870970162673966D-4, &.12201867415273722D-4,.12201867415273722D-4,.16489249984395503D-4, &.12201867415273722D-4,.16489249984395503D-4,.70520841769397043D-5, &.12528455501298777D-2,.12528455501298777D-2,.98415547727577021D-5, &.12528455501298777D-2,.98415547727577021D-5,.24396027918411155D-4, &.73616190579116149D-5,.73616190579116149D-5,.14791680756823022D-3, &.73616190579116149D-5,.14791680756823022D-3,.10470744926788519D-4, &.21233202323452044D-4,.21233202323452044D-4,.72270639537941898D-5, &.21233202323452044D-4,.72270639537941898D-5,.23632490318042525D-3, &.83802749228177005D-5,.83802749228177005D-5,.42969329942629927D-4, &.83802749228177005D-5,.42969329942629927D-4,.80644521684489196D-5, &.53765588308960295D-4,.53765588308960295D-4,.13404009874506138D-4, &.53765588308960295D-4,.13404009874506138D-4,.14706812328188258D-4, &.70162649105540367D-5,.70162649105540367D-5,.13478768207244246D-1, &.70162649105540367D-5,.13478768207244246D-1,.95162661094599139D-5, &.26654577119062012D-4,.26654577119062012D-4,.74576006200979565D-5, &.26654577119062012D-4,.74576006200979565D-5,.11752463037037148D-3, &.10891824505456188D-4,.10891824505456188D-4,.19689590333047166D-4, &.10891824505456188D-4,.19689590333047166D-4,.71639130531249136D-5, &.33203596934200402D-3,.33203596934200402D-3,.85922871127669669D-5, &.33203596934200402D-3,.85922871127669669D-5,.38143491148269066D-4, &.79025881036570546D-5,.79025881036570546D-5,.62268757844612755D-4, &.79025881036570546D-5,.62268757844612755D-4,.12741699052098109D-4, &.15596299928649304D-4,.15596299928649304D-4,.70303750478911679D-5, &.15596299928649304D-4,.70303750478911679D-5,.27754390877647839D-2, &.10161915328548843D-4,.10161915328548843D-4,.26413113002045082D-4, &.10161915328548843D-4,.26413113002045082D-4,.74474188596923435D-5, &.12011245561186619D-3,.12011245561186619D-3,.10847854073726147D-4, &.12011245561186619D-3,.10847854073726147D-4,.19834929815447473D-4, &.71697114440888450D-5,.71697114440888450D-5,.32003973873601668D-3, &.71697114440888450D-5,.32003973873601668D-3,.85701503562527331D-5, &.38585943320579917D-4,.38585943320579917D-4,.79180255770751269D-5, &.38585943320579917D-4,.79180255770751269D-5,.61326631179228026D-4, &.12804843726408045D-4,.12804843726408045D-4,.15502723882241799D-4, &.12804843726408045D-4,.15502723882241799D-4,.70284915113695867D-5, &.31038237254268014D-2,.31038237254268014D-2,.90291966870403351D-5, &.31038237254268014D-2,.90291966870403351D-5,.31398734883402714D-4/ DATA ((COEF(I,J),J=1,6),I=725,736)/ &.76514745710831673D-5,.76514745710831673D-5,.83987893343168382D-4, &.76514745710831673D-5,.83987893343168382D-4,.11709437272952442D-4, &.17482810818312803D-4,.17482810818312803D-4,.70824371537868323D-5, &.17482810818312803D-4,.70824371537868323D-5,.71128686246649871D-3, &.90029115721813132D-5,.90029115721813132D-5,.31720792790347823D-4, &.90029115721813132D-5,.31720792790347823D-4,.76641446464951208D-5, &.82490990576319028D-4,.82490990576319028D-4,.11761968965338093D-4, &.82490990576319028D-4,.11761968965338093D-4,.17367002047293309D-4, &.70786335148686870D-5,.70786335148686870D-5,.75186114338431098D-3, &.70786335148686870D-5,.75186114338431098D-3,.10199237216518730D-4, &.22627713278614384D-4,.22627713278614384D-4,.72860182531795594D-5, &.22627713278614384D-4,.72860182531795594D-5,.18687871989425724D-3, &.82437233965993013D-5,.82437233965993013D-5,.47575220605350039D-4, &.82437233965993013D-5,.47575220605350039D-4,.13911406857838684D-4, &.10124948250071571D-4,.10124948250071571D-4,.22813183119737957D-4, &.10124948250071571D-4,.22813183119737957D-4,.72782302873419314D-5, &.17697022092294363D-3,.17697022092294363D-3,.70126128917403761D-5, &.17697022092294363D-3,.70126128917403761D-5,.22444828609212017D-4, &.14140914363232827D-4,.14140914363232827D-4,.81880385993792202D-5, &.14140914363232827D-4,.81880385993792202D-5,.46957547824856738D-4, &.13987076721131579D-4,.13987076721131579D-4,.70125085855924287D-5, &.13987076721131579D-4,.70125085855924287D-5,.95474273078995595D-5, &-.47318385459766968D0,.49616903185419430D0,-.25938664641601587D0, &.49616903185419430D0,-.25938664641601587D0,.71215820993688990D-5/ C ---- INPUT VALIDATION ZEROVALUE=1.D-14 IF(H0.LT.ZEROVALUE) THEN WRITE(*,*) "INITIAL STEP (H0) IS TOO SMALL" RETURN ENDIF IF (TEND.LE.T) THEN WRITE(*,*) "END-TIME (TEND) IS LESS THAN INITIAL TIME (T0)" RETURN END IF IF (IWORK(4).EQ.0) THEN IF (ATOL(1).LE.ZEROVALUE.OR.RTOL(1).LE.10.D0*ZEROVALUE) THEN WRITE(6,*) 'TOLERANCES ARE TOO SMALL', ATOL(1)," ", RTOL(1) RETURN END IF ELSE DO I=1,N IF (ATOL(I).LE.ZEROVALUE.OR.RTOL(I).LE.10.D0*ZEROVALUE) THEN WRITE(6,*) 'TOLERANCES ARE TOO SMALL' RETURN END IF END DO END IF IF ( IWORK(1).EQ.0 ) THEN WRITE(*,*) 'IWORK(1)=0 - Compute the spectral radius internally' ELSE IF (IWORK(1).EQ.1) THEN WRITE(*,*) 'IWORK(1)=1 - RHO is used to obtain spectral radius' ELSE WRITE(*,*) 'INCORRECT IWORK(1) PARAMETER - IT CAN BE 0 OR 1' END IF IF ( IWORK(2).EQ.0 ) THEN WRITE(*,*) 'IWORK(2)=0 - The Jacobian is not constant' ELSE IF (IWORK(2).EQ.1) THEN WRITE(*,*) 'IWORK(2)=1 - The Jacobian is constant' ELSE WRITE(*,*) 'INCORRECT IWORK(2) PARAMETER - IT CAN BE 0 OR 1' END IF IF ( IWORK(3).EQ.0 ) THEN WRITE(*,*) 'IWORK(3)=0 - DUMKA3 compute until t=tend' ELSE IF (IWORK(3).EQ.1) THEN WRITE(*,*) 'IWORK(3)=1 - DUMKA3 returns after each step' ELSE WRITE(*,*) 'INCORRECT IWORK(3) PARAMETER - IT CAN BE 0 OR 1' END IF IF ( IWORK(4).EQ.0 ) THEN WRITE(*,*) 'IWORK(4)=0 - Atol and rtol are scalar' ELSE IF (IWORK(4).EQ.1) THEN WRITE(*,*) 'IWORK(4)=1 - Atol and rtol are arrays of length n' ELSE WRITE(*,*) 'INCORRECT IWORK(4) PARAMETER - IT CAN BE 0 OR 1' END IF C ---- INITIALIZATION ------------------- IWORK(5)=0 IWORK(6)=0 IWORK(7)=0 IWORK(8)=0 IWORK(9)=0 IWORK(10)=0 IWORK(11)=0 IWORK(12)=0 I_ERR=0 ERR_N = 0.D0 H_N = 0.D0 H_NEW=H0 CALL FUN(N,T,Y,Z) IWORK(5)=IWORK(5)+1 IF ( IWORK(1).EQ.0 ) THEN CALL FUN(N,T,Z,OLDEIGENVECTOR) IWORK(9)=IWORK(9)+1 ENDIF 200 CONTINUE C ---- CALCULATE 2/MAX_EIGENVALUE (WHICH IS MAXIMUM STEP C ---- SIZE IN EULER METHOD) C ---- VIA INTERNAL(EIGENVALUE) OR EXTERNAL(RHO) PROCEDURE IWORK(6)=IWORK(6)+1 IF ( IWORK(2).EQ.0.OR.IWORK(6).EQ.1) THEN IF ( IWORK(1).EQ.0 ) THEN IF ((IWORK(6).EQ.1).OR. & (MOD(IWORK(6),20).EQ.0).OR.(I_ERR.EQ.1)) THEN EV = DABS( EIGENVALUE(N, Y, Z, OLDEIGENVECTOR, Z1, & EIGENVALUEOLD, T, FUN, ZEROVALUE, IWORK(9)) ) IF ( EV.LT.ZEROVALUE ) THEN EV = ZEROVALUE END IF COU = 2.D0/ (EV*1.2D0) END IF ELSE EV = DABS ( RHO(N,T,Y) ) IF (EV.EQ.0.D0) THEN EV = ZEROVALUE WRITE(*,*) "RHO RETURNS 0: ", EV," WILL BE USED" END IF COU = 2.D0/EV ENDIF ENDIF IF ( IWORK(11).LT.EV ) IWORK(11)=EV IF ( IWORK(12).GT.EV.OR.IWORK(6).EQ.1 ) IWORK(12)=EV I_ERR=0 H=DMIN1(H_NEW,(TEND-T)) INDEX=1 5 IF((2.D0*H/COU).GT.STAB_REG(INDEX)) THEN IF(INDEX.EQ.N_P) GOTO 6 INDEX=INDEX+1 GOTO 5 END IF 6 CONTINUE IF(INDEX.GT.1) THEN IF(((STAB_REG(INDEX-1)/(N_DEG(INDEX-1)*1.D0))*N_DEG(INDEX)*1.D0) & .GT.(2.d0*H/COU))THEN INDEX=INDEX-1 END IF END IF N_POL_DEGREE=N_DEG(INDEX) IF(IWORK(10).LT.N_POL_DEGREE) IWORK(10)=N_POL_DEGREE H=DMIN1((STAB_REG(INDEX)*COU*0.5d0),H) C WRITE(*,300) N_POL_DEGREE,H/(N_POL_DEGREE*COU),T 300 FORMAT (1X,'SPEED UP=',I3,' TIMES',' H/COU=',F14.7,' TIME',F14.7) N_POL=0 DO K=INDEX_FIRST(INDEX),INDEX_LAST(INDEX) IF(K.EQ.INDEX_FIRST(INDEX)) THEN DO I=1,N Y_EN(I)=Y(I) END DO DO I=1,N Z_EN(I)=Z(I) END DO T_EN=T END IF N_POL=N_POL+3 A_21=H*COEF(K,1) C_2=A_21 A_31=H*COEF(K,2) A_32=H*COEF(K,3) C_3=A_31+A_32 A_41=H*COEF(K,4) A_42=H*COEF(K,5) A_43=H*COEF(K,6) C_4=A_41+A_42+A_43 DO I=1,N Y(I)=Y(I)+A_21*Z(I) END DO T2=T+C_2 CALL FUN(N,T2,Y,Z1) IWORK(5)=IWORK(5)+1 IF(N_POL.EQ.N_DEG(INDEX)) THEN R=H*(COEF(K,2)-COEF(K,1)) DO I=1,N Y(I)=Y(I)+R*Z(I)+A_32*Z1(I) END DO ELSE DO I=1,N Y(I)=Y(I)+A_32*Z1(I) END DO END IF T3=T+C_3 IF(N_POL.EQ.N_DEG(INDEX)) THEN TMP_1=C_3/2.D0 TMP_2=(C_4-C_2)/2.D0 DO I=1,N Z1(I)=TMP_1*Z1(I)-TMP_2*Z(I) END DO END IF CALL FUN(N,T3,Y,Z) IWORK(5)=IWORK(5)+1 IF(N_POL.EQ.N_DEG(INDEX)) THEN DO I=1,N Z1(I)=TMP_2*Z(I)+Z1(I) END DO END IF DO I=1,N Y(I)=Y(I)+A_43*Z(I) END DO T4=T+C_4 T=T4 CALL FUN(N,T,Y,Z) IWORK(5)=IWORK(5)+1 END DO IF(IWORK(4).EQ.0) THEN ATOL_SCALAR = ATOL(1) RTOL_SCALAR = RTOL(1) DO I=1,N Z1(I)=((-TMP_1)*Z(I)+Z1(I))/ & (RTOL_SCALAR*DMAX1(DABS(Y(I)),DABS(Y_EN(I)))+ATOL_SCALAR) END DO ELSE DO I=1,N Z1(I)=((-TMP_1)*Z(I)+Z1(I))/ & (RTOL(I)*DMAX1(DABS(Y(I)),DABS(Y_EN(I)))+ATOL(I)) END DO END IF C IF(N_POL_DEGREE.EQ.3) THEN C ERROR_CORRECTOR = -3.4785D0/ 8.3338D0 C DO I=1,N C Z1(I)=ERROR_CORRECTOR*Z1(I) C END DO C END IF H_NEW=H CALL ST(N,T,Y,Z1,H_NEW,I_ERR, & N_POL_DEGREE,ERR_N,H_N,IWORK(6)) IF(I_ERR.EQ.1) THEN DO I=1,N Y(I)=Y_EN(I) END DO DO I=1,N Z(I)=Z_EN(I) END DO T=T_EN IWORK(8)=IWORK(8)+1 GOTO 200 ELSE IWORK(7)=IWORK(7)+1 IF(T.GE.(TEND).OR.IWORK(3).EQ.1) THEN H0=H_NEW RETURN ELSE GOTO 200 END IF END IF RETURN END c c c SUBROUTINE ST(N,T,Y,Z1,H_NEW,I_ERR, & N_POL_DEGREE,ERR_N,H_N,ISTEP_ID) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DOUBLE PRECISION Y(N),Z1(N) EPS1=0.D0 DO I=1,N EPS1=EPS1+Z1(I)*Z1(I) END DO EPS=DSQRT( EPS1/(1.D0*N) ) FRACMIN=0.1D0 IF(EPS.EQ.0.D0) EPS=1.D-14 FRAC=(1.D0/EPS)**(1.D0/3.D0) IF(EPS.LE.1.D0) THEN IF (ERR_N.GT.0.D0.AND.H_N.GT.0.D0) THEN FRAC2=ERR_N**(1.D0/3.D0)*FRAC**2.d0*(H_NEW/H_N) FRAC=DMIN1(FRAC,FRAC2) END IF I_ERR=0 FRACMAX=2.D0 FRAC=DMIN1(FRACMAX,DMAX1(FRACMIN,0.8D0*FRAC)) H_OLD = H_NEW H_NEW=FRAC*H_NEW H_N=H_OLD ERR_N=EPS C WRITE(*,102) EPS, T C WRITE(*,101) N_POL_DEGREE, H_NEW, H_OLD ELSE I_ERR=1 FRACMAX=1.D0 FRAC=0.8D0*DMIN1(FRACMAX,DMAX1(FRACMIN,0.8D0*FRAC)) IF(ISTEP_ID.EQ.1) THEN H_NEW=FRACMIN*H_NEW ELSE H_NEW=FRAC*H_NEW ENDIF WRITE(*,100) EPS, T WRITE(*,101) N_POL_DEGREE, H_NEW, H_OLD END IF 100 FORMAT('STEP REJECTED, BECAUSE ERROR=',1E18.10, & ' TIME=',1E18.10) 102 FORMAT('STEP ACCEPTED, BECAUSE ERROR=',1E18.10, & ' TIME=',1E18.10) 101 FORMAT(' POLYNOMIAL DEGREE =', I3, & ' NEW STEP =',1E18.10,' OLD STEP',1E18.10) RETURN END DOUBLE PRECISION FUNCTION EIGENVALUE(N,Y,FY,V, & TMP,EIGENVALUEOLD,T,FUN,ZEROVALUE,NUM_FUN_EVAL) IMPLICIT DOUBLE PRECISION(A-H,O-Z) EXTERNAL FUN DOUBLE PRECISION Y(N),FY(N),V(N),TMP(N) WRITE(*,*) "COMPUTE EIGENVALUE ..." MAXIT = 30 TOL = 0.01 IDX = 0 RADIUS = 1.D-10 ALMBDA = EIGENVALUEOLD DNRMV = 0.D0 DNRMV = DNORM(N,V) DO WHILE ( IDX.EQ.0.OR.(DABS(ALMBDA - EIGENVALUEOLD) & .GT.DABS(ALMBDA)*TOL.AND.IDX.LE.MAXIT ) ) IDX = IDX+1 C ALL VECTORS SHOULD BE ||Y-V|| < RADIUS, BECAUSE PROBLEM CAN BE NON-LINEAR IF ( DNRMV.LT.ZEROVALUE ) THEN DO I=1,N V(I) = RAND(1) END DO DNRMV = DNORM(N,V) ENDIF DNORMINV= RADIUS/DNRMV DO I=1,N V(I) = Y(I) + DNORMINV * V(I) END DO CALL FUN(N,T,V,TMP) NUM_FUN_EVAL=NUM_FUN_EVAL+1 C APPROXIMATION OF EIGENVECTOR DO I=1,N V(I) = TMP(I) - FY(I) END DO C // APPROXIMATION OF EIGENVALUE DNRMV = DNORM(N,V) EIGENVALUEOLD = ALMBDA ALMBDA = DNRMV/RADIUS WRITE(*,*) "ITERATION = ", IDX, " EIGENVALUE = ", ALMBDA END DO EIGENVALUEOLD = ALMBDA EIGENVALUE = ALMBDA RETURN END DOUBLE PRECISION FUNCTION DNORM ( N, Y ) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DOUBLE PRECISION Y(N) DNORM = 0.D0 DO I=1,N DNORM = DNORM + Y(I)*Y(I) END DO DNORM = DSQRT(DNORM) RETURN END