May 15 2018 7:37:07.198 PM SPECIAL_FUNCTIONS_TEST: FORTRAN90 version Test the SPECIAL_FUNCTIONS library. AIRYA_TEST Test AIRYA x Ai(x) Bi(x) Ai'(x) Bi'(x) 0.0 0.35502805 0.61492663 -0.25881940 0.44828836 10.0 0.11047533E-09 0.45564115E+09 -0.35206337E-09 0.14292361E+10 20.0 0.16916729E-26 0.21037650E+26 -0.75863916E-26 0.93818393E+26 30.0 0.32082176E-48 0.90572885E+47 -0.17598766E-47 0.49533045E+48 x Ai(-x) Bi(-x) Ai'(-x) Bi'(-x) 0.0 0.35502805 0.61492663 -0.25881940 0.44828836 10.0 0.40241238E-01 -0.31467983 0.99626504 0.11941411 20.0 -0.17640613 -0.20013931 0.89286286 -0.79142903 30.0 -0.87968188E-01 -0.22444694 1.2286206 -0.48369473 BETA_TEST: Test BETA. p q B(p,q) --------------------------------- 1.5 2.0 0.2666666667D+00 2.5 2.0 0.1142857143D+00 1.5 3.0 0.1523809524D+00 CISIA_TEST CISIA computes the cosine and sine integrals. x ci(x) si(x) ------------------------------------ 0.0 -0.10000000+301 0.0000000 5.0 -0.19002975 1.5499312 10.0 -0.45456433E-01 1.6583476 20.0 0.44419821E-01 1.5482417 30.0 -0.33032417E-01 1.5667565 40.0 0.19020008E-01 1.5869851 CISIB_TEST CISIB computes the cosine and sine integrals. x ci(x) si(x) ------------------------------------ 0.0 -0.10000000+301 0.0000000 5.0 -0.19002984 1.5499311 10.0 -0.45456288E-01 1.6583479 20.0 0.44420091E-01 1.5482414 30.0 -0.33032580E-01 1.5667566 40.0 0.19020075E-01 1.5869852 CJY01_TEST Test CJY01 Z = 4.0000000 + i * 2.0000000 n Re[Jn(z)] Im[Jn(z)] Re[Jn'(z)] Im[Jn'(z)] -------------------------------------------------------------------- -0.13787022D+01 0.39054236D+00 0.50735255D+00 0.12263041D+01 -0.50735255D+00 -0.12263041D+01 -0.11546013D+01 0.58506793D+00 n Re[Yn(z)] Im[Yn(z)] Re[Yn'(z)] Im[Yn'(z)] -------------------------------------------------------------------- -0.38145893D+00 -0.13291649D+01 -0.12793101D+01 0.51220420D+00 0.12793101D+01 -0.51220420D+00 -0.58610052D+00 -0.10987930D+01 COMELP_TEST COMELP computes complete elliptic integrals H(K), E(K). k K(k) E(K) --------------------------------- 0.00 1.570796 1.570796 0.25 1.596242 1.545957 0.50 1.685750 1.467462 0.75 1.910990 1.318472 1.00 Infinity 1.000000 HYGFX_TEST: HYGFX evaluates the hypergeometric function 2F1. A B C X F(A,B,C,X) -2.50000 3.30000 6.70000 0.250000 0.7235612934899781 -0.500000 3.30000 6.70000 0.250000 0.9361014508109931 0.500000 3.30000 6.70000 0.250000 1.068969481117925 2.50000 3.30000 6.70000 0.250000 1.405156320011212 -2.50000 3.30000 6.70000 0.550000 0.4696143163982162 -0.500000 3.30000 6.70000 0.550000 0.8518739037527251 0.500000 3.30000 6.70000 0.550000 1.179535753216077 2.50000 3.30000 6.70000 0.550000 2.399906290477784 -2.50000 3.30000 6.70000 0.850000 0.2910609592841473 -0.500000 3.30000 6.70000 0.850000 0.7554318669735642 0.500000 3.30000 6.70000 0.850000 1.351049685305985 2.50000 3.30000 6.70000 0.850000 5.738156552619188 3.30000 6.70000 -5.50000 0.250000 15090.66974870460 3.30000 6.70000 -0.500000 0.250000 -2163.147887081946 3.30000 6.70000 0.500000 0.250000 264.5167677380311 3.30000 6.70000 4.50000 0.250000 4.194691581903191 3.30000 6.70000 -5.50000 0.550000 10170777974.04883 3.30000 6.70000 -0.500000 0.550000 -3085477.230745004 3.30000 6.70000 0.500000 0.550000 119678.5951809830 3.30000 6.70000 4.50000 0.550000 58.09272870639462 3.30000 6.70000 -5.50000 0.850000 0.5868208761512401E+19 3.30000 6.70000 -0.500000 0.850000 -1021736954816.628 3.30000 6.70000 0.500000 0.850000 9237064769.046463 3.30000 6.70000 4.50000 0.850000 20396.91377601964 MSPHJ SPHJ evaluates spherical Bessel J functions n x jn(x) jn'(x) 0 0.9050000000D+00 0.8689780718D+00 -0.2776712617D+00 1 0.9050000000D+00 0.2776712617D+00 0.2553399244D+00 2 0.9050000000D+00 0.5147914933D-01 0.1070221479D+00 3 0.9050000000D+00 0.6743927972D-02 0.2167173288D-01 4 0.9050000000D+00 0.6838294584D-03 0.2965864666D-02 5 0.9050000000D+00 0.5658597917D-04 0.3086737954D-03 n x jn(x) jn'(x) 0 0.1000000000D+02 -0.5440211109D-01 -0.7846694180D-01 1 0.1000000000D+02 0.7846694180D-01 -0.7009549945D-01 2 0.1000000000D+02 0.7794219363D-01 0.5508428371D-01 3 0.1000000000D+02 -0.3949584498D-01 0.9374053162D-01 4 0.1000000000D+02 -0.1055892851D+00 0.1329879757D-01 5 0.1000000000D+02 -0.5553451162D-01 -0.7226857814D-01 SPECIAL_FUNCTIONS_TEST: Normal end of execution. May 15 2018 7:37:07.216 PM