30 May 2018 11:01:58.475 AM TOMS577_TEST FORTRAN90 version TOMS577 evaluates Carlson's elliptic functions. RC_TEST RC evaluates the elementary integral RC(X,Y) X Y RC(X,Y) 0.1510000000000000D-77 0.1510000000000000D-77 0.8137884587711593D+39 0.3010000000000000D-77 0.5500000000000000D-78 0.9524521540989471D+39 0.0000000000000000D+00 0.3010000000000000D-77 0.9053919509295885D+39 0.9900000000000000D+75 0.5500000000000000D-78 0.5629713505821428D-35 0.0000000000000000D+00 0.9900000000000000D+75 0.4992318443077758D-37 0.9900000000000000D+75 0.9900000000000000D+75 0.3178208630818641D-37 RC - Error! Invalid input arguments. X = 0.0000000000000000D+00 Y = 0.2000000000000000D-77 0.0000000000000000D+00 0.2000000000000000D-77 ***Error*** 0.0000000000000000D+00 0.2500000000000000D+00 0.3141592653589793D+01 0.6250000000000000D-01 0.1250000000000000D+00 0.3141592653589793D+01 0.2250000000000000D+01 0.2000000000000000D+01 0.6931471805599454D+00 0.1000000000000000D-01 0.1000000000000000D+01 0.1478037662374775D+01 0.2000000000000000D-01 0.1000000000000000D+01 0.1443406235712271D+01 0.5000000000000000D-01 0.1000000000000000D+01 0.1380231154269966D+01 0.1000000000000000D+00 0.1000000000000000D+01 0.1316609847527586D+01 0.2000000000000000D+00 0.1000000000000000D+01 0.1237829897094659D+01 0.4000000000000000D+00 0.1000000000000000D+01 0.1143920647968054D+01 0.6000000000000000D+00 0.1000000000000000D+01 0.1082636119571208D+01 0.8000000000000000D+00 0.1000000000000000D+01 0.1036747571331046D+01 0.1000000000000000D+01 0.1000000000000000D+01 0.1000000000000000D+01 0.1200000000000000D+01 0.1000000000000000D+01 0.9693519329631458D+00 0.1500000000000000D+01 0.1000000000000000D+01 0.9312298594527122D+00 0.2000000000000000D+01 0.1000000000000000D+01 0.8813735870195432D+00 0.3000000000000000D+01 0.1000000000000000D+01 0.8104969894767539D+00 0.4000000000000000D+01 0.1000000000000000D+01 0.7603459963009465D+00 0.5000000000000000D+01 0.1000000000000000D+01 0.7218177375894052D+00 0.1000000000000000D+02 0.1000000000000000D+01 0.6061488197440222D+00 0.2000000000000000D+02 0.1000000000000000D+01 0.4997299176960261D+00 0.5000000000000000D+02 0.1000000000000000D+01 0.3777315372940898D+00 0.1000000000000000D+03 0.1000000000000000D+01 0.3008302149854818D+00 0.1000000000000000D+04 0.1000000000000000D+01 0.1311981462667840D+00 0.1000000000000000D+05 0.1000000000000000D+01 0.5298557300099017D-01 0.1000000000000000D+06 0.1000000000000000D+01 0.2039555141711115D-01 0.1000000000000000D+07 0.1000000000000000D+01 0.7600906009995943D-02 0.1000000000000000D+08 0.1000000000000000D+01 0.2767687204982141D-02 0.1000000000000000D+09 0.1000000000000000D+01 0.9903487599553565D-03 0.1000000000000000D+10 0.1000000000000000D+01 0.3495828416193509D-03 0.1000000000000000D+11 0.1000000000000000D+01 0.1220607264611548D-03 0.1000000000000000D+13 0.1000000000000000D+01 0.1450865773853122D-04 0.1000000000000000D+16 0.1000000000000000D+01 0.5680252434588194D-06 0.1000000000000000D+21 0.1000000000000000D+01 0.2371899811050040D-08 0.1000000000000000D+31 0.1000000000000000D+01 0.3523192357547063D-13 0.1000000000000000D+41 0.1000000000000000D+01 0.4674484904044085D-18 0.1000000000000000D+51 0.1000000000000000D+01 0.5825777450541109D-23 RC_TEST2 Compare LOG(X)/(X-1) and ARCTAN(X) with values based on RC. X From LOG From RC 0.2D+00 0.2011797390542625D+01 0.2011797390542625D+01 0.4D+00 0.1527151219790258D+01 0.1527151219790258D+01 0.6D+00 0.1277064059414976D+01 0.1277064059414977D+01 0.8D+00 0.1115717756571049D+01 0.1115717756571049D+01 0.1D+01 **** ZERO DIVIDE ***** 0.1000000000000000D+01 0.1D+01 0.9116077839697732D+00 0.9116077839697730D+00 0.1D+01 0.8411805915530323D+00 0.8411805915530324D+00 0.2D+01 0.7833393820762259D+00 0.7833393820762258D+00 0.2D+01 0.7347333311276488D+00 0.7347333311276486D+00 0.2D+01 0.6931471805599453D+00 0.6931471805599454D+00 Extreme values of X X From LOG From RC 0.1D-74 0.1726938819745534D+03 0.1726938819745534D+03 0.1D-64 0.1496680310446130D+03 0.1496680310446130D+03 0.1D-54 0.1266421801146725D+03 0.1266421801146725D+03 0.1D-44 0.1036163291847321D+03 0.1036163291847320D+03 0.1D-34 0.8059047825479161D+02 0.8059047825479161D+02 0.1D-24 0.5756462732485114D+02 0.5756462732485113D+02 0.1D-14 0.3453877639491072D+02 0.3453877639491072D+02 0.1D-04 0.1151304059537618D+02 0.1151304059537618D+02 0.1D+06 0.1151304059537618D-03 0.1151304059537618D-03 0.1D+16 0.3453877639491072D-13 0.3453877639491072D-13 0.1D+26 0.5756462732485114D-23 0.5756462732485114D-23 0.1D+36 0.8059047825479162D-33 0.8059047825479158D-33 0.1D+46 0.1036163291847320D-42 0.1036163291847321D-42 0.1D+56 0.1266421801146725D-52 0.1266421801146725D-52 0.1D+66 0.1496680310446129D-62 0.1496680310446130D-62 0.1D+76 0.1726938819745534D-72 0.1726938819745534D-72 X From ARCTAN From RC RC - Error! Invalid input arguments. X = 0.1000000000000000D+76 Y = 0.1000000000000000D+76 0.1D-74 0.1000000000000000D-74 0.0000000000000000D+00 0.1D-14 0.1000000000000000D-14 0.9999999999999999D-15 0.1D-02 0.9999996666668668D-03 0.9999996666668666D-03 0.1D+00 0.9966865249116204D-01 0.9966865249116202D-01 0.2D+00 0.1973955598498808D+00 0.1973955598498808D+00 0.5D+00 0.4636476090008061D+00 0.4636476090008061D+00 0.1D+01 0.7853981633974483D+00 0.7853981633974482D+00 0.2D+01 0.1107148717794090D+01 0.1107148717794091D+01 0.5D+01 0.1373400766945016D+01 0.1373400766945016D+01 0.1D+02 0.1471127674303735D+01 0.1471127674303734D+01 0.1D+04 0.1569796327128230D+01 0.1569796327128230D+01 0.1D+16 0.1570796326794896D+01 0.1570796326794896D+01 0.1D+76 0.1570796326794897D+01 0.1570796326794897D+01 RD_TEST RD evaluates the Carlson elliptic integral of the second kind, RD(X,Y,Z) X Y Z RD(X,Y,Z) 0.0000000000000000D+00 0.6010000000000001D-50 0.6010000000000001D-50 0.5057075390167092D+76 0.5500000000000000D-78 0.6010000000000001D-50 0.6010000000000001D-50 0.5057075390167032D+76 0.0000000000000000D+00 0.6010000000000001D-50 0.9900000000000000D+48 0.3455579458524223D-69 0.5500000000000000D-78 0.6010000000000001D-50 0.9900000000000000D+48 0.3455579458524223D-69 0.0000000000000000D+00 0.9900000000000000D+48 0.6010000000000001D-50 0.5016827696801392D+27 0.5500000000000000D-78 0.9900000000000000D+48 0.6010000000000001D-50 0.5016827696801343D+27 0.0000000000000000D+00 0.9900000000000000D+48 0.9900000000000000D+48 0.2391984406816877D-71 0.5500000000000000D-78 0.9900000000000000D+48 0.9900000000000000D+48 0.2391984406816877D-71 0.3010000000000000D-50 0.3010000000000000D-50 0.6010000000000001D-50 0.3170935022408316D+76 0.3010000000000000D-50 0.3010000000000000D-50 0.9900000000000000D+48 0.3444999003763636D-69 0.9900000000000000D+48 0.9900000000000000D+48 0.6010000000000001D-50 0.3908848372573853D-22 0.9900000000000000D+48 0.9900000000000000D+48 0.9900000000000000D+48 0.1015189712383042D-71 RD - Error! Invalid input arguments. X = 0.0000000000000000D+00 Y = 0.3010000000000000D-50 Z = 0.1000000000000000D+01 0.0000000000000000D+00 0.3010000000000000D-50 0.1000000000000000D+01 ***Error*** 0.0000000000000000D+00 0.2000000000000000D+01 0.1000000000000000D+01 0.1797210352103389D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D-09 0.4242559649759194D+06 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D-04 0.1333576985789446D+04 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D-01 0.3543137955816859D+02 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+00 0.8049184457709792D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+00 0.4795496764410121D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+00 0.2260648559158617D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.1210841860059132D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.6176739675072550D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.2361769958541499D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.1086719981146340D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.6540379849123936D-02 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+06 0.5320236560151541D-06 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+11 0.3409325949193374D-13 RF_TEST RF evaluates the Carlson elliptic integral of the first kind, RF(X,Y,Z) X Y Z RF(X,Y,Z) 0.1510000000000000D-77 0.1510000000000000D-77 0.1510000000000000D-77 0.8137884587711593D+39 0.1510000000000000D-77 0.1510000000000000D-77 0.9900000000000000D+75 0.5613664399997798D-35 0.0000000000000000D+00 0.3010000000000000D-77 0.3010000000000000D-77 0.9053919509295885D+39 0.0000000000000000D+00 0.3010000000000000D-77 0.9900000000000000D+75 0.5624731938415326D-35 0.0000000000000000D+00 0.9900000000000000D+75 0.9900000000000000D+75 0.4992318443077759D-37 0.9900000000000000D+75 0.9900000000000000D+75 0.9900000000000000D+75 0.3178208630818641D-37 0.5500000000000000D-78 0.3010000000000000D-77 0.3010000000000000D-77 0.7198956468812792D+39 0.5500000000000000D-78 0.3010000000000000D-77 0.9900000000000000D+75 0.5613420743331975D-35 0.5500000000000000D-78 0.9900000000000000D+75 0.9900000000000000D+75 0.4992318443077759D-37 RF - Error! Invalid input arguments. X = 0.0000000000000000D+00 Y = 0.2000000000000000D-77 Z = 0.1000000000000000D+01 0.0000000000000000D+00 0.2000000000000000D-77 0.1000000000000000D+01 ***Error*** 0.0000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.1311028777146060D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.1110720734539592D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1100000000000000D+01 0.1091430716268119D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1200000000000000D+01 0.1073732561377385D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1300000000000000D+01 0.1057389655801202D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1400000000000000D+01 0.1042214695930046D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1500000000000000D+01 0.1028056801052127D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1600000000000000D+01 0.1014792597551546D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1700000000000000D+01 0.1002319885847406D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1800000000000000D+01 0.9905530382414446D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1900000000000000D+01 0.9794195874708435D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.9688576532724523D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2200000000000000D+01 0.9492423586314841D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2400000000000000D+01 0.9313590566024380D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2600000000000000D+01 0.9149398964447096D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2800000000000000D+01 0.8997743463505848D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.3000000000000000D+01 0.8856939232035483D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.3500000000000000D+01 0.8543997654075564D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.4000000000000000D+01 0.8275182379903345D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.4500000000000000D+01 0.8040168167992912D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.7831840398091505D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.6000000000000000D+01 0.7476135036481746D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.7000000000000000D+01 0.7180694683312081D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.8000000000000000D+01 0.6929108826625124D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.9000000000000000D+01 0.6710782484938838D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.6518490214136576D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+02 0.5331872892861481D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.3000000000000000D+02 0.4706527717638377D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.4000000000000000D+02 0.4295233682160168D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+02 0.3994861473220863D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.3163981445460594D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+03 0.2479287246742144D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+03 0.1771580322263678D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+04 0.1361926299714245D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+05 0.5456849363690315D-01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+06 0.2089626719954944D-01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+07 0.7759252354859355D-02 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+09 0.1006183477210615D-02 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+11 0.1236441982979439D-03 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+13 0.1466700492234990D-04 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+16 0.5730326210782770D-06 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+21 0.2387734529432078D-08 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+31 0.3539027075929101D-13 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+41 0.4690319622426123D-18 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+51 0.5841612168923146D-23 RJ_TEST RJ evaluates the Carlson elliptic integral of the third kind, RJ(X,Y,Z,P) X Y Z P RJ(X,Y,Z,P) 0.1010000000000000D-25 0.1010000000000000D-25 0.1010000000000000D-25 0.2010000000000000D-25 0.6363797977910284D+39 0.1010000000000000D-25 0.1010000000000000D-25 0.2990000000000000D+25 0.2010000000000000D-25 0.5969809640810398D+14 0.0000000000000000D+00 0.2010000000000000D-25 0.2010000000000000D-25 0.2010000000000000D-25 0.8268315690384853D+39 0.0000000000000000D+00 0.2010000000000000D-25 0.2990000000000000D+25 0.2010000000000000D-25 0.8631566148749853D+14 0.0000000000000000D+00 0.2990000000000000D+25 0.2990000000000000D+25 0.2010000000000000D-25 0.1111659835295890D-10 0.2990000000000000D+25 0.2990000000000000D+25 0.2990000000000000D+25 0.2010000000000000D-25 0.3333899154957826D-34 0.5500000000000000D-78 0.2010000000000000D-25 0.2010000000000000D-25 0.2010000000000000D-25 0.8268315690384853D+39 0.5500000000000000D-78 0.2010000000000000D-25 0.2990000000000000D+25 0.2010000000000000D-25 0.8631566148749853D+14 0.5500000000000000D-78 0.2990000000000000D+25 0.2990000000000000D+25 0.2010000000000000D-25 0.1111659835295890D-10 0.2010000000000000D-25 0.2010000000000000D-25 0.2010000000000000D-25 0.2010000000000000D-25 0.3509182168450738D+39 0.1010000000000000D-25 0.1010000000000000D-25 0.1010000000000000D-25 0.2990000000000000D+25 0.9983650737892871D-11 0.1010000000000000D-25 0.1010000000000000D-25 0.2990000000000000D+25 0.2990000000000000D+25 0.3353865073488965D-34 0.0000000000000000D+00 0.2010000000000000D-25 0.2010000000000000D-25 0.2990000000000000D+25 0.1111659835295890D-10 0.0000000000000000D+00 0.2010000000000000D-25 0.2990000000000000D+25 0.2990000000000000D+25 0.3374118957422856D-34 0.0000000000000000D+00 0.2990000000000000D+25 0.2990000000000000D+25 0.2990000000000000D+25 0.4557265740116441D-36 0.2990000000000000D+25 0.2990000000000000D+25 0.2990000000000000D+25 0.2990000000000000D+25 0.1934163652061004D-36 0.5500000000000000D-78 0.2010000000000000D-25 0.2010000000000000D-25 0.2990000000000000D+25 0.1111659835295890D-10 0.5500000000000000D-78 0.2010000000000000D-25 0.2990000000000000D+25 0.2990000000000000D+25 0.3374118957422856D-34 0.5500000000000000D-78 0.2990000000000000D+25 0.2990000000000000D+25 0.2990000000000000D+25 0.4557265740116441D-36 0.2010000000000000D-25 0.2010000000000000D-25 0.2010000000000000D-25 0.2990000000000000D+25 0.7077046312962523D-11 RJ - Error! Invalid input arguments. X = 0.0000000000000000D+00 Y = 0.1900000000000000D-25 Z = 0.1900000000000000D-25 P = 0.1000000000000000D+01 0.0000000000000000D+00 0.1900000000000000D-25 0.1900000000000000D-25 0.1000000000000000D+01 ***Error*** 0.2000000000000000D+01 0.3000000000000000D+01 0.4000000000000000D+01 0.5000000000000000D+01 0.1429757966715675D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.2500000000000000D+00 0.2608105757997978D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.7500000000000000D+00 0.1442864193614197D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.1000000000000000D+01 0.1210841860059132D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+01 0.2000000000000000D+01 0.7670625432950464D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.2500000000000000D+00 0.2008380291058661D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.7500000000000000D+00 0.1136845217951276D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.1500000000000000D+01 0.7460643441600711D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.2000000000000000D+01 0.4000000000000000D+01 0.3794195060366037D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.2500000000000000D+00 0.1370892718222587D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.7500000000000000D+00 0.7949204422899481D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.3000000000000000D+01 0.3385294211417840D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.5000000000000000D+01 0.1000000000000000D+02 0.1398476510985250D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.2500000000000000D+00 0.1005572682199914D+01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.7500000000000000D+00 0.5908850286211613D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.5000000000000000D+01 0.1811791604390096D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+02 0.2000000000000000D+02 0.6283564851298871D-01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.2500000000000000D+00 0.3343272158567165D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.7500000000000000D+00 0.2006464395741285D+00 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.5000000000000000D+02 0.1167793209710091D-01 0.5000000000000000D+00 0.1000000000000000D+01 0.1000000000000000D+03 0.2000000000000000D+03 0.3581148146400939D-02 TOMS577_TEST Normal end of execution. 30 May 2018 11:01:58.476 AM