PROGRAM main c*********************************************************************72 c cc MAIN is the main program for TOMS757_PRB. c c Discussion: c c TOMS757_PRB tests the TOMS757 library. C C This program tests the 37 functions in the file MISCFUN. C It is a fairly simple code with each function being tested C at 20 different arguments. The code compares the value C from the function with a pre-computed value, and produces C the absolute and relative errors. C C The function codes are contained in he file MISCFUN.F. These C MUST be changed to double precision form, appended to this C file, and the complete program compiled. C C The program is written in double-precision Fortran 77. For C single-precision tests use the file MSCTSTSP.F C C The variable IOUT is used to direct the program output C and is set to the standard value 6. For other output streams C change the data statement. C C Author: Allan MacLeod, C Dept. of Mathematics and Statistics, C University of Paisley, C High St., C Paisley, C SCOTLAND C PA1 2BE C C (e-mail: macl_ms0@paisley.ac.uk ) C C INTEGER I,IOUT DOUBLE PRECISION ABSERR,COMP,DEN,NUM,PT,RELERR,RES DOUBLE PRECISION AB0DAT(3,20),AB1DAT(3,20), & AB2DAT(3,20),AIIDAT(3,20),AGIDAT(3,20), & AHIDAT(3,20),ATNDAT(3,20),BIIDAT(3,20),CLNDAT(3,20), & DB1DAT(3,20),DB2DAT(3,20),DB3DAT(3,20),DB4DAT(3,20), & EX3DAT(3,20),GSTDAT(3,20),I0IDAT(3,20),IL0DAT(3,20), & IL1DAT(3,20),J0IDAT(3,20),K0IDAT(3,20),LOBDAT(3,20), & STRDAT(3,20),SH0DAT(3,20),SH1DAT(3,20),SL0DAT(3,20), & SL1DAT(3,20),SY1DAT(3,20),SY2DAT(3,20),TR2DAT(3,20), & TR3DAT(3,20),TR4DAT(3,20),TR5DAT(3,20),TR6DAT(3,20), & TR7DAT(3,20),TR8DAT(3,20),TR9DAT(3,20),Y0IDAT(3,20) DOUBLE PRECISION ABRAM0,ABRAM1,ABRAM2,AIRINT,AIRYGI,AIRYHI, & ATNINT,BIRINT,CLAUSN,DEBYE1,DEBYE2,DEBYE3,DEBYE4, & EXP3,GOODST,I0INT,I0ML0,I1ML1,J0INT,K0INT,LOBACH, & STROM,STRVH0,STRVH1,STRVL0,STRVL1,SYNCH1,SYNCH2, & TRAN02,TRAN03,TRAN04,TRAN05,TRAN06,TRAN07,TRAN08, & TRAN09,Y0INT CHARACTER*6 FNAME DATA IOUT/6/ DATA AB0DAT/1.0D0, 512.0D0, 0.87377726306985360531D0, & 1.0D0, 128.0D0, 0.84721859650456925922D0, & 1.0D0, 32.0D0, 0.77288934483988301615D0, & 1.0D0, 8.0D0, 0.59684345853450151603D0, & 1.0D0, 2.0D0, 0.29871735283675888392D0, & 1.0D0, 1.0D0, 0.15004596450516388138D0, & 5.0D0, 4.0D0, 0.11114662419157955096D0, & 3.0D0, 2.0D0, 0.83909567153151897766D-1, & 15.0D0, 8.0D0, 0.56552321717943417515D-1, & 2.0D0, 1.0D0, 0.49876496603033790206D-1, & 17.0D0, 8.0D0, 0.44100889219762791328D-1, & 3.0D0, 1.0D0, 0.19738535180254062496D-1, & 4.0D0, 1.0D0, 0.86193088287161479900D-2, & 5.0D0, 1.0D0, 0.40224788162540127227D-2, & 6.0D0, 1.0D0, 0.19718658458164884826D-2, & 7.0D0, 1.0D0, 0.10045868340133538505D-2, & 10.0D0, 1.0D0, 0.15726917263304498649D-3, & 15.0D0, 1.0D0, 0.10352666912350263437D-4, & 20.0D0, 1.0D0, 0.91229759190956745069D-6, & 40.0D0, 1.0D0, 0.25628287737952698742D-9/ DATA AB1DAT/1.0D0, 512.0D0, 0.49828219848799921792D0, & 1.0D0, 128.0D0, 0.49324391773047288556D0, & 1.0D0, 32.0D0, 0.47431612784691234649D0, & 1.0D0, 8.0D0, 0.41095983258760410149D0, & 1.0D0, 2.0D0, 0.25317617388227035867D0, & 1.0D0, 1.0D0, 0.14656338138597777543D0, & 5.0D0, 4.0D0, 0.11421547056018366587D0, & 3.0D0, 2.0D0, 0.90026307383483764795D-1, & 15.0D0, 8.0D0, 0.64088214170742303375D-1, & 2.0D0, 1.0D0, 0.57446614314166191085D-1, & 17.0D0, 8.0D0, 0.51581624564800730959D-1, & 3.0D0, 1.0D0, 0.25263719555776416016D-1, & 4.0D0, 1.0D0, 0.11930803330196594536D-1, & 5.0D0, 1.0D0, 0.59270542280915272465D-2, & 6.0D0, 1.0D0, 0.30609215358017829567D-2, & 7.0D0, 1.0D0, 0.16307382136979552833D-2, & 10.0D0, 1.0D0, 0.28371851916959455295D-3, & 15.0D0, 1.0D0, 0.21122150121323238154D-4, & 20.0D0, 1.0D0, 0.20344578892601627337D-5, & 40.0D0, 1.0D0, 0.71116517236209642290D-9/ DATA AB2DAT/1.0D0, 512.0D0, 0.44213858162107913430D0, & 1.0D0, 128.0D0, 0.43923379545684026308D0, & 1.0D0, 32.0D0, 0.42789857297092602234D0, & 1.0D0, 8.0D0, 0.38652825661854504406D0, & 1.0D0, 2.0D0, 0.26538204413231368110D0, & 1.0D0, 1.0D0, 0.16848734838334595000D0, & 5.0D0, 4.0D0, 0.13609200032513227112D0, & 3.0D0, 2.0D0, 0.11070330027727917352D0, & 15.0D0, 8.0D0, 0.82126019995530382267D-1, & 2.0D0, 1.0D0, 0.74538781999594581763D-1, & 17.0D0, 8.0D0, 0.67732034377612811390D-1, & 3.0D0, 1.0D0, 0.35641808698811851022D-1, & 4.0D0, 1.0D0, 0.17956589956618269083D-1, & 5.0D0, 1.0D0, 0.94058737143575370625D-2, & 6.0D0, 1.0D0, 0.50809356204299213556D-2, & 7.0D0, 1.0D0, 0.28149565414209719359D-2, & 10.0D0, 1.0D0, 0.53808696422559303431D-3, & 15.0D0, 1.0D0, 0.44821756380146327259D-4, & 20.0D0, 1.0D0, 0.46890678427324100410D-5, & 40.0D0, 1.0D0, 0.20161544850996420504D-8/ DATA AIIDAT/-12.0D0, 1.0D0, -0.75228838916610124300D0, & -11.0D0, 1.0D0, -0.57348350185854889466D0, & -10.0D0, 1.0D0, -0.76569840313421291743D0, & -19.0D0, 2.0D0, -0.65181015505382467421D0, & -9.0D0, 1.0D0, -0.55881974894471876922D0, & -13.0D0, 2.0D0, -0.56902352870716815309D0, & -4.0D0, 1.0D0, -0.47800749642926168100D0, & -1.0D0, 1.0D0, -0.46567398346706861416D0, & -1.0D0, 4.0D0, -0.96783140945618013679D-1, & -1.0D0, 1024.0D0, -0.34683049857035607494D-3, & 1.0D0, 1024.0D0, 0.34658366917927930790D-3, & 1.0D0, 128.0D0, 0.27657581846051227124D-2, & 1.0D0, 2.0D0, 0.14595330491185717833D0, & 1.0D0, 1.0D0, 0.23631734191710977960D0, & 4.0D0, 1.0D0, 0.33289264538612212697D0, & 9.0D0, 2.0D0, 0.33318759129779422976D0, & 6.0D0, 1.0D0, 0.33332945170523851439D0, & 8.0D0, 1.0D0, 0.33333331724248357420D0, & 10.0D0, 1.0D0, 0.33333333329916901594D0, & 12.0D0, 1.0D0, 0.33333333333329380187D0/ DATA AGIDAT/-1.0D0, 512.0D0, 0.20468308070040542435D0, & -1.0D0, 8.0D0, 0.18374662832557904078D0, & -1.0D0, 1.0D0, -0.11667221729601528265D0, & -4.0D0, 1.0D0, 0.31466934902729557596D0, & -8.0D0, 1.0D0, -0.37089040722426257729D0, & -33.0D0, 4.0D0, -0.25293059772424019694D0, & -9.0D0, 1.0D0, 0.28967410658692701936D0, & -10.0D0, 1.0D0, -0.34644836492634090590D0, & -11.0D0, 1.0D0, 0.28076035913873049496D0, & -13.0D0, 1.0D0, 0.21814994508094865815D0, & 1.0D0, 512.0D0, 0.20526679000810503329D0, & 1.0D0, 8.0D0, 0.22123695363784773258D0, & 1.0D0, 1.0D0, 0.23521843981043793760D0, & 4.0D0, 1.0D0, 0.82834303363768729338D-1, & 7.0D0, 1.0D0, 0.45757385490989281893D-1, & 29.0D0, 4.0D0, 0.44150012014605159922D-1, & 8.0D0, 1.0D0, 0.39951133719508907541D-1, & 9.0D0, 1.0D0, 0.35467706833949671483D-1, & 10.0D0, 1.0D0, 0.31896005100679587981D-1, & 12.0D0, 1.0D0, 0.26556892713512410405D-1/ DATA AHIDAT/-1.0D0, 512.0D0, 0.40936798278458884024D0, & -1.0D0, 8.0D0, 0.37495291608048868619D0, & -1.0D0, 1.0D0, 0.22066960679295989454D0, & -4.0D0, 1.0D0, 0.77565356679703713590D-1, & -8.0D0, 1.0D0, 0.39638826473124717315D-1, & -33.0D0, 4.0D0, 0.38450072575004151871D-1, & -9.0D0, 1.0D0, 0.35273216868317898556D-1, & -10.0D0, 1.0D0, 0.31768535282502272742D-1, & -11.0D0, 1.0D0, 0.28894408288051391369D-1, & -13.0D0, 1.0D0, 0.24463284011678541180D-1, & 1.0D0, 512.0D0, 0.41053540139998941517D0, & 1.0D0, 8.0D0, 0.44993502381204990817D0, & 1.0D0, 1.0D0, 0.97220515514243332184D0, & 4.0D0, 1.0D0, 0.83764237105104371193D2, & 7.0D0, 1.0D0, 0.80327744952044756016D5, & 29.0D0, 4.0D0, 0.15514138847749108298D6, & 8.0D0, 1.0D0, 0.11995859641733262114D7, & 9.0D0, 1.0D0, 0.21472868855967642259D8, & 10.0D0, 1.0D0, 0.45564115351632913590D9, & 12.0D0, 1.0D0, 0.32980722582904761929D12/ DATA ATNDAT/1.0D0, 512.0D0, 0.19531241721588483191D-2, & -1.0D0, 256.0D0, -0.39062433772980711281D-2, & 1.0D0, 128.0D0, 0.78124470192576499535D-2, & 1.0D0, 64.0D0, 0.15624576181996527280D-1, & -1.0D0, 32.0D0, -0.31246610349485401551D-1, & 1.0D0, 16.0D0, 0.62472911335014397321D-1, & 1.0D0, 8.0D0, 0.12478419717389654039D0, & -1.0D0, 4.0D0, -0.24830175098230686908D0, & 1.0D0, 2.0D0, 0.48722235829452235711D0, & 1.0D0, 1.0D0, 0.91596559417721901505D0, & 3.0D0, 2.0D0, 0.12749694484943800618D1, & -2.0D0, 1.0D0, -0.15760154034463234224D1, & 4.0D0, 1.0D0, 0.24258878412859089996D1, & 8.0D0, 1.0D0, 0.33911633326292997361D1, & 16.0D0, 1.0D0, 0.44176450919422186583D1, & -20.0D0, 1.0D0, -0.47556713749547247774D1, & 25.0D0, 1.0D0, 0.50961912150934111303D1, & 30.0D0, 1.0D0, 0.53759175735714876256D1, & -50.0D0, 1.0D0, -0.61649904785027487422D1, & 100.0D0, 1.0D0, 0.72437843013083534973D1/ DATA BIIDAT/-12.0D0, 1.0D0, 0.17660819031554631869D-1, & -10.0D0, 1.0D0, -0.15040424806140020451D-1, & -8.0D0, 1.0D0, 0.14756446293227661920D-1, & -15.0D0, 2.0D0, -0.11847304264848446271D0, & -7.0D0, 1.0D0, -0.64916741266165856037D-1, & -13.0D0, 2.0D0, 0.97260832464381044540D-1, & -4.0D0, 1.0D0, 0.50760058495287539119D-1, & -1.0D0, 1.0D0, -0.37300500963429492179D0, & -1.0D0, 4.0D0, -0.13962988442666578531D0, & -1.0D0, 512.0D0, -0.12001735266723296160D-2, & 1.0D0, 512.0D0, 0.12018836117890354598D-2, & 1.0D0, 2.0D0, 0.36533846550952011043D0, & 1.0D0, 1.0D0, 0.87276911673800812196D0, & 4.0D0, 1.0D0, 0.48219475263803429675D2, & 8.0D0, 1.0D0, 0.44006525804904178439D6, & 17.0D0, 2.0D0, 0.17608153976228301458D7, & 9.0D0, 1.0D0, 0.73779211705220007228D7, & 10.0D0, 1.0D0, 0.14780980310740671617D9, & 12.0D0, 1.0D0, 0.97037614223613433849D11, & 14.0D0, 1.0D0, 0.11632737638809878460D15/ DATA CLNDAT/1.0D0, 512.0D0, 0.14137352886760576684D-1, & 1.0D0, 32.0D0, 0.13955467081981281934D0, & -1.0D0, 8.0D0, -0.38495732156574238507D0, & 1.0D0, 2.0D0, 0.84831187770367927099D0, & 1.0D0, 1.0D0, 0.10139591323607685043D1, & -3.0D0, 2.0D0, -0.93921859275409211003D0, & 2.0D0, 1.0D0, 0.72714605086327924743D0, & 5.0D0, 2.0D0, 0.43359820323553277936D0, & -3.0D0, 1.0D0, -0.98026209391301421161D-1, & 4.0D0, 1.0D0, -0.56814394442986978080D0, & 17.0D0, 4.0D0, -0.70969701784448921625D0, & -5.0D0, 1.0D0, 0.99282013254695671871D0, & 11.0D0, 2.0D0, -0.98127747477447367875D0, & 6.0D0, 1.0D0, -0.64078266570172320959D0, & 8.0D0, 1.0D0, 0.86027963733231192456D0, & -10.0D0, 1.0D0, 0.39071647608680211043D0, & 15.0D0, 1.0D0, 0.47574793926539191502D0, & 20.0D0, 1.0D0, 0.10105014481412878253D1, & -30.0D0, 1.0D0, 0.96332089044363075154D0, & 50.0D0, 1.0D0, -0.61782699481929311757D0/ DATA DB1DAT/1.0D0, 512.0D0, 0.99951182471380889183D0, & 1.0D0, 32.0D0, 0.99221462647120597836D0, & 1.0D0, 8.0D0, 0.96918395997895308324D0, & 1.0D0, 2.0D0, 0.88192715679060552968D0, & 1.0D0, 1.0D0, 0.77750463411224827642D0, & 3.0D0, 2.0D0, 0.68614531078940204342D0, & 2.0D0, 1.0D0, 0.60694728460981007205D0, & 5.0D0, 2.0D0, 0.53878956907785587703D0, & 3.0D0, 1.0D0, 0.48043521957304283829D0, & 4.0D0, 1.0D0, 0.38814802129793784501D0, & 17.0D0, 4.0D0, 0.36930802829242526815D0, & 5.0D0, 1.0D0, 0.32087619770014612104D0, & 11.0D0, 2.0D0, 0.29423996623154246701D0, & 6.0D0, 1.0D0, 0.27126046678502189985D0, & 8.0D0, 1.0D0, 0.20523930310221503723D0, & 10.0D0, 1.0D0, 0.16444346567994602563D0, & 15.0D0, 1.0D0, 0.10966194482735821276D0, & 20.0D0, 1.0D0, 0.82246701178200016086D-1, & 30.0D0, 1.0D0, 0.54831135561510852445D-1, & 50.0D0, 1.0D0, 0.32898681336964528729D-1/ DATA DB2DAT/1.0D0, 512.0D0, 0.99934911727904599738D0, & 1.0D0, 32.0D0, 0.98962402299599181205D0, & 1.0D0, 8.0D0, 0.95898426200345986743D0, & 1.0D0, 2.0D0, 0.84372119334725358934D0, & 1.0D0, 1.0D0, 0.70787847562782928288D0, & 3.0D0, 2.0D0, 0.59149637225671282917D0, & 2.0D0, 1.0D0, 0.49308264399053185014D0, & 5.0D0, 2.0D0, 0.41079413579749669069D0, & 3.0D0, 1.0D0, 0.34261396060786351671D0, & 4.0D0, 1.0D0, 0.24055368752127897660D0, & 17.0D0, 4.0D0, 0.22082770061202308232D0, & 5.0D0, 1.0D0, 0.17232915939014138975D0, & 11.0D0, 2.0D0, 0.14724346738730182894D0, & 6.0D0, 1.0D0, 0.12666919046715789982D0, & 8.0D0, 1.0D0, 0.74268805954862854626D-1, & 10.0D0, 1.0D0, 0.47971498020121871622D-1, & 15.0D0, 1.0D0, 0.21369201683658373846D-1, & 20.0D0, 1.0D0, 0.12020564476446432799D-1, & 30.0D0, 1.0D0, 0.53424751249537071952D-2, & 50.0D0, 1.0D0, 0.19232910450553508562D-2/ DATA DB3DAT/1.0D0, 512.0D0, 0.99926776885985461940D0, & 1.0D0, 32.0D0, 0.98833007755734698212D0, & 1.0D0, 8.0D0, 0.95390610472023510237D0, & 1.0D0, 2.0D0, 0.82496296897623372315D0, & 1.0D0, 1.0D0, 0.67441556407781468010D0, & 3.0D0, 2.0D0, 0.54710665141286285468D0, & 2.0D0, 1.0D0, 0.44112847372762418113D0, & 5.0D0, 2.0D0, 0.35413603481042394211D0, & 3.0D0, 1.0D0, 0.28357982814342246206D0, & 4.0D0, 1.0D0, 0.18173691382177474795D0, & 17.0D0, 4.0D0, 0.16277924385112436877D0, & 5.0D0, 1.0D0, 0.11759741179993396450D0, & 11.0D0, 2.0D0, 0.95240802723158889887D-1, & 6.0D0, 1.0D0, 0.77581324733763020269D-1, & 8.0D0, 1.0D0, 0.36560295673194845002D-1, & 10.0D0, 1.0D0, 0.19295765690345489563D-1, & 15.0D0, 1.0D0, 0.57712632276188798621D-2, & 20.0D0, 1.0D0, 0.24352200674805479827D-2, & 30.0D0, 1.0D0, 0.72154882216335666096D-3, & 50.0D0, 1.0D0, 0.15585454565440389896D-3/ DATA DB4DAT/1.0D0, 512.0D0, 0.99921896192761576256D0, & 1.0D0, 32.0D0, 0.98755425280996071022D0, & 1.0D0, 8.0D0, 0.95086788606389739976D0, & 1.0D0, 2.0D0, 0.81384569172034042516D0, & 1.0D0, 1.0D0, 0.65487406888673697092D0, & 3.0D0, 2.0D0, 0.52162830964878715188D0, & 2.0D0, 1.0D0, 0.41189273671788528876D0, & 5.0D0, 2.0D0, 0.32295434858707304628D0, & 3.0D0, 1.0D0, 0.25187863642883314410D0, & 4.0D0, 1.0D0, 0.15185461258672022043D0, & 17.0D0, 4.0D0, 0.13372661145921413299D0, & 5.0D0, 1.0D0, 0.91471377664481164749D-1, & 11.0D0, 2.0D0, 0.71227828197462523663D-1, & 6.0D0, 1.0D0, 0.55676547822738862783D-1, & 8.0D0, 1.0D0, 0.21967566525574960096D-1, & 10.0D0, 1.0D0, 0.96736755602711590082D-2, & 15.0D0, 1.0D0, 0.19646978158351837850D-2, & 20.0D0, 1.0D0, 0.62214648623965450200D-3, & 30.0D0, 1.0D0, 0.12289514092077854510D-3, & 50.0D0, 1.0D0, 0.15927210319002161231D-4/ DATA EX3DAT/1.0D0, 512.0D0, 0.19531249963620212007D-2, & 1.0D0, 128.0D0, 0.78124990686775522671D-2, & 1.0D0, 32.0D0, 0.31249761583499728667D-1, & 1.0D0, 8.0D0, 0.12493899888803079984D0, & 1.0D0, 2.0D0, 0.48491714311363971332D0, & 1.0D0, 1.0D0, 0.80751118213967145286D0, & 5.0D0, 4.0D0, 0.86889265412623270696D0, & 3.0D0, 2.0D0, 0.88861722235357162648D0, & 15.0D0, 8.0D0, 0.89286018500218176869D0, & 2.0D0, 1.0D0, 0.89295351429387631138D0, & 17.0D0, 8.0D0, 0.89297479112737843939D0, & 18.0D0, 8.0D0, 0.89297880579798112220D0, & 5.0D0, 2.0D0, 0.89297950317496621294D0, & 11.0D0, 4.0D0, 0.89297951152951902903D0, & 3.0D0, 1.0D0, 0.89297951156918122102D0, & 25.0D0, 8.0D0, 0.89297951156924734716D0, & 13.0D0, 4.0D0, 0.89297951156924917298D0, & 7.0D0, 2.0D0, 0.89297951156924921121D0, & 15.0D0, 4.0D0, 0.89297951156924921122D0, & 4.0D0, 1.0D0, 0.89297951156924921122D0/ DATA GSTDAT/1.0D0, 512.0D0, 0.59531540040441651584D1, & 1.0D0, 128.0D0, 0.45769601268624494109D1, & 1.0D0, 32.0D0, 0.32288921331902217638D1, & 1.0D0, 8.0D0, 0.19746110873568719362D1, & 1.0D0, 2.0D0, 0.96356046208697728563D0, & 1.0D0, 1.0D0, 0.60513365250334458174D0, & 5.0D0, 4.0D0, 0.51305506459532198016D0, & 3.0D0, 2.0D0, 0.44598602820946133091D0, & 15.0D0, 8.0D0, 0.37344458206879749357D0, & 2.0D0, 1.0D0, 0.35433592884953063055D0, & 17.0D0, 8.0D0, 0.33712156518881920994D0, & 5.0D0, 2.0D0, 0.29436170729362979176D0, & 3.0D0, 1.0D0, 0.25193499644897222840D0, & 7.0D0, 2.0D0, 0.22028778222123939276D0, & 4.0D0, 1.0D0, 0.19575258237698917033D0, & 9.0D0, 2.0D0, 0.17616303166670699424D0, & 5.0D0, 1.0D0, 0.16015469479664778673D0, & 23.0D0, 4.0D0, 0.14096116876193391066D0, & 6.0D0, 1.0D0, 0.13554987191049066274D0, & 7.0D0, 1.0D0, 0.11751605060085098084D0/ DATA I0IDAT/1.0D0, 512.0D0, 0.19531256208818052282D-2, & -1.0D0, 256.0D0, -0.39062549670565734544D-2, & 1.0D0, 16.0D0, 0.62520348032546565850D-1, & 1.0D0, 8.0D0, 0.12516285581366971819D0, & -1.0D0, 2.0D0, -0.51051480879740303760D0, & 1.0D0, 1.0D0, 0.10865210970235898158D1, & 2.0D0, 1.0D0, 0.27750019054282535299D1, & -4.0D0, 1.0D0, -0.13775208868039716639D2, & 8.0D0, 1.0D0, 0.46424372058106108576D3, & 18.0D0, 1.0D0, 0.64111867658021584522D7, & -37.0D0, 2.0D0, -0.10414860803175857953D8, & 20.0D0, 1.0D0, 0.44758598913855743089D8, & -21.0D0, 1.0D0, -0.11852985311558287888D9, & 22.0D0, 1.0D0, 0.31430078220715992752D9, & -23.0D0, 1.0D0, -0.83440212900794309620D9, & 24.0D0, 1.0D0, 0.22175367579074298261D10, & 25.0D0, 1.0D0, 0.58991731842803636487D10, & -27.0D0, 1.0D0, -0.41857073244691522147D11, & 30.0D0, 1.0D0, 0.79553885818472357663D12, & 40.0D0, 1.0D0, 0.15089715082719201025D17/ DATA IL0DAT/1.0D0, 512.0D0, 0.99875755515461749793D0, & 1.0D0, 64.0D0, 0.99011358230706643807D0, & 1.0D0, 8.0D0, 0.92419435310023947018D0, & 1.0D0, 2.0D0, 0.73624267134714273902D0, & 1.0D0, 1.0D0, 0.55582269181411744686D0, & 2.0D0, 1.0D0, 0.34215154434462160628D0, & 4.0D0, 1.0D0, 0.17087174888774706539D0, & 8.0D0, 1.0D0, 0.81081008709219208918D-1, & 12.0D0, 1.0D0, 0.53449421441089580702D-1, & 16.0D0, 1.0D0, 0.39950321008923244846D-1, & 65.0D0, 4.0D0, 0.39330637437584921392D-1, & 17.0D0, 1.0D0, 0.37582274342808670750D-1, & 20.0D0, 1.0D0, 0.31912486554480390343D-1, & 25.0D0, 1.0D0, 0.25506146883504738403D-1, & 30.0D0, 1.0D0, 0.21244480317825292412D-1, & 40.0D0, 1.0D0, 0.15925498348551684335D-1, & 50.0D0, 1.0D0, 0.12737506927242585015D-1, & 75.0D0, 1.0D0, 0.84897750814784916847D-2, & 100.0D0, 1.0D0, 0.63668349178454469153D-2, & 125.0D0, 1.0D0, 0.50932843163122551114D-2/ DATA IL1DAT/1.0D0, 512.0D0, 0.97575346155386267134D-3, & 1.0D0, 64.0D0, 0.77609293280609272733D-2, & 1.0D0, 8.0D0, 0.59302966404545373770D-1, & 1.0D0, 2.0D0, 0.20395212276737365307D0, & 1.0D0, 1.0D0, 0.33839472293667639038D0, & 2.0D0, 1.0D0, 0.48787706726961324579D0, & 4.0D0, 1.0D0, 0.59018734196576517506D0, & 8.0D0, 1.0D0, 0.62604539530312149476D0, & 12.0D0, 1.0D0, 0.63209315274909764698D0, & 16.0D0, 1.0D0, 0.63410179313235359215D0, & 65.0D0, 4.0D0, 0.63417966797578128188D0, & 17.0D0, 1.0D0, 0.63439268632392089434D0, & 20.0D0, 1.0D0, 0.63501579073257770690D0, & 25.0D0, 1.0D0, 0.63559616677359459337D0, & 30.0D0, 1.0D0, 0.63591001826697110312D0, & 40.0D0, 1.0D0, 0.63622113181751073643D0, & 50.0D0, 1.0D0, 0.63636481702133606597D0, & 75.0D0, 1.0D0, 0.63650653499619902120D0, & 100.0D0, 1.0D0, 0.63655609126300261851D0, & 125.0D0, 1.0D0, 0.63657902087183929223D0/ DATA J0IDAT/1.0D0, 1024.0D0, 0.97656242238978822427D-3, & 1.0D0, 256.0D0, 0.39062450329491108875D-2, & -1.0D0, 16.0D0, -0.62479657927917933620D-1, & 1.0D0, 8.0D0, 0.12483733492120479139D0, & -1.0D0, 2.0D0, -0.48968050664604505505D0, & 1.0D0, 1.0D0, 0.91973041008976023931D0, & -2.0D0, 1.0D0, -0.14257702931970265690D1, & 4.0D0, 1.0D0, 0.10247341594606064818D1, & -8.0D0, 1.0D0, -0.12107468348304501655D1, & 16.0D0, 1.0D0, 0.11008652032736190799D1, & -33.0D0, 2.0D0, -0.10060334829904124192D1, & 18.0D0, 1.0D0, 0.81330572662485953519D0, & -20.0D0, 1.0D0, -0.10583788214211277585D1, & 25.0D0, 1.0D0, 0.87101492116545875169D0, & -30.0D0, 1.0D0, -0.88424908882547488420D0, & 40.0D0, 1.0D0, 0.11257761503599914603D1, & -50.0D0, 1.0D0, -0.90141212258183461184D0, & 75.0D0, 1.0D0, 0.91441344369647797803D0, & -80.0D0, 1.0D0, -0.94482281938334394886D0, & 100.0D0, 1.0D0, 0.92266255696016607257D0/ DATA K0IDAT/1.0D0, 1024.0D0, 0.78587929563466784589D-2, & 1.0D0, 256.0D0, 0.26019991617330578111D-1, & 1.0D0, 16.0D0, 0.24311842237541167904D0, & 1.0D0, 8.0D0, 0.39999633750480508861D0, & 1.0D0, 2.0D0, 0.92710252093114907345D0, & 1.0D0, 1.0D0, 0.12425098486237782662D1, & 2.0D0, 1.0D0, 0.14736757343168286825D1, & 4.0D0, 1.0D0, 0.15606495706051741364D1, & 5.0D0, 1.0D0, 0.15673873907283660493D1, & 6.0D0, 1.0D0, 0.15696345532693743714D1, & 13.0D0, 2.0D0, 0.15701153443250786355D1, & 8.0D0, 1.0D0, 0.15706574852894436220D1, & 10.0D0, 1.0D0, 0.15707793116159788598D1, & 12.0D0, 1.0D0, 0.15707942066465767196D1, & 15.0D0, 1.0D0, 0.15707962315469192247D1, & 20.0D0, 1.0D0, 0.15707963262340149876D1, & 30.0D0, 1.0D0, 0.15707963267948756308D1, & 50.0D0, 1.0D0, 0.15707963267948966192D1, & 80.0D0, 1.0D0, 0.15707963267948966192D1, & 100.0D0, 1.0D0, 0.15707963267948966192D1/ DATA LOBDAT/1.0D0, 512.0D0, 0.12417639065161393857D-8, & 1.0D0, 128.0D0, 0.79473344770001088225D-7, & 1.0D0, 32.0D0, 0.50867598186208834198D-5, & 1.0D0, 8.0D0, 0.32603097901207200319D-3, & 1.0D0, 2.0D0, 0.21380536815408214419D-1, & 1.0D0, 1.0D0, 0.18753816902083824050D0, & 3.0D0, 2.0D0, 0.83051199971883645115D0, & 2.0D0, 1.0D0, 0.18854362426679034904D1, & 5.0D0, 2.0D0, 0.21315988986516411053D1, & 3.0D0, 1.0D0, 0.21771120185613427221D1, & 4.0D0, 1.0D0, 0.22921027921896650849D1, & 5.0D0, 1.0D0, 0.39137195028784495586D1, & 6.0D0, 1.0D0, 0.43513563983836427904D1, & 7.0D0, 1.0D0, 0.44200644968478185898D1, & 10.0D0, 1.0D0, 0.65656013133623829156D1, & 15.0D0, 1.0D0, 0.10825504661504599479D2, & 20.0D0, 1.0D0, 0.13365512855474227325D2, & 30.0D0, 1.0D0, 0.21131002685639959927D2, & 50.0D0, 1.0D0, 0.34838236589449117389D2, & 100.0D0, 1.0D0, 0.69657062437837394278D2/ DATA STRDAT/1.0D0, 512.0D0, 0.21901065985698662316D-15, & 1.0D0, 128.0D0, 0.22481399438625244761D-12, & 1.0D0, 32.0D0, 0.23245019579558857124D-9, & 1.0D0, 8.0D0, 0.24719561475975007037D-6, & 1.0D0, 2.0D0, 0.28992610989833245669D-3, & 1.0D0, 1.0D0, 0.10698146390809715091D-1, & 3.0D0, 2.0D0, 0.89707650964424730705D-1, & 2.0D0, 1.0D0, 0.40049605719592888440D0, & 3.0D0, 1.0D0, 0.30504104398079096598D1, & 4.0D0, 1.0D0, 0.11367704858439426431D2, & 33.0D0, 8.0D0, 0.12960679405324786954D2, & 9.0D0, 2.0D0, 0.18548713944748505675D2, & 5.0D0, 1.0D0, 0.27866273821903121400D2, & 6.0D0, 1.0D0, 0.51963334071699323351D2, & 8.0D0, 1.0D0, 0.10861016747891228129D3, & 10.0D0, 1.0D0, 0.15378903316556621624D3, & 15.0D0, 1.0D0, 0.19302665532558721516D3, & 20.0D0, 1.0D0, 0.19636850166006541482D3, & 30.0D0, 1.0D0, 0.19651946766008214217D3, & 50.0D0, 1.0D0, 0.19651956920868316152D3/ DATA SH0DAT/1.0D0, 512.0D0, 0.12433974658847434366D-2, & -1.0D0, 128.0D0, -0.49735582423748415045D-2, & 1.0D0, 16.0D0, 0.39771469054536941564D-1, & -1.0D0, 4.0D0, -0.15805246001653314198D0, & 1.0D0, 1.0D0, 0.56865662704828795099D0, & 5.0D0, 4.0D0, 0.66598399314899916605D0, & 2.0D0, 1.0D0, 0.79085884950809589255D0, & -4.0D0, 1.0D0, -0.13501457342248639716D0, & 15.0D0, 2.0D0, 0.20086479668164503137D0, & 11.0D0, 1.0D0, -0.11142097800261991552D0, & 23.0D0, 2.0D0, -0.17026804865989885869D0, & -16.0D0, 1.0D0, -0.13544931808186467594D0, & 20.0D0, 1.0D0, 0.94393698081323450897D-1, & 25.0D0, 1.0D0, -0.10182482016001510271D0, & -30.0D0, 1.0D0, 0.96098421554162110012D-1, & 50.0D0, 1.0D0, -0.85337674826118998952D-1, & 75.0D0, 1.0D0, -0.76882290637052720045D-1, & -80.0D0, 1.0D0, 0.47663833591418256339D-1, & 100.0D0, 1.0D0, -0.70878751689647343204D-1, & -125.0D0, 1.0D0, 0.65752908073352785368D-1/ DATA SH1DAT/1.0D0, 512.0D0, 0.80950369576367526071D-6, & -1.0D0, 128.0D0, 0.12952009724113229165D-4, & 1.0D0, 16.0D0, 0.82871615165407083021D-3, & -1.0D0, 4.0D0, 0.13207748375849572564D-1, & 1.0D0, 1.0D0, 0.19845733620194439894D0, & 5.0D0, 4.0D0, 0.29853823231804706294D0, & 2.0D0, 1.0D0, 0.64676372828356211712D0, & -4.0D0, 1.0D0, 0.10697266613089193593D1, & 15.0D0, 2.0D0, 0.38831308000420560970D0, & 9.0D0, 1.0D0, 0.74854243745107710333D0, & 19.0D0, 2.0D0, 0.84664854642567359993D0, & -12.0D0, 1.0D0, 0.58385732464244384564D0, & 17.0D0, 1.0D0, 0.80600584524215772824D0, & 25.0D0, 1.0D0, 0.53880362132692947616D0, & -30.0D0, 1.0D0, 0.72175037834698998506D0, & 50.0D0, 1.0D0, 0.58007844794544189900D0, & 75.0D0, 1.0D0, 0.60151910385440804463D0, & -80.0D0, 1.0D0, 0.70611511147286827018D0, & 100.0D0, 1.0D0, 0.61631110327201338454D0, & -125.0D0, 1.0D0, 0.62778480765443656489D0/ DATA SL0DAT/1.0D0, 512.0D0, 0.12433985199262820188D-2, & -1.0D0, 32.0D0, -0.19896526647882937004D-1, & 1.0D0, 8.0D0, 0.79715713253115014945D-1, & -1.0D0, 2.0D0, -0.32724069939418078025D0, & 1.0D0, 1.0D0, 0.71024318593789088874D0, & 2.0D0, 1.0D0, 0.19374337579914456612D1, & -4.0D0, 1.0D0, -0.11131050203248583431D2, & 7.0D0, 1.0D0, 0.16850062034703267148D3, & -10.0D0, 1.0D0, -0.28156522493745948555D4, & 16.0D0, 1.0D0, 0.89344618796978400815D6, & 65.0D0, 4.0D0, 0.11382025002851451057D7, & -17.0D0, 1.0D0, -0.23549701855860190304D7, & 20.0D0, 1.0D0, 0.43558282527641046718D8, & 45.0D0, 2.0D0, 0.49993516476037957165D9, & -25.0D0, 1.0D0, -0.57745606064408041689D10, & 30.0D0, 1.0D0, 0.78167229782395624524D12, & -40.0D0, 1.0D0, -0.14894774793419899908D17, & 50.0D0, 1.0D0, 0.29325537838493363267D21, & 60.0D0, 1.0D0, 0.58940770556098011683D25, & -70.0D0, 1.0D0, -0.12015889579125463605D30/ DATA SL1DAT/1.0D0, 512.0D0, 0.80950410749865126939D-6, & -1.0D0, 32.0D0, 0.20724649092571514607D-3, & 1.0D0, 8.0D0, 0.33191834066894516744D-2, & -1.0D0, 2.0D0, 0.53942182623522663292D-1, & 1.0D0, 1.0D0, 0.22676438105580863683D0, & 2.0D0, 1.0D0, 0.11027597873677158176D1, & -4.0D0, 1.0D0, 0.91692778117386847344D1, & 7.0D0, 1.0D0, 0.15541656652426660966D3, & -10.0D0, 1.0D0, 0.26703582852084829694D4, & 16.0D0, 1.0D0, 0.86505880175304633906D6, & 65.0D0, 4.0D0, 0.11026046613094942620D7, & -17.0D0, 1.0D0, 0.22846209494153934787D7, & 20.0D0, 1.0D0, 0.42454972750111979449D8, & 45.0D0, 2.0D0, 0.48869614587997695539D9, & -25.0D0, 1.0D0, 0.56578651292431051863D10, & 30.0D0, 1.0D0, 0.76853203893832108948D12, & -40.0D0, 1.0D0, 0.14707396163259352103D17, & 50.0D0, 1.0D0, 0.29030785901035567967D21, & 60.0D0, 1.0D0, 0.58447515883904682813D25, & -70.0D0, 1.0D0, 0.11929750788892311875D30/ DATA SY1DAT/1.0D0, 512.0D0, 0.26514864547487397044D0, & 1.0D0, 32.0D0, 0.62050129979079045645D0, & 1.0D0, 8.0D0, 0.85112572132368011206D0, & 1.0D0, 2.0D0, 0.87081914687546885094D0, & 1.0D0, 1.0D0, 0.65142281535536396975D0, & 3.0D0, 2.0D0, 0.45064040920322354579D0, & 2.0D0, 1.0D0, 0.30163590285073940285D0, & 5.0D0, 2.0D0, 0.19814490804441305867D0, & 3.0D0, 1.0D0, 0.12856571000906381300D0, & 4.0D0, 1.0D0, 0.52827396697866818297D-1, & 17.0D0, 4.0D0, 0.42139298471720305542D-1, & 5.0D0, 1.0D0, 0.21248129774981984268D-1, & 11.0D0, 2.0D0, 0.13400258907505536491D-1, & 6.0D0, 1.0D0, 0.84260797314108699935D-2, & 8.0D0, 1.0D0, 0.12884516186754671469D-2, & 10.0D0, 1.0D0, 0.19223826430086897418D-3, & 12.0D0, 1.0D0, 0.28221070834007689394D-4, & 15.0D0, 1.0D0, 0.15548757973038189372D-5, & 20.0D0, 1.0D0, 0.11968634456097453636D-7, & 25.0D0, 1.0D0, 0.89564246772237127742D-10/ DATA SY2DAT/1.0D0, 512.0D0, 0.13430727275667378338D0, & 1.0D0, 32.0D0, 0.33485265272424176976D0, & 1.0D0, 8.0D0, 0.50404224110911078651D0, & 1.0D0, 2.0D0, 0.60296523236016785113D0, & 1.0D0, 1.0D0, 0.49447506210420826699D0, & 3.0D0, 2.0D0, 0.36036067860473360389D0, & 2.0D0, 1.0D0, 0.24967785497625662113D0, & 5.0D0, 2.0D0, 0.16813830542905833533D0, & 3.0D0, 1.0D0, 0.11117122348556549832D0, & 4.0D0, 1.0D0, 0.46923205826101330711D-1, & 17.0D0, 4.0D0, 0.37624545861980001482D-1, & 5.0D0, 1.0D0, 0.19222123172484106436D-1, & 11.0D0, 2.0D0, 0.12209535343654701398D-1, & 6.0D0, 1.0D0, 0.77249644268525771866D-2, & 8.0D0, 1.0D0, 0.12029044213679269639D-2, & 10.0D0, 1.0D0, 0.18161187569530204281D-3, & 12.0D0, 1.0D0, 0.26884338006629353506D-4, & 15.0D0, 1.0D0, 0.14942212731345828759D-5, & 20.0D0, 1.0D0, 0.11607696854385161390D-7, & 25.0D0, 1.0D0, 0.87362343746221526073D-10/ DATA TR2DAT/1.0D0, 512.0D0, 0.19531247930394515480D-2, & 1.0D0, 32.0D0, 0.31249152314331109004D-1, & 1.0D0, 8.0D0, 0.12494577194783451032D0, & 1.0D0, 2.0D0, 0.49655363615640595865D0, & 1.0D0, 1.0D0, 0.97303256135517012845D0, & 3.0D0, 2.0D0, 0.14121978695932525805D1, & 2.0D0, 1.0D0, 0.18017185674405776809D1, & 5.0D0, 2.0D0, 0.21350385339277043015D1, & 3.0D0, 1.0D0, 0.24110500490169534620D1, & 4.0D0, 1.0D0, 0.28066664045631179931D1, & 17.0D0, 4.0D0, 0.28777421863296234131D1, & 5.0D0, 1.0D0, 0.30391706043438554330D1, & 11.0D0, 2.0D0, 0.31125074928667355940D1, & 6.0D0, 1.0D0, 0.31656687817738577185D1, & 8.0D0, 1.0D0, 0.32623520367816009184D1, & 10.0D0, 1.0D0, 0.32843291144979517358D1, & 15.0D0, 1.0D0, 0.32897895167775788137D1, & 20.0D0, 1.0D0, 0.32898672226665499687D1, & 30.0D0, 1.0D0, 0.32898681336064325400D1, & 50.0D0, 1.0D0, 0.32898681336964528724D1/ DATA TR3DAT/1.0D0, 512.0D0, 0.19073483296476379584D-5, & 1.0D0, 32.0D0, 0.48826138243180786081D-3, & 1.0D0, 8.0D0, 0.78074163848431205820D-2, & 1.0D0, 2.0D0, 0.12370868718812031049D0, & 1.0D0, 1.0D0, 0.47984100657241749994D0, & 3.0D0, 2.0D0, 0.10269431622039754738D1, & 2.0D0, 1.0D0, 0.17063547219458658863D1, & 5.0D0, 2.0D0, 0.24539217444475937661D1, & 3.0D0, 1.0D0, 0.32106046629422467723D1, & 4.0D0, 1.0D0, 0.45792174372291563703D1, & 17.0D0, 4.0D0, 0.48722022832940370805D1, & 5.0D0, 1.0D0, 0.56143866138422732286D1, & 11.0D0, 2.0D0, 0.59984455864575470009D1, & 6.0D0, 1.0D0, 0.63033953673480961120D1, & 8.0D0, 1.0D0, 0.69579908688361166266D1, & 10.0D0, 1.0D0, 0.71503227120085929750D1, & 15.0D0, 1.0D0, 0.72110731475871876393D1, & 20.0D0, 1.0D0, 0.72123221966388461839D1, & 30.0D0, 1.0D0, 0.72123414161609465119D1, & 50.0D0, 1.0D0, 0.72123414189575656868D1/ DATA TR4DAT/1.0D0, 512.0D0, 0.24835263919461834041D-8, & 1.0D0, 32.0D0, 0.10172029353616724881D-4, & 1.0D0, 8.0D0, 0.65053332405940765479D-3, & 1.0D0, 2.0D0, 0.41150448004155727767D-1, & 1.0D0, 1.0D0, 0.31724404523442648241D0, & 3.0D0, 2.0D0, 0.10079442901142373591D1, & 2.0D0, 1.0D0, 0.22010881024333408363D1, & 5.0D0, 2.0D0, 0.38846508619156545210D1, & 3.0D0, 1.0D0, 0.59648223973714765245D1, & 4.0D0, 1.0D0, 0.10731932392998622219D2, & 17.0D0, 4.0D0, 0.11940028876819364777D2, & 5.0D0, 1.0D0, 0.15359784316882182982D2, & 11.0D0, 2.0D0, 0.17372587633093742893D2, & 6.0D0, 1.0D0, 0.19122976016053166969D2, & 8.0D0, 1.0D0, 0.23583979156921941515D2, & 10.0D0, 1.0D0, 0.25273667677030441733D2, & 15.0D0, 1.0D0, 0.25955198214572256372D2, & 20.0D0, 1.0D0, 0.25975350935212241910D2, & 30.0D0, 1.0D0, 0.25975757522084093747D2, & 50.0D0, 1.0D0, 0.25975757609067315288D2/ DATA TR5DAT/1.0D0, 512.0D0, 0.36379780361036116971D-11, & 1.0D0, 32.0D0, 0.23840564453948442379D-6, & 1.0D0, 8.0D0, 0.60982205372226969189D-4, & 1.0D0, 2.0D0, 0.15410004586376649337D-1, & 1.0D0, 1.0D0, 0.23661587923909478926D0, & 3.0D0, 2.0D0, 0.11198756851307629651D1, & 2.0D0, 1.0D0, 0.32292901663684049171D1, & 5.0D0, 2.0D0, 0.70362973105160654056D1, & 3.0D0, 1.0D0, 0.12770557691044159511D2, & 4.0D0, 1.0D0, 0.29488339015245845447D2, & 17.0D0, 4.0D0, 0.34471340540362254586D2, & 5.0D0, 1.0D0, 0.50263092218175187785D2, & 11.0D0, 2.0D0, 0.60819909101127165207D2, & 6.0D0, 1.0D0, 0.70873334429213460498D2, & 8.0D0, 1.0D0, 0.10147781242977788097D3, & 10.0D0, 1.0D0, 0.11638074540242071077D3, & 15.0D0, 1.0D0, 0.12409623901262967878D3, & 20.0D0, 1.0D0, 0.12442270155632550228D3, & 30.0D0, 1.0D0, 0.12443132790838589548D3, & 50.0D0, 1.0D0, 0.12443133061720432435D3/ DATA TR6DAT/1.0D0, 512.0D0, 0.56843405953641209574D-14, & 1.0D0, 32.0D0, 0.59601180165247401484D-8, & 1.0D0, 8.0D0, 0.60978424397580572815D-5, & 1.0D0, 2.0D0, 0.61578909866319494394D-2, & 1.0D0, 1.0D0, 0.18854360275680840514D0, & 3.0D0, 2.0D0, 0.13319251347921659134D1, & 2.0D0, 1.0D0, 0.50857202271697616755D1, & 5.0D0, 2.0D0, 0.13729222365466557122D2, & 3.0D0, 1.0D0, 0.29579592481641441292D2, & 4.0D0, 1.0D0, 0.88600835706899853768D2, & 17.0D0, 4.0D0, 0.10916037113373004909D3, & 5.0D0, 1.0D0, 0.18224323749575359518D3, & 11.0D0, 2.0D0, 0.23765383125586756031D3, & 6.0D0, 1.0D0, 0.29543246745959381136D3, & 8.0D0, 1.0D0, 0.50681244381280455592D3, & 10.0D0, 1.0D0, 0.63878231134946125623D3, & 15.0D0, 1.0D0, 0.72699203556994876111D3, & 20.0D0, 1.0D0, 0.73230331643146851717D3, & 30.0D0, 1.0D0, 0.73248692015882096369D3, & 50.0D0, 1.0D0, 0.73248700462879996604D3/ DATA TR7DAT/1.0D0, 512.0D0, 0.92518563327283409427D-17, & 1.0D0, 32.0D0, 0.15521095556949867541D-9, & 1.0D0, 8.0D0, 0.63516238373841716290D-6, & 1.0D0, 2.0D0, 0.25638801246626135714D-2, & 1.0D0, 1.0D0, 0.15665328993811649746D0, & 3.0D0, 2.0D0, 0.16538225039181097423D1, & 2.0D0, 1.0D0, 0.83763085709508211054D1, & 5.0D0, 2.0D0, 0.28078570717830763747D2, & 3.0D0, 1.0D0, 0.72009676046751991365D2, & 4.0D0, 1.0D0, 0.28174905701691911450D3, & 17.0D0, 4.0D0, 0.36660227975327792529D3, & 5.0D0, 1.0D0, 0.70556067982603601123D3, & 11.0D0, 2.0D0, 0.99661927562755629434D3, & 6.0D0, 1.0D0, 0.13288914430417403901D4, & 8.0D0, 1.0D0, 0.27987640273169129925D4, & 10.0D0, 1.0D0, 0.39721376409416504325D4, & 15.0D0, 1.0D0, 0.49913492839319899726D4, & 20.0D0, 1.0D0, 0.50781562639825019000D4, & 30.0D0, 1.0D0, 0.50820777202028708434D4, & 50.0D0, 1.0D0, 0.50820803580047164618D4/ DATA TR8DAT/1.0D0, 512.0D0, 0.15488598634539359463D-19, & 1.0D0, 32.0D0, 0.41574269117845953797D-11, & 1.0D0, 8.0D0, 0.68050651245227411689D-7, & 1.0D0, 2.0D0, 0.10981703519563009836D-2, & 1.0D0, 1.0D0, 0.13396432776187883834D0, & 3.0D0, 2.0D0, 0.21153387806998617182D1, & 2.0D0, 1.0D0, 0.14227877028750735641D2, & 5.0D0, 2.0D0, 0.59312061431647843226D2, & 3.0D0, 1.0D0, 0.18139614577043147745D3, & 4.0D0, 1.0D0, 0.93148001928992220863D3, & 17.0D0, 4.0D0, 0.12817928112604611804D4, & 5.0D0, 1.0D0, 0.28572838386329242218D4, & 11.0D0, 2.0D0, 0.43872971687877730010D4, & 6.0D0, 1.0D0, 0.62993229139406657611D4, & 8.0D0, 1.0D0, 0.16589426277154888511D5, & 10.0D0, 1.0D0, 0.27064780798797398935D5, & 15.0D0, 1.0D0, 0.38974556062543661284D5, & 20.0D0, 1.0D0, 0.40400240716905025786D5, & 30.0D0, 1.0D0, 0.40484316504120655568D5, & 50.0D0, 1.0D0, 0.40484399001892184901D5/ DATA TR9DAT/1.0D0, 512.0D0, 0.26469772870084897671D-22, & 1.0D0, 32.0D0, 0.11367943653594246210D-12, & 1.0D0, 8.0D0, 0.74428246255329800255D-8, & 1.0D0, 2.0D0, 0.48022728485415366194D-3, & 1.0D0, 1.0D0, 0.11700243014358676725D0, & 3.0D0, 2.0D0, 0.27648973910899914391D1, & 2.0D0, 1.0D0, 0.24716631405829192997D2, & 5.0D0, 2.0D0, 0.12827119828849828583D3, & 3.0D0, 1.0D0, 0.46842894800662208986D3, & 4.0D0, 1.0D0, 0.31673967371627895718D4, & 17.0D0, 4.0D0, 0.46140886546630195390D4, & 5.0D0, 1.0D0, 0.11952718545392302185D5, & 11.0D0, 2.0D0, 0.20001612666477027728D5, & 6.0D0, 1.0D0, 0.31011073271851366554D5, & 8.0D0, 1.0D0, 0.10352949905541130133D6, & 10.0D0, 1.0D0, 0.19743173017140591390D6, & 15.0D0, 1.0D0, 0.33826030414658460679D6, & 20.0D0, 1.0D0, 0.36179607036750755227D6, & 30.0D0, 1.0D0, 0.36360622124777561525D6, & 50.0D0, 1.0D0, 0.36360880558827162725D6/ DATA Y0IDAT/1.0D0, 512.0D0, -0.91442642860172110926D-2, & 1.0D0, 128.0D0, -0.29682047390397591290D-1, & 1.0D0, 8.0D0, -0.25391431276585388961D0, & 1.0D0, 2.0D0, -0.56179545591464028187D0, & 1.0D0, 1.0D0, -0.63706937660742309754D0, & 2.0D0, 1.0D0, -0.28219285008510084123D0, & 4.0D0, 1.0D0, 0.38366964785312561103D0, & 6.0D0, 1.0D0, -0.12595061285798929390D0, & 10.0D0, 1.0D0, 0.24129031832266684828D0, & 16.0D0, 1.0D0, 0.17138069757627037938D0, & 65.0D0, 4.0D0, 0.18958142627134083732D0, & 17.0D0, 1.0D0, 0.17203846136449706946D0, & 20.0D0, 1.0D0, -0.16821597677215029611D0, & 25.0D0, 1.0D0, -0.93607927351428988679D-1, & 30.0D0, 1.0D0, 0.88229711948036648408D-1, & 40.0D0, 1.0D0, -0.89324662736274161841D-2, & 50.0D0, 1.0D0, -0.54814071000063488284D-1, & 70.0D0, 1.0D0, -0.94958246003466381588D-1, & 100.0D0, 1.0D0, -0.19598064853404969850D-1, & 125.0D0, 1.0D0, -0.83084772357154773468D-2/ C C TEST ABRAM0 C FNAME = 'ABRAM0' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 100 I = 1 , 20 NUM = AB0DAT ( 1 , I ) DEN = AB0DAT ( 2 , I ) RES = AB0DAT ( 3 , I ) PT = NUM / DEN COMP = ABRAM0(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 100 CONTINUE C C TEST ABRAM1 C FNAME = 'ABRAM1' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 110 I = 1 , 20 NUM = AB1DAT ( 1 , I ) DEN = AB1DAT ( 2 , I ) RES = AB1DAT ( 3 , I ) PT = NUM / DEN COMP = ABRAM1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 110 CONTINUE C C TEST ABRAM2 C FNAME = 'ABRAM2' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 120 I = 1 , 20 NUM = AB2DAT ( 1 , I ) DEN = AB2DAT ( 2 , I ) RES = AB2DAT ( 3 , I ) PT = NUM / DEN COMP = ABRAM2(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 120 CONTINUE C C TEST AIRINT C FNAME = 'AIRINT' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 130 I = 1 , 20 NUM = AIIDAT ( 1 , I ) DEN = AIIDAT ( 2 , I ) RES = AIIDAT ( 3 , I ) PT = NUM / DEN COMP = AIRINT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 130 CONTINUE C C TEST AIRYGI C FNAME = 'AIRYGI' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 140 I = 1 , 20 NUM = AGIDAT ( 1 , I ) DEN = AGIDAT ( 2 , I ) RES = AGIDAT ( 3 , I ) PT = NUM / DEN COMP = AIRYGI(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 140 CONTINUE C C TEST AIRYHI C FNAME = 'AIRYHI' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 150 I = 1 , 20 NUM = AHIDAT ( 1 , I ) DEN = AHIDAT ( 2 , I ) RES = AHIDAT ( 3 , I ) PT = NUM / DEN COMP = AIRYHI(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 150 CONTINUE C C TEST ATNINT C FNAME = 'ATNINT' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 160 I = 1 , 20 NUM = ATNDAT ( 1 , I ) DEN = ATNDAT ( 2 , I ) RES = ATNDAT ( 3 , I ) PT = NUM / DEN COMP = ATNINT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 160 CONTINUE C C TEST BIRINT C FNAME = 'BIRINT' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 170 I = 1 , 20 NUM = BIIDAT ( 1 , I ) DEN = BIIDAT ( 2 , I ) RES = BIIDAT ( 3 , I ) PT = NUM / DEN COMP = BIRINT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 170 CONTINUE C C TEST CLAUSN C FNAME = 'CLAUSN' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 180 I = 1 , 20 NUM = CLNDAT ( 1 , I ) DEN = CLNDAT ( 2 , I ) RES = CLNDAT ( 3 , I ) PT = NUM / DEN COMP = CLAUSN(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 180 CONTINUE C C TEST DEBYE1 C FNAME = 'DEBYE1' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 190 I = 1 , 20 NUM = DB1DAT ( 1 , I ) DEN = DB1DAT ( 2 , I ) RES = DB1DAT ( 3 , I ) PT = NUM / DEN COMP = DEBYE1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 190 CONTINUE C C TEST DEBYE2 C FNAME = 'DEBYE2' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 200 I = 1 , 20 NUM = DB2DAT ( 1 , I ) DEN = DB2DAT ( 2 , I ) RES = DB2DAT ( 3 , I ) PT = NUM / DEN COMP = DEBYE2(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 200 CONTINUE C C TEST DEBYE3 C FNAME = 'DEBYE3' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 210 I = 1 , 20 NUM = DB3DAT ( 1 , I ) DEN = DB3DAT ( 2 , I ) RES = DB3DAT ( 3 , I ) PT = NUM / DEN COMP = DEBYE3(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 210 CONTINUE C C TEST DEBYE4 C FNAME = 'DEBYE4' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 220 I = 1 , 20 NUM = DB4DAT ( 1 , I ) DEN = DB4DAT ( 2 , I ) RES = DB4DAT ( 3 , I ) PT = NUM / DEN COMP = DEBYE4(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 220 CONTINUE C C TEST EXP3 C FNAME = 'EXP3 ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 230 I = 1 , 20 NUM = EX3DAT ( 1 , I ) DEN = EX3DAT ( 2 , I ) RES = EX3DAT ( 3 , I ) PT = NUM / DEN COMP = EXP3(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 230 CONTINUE C C TEST GOODST C FNAME = 'GOODST' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 240 I = 1 , 20 NUM = GSTDAT ( 1 , I ) DEN = GSTDAT ( 2 , I ) RES = GSTDAT ( 3 , I ) PT = NUM / DEN COMP = GOODST(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 240 CONTINUE C C TEST I0INT C FNAME = 'I0INT ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 250 I = 1 , 20 NUM = I0IDAT ( 1 , I ) DEN = I0IDAT ( 2 , I ) RES = I0IDAT ( 3 , I ) PT = NUM / DEN COMP = I0INT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 250 CONTINUE C C TEST I0ML0 C FNAME = 'I0ML0 ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 260 I = 1 , 20 NUM = IL0DAT ( 1 , I ) DEN = IL0DAT ( 2 , I ) RES = IL0DAT ( 3 , I ) PT = NUM / DEN COMP = I0ML0(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 260 CONTINUE C C TEST I1ML1 C FNAME = 'I1ML1 ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 270 I = 1 , 20 NUM = IL1DAT ( 1 , I ) DEN = IL1DAT ( 2 , I ) RES = IL1DAT ( 3 , I ) PT = NUM / DEN COMP = I1ML1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 270 CONTINUE C C TEST J0INT C FNAME = 'J0INT ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 280 I = 1 , 20 NUM = J0IDAT ( 1 , I ) DEN = J0IDAT ( 2 , I ) RES = J0IDAT ( 3 , I ) PT = NUM / DEN COMP = J0INT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 280 CONTINUE C C TEST K0INT C FNAME = 'K0INT ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 290 I = 1 , 20 NUM = K0IDAT ( 1 , I ) DEN = K0IDAT ( 2 , I ) RES = K0IDAT ( 3 , I ) PT = NUM / DEN COMP = K0INT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 290 CONTINUE C C TEST LOBACH C FNAME = 'LOBACH' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 300 I = 1 , 20 NUM = LOBDAT ( 1 , I ) DEN = LOBDAT ( 2 , I ) RES = LOBDAT ( 3 , I ) PT = NUM / DEN COMP = LOBACH(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 300 CONTINUE C C TEST STROM C FNAME = 'STROM ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 310 I = 1 , 20 NUM = STRDAT ( 1 , I ) DEN = STRDAT ( 2 , I ) RES = STRDAT ( 3 , I ) PT = NUM / DEN COMP = STROM(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 310 CONTINUE C C TEST STRVH0 C FNAME = 'STRVH0' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 320 I = 1 , 20 NUM = SH0DAT ( 1 , I ) DEN = SH0DAT ( 2 , I ) RES = SH0DAT ( 3 , I ) PT = NUM / DEN COMP = STRVH0(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 320 CONTINUE C C TEST STRVH1 C FNAME = 'STRVH1' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 330 I = 1 , 20 NUM = SH1DAT ( 1 , I ) DEN = SH1DAT ( 2 , I ) RES = SH1DAT ( 3 , I ) PT = NUM / DEN COMP = STRVH1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 330 CONTINUE C C TEST STRVL0 C FNAME = 'STRVL0' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 340 I = 1 , 20 NUM = SL0DAT ( 1 , I ) DEN = SL0DAT ( 2 , I ) RES = SL0DAT ( 3 , I ) PT = NUM / DEN COMP = STRVL0(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 340 CONTINUE C C TEST STRVL1 C FNAME = 'STRVL1' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 350 I = 1 , 20 NUM = SL1DAT ( 1 , I ) DEN = SL1DAT ( 2 , I ) RES = SL1DAT ( 3 , I ) PT = NUM / DEN COMP = STRVL1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 350 CONTINUE C C TEST SYNCH1 C FNAME = 'SYNCH1' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 360 I = 1 , 20 NUM = SY1DAT ( 1 , I ) DEN = SY1DAT ( 2 , I ) RES = SY1DAT ( 3 , I ) PT = NUM / DEN COMP = SYNCH1(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 360 CONTINUE C C TEST SYNCH2 C FNAME = 'SYNCH2' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 370 I = 1 , 20 NUM = SY2DAT ( 1 , I ) DEN = SY2DAT ( 2 , I ) RES = SY2DAT ( 3 , I ) PT = NUM / DEN COMP = SYNCH2(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 370 CONTINUE C C TEST TRAN02 C FNAME = 'TRAN02' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 380 I = 1 , 20 NUM = TR2DAT ( 1 , I ) DEN = TR2DAT ( 2 , I ) RES = TR2DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN02(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 380 CONTINUE C C TEST TRAN03 C FNAME = 'TRAN03' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 390 I = 1 , 20 NUM = TR3DAT ( 1 , I ) DEN = TR3DAT ( 2 , I ) RES = TR3DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN03(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 390 CONTINUE C C TEST TRAN04 C FNAME = 'TRAN04' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 400 I = 1 , 20 NUM = TR4DAT ( 1 , I ) DEN = TR4DAT ( 2 , I ) RES = TR4DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN04(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 400 CONTINUE C C TEST TRAN05 C FNAME = 'TRAN05' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 410 I = 1 , 20 NUM = TR5DAT ( 1 , I ) DEN = TR5DAT ( 2 , I ) RES = TR5DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN05(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 410 CONTINUE C C TEST TRAN06 C FNAME = 'TRAN06' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 420 I = 1 , 20 NUM = TR6DAT ( 1 , I ) DEN = TR6DAT ( 2 , I ) RES = TR6DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN06(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 420 CONTINUE C C TEST TRAN07 C FNAME = 'TRAN07' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 430 I = 1 , 20 NUM = TR7DAT ( 1 , I ) DEN = TR7DAT ( 2 , I ) RES = TR7DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN07(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 430 CONTINUE C C TEST TRAN08 C FNAME = 'TRAN08' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 440 I = 1 , 20 NUM = TR8DAT ( 1 , I ) DEN = TR8DAT ( 2 , I ) RES = TR8DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN08(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 440 CONTINUE C C TEST TRAN09 C FNAME = 'TRAN09' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 450 I = 1 , 20 NUM = TR9DAT ( 1 , I ) DEN = TR9DAT ( 2 , I ) RES = TR9DAT ( 3 , I ) PT = NUM / DEN COMP = TRAN09(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 450 CONTINUE C C TEST Y0INT C FNAME = 'Y0INT ' WRITE(IOUT,1000)FNAME WRITE(IOUT,1010) DO 460 I = 1 , 20 NUM = Y0IDAT ( 1 , I ) DEN = Y0IDAT ( 2 , I ) RES = Y0IDAT ( 3 , I ) PT = NUM / DEN COMP = Y0INT(PT) ABSERR = ABS ( RES - COMP ) RELERR = ABSERR / ABS(RES) WRITE(IOUT,1020)PT,ABSERR,RELERR 460 CONTINUE C C PRINT STATEMENTS C 1000 FORMAT(////15X,'TESTING THE FUNCTION ',A6) 1010 FORMAT(/5X,'ARGUMENT',13X,'ABS. ERROR',13X,'REL. ERROR') 1020 FORMAT(3X,F10.4,8X,D15.5,8X,D15.5) END