Tue Oct 19 11:33:55 2021 elliptic_integral_test() Python version: 3.6.9 Test elliptic_integral(). rc_test rc evaluates the elementary integral RC(X,Y) X Y RC(X,Y) 1.51e-78 1.51e-78 8.137884587711593e+38 3.01e-78 5.5e-79 9.524521540989471e+38 0 3.01e-78 9.053919509295885e+38 9.9e+74 5.5e-79 5.629713505821428e-36 0 9.9e+74 4.992318443077758e-38 9.9e+74 9.9e+74 3.178208630818641e-38 RC - Error! Invalid input arguments. X = 0 Y = 2e-78 0 2e-78 **Error*** 0 0.25 3.141592653589793 0.0625 0.125 3.141592653589793 2.25 2 0.6931471805599454 0.01 1 1.478037662374775 0.02 1 1.443406235712271 0.05 1 1.380231154269966 0.1 1 1.316609847527586 0.2 1 1.237829897094659 0.4 1 1.143920647968054 0.6 1 1.082636119571208 0.8 1 1.036747571331046 1 1 1 1.2 1 0.9693519329631458 1.5 1 0.9312298594527122 2 1 0.8813735870195432 3 1 0.8104969894767539 4 1 0.7603459963009465 5 1 0.7218177375894052 10 1 0.6061488197440222 20 1 0.4997299176960261 50 1 0.3777315372940898 100 1 0.3008302149854818 1000 1 0.131198146266784 10000 1 0.05298557300099017 100000 1 0.02039555141711115 1000000 1 0.007600906009995943 10000000 1 0.002767687204982141 100000000 1 0.0009903487599553565 1000000000 1 0.0003495828416193509 10000000000 1 0.0001220607264611548 1000000000000 1 1.450865773853122e-05 1000000000000000 1 5.680252434588194e-07 1e+20 1 2.37189981105004e-09 1e+30 1 3.523192357547063e-14 1e+40 1 4.674484904044085e-19 1e+50 1 5.825777450541109e-24 rc_test2 Compare LOG(X)/(X-1) and ARCTAN(X) with values based on rc(). X From LOG From RC 0.2 2.011797390542625 2.011797390542625 0.4 1.527151219790258 1.527151219790258 0.6 1.277064059414976 1.277064059414977 0.8 1.115717756571049 1.115717756571049 1 ***** ZERO DIVIDE ***** 1 1 0.9116077839697732 0.911607783969773 1 0.8411805915530323 0.8411805915530324 2 0.7833393820762259 0.7833393820762258 2 0.7347333311276488 0.7347333311276486 2 0.6931471805599453 0.6931471805599454 Extreme values of X X From LOG From RC 1e-75 172.6938819745534 172.6938819745534 1e-65 149.668031044613 149.668031044613 1e-55 126.6421801146725 126.6421801146725 1e-45 103.6163291847321 103.616329184732 1e-35 80.59047825479161 80.59047825479161 1e-25 57.56462732485114 57.56462732485113 1e-15 34.53877639491072 34.53877639491072 1e-05 11.51304059537618 11.51304059537618 1e+05 0.0001151304059537618 0.0001151304059537618 1e+15 3.453877639491072e-14 3.453877639491072e-14 1e+25 5.756462732485114e-24 5.756462732485114e-24 1e+35 8.059047825479162e-34 8.059047825479158e-34 1e+45 1.036163291847321e-43 1.036163291847321e-43 1e+55 1.266421801146725e-53 1.266421801146725e-53 1e+65 1.49668031044613e-63 1.49668031044613e-63 1e+75 1.726938819745534e-73 1.726938819745534e-73 X From ARCTAN From RC RC - Error! Invalid input arguments. X = 1e+75 Y = 1e+75 1e-75 1e-75 0 1e-15 1e-15 9.999999999999999e-16 0.001 0.0009999996666668668 0.0009999996666668666 0.1 0.09966865249116204 0.09966865249116202 0.2 0.1973955598498808 0.1973955598498808 0.5 0.4636476090008061 0.4636476090008061 1 0.7853981633974483 0.7853981633974482 2 1.10714871779409 1.107148717794091 5 1.373400766945016 1.373400766945016 1e+01 1.471127674303735 1.471127674303734 1e+03 1.56979632712823 1.56979632712823 1e+15 1.570796326794896 1.570796326794896 1e+75 1.570796326794897 1.570796326794897 rd_test rd() evaluates the Carlson elliptic integral of the second kind, RD(X,Y,Z) X Y Z RD(X,Y,Z) 0 6.010000000000001e-51 6.010000000000001e-51 5.057075390167092e+75 5.5e-79 6.010000000000001e-51 6.010000000000001e-51 5.057075390167032e+75 0 6.010000000000001e-51 9.9e+47 3.455579458524223e-70 5.5e-79 6.010000000000001e-51 9.9e+47 3.455579458524223e-70 0 9.9e+47 6.010000000000001e-51 5.016827696801392e+26 5.5e-79 9.9e+47 6.010000000000001e-51 5.016827696801343e+26 0 9.9e+47 9.9e+47 2.391984406816877e-72 5.5e-79 9.9e+47 9.9e+47 2.391984406816877e-72 3.01e-51 3.01e-51 6.010000000000001e-51 3.170935022408316e+75 3.01e-51 3.01e-51 9.9e+47 3.444999003763636e-70 9.9e+47 9.9e+47 6.010000000000001e-51 3.908848372573853e-23 9.9e+47 9.9e+47 9.9e+47 1.015189712383042e-72 rd() - Error! Invalid input arguments. X = 0 Y = 3.01e-51 Z = 1 0 3.01e-51 1 ***Error*** 0 2 1 1.797210352103389 0.5 1 1e-10 424255.9649759194 0.5 1 1e-05 1333.576985789446 0.5 1 0.01 35.43137955816859 0.5 1 0.1 8.049184457709792 0.5 1 0.2 4.795496764410121 0.5 1 0.5 2.260648559158617 0.5 1 1 1.210841860059132 0.5 1 2 0.617673967507255 0.5 1 5 0.2361769958541499 0.5 1 10 0.108671998114634 0.5 1 100 0.006540379849123936 0.5 1 100000 5.320236560151541e-07 0.5 1 10000000000 3.409325949193374e-14 rf_test rf() evaluates the Carlson elliptic integral of the first kind, RF(X,Y,Z) X Y Z RF(X,Y,Z) 1.51e-78 1.51e-78 1.51e-78 8.137884587711593e+38 1.51e-78 1.51e-78 9.9e+74 5.613664399997798e-36 0 3.01e-78 3.01e-78 9.053919509295885e+38 0 3.01e-78 9.9e+74 5.624731938415326e-36 0 9.9e+74 9.9e+74 4.992318443077759e-38 9.9e+74 9.9e+74 9.9e+74 3.178208630818641e-38 5.5e-79 3.01e-78 3.01e-78 7.198956468812792e+38 5.5e-79 3.01e-78 9.9e+74 5.613420743331975e-36 5.5e-79 9.9e+74 9.9e+74 4.992318443077759e-38 rd() - Error! Invalid input arguments. X = 0 Y = 2e-78 Z = 1 0 2e-78 1 ***Error*** 0 1 2 1.31102877714606 0.5 1 1 1.110720734539592 0.5 1 1.1 1.091430716268119 0.5 1 1.2 1.073732561377385 0.5 1 1.3 1.057389655801202 0.5 1 1.4 1.042214695930046 0.5 1 1.5 1.028056801052127 0.5 1 1.6 1.014792597551546 0.5 1 1.7 1.002319885847406 0.5 1 1.8 0.9905530382414446 0.5 1 1.9 0.9794195874708435 0.5 1 2 0.9688576532724523 0.5 1 2.2 0.9492423586314841 0.5 1 2.4 0.931359056602438 0.5 1 2.6 0.9149398964447096 0.5 1 2.8 0.8997743463505848 0.5 1 3 0.8856939232035483 0.5 1 3.5 0.8543997654075564 0.5 1 4 0.8275182379903345 0.5 1 4.5 0.8040168167992912 0.5 1 5 0.7831840398091505 0.5 1 6 0.7476135036481746 0.5 1 7 0.7180694683312081 0.5 1 8 0.6929108826625124 0.5 1 9 0.6710782484938838 0.5 1 10 0.6518490214136576 0.5 1 20 0.5331872892861481 0.5 1 30 0.4706527717638377 0.5 1 40 0.4295233682160168 0.5 1 50 0.3994861473220863 0.5 1 100 0.3163981445460594 0.5 1 200 0.2479287246742144 0.5 1 500 0.1771580322263678 0.5 1 1000 0.1361926299714245 0.5 1 10000 0.05456849363690315 0.5 1 100000 0.02089626719954944 0.5 1 1000000 0.007759252354859355 0.5 1 100000000 0.001006183477210615 0.5 1 10000000000 0.0001236441982979439 0.5 1 1000000000000 1.46670049223499e-05 0.5 1 1000000000000000 5.73032621078277e-07 0.5 1 1e+20 2.387734529432078e-09 0.5 1 1e+30 3.539027075929101e-14 0.5 1 1e+40 4.690319622426123e-19 0.5 1 1e+50 5.841612168923146e-24 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) 1.01e-26 1.01e-26 1.01e-26 2.01e-26 6.363797977910284e+38 1.01e-26 1.01e-26 2.99e+24 2.01e-26 59698096408103.98 0 2.01e-26 2.01e-26 2.01e-26 8.268315690384853e+38 0 2.01e-26 2.99e+24 2.01e-26 86315661487498.53 0 2.99e+24 2.99e+24 2.01e-26 1.11165983529589e-11 2.99e+24 2.99e+24 2.99e+24 2.01e-26 3.333899154957826e-35 5.5e-79 2.01e-26 2.01e-26 2.01e-26 8.268315690384853e+38 5.5e-79 2.01e-26 2.99e+24 2.01e-26 86315661487498.53 5.5e-79 2.99e+24 2.99e+24 2.01e-26 1.11165983529589e-11 2.01e-26 2.01e-26 2.01e-26 2.01e-26 3.509182168450738e+38 1.01e-26 1.01e-26 1.01e-26 2.99e+24 9.983650737892871e-12 1.01e-26 1.01e-26 2.99e+24 2.99e+24 3.353865073488965e-35 0 2.01e-26 2.01e-26 2.99e+24 1.11165983529589e-11 0 2.01e-26 2.99e+24 2.99e+24 3.374118957422856e-35 0 2.99e+24 2.99e+24 2.99e+24 4.557265740116441e-37 2.99e+24 2.99e+24 2.99e+24 2.99e+24 1.934163652061004e-37 5.5e-79 2.01e-26 2.01e-26 2.99e+24 1.11165983529589e-11 5.5e-79 2.01e-26 2.99e+24 2.99e+24 3.374118957422856e-35 5.5e-79 2.99e+24 2.99e+24 2.99e+24 4.557265740116441e-37 2.01e-26 2.01e-26 2.01e-26 2.99e+24 7.077046312962523e-12 rj() - Error! Invalid input arguments. X = 0 Y = 1.9e-26 Z = 1.9e-26 P = 1 0 1.9e-26 1.9e-26 1 ***Error*** 2 3 4 5 0.1429757966715675 0.5 1 1 0.25 2.608105757997978 0.5 1 1 0.75 1.442864193614197 0.5 1 1 1 1.210841860059132 0.5 1 1 2 0.7670625432950464 0.5 1 2 0.25 2.008380291058661 0.5 1 2 0.75 1.136845217951276 0.5 1 2 1.5 0.7460643441600711 0.5 1 2 4 0.3794195060366037 0.5 1 5 0.25 1.370892718222587 0.5 1 5 0.75 0.7949204422899481 0.5 1 5 3 0.338529421141784 0.5 1 5 10 0.139847651098525 0.5 1 10 0.25 1.005572682199914 0.5 1 10 0.75 0.5908850286211613 0.5 1 10 5 0.1811791604390096 0.5 1 10 20 0.06283564851298871 0.5 1 100 0.25 0.3343272158567165 0.5 1 100 0.75 0.2006464395741285 0.5 1 100 50 0.01167793209710091 0.5 1 100 200 0.003581148146400939 elliptic_fa_test: elliptic_fa returns values of the complete elliptic integral of the first kind, with parameter angle A. A F(A) F(A) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 5.000000 1.573792130924768 1.573792130924768 10.000000 1.582842804338351 1.582842804338351 15.000000 1.59814200211254 1.59814200211254 20.000000 1.620025899124204 1.620025899124204 25.000000 1.64899521847853 1.648995218478531 30.000000 1.685750354812596 1.685750354812596 35.000000 1.731245175657058 1.731245175657058 40.000000 1.786769134885021 1.786769134885021 45.000000 1.854074677301372 1.854074677301372 50.000000 1.935581096004722 1.935581096004722 55.000000 2.034715312185791 2.03471531218579 60.000000 2.156515647499643 2.156515647499643 65.000000 2.308786798167196 2.308786798167196 70.000000 2.504550079001634 2.504550079001633 75.000000 2.768063145368768 2.768063145368768 80.000000 3.153385251887839 3.153385251887837 85.000000 3.831741999784146 3.831741999784147 elliptic_fk_test: elliptic_fk returns values of the complete elliptic integral of the first kind, with parameter K. K F(K) F(K) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 0.223607 1.591003453790792 1.591003453790792 0.316228 1.612441348720219 1.612441348720219 0.387298 1.63525673226458 1.63525673226458 0.447214 1.659623598610528 1.659623598610528 0.500000 1.685750354812596 1.685750354812596 0.547723 1.713889448178791 1.713889448178791 0.591608 1.744350597225613 1.744350597225613 0.632456 1.777519371491253 1.777519371491253 0.670820 1.813883936816983 1.813883936816982 0.707107 1.854074677301372 1.854074677301372 0.741620 1.898924910271554 1.898924910271553 0.774597 1.949567749806026 1.949567749806026 0.806226 2.007598398424376 2.007598398424376 0.836660 2.075363135292469 2.07536313529247 0.866025 2.156515647499643 2.156515647499643 0.894427 2.257205326820854 2.257205326820854 0.921954 2.38901648632558 2.389016486325581 0.948683 2.578092113348173 2.578092113348173 0.974679 2.908337248444552 2.908337248444551 elliptic_fm_test: elliptic_fm returns values of the complete elliptic integral of the first kind, with parameter M. M F(M) F(M) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 0.050000 1.591003453790792 1.591003453790793 0.100000 1.612441348720219 1.61244134872022 0.150000 1.63525673226458 1.63525673226458 0.200000 1.659623598610528 1.659623598610528 0.250000 1.685750354812596 1.685750354812596 0.300000 1.713889448178791 1.713889448178791 0.350000 1.744350597225613 1.744350597225613 0.400000 1.777519371491253 1.777519371491253 0.450000 1.813883936816983 1.813883936816982 0.500000 1.854074677301372 1.854074677301372 0.550000 1.898924910271554 1.898924910271553 0.600000 1.949567749806026 1.949567749806026 0.650000 2.007598398424376 2.007598398424376 0.700000 2.075363135292469 2.075363135292469 0.750000 2.156515647499643 2.156515647499643 0.800000 2.257205326820854 2.257205326820854 0.850000 2.38901648632558 2.38901648632558 0.900000 2.578092113348173 2.578092113348173 0.950000 2.908337248444552 2.908337248444551 elliptic_ea_test: elliptic_ea returns values of the complete elliptic integral of the second kind, with parameter angle A. A E(A) E(A) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 5.000000 1.567809073977622 1.567809073977621 10.000000 1.558887196601596 1.558887196601596 15.000000 1.544150496914673 1.544150496914673 20.000000 1.523799205259774 1.523799205259774 25.000000 1.498114928422116 1.498114928422115 30.000000 1.467462209339427 1.467462209339427 35.000000 1.432290969306756 1.432290969306756 40.000000 1.393140248523812 1.393140248523812 45.000000 1.350643881047676 1.350643881047676 50.000000 1.305539094297794 1.305539094297794 55.000000 1.258679624779997 1.258679624779997 60.000000 1.211056027568459 1.21105602756846 65.000000 1.163827964493139 1.163827964493139 70.000000 1.118377737969864 1.118377737969863 75.000000 1.076405113076403 1.076405113076403 80.000000 1.04011439570601 1.04011439570601 85.000000 1.012663506234396 1.012663506234396 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 90.000000 1 0 elliptic_ek_test: elliptic_ek returns values of the complete elliptic integral of the second kind, with parameter K. K E(K) E(K) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 0.223607 1.550973351780472 1.550973351780472 0.316228 1.530757636897763 1.530757636897763 0.387298 1.510121832092819 1.510121832092819 0.447214 1.489035058095853 1.489035058095853 0.500000 1.467462209339427 1.467462209339427 0.547723 1.445363064412665 1.445363064412665 0.591608 1.422691133490879 1.422691133490879 0.632456 1.399392138897432 1.399392138897432 0.670820 1.375401971871116 1.375401971871116 0.707107 1.350643881047676 1.350643881047676 0.741620 1.32502449795823 1.32502449795823 0.774597 1.298428035046913 1.298428035046913 0.806226 1.270707479650149 1.27070747965015 0.836660 1.241670567945823 1.241670567945823 0.866025 1.211056027568459 1.21105602756846 0.894427 1.178489924327839 1.178489924327839 0.921954 1.143395791883166 1.143395791883166 0.948683 1.104774732704073 1.104774732704073 0.974679 1.060473727766278 1.060473727766278 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 1.000000 1 0 elliptic_em_test: elliptic_em returns values of the complete elliptic integral of the second kind, with parameter M. M E(M) E(M) Tabulated Calculated 0.000000 1.570796326794897 1.570796326794897 0.050000 1.550973351780472 1.550973351780473 0.100000 1.530757636897763 1.530757636897763 0.150000 1.510121832092819 1.510121832092819 0.200000 1.489035058095853 1.489035058095852 0.250000 1.467462209339427 1.467462209339427 0.300000 1.445363064412665 1.445363064412665 0.350000 1.422691133490879 1.422691133490879 0.400000 1.399392138897432 1.399392138897432 0.450000 1.375401971871116 1.375401971871116 0.500000 1.350643881047676 1.350643881047676 0.550000 1.32502449795823 1.32502449795823 0.600000 1.298428035046913 1.298428035046913 0.650000 1.270707479650149 1.270707479650149 0.700000 1.241670567945823 1.241670567945823 0.750000 1.211056027568459 1.21105602756846 0.800000 1.178489924327839 1.178489924327839 0.850000 1.143395791883166 1.143395791883166 0.900000 1.104774732704073 1.104774732704073 0.950000 1.060473727766278 1.060473727766278 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 rd() - Error! Invalid input arguments. X = 0 Y = 0 Z = 1 1.000000 1 0 elliptic_pia_test: elliptic_pia returns values of the complete elliptic integral of the third kind, with parameter angle A. N A Pi(N,A) Pi(N,A) Tabulated Calculated -10.000000 30.000000 0.4892245275965397 0.4892245275965397 -10.000000 45.000000 0.5106765677902629 0.5106765677902629 -10.000000 60.000000 0.5460409271920561 0.5460409271920565 -10.000000 77.079034 0.6237325893535237 0.6237325893535237 -3.000000 30.000000 0.823045542660675 0.8230455426606752 -3.000000 45.000000 0.8760028274011437 0.8760028274011439 -3.000000 60.000000 0.9660073560143946 0.9660073560143951 -3.000000 77.079034 1.171952391481798 1.171952391481798 -1.000000 30.000000 1.177446843000566 1.177446843000566 -1.000000 45.000000 1.273127366749682 1.273127366749682 -1.000000 60.000000 1.440034318657551 1.440034318657551 -1.000000 77.079034 1.836472172302591 1.83647217230259 0.000000 30.000000 1.685750354812596 1.685750354812596 0.000000 45.000000 1.854074677301372 1.854074677301372 0.000000 60.000000 2.156515647499643 2.156515647499643 0.000000 77.079034 2.908337248444552 2.908337248444551 0.500000 30.000000 2.413671504201195 2.413671504201195 0.500000 45.000000 2.701287762095351 2.701287762095351 0.500000 60.000000 3.234773471249465 3.234773471249465 0.500000 77.079034 4.633308147279891 4.633308147279889 elliptic_pik_test: elliptic_pik returns values of the complete elliptic integral of the third kind, with parameter K. N K Pi(N,K) Pi(N,K) Tabulated Calculated -10.000000 0.500000 0.4892245275965397 0.4892245275965397 -10.000000 0.707107 0.5106765677902629 0.5106765677902632 -10.000000 0.866025 0.5460409271920561 0.5460409271920565 -10.000000 0.974679 0.6237325893535237 0.6237325893535237 -3.000000 0.500000 0.823045542660675 0.8230455426606752 -3.000000 0.707107 0.8760028274011437 0.8760028274011439 -3.000000 0.866025 0.9660073560143946 0.9660073560143951 -3.000000 0.974679 1.171952391481798 1.171952391481798 -1.000000 0.500000 1.177446843000566 1.177446843000566 -1.000000 0.707107 1.273127366749682 1.273127366749683 -1.000000 0.866025 1.440034318657551 1.440034318657551 -1.000000 0.974679 1.836472172302591 1.83647217230259 0.000000 0.500000 1.685750354812596 1.685750354812596 0.000000 0.707107 1.854074677301372 1.854074677301372 0.000000 0.866025 2.156515647499643 2.156515647499643 0.000000 0.974679 2.908337248444552 2.908337248444551 0.500000 0.500000 2.413671504201195 2.413671504201195 0.500000 0.707107 2.701287762095351 2.701287762095351 0.500000 0.866025 3.234773471249465 3.234773471249465 0.500000 0.974679 4.633308147279891 4.633308147279889 elliptic_pim_test: elliptic_pim returns values of the complete elliptic integral of the third kind, with parameter M. N M Pi(N,M) Pi(N,M) Tabulated Calculated -10.000000 0.250000 0.4892245275965397 0.4892245275965397 -10.000000 0.500000 0.5106765677902629 0.5106765677902629 -10.000000 0.750000 0.5460409271920561 0.5460409271920565 -10.000000 0.950000 0.6237325893535237 0.6237325893535228 -3.000000 0.250000 0.823045542660675 0.8230455426606752 -3.000000 0.500000 0.8760028274011437 0.8760028274011439 -3.000000 0.750000 0.9660073560143946 0.9660073560143951 -3.000000 0.950000 1.171952391481798 1.171952391481798 -1.000000 0.250000 1.177446843000566 1.177446843000566 -1.000000 0.500000 1.273127366749682 1.273127366749682 -1.000000 0.750000 1.440034318657551 1.440034318657551 -1.000000 0.950000 1.836472172302591 1.83647217230259 0.000000 0.250000 1.685750354812596 1.685750354812596 0.000000 0.500000 1.854074677301372 1.854074677301372 0.000000 0.750000 2.156515647499643 2.156515647499643 0.000000 0.950000 2.908337248444552 2.908337248444551 0.500000 0.250000 2.413671504201195 2.413671504201195 0.500000 0.500000 2.701287762095351 2.701287762095351 0.500000 0.750000 3.234773471249465 3.234773471249465 0.500000 0.950000 4.633308147279891 4.633308147279889 elliptic_inc_fa_test: elliptic_inc_fa returns values of the incomplete elliptic integral of the first kind, with parameters PHI, A. Compare with tabulated value. Phi A Tabulated elliptic_inc_fa(phi,a) 0.343091 123.082123 0.3478806460316299 0.3478806460316299 1.302990 11.269317 1.3131805770095839 1.3131805770095843 0.652363 -94.888065 0.7037956689264326 0.7037956689264325 0.404602 -99.714079 0.4157626844675118 0.4157626844675119 0.068846 57.058810 0.0688847548328514 0.0688847548328514 0.096961 -19.713633 0.0969781675484583 0.0969781675484583 0.630370 56.312303 0.6605394722518515 0.6605394722518516 1.252375 -91.556053 1.8275834603675101 1.8275834603675094 1.409796 -27.006546 1.4822587833924870 1.4822587833924878 0.148511 -169.229373 0.1485295339221232 0.1485295339221233 1.349466 61.968596 1.7538000627014940 1.7538000627014940 0.193371 -158.732440 0.1935288964653510 0.1935288964653510 0.408883 105.088396 0.4199100508706138 0.4199100508706138 0.178543 -48.958839 0.1790836490491233 0.1790836490491232 1.292588 -42.585688 1.4460488322797631 1.4460488322797631 1.087096 11.656033 1.0940976521009840 1.0940976521009838 1.352795 -8.398114 1.3589479084270351 1.3589479084270346 1.432530 17.693622 1.4640007823153800 1.4640007823153802 0.296809 73.880342 0.3009092014525799 0.3009092014525800 0.623588 -69.824923 0.6621341112075102 0.6621341112075103 elliptic_inc_fk_test: elliptic_inc_fk returns values of the incomplete elliptic integral of the first kind, with parameters PHI, K. Compare with tabulated value. Phi K Tabulated elliptic_inc_fk(phi,k) 0.343091 2.712953 0.4340870330108736 0.4340870330108734 1.302990 0.127952 1.3073125113981141 1.3073125113981143 0.652363 -1.429438 0.8005154258533936 0.8005154258533933 0.404602 -1.981659 0.4656721451084328 0.4656721451084328 0.068846 3.894802 0.0696984961344177 0.0696984961344177 0.096961 -1.042486 0.0971264670875049 0.0971264670875049 0.630370 0.864114 0.6632598061016007 0.6632598061016008 1.252375 -1.049058 2.2308677858579000 2.2308677858579107 1.409796 -0.302406 1.4398462828880190 1.4398462828880187 0.148511 -6.574289 0.2043389243773096 0.2043389243773098 1.349466 0.698740 1.5371835748817710 1.5371835748817699 0.193371 -5.125586 0.2749229901565622 0.2749229901565621 0.408883 2.074948 0.4828388342828284 0.4828388342828284 0.178543 -1.670886 0.1812848567886627 0.1812848567886627 1.292588 -0.484360 1.3607295223418410 1.3607295223418410 1.087096 0.139306 1.0903968091202700 1.0903968091202703 1.352795 -0.094653 1.3553630515818080 1.3553630515818083 1.432530 0.197721 1.4454628197324411 1.4454628197324413 0.296809 1.788160 0.3125355489354676 0.3125355489354676 0.623588 -1.077781 0.6775731623807174 0.6775731623807173 elliptic_inc_fm_test: elliptic_inc_fm returns values of the incomplete elliptic integral of the first kind, with parameters PHI, M. Compare with tabulated value. Phi M Tabulated elliptic_inc_fm(phi,m) 0.343091 8.450690 0.4804314075855023 0.4804314075855022 1.302990 0.603988 1.5356349810920249 1.5356349810920260 0.652363 0.179413 0.6602285297476601 0.6602285297476599 0.404602 0.709569 0.4125884303785135 0.4125884303785134 0.068846 133.964339 0.0796456600715538 0.0796456600715537 0.096961 47.966214 0.1062834070535258 0.1062834070535258 0.630370 2.172071 0.7733990864393913 0.7733990864393913 1.252375 0.002038 1.2528624998922280 1.2528624998922278 1.409796 0.360004 1.5499886866115320 1.5499886866115324 0.148511 0.621954 0.1488506735822822 0.1488506735822822 1.349466 0.883422 1.8922299007996619 1.8922299007996617 0.193371 0.203429 0.1936153327753556 0.1936153327753556 0.408883 5.772526 0.5481932935424454 0.5481932935424454 0.178543 11.148539 0.1911795073571756 0.1911795073571756 1.292588 0.288924 1.3792250693497561 1.3792250693497552 1.087096 0.716662 1.2612824533314020 1.2612824533314022 1.352795 0.476062 1.5352398385253780 1.5352398385253774 1.432530 0.609495 1.7397824181560710 1.7397824181560722 0.296809 8.902277 0.3616930047198503 0.3616930047198503 0.623588 0.543444 0.6458627645916422 0.6458627645916422 elliptic_inc_ea_test: elliptic_inc_ea returns values of the incomplete elliptic integral of the second kind, with parameters PHI, A. Compare with tabulated value. Phi A Tabulated elliptic_inc_ea(phi,a) 0.343091 123.082123 0.3384181367348019 0.3384181367348018 1.302990 11.269317 1.2929246245095061 1.2929246245095070 0.652363 -94.888065 0.6074183768796306 0.6074183768796306 0.404602 -99.714079 0.3939726730783567 0.3939726730783568 0.068846 57.058810 0.0688081409708980 0.0688081409708980 0.096961 -19.713633 0.0969436473376824 0.0969436473376824 0.630370 56.312303 0.6025937791452033 0.6025937791452034 1.252375 -91.556053 0.9500549494837583 0.9500549494837583 1.409796 -27.006546 1.3427833721404860 1.3427833721404863 0.148511 -169.229373 0.1484915631401388 0.1484915631401388 1.349466 61.968596 1.0854328870509260 1.0854328870509264 0.193371 -158.732440 0.1932136916085597 0.1932136916085597 0.408883 105.088396 0.3983689593057807 0.3983689593057807 0.178543 -48.958839 0.1780054133336934 0.1780054133336934 1.292588 -42.585688 1.1645252702735360 1.1645252702735363 1.087096 11.656033 1.0801670475418450 1.0801670475418454 1.352795 -8.398114 1.3466849638303120 1.3466849638303116 1.432530 17.693622 1.4021002726855041 1.4021002726855041 0.296809 73.880342 0.2928091845544553 0.2928091845544553 0.623588 -69.824923 0.5889342583405707 0.5889342583405708 elliptic_inc_ek_test: elliptic_inc_ek returns values of the incomplete elliptic integral of the second kind, with parameters PHI, K. Compare with tabulated value. Phi K Tabulated elliptic_inc_ek(phi,k) 0.343091 2.712953 0.2852345328295404 0.2852345328295403 1.302990 0.127952 1.2986902255679209 1.2986902255679211 0.652363 -1.429438 0.5508100202571943 0.5508100202571941 0.404602 -1.981659 0.3575401358115371 0.3575401358115370 0.068846 3.894802 0.0680130780550745 0.0680130780550745 0.096961 -1.042486 0.0967958498023184 0.0967958498023184 0.630370 0.864114 0.6003112504412838 0.6003112504412839 1.252375 -1.049058 0.8996717721794724 0.8996717721794723 1.409796 -0.302406 1.3807152614538749 1.3807152614538749 0.148511 -6.574289 0.1191644625202453 0.1191644625202453 1.349466 0.698740 1.1969948381715569 1.1969948381715565 0.193371 -5.125586 0.1536260979667945 0.1536260979667945 0.408883 2.074948 0.3546768920544152 0.3546768920544152 0.178543 -1.670886 0.1758756066650882 0.1758756066650882 1.292588 -0.484360 1.2298191094105690 1.2298191094105686 1.087096 0.139306 1.0838106611433700 1.0838106611433700 1.352795 -0.094653 1.3502337815737799 1.3502337815737802 1.432530 0.197721 1.4197758847092179 1.4197758847092175 0.296809 1.788160 0.2824895528020034 0.2824895528020034 0.623588 -1.077781 0.5770427720982867 0.5770427720982867 elliptic_inc_em_test: elliptic_inc_em returns values of the incomplete elliptic integral of the second kind, with parameters PHI, M. Compare with tabulated value. Phi M Tabulated elliptic_inc_em(phi,m) 0.343091 8.450690 0.2732317284159052 0.2732317284159052 1.302990 0.603988 1.1247497250997811 1.1247497250997813 0.652363 0.179413 0.6446601913679151 0.6446601913679149 0.404602 0.709569 0.3968902354370061 0.3968902354370061 0.068846 133.964339 0.0606396079994467 0.0606396079994467 0.096961 47.966214 0.0890941157794873 0.0890941157794873 0.630370 2.172071 0.5324020148020150 0.5324020148020149 1.252375 0.002038 1.2518886406602650 1.2518886406602647 1.409796 0.360004 1.2889711619162600 1.2889711619162605 0.148511 0.621954 0.1481718153599732 0.1481718153599732 1.349466 0.883422 1.0380901856399130 1.0380901856399132 0.193371 0.203429 0.1931275771541276 0.1931275771541275 0.408883 5.772526 0.3304419611986801 0.3304419611986802 0.178543 11.148539 0.1673947960639630 0.1673947960639630 1.292588 0.288924 1.2145011753247359 1.2145011753247354 1.087096 0.716662 0.9516560179840655 0.9516560179840655 1.352795 0.476062 1.2036829595261760 1.2036829595261755 1.432530 0.609495 1.2064263261854189 1.2064263261854191 0.296809 8.902277 0.2522791382096692 0.2522791382096691 0.623588 0.543444 0.6026499038720986 0.6026499038720985 elliptic_inc_pia_test: elliptic_inc_pia returns values of the incomplete elliptic integral of the third kind, with parameters PHI, N, A. Compare with tabulated value. Phi N A Tabulated elliptic_inc_pia(phi,n,a) 0.343091 8.064681 88.878225 0.7099335174334724 0.7099335174334722 0.882309 -0.284059 -86.552087 0.9601963779142505 0.9601963779142507 0.404602 -5.034023 -116.619570 0.3362852532098376 0.3362852532098377 0.995831 -1.244606 -9.742878 0.7785343427543768 0.7785343427543769 0.630370 1.465982 65.734809 0.8578897552144780 0.8578897552144780 0.002888 95338.128573 -115.038772 0.0046307723449318 0.0046307723449318 0.148511 -44.431306 124.942118 0.1173842687902911 0.1173842687902912 1.320800 -0.802937 -89.787044 1.5057880706602671 1.5057880706602669 0.408883 5.218883 -98.426738 0.7213264194624553 0.7213264194624552 0.552337 2.345822 -53.749362 0.8073261799642218 0.8073261799642221 1.087096 0.157358 68.280476 1.4028538111108380 1.4028538111108377 0.712818 1.926593 20.821747 1.2592453314745129 1.2592453314745127 0.296809 6.113983 -29.104236 0.3779079263971614 0.3779079263971614 0.291091 1.805711 -37.801767 0.3088493910496766 0.3088493910496766 0.969503 -0.407285 -55.811734 0.9782829177005183 0.9782829177005182 1.122289 -0.941640 -37.665946 0.9430491574504173 0.9430491574504173 1.295912 0.700966 -80.094082 3.3207962773841548 3.3207962773841562 1.116491 -1.019831 52.238065 0.9730988737054799 0.9730988737054804 1.170719 -0.451080 74.309452 1.3019880949537890 1.3019880949537890 1.199361 0.602882 -17.229207 1.6455836044525900 1.6455836044525907 elliptic_inc_pik_test: elliptic_inc_pik returns values of the incomplete elliptic integral of the third kind, with parameters PHI, N, K. Compare with tabulated value. Phi N K Tabulated elliptic_inc_pik(phi,n,k) 0.343091 8.064681 1.959037 0.7982975462595892 0.7982975462595890 0.882309 -0.284059 -1.123742 1.0240221347260361 1.0240221347260359 0.404602 -5.034023 -2.317629 0.4015812085264200 0.4015812085264200 0.995831 -1.244606 -0.120258 0.7772649487439858 0.7772649487439859 0.630370 1.465982 1.008703 0.8737159913132074 0.8737159913132075 0.002888 95338.128573 -103.367749 0.0047333342976913 0.0047333342976913 0.148511 -44.431306 4.853800 0.1280656893638068 0.1280656893638067 1.320800 -0.802937 -1.016577 1.5943760375125640 1.5943760375125646 0.408883 5.218883 -1.943415 0.8521145133671923 0.8521145133671921 0.552337 2.345822 -0.887659 0.8154325229803082 0.8154325229803081 1.087096 0.157358 0.816049 1.3159451407542699 1.3159451407542697 0.712818 1.926593 0.299455 1.2539462314842400 1.2539462314842400 0.296809 6.113983 -0.704423 0.3796503567258643 0.3796503567258642 0.291091 1.805711 -0.926652 0.3111034454739552 0.3111034454739551 0.969503 -0.407285 -0.696261 0.9442477901112342 0.9442477901112344 1.122289 -0.941640 -0.445393 0.9153111661980959 0.9153111661980956 1.295912 0.700966 -0.910458 2.8420806443283930 2.8420806443283930 1.116491 -1.019831 0.618750 0.9263253777034376 0.9263253777034379 1.170719 -0.451080 0.867231 1.2123960187576239 1.2123960187576241 1.199361 0.602882 -0.199677 1.6280835727104710 1.6280835727104708 elliptic_inc_pim_test: elliptic_inc_pim returns values of the incomplete elliptic integral of the third kind, with parameters PHI, N, M. Compare with tabulated value. Phi N M Tabulated elliptic_inc_pim(phi,n,m) 0.343091 8.064681 7.330123 1.0469349800784999 1.0469349800784995 0.882309 -0.284059 0.110881 0.8421144481406690 0.8421144481406688 0.404602 -5.034023 0.282836 0.3321642201520043 0.3321642201520044 0.995831 -1.244606 0.638300 0.8483033529960849 0.8483033529960846 0.630370 1.465982 2.294719 1.0557538176567720 1.0557538176567718 0.002888 95338.128573 42062.553298 0.0051088961442656 0.0051088961442656 0.148511 -44.431306 39.239434 0.1426848042785896 0.1426848042785896 1.320800 -0.802937 0.008002 1.0313509582064240 1.0313509582064242 0.408883 5.218883 0.719058 0.7131013701418496 0.7131013701418496 0.552337 2.345822 0.970377 0.8268044665355507 0.8268044665355507 1.087096 0.157358 1.098881 1.5763286789601501 1.5763286789601489 0.712818 1.926593 1.398067 1.5428171208572110 1.5428171208572108 0.296809 6.113983 4.641022 0.4144629799126912 0.4144629799126910 0.291091 1.805711 4.455969 0.3313231611366746 0.3313231611366746 0.969503 -0.407285 0.313145 0.9195822851915201 0.9195822851915201 1.122289 -0.941640 0.368644 0.9422320754002217 0.9422320754002216 1.295912 0.700966 0.066782 2.0365990028158589 2.0365990028158603 1.116491 -1.019831 0.963554 1.0767992314998820 1.0767992314998829 1.170719 -0.451080 1.060209 1.4160844629578520 1.4160844629578537 1.199361 0.602882 0.468716 1.8241249223108911 1.8241249223108922 jacobi_cn_test: jacobi_cn() evaluates the Jacobi elliptic function CN. U M Exact CN CN(U,M) 0.1000 0.0000 0.9950041652780258 0.9950041652780258 0.2000 0.0000 0.9800665778412416 0.9800665778412416 0.5000 0.0000 0.8775825618903726 0.8775825618903728 1.0000 0.0000 0.5403023058681397 0.5403023058681398 2.0000 0.0000 -0.4161468365471424 -0.4161468365471424 0.1000 0.5000 0.9950124626090582 0.9950124626090583 0.2000 0.5000 0.9801976276784098 0.9801976276784097 0.5000 0.5000 0.8822663948904403 0.8822663948904402 1.0000 0.5000 0.5959765676721407 0.5959765676721407 2.0000 0.5000 -0.1031836155277618 -0.103183615527762 0.1000 1.0000 0.9950207489532265 0.9950207489532266 0.2000 1.0000 0.9803279976447253 0.9803279976447253 0.5000 1.0000 0.8868188839700739 0.886818883970074 1.0000 1.0000 0.6480542736638853 0.6480542736638855 2.0000 1.0000 0.2658022288340797 0.2658022288340797 4.0000 1.0000 0.03661899347368653 0.03661899347368653 -0.2000 1.0000 0.9803279976447253 0.9803279976447253 -0.5000 1.0000 0.8868188839700739 0.886818883970074 -1.0000 1.0000 0.6480542736638853 0.6480542736638855 -2.0000 1.0000 0.2658022288340797 0.2658022288340797 jacobi_dn_test: jacobi_dn() evaluates the Jacobi elliptic function DN. U M Exact DN DN(U,M) 0.1000 0.0000 1 1 0.2000 0.0000 1 1 0.5000 0.0000 1 1 1.0000 0.0000 1 1 2.0000 0.0000 1 1 0.1000 0.5000 0.9975093485144243 0.9975093485144244 0.2000 0.5000 0.9901483195224799 0.9901483195224799 0.5000 0.5000 0.9429724257773857 0.9429724257773857 1.0000 0.5000 0.8231610016315963 0.8231610016315963 2.0000 0.5000 0.7108610477840873 0.7108610477840874 0.1000 1.0000 0.9950207489532265 0.9950207489532266 0.2000 1.0000 0.9803279976447253 0.9803279976447253 0.5000 1.0000 0.8868188839700739 0.886818883970074 1.0000 1.0000 0.6480542736638853 0.6480542736638855 2.0000 1.0000 0.2658022288340797 0.2658022288340797 4.0000 1.0000 0.03661899347368653 0.03661899347368653 -0.2000 1.0000 0.9803279976447253 0.9803279976447253 -0.5000 1.0000 0.8868188839700739 0.886818883970074 -1.0000 1.0000 0.6480542736638853 0.6480542736638855 -2.0000 1.0000 0.2658022288340797 0.2658022288340797 jacobi_sn_test: jacobi_sn() evaluates the Jacobi elliptic function SN. U M Exact SN SN(U,M) 0.1000 0.0000 0.09983341664682815 0.09983341664682815 0.2000 0.0000 0.1986693307950612 0.1986693307950612 0.5000 0.0000 0.479425538604203 0.479425538604203 1.0000 0.0000 0.8414709848078965 0.8414709848078965 2.0000 0.0000 0.9092974268256817 0.9092974268256817 0.1000 0.5000 0.09975068547462485 0.09975068547462485 0.2000 0.5000 0.1980217429819704 0.1980217429819704 0.5000 0.5000 0.4707504736556573 0.4707504736556573 1.0000 0.5000 0.8030018248956439 0.8030018248956439 2.0000 0.5000 0.9946623253580177 0.9946623253580177 0.1000 1.0000 0.09966799462495582 0.09966799462495582 0.2000 1.0000 0.197375320224904 0.197375320224904 0.5000 1.0000 0.4621171572600098 0.4621171572600097 1.0000 1.0000 0.7615941559557649 0.7615941559557649 2.0000 1.0000 0.9640275800758169 0.9640275800758169 4.0000 1.0000 0.999329299739067 0.999329299739067 -0.2000 1.0000 -0.197375320224904 -0.197375320224904 -0.5000 1.0000 -0.4621171572600098 -0.4621171572600097 -1.0000 1.0000 -0.7615941559557649 -0.7615941559557649 -2.0000 1.0000 -0.9640275800758169 -0.9640275800758169 sncndn_test(): sncndn() evaluates the Jacobi elliptic functions SN, DN, CN. U M SN(U,M) CN(U,M) DN(U,M) Tab: 0.1000 0.0000 0.0998334166468282 0.9950041652780258 1.0000000000000000 Calc: 0.1000 0.0000 0.0998334166468282 0.9950041652780258 1.0000000000000000 Tab: 0.2000 0.0000 0.1986693307950612 0.9800665778412416 1.0000000000000000 Calc: 0.2000 0.0000 0.1986693307950612 0.9800665778412416 1.0000000000000000 Tab: 0.5000 0.0000 0.4794255386042030 0.8775825618903726 1.0000000000000000 Calc: 0.5000 0.0000 0.4794255386042030 0.8775825618903728 1.0000000000000000 Tab: 1.0000 0.0000 0.8414709848078965 0.5403023058681397 1.0000000000000000 Calc: 1.0000 0.0000 0.8414709848078965 0.5403023058681398 1.0000000000000000 Tab: 2.0000 0.0000 0.9092974268256817 -0.4161468365471424 1.0000000000000000 Calc: 2.0000 0.0000 0.9092974268256817 -0.4161468365471424 1.0000000000000000 Tab: 0.1000 0.5000 0.0997506854746248 0.9950124626090582 0.9975093485144243 Calc: 0.1000 0.5000 0.0997506854746248 0.9950124626090583 0.9975093485144244 Tab: 0.2000 0.5000 0.1980217429819704 0.9801976276784098 0.9901483195224799 Calc: 0.2000 0.5000 0.1980217429819704 0.9801976276784097 0.9901483195224799 Tab: 0.5000 0.5000 0.4707504736556573 0.8822663948904403 0.9429724257773857 Calc: 0.5000 0.5000 0.4707504736556573 0.8822663948904402 0.9429724257773857 Tab: 1.0000 0.5000 0.8030018248956439 0.5959765676721407 0.8231610016315963 Calc: 1.0000 0.5000 0.8030018248956439 0.5959765676721407 0.8231610016315963 Tab: 2.0000 0.5000 0.9946623253580177 -0.1031836155277618 0.7108610477840873 Calc: 2.0000 0.5000 0.9946623253580177 -0.1031836155277620 0.7108610477840874 Tab: 0.1000 1.0000 0.0996679946249558 0.9950207489532265 0.9950207489532265 Calc: 0.1000 1.0000 0.0996679946249558 0.9950207489532266 0.9950207489532266 Tab: 0.2000 1.0000 0.1973753202249040 0.9803279976447253 0.9803279976447253 Calc: 0.2000 1.0000 0.1973753202249040 0.9803279976447253 0.9803279976447253 Tab: 0.5000 1.0000 0.4621171572600098 0.8868188839700739 0.8868188839700739 Calc: 0.5000 1.0000 0.4621171572600097 0.8868188839700740 0.8868188839700740 Tab: 1.0000 1.0000 0.7615941559557649 0.6480542736638853 0.6480542736638853 Calc: 1.0000 1.0000 0.7615941559557649 0.6480542736638855 0.6480542736638855 Tab: 2.0000 1.0000 0.9640275800758169 0.2658022288340797 0.2658022288340797 Calc: 2.0000 1.0000 0.9640275800758169 0.2658022288340797 0.2658022288340797 Tab: 4.0000 1.0000 0.9993292997390670 0.0366189934736865 0.0366189934736865 Calc: 4.0000 1.0000 0.9993292997390670 0.0366189934736865 0.0366189934736865 Tab: -0.2000 1.0000 -0.1973753202249040 0.9803279976447253 0.9803279976447253 Calc: -0.2000 1.0000 -0.1973753202249040 0.9803279976447253 0.9803279976447253 Tab: -0.5000 1.0000 -0.4621171572600098 0.8868188839700739 0.8868188839700739 Calc: -0.5000 1.0000 -0.4621171572600097 0.8868188839700740 0.8868188839700740 Tab: -1.0000 1.0000 -0.7615941559557649 0.6480542736638853 0.6480542736638853 Calc: -1.0000 1.0000 -0.7615941559557649 0.6480542736638855 0.6480542736638855 Tab: -2.0000 1.0000 -0.9640275800758169 0.2658022288340797 0.2658022288340797 Calc: -2.0000 1.0000 -0.9640275800758169 0.2658022288340797 0.2658022288340797 elliptic_integral_test Normal end of execution. Tue Oct 19 11:33:55 2021