May 17 2007 2:48:12.368 PM NSWC_PRB FORTRAN90 version Tests for the NSWC library. TEST01: AI evaluates the Airy AI function. BI evaluates the Airy BI function. X Exact F AI(X) 0.0000 .355028 .355028 0.1000 .329203 .329203 0.2000 .303703 .303703 0.3000 .278806 .278807 0.4000 .254742 .254742 0.5000 .231694 .231694 0.6000 .209800 .209800 0.7000 .189162 .189162 0.8000 .169846 .169846 0.9000 .151887 .151887 1.0000 .135292 .135292 X Exact F BI(X) 0.0000 .614927 .614927 0.1000 .659862 .659862 0.2000 .705464 .705464 0.3000 .752486 .752486 0.4000 .801773 .801773 0.5000 .854277 .854277 0.6000 .911063 .911063 0.7000 .973329 .973329 0.8000 1.04242 1.04242 0.9000 1.11987 1.11987 1.0000 1.20742 1.20742 TEST02: BESI evaluates the Bessel I function. X Exact F BESI(0)(X) 0.0000 1.00000 1.00000 0.2000 1.01003 1.01003 0.4000 1.04040 1.04040 0.6000 1.09205 1.09205 0.8000 1.16651 1.16651 1.0000 1.26607 1.26607 1.2000 1.39373 1.39373 1.4000 1.55340 1.55340 1.6000 1.74998 1.74998 1.8000 1.98956 1.98956 2.0000 2.27959 2.27959 2.5000 3.28984 3.28984 3.0000 4.88079 4.88079 3.5000 7.37820 7.37820 4.0000 11.3019 11.3019 4.5000 17.4812 17.4812 5.0000 27.2399 27.2399 6.0000 67.2344 67.2344 8.0000 427.564 427.564 10.0000 2815.72 2815.73 TEST03: BESI evaluates the Bessel I function. X Exact F BESI(1)(X) 0.0000 .000000 .000000 0.2000 .100501 .100501 0.4000 .204027 .204027 0.6000 .313704 .313704 0.8000 .432865 .432865 1.0000 .565159 .565159 1.2000 .714678 .714678 1.4000 .886092 .886092 1.6000 1.08481 1.08481 1.8000 1.31717 1.31717 2.0000 1.59064 1.59064 2.5000 2.51672 2.51672 3.0000 3.95337 3.95337 3.5000 6.20583 6.20583 4.0000 9.75947 9.75946 4.5000 15.3892 15.3892 5.0000 24.3356 24.3356 6.0000 61.3419 61.3419 8.0000 399.873 399.873 10.0000 2670.99 2671.00 TEST04: BESI evaluates the Bessel I function. NU X Exact F BESI(NU)(X) 2 0.2000 0.501669E-02 0.501669E-02 2 1.0000 .135748 .135748 2 2.0000 .688948 .688949 2 2.5000 1.27647 1.27647 2 3.0000 2.24521 2.24521 2 5.0000 17.5056 17.5056 2 10.0000 2281.52 2281.53 2 20.0000 0.393128E+08 0.393128E+08 3 1.0000 0.221684E-01 0.221684E-01 3 2.0000 .212740 .212740 3 5.0000 10.3312 10.3312 3 10.0000 17.5838 1758.39 3 50.0000 0.267776E+21 0.267776E+21 5 1.0000 0.271463E-03 0.271463E-03 5 2.0000 0.982568E-02 0.982568E-02 5 5.0000 2.15797 2.15797 5 10.0000 777.188 777.191 5 50.0000 0.227855E+21 0.227855E+21 10 1.0000 0.275295E-09 0.275295E-09 10 2.0000 0.301696E-06 0.301696E-06 10 5.0000 0.458004E-02 0.458004E-02 10 10.0000 21.8917 21.8918 10 50.0000 0.107160E+21 0.107160E+21 20 1.0000 0.396684E-24 0.396683E-24 20 2.0000 0.431056E-18 0.431056E-18 20 5.0000 0.502424E-10 0.502424E-10 20 10.0000 0.125080E-03 0.125080E-03 20 50.0000 0.544201E+19 0.544200E+19 TEST05: BESJ evaluates the Bessel J function. X Exact F BESJ(0)(X) 0.0000 1.00000 1.00000 1.0000 .765198 .765198 2.0000 .223891 .223891 3.0000 -.260052 -.260052 4.0000 -.397150 -.397150 5.0000 -.177597 -.177597 6.0000 .150645 .150645 7.0000 .300079 .300079 8.0000 .171651 .171651 9.0000 -0.903336E-01 -0.903338E-01 10.0000 -.245936 -.245936 11.0000 -.171190 -.171190 12.0000 0.476893E-01 0.476895E-01 13.0000 .206926 .206926 14.0000 .171073 .171073 15.0000 -0.142245E-01 -0.142243E-01 TEST06: BESJ evaluates the Bessel J function. X Exact F BESJ(1)(X) 0.0000 .000000 .000000 1.0000 .440051 .440051 2.0000 .576725 .576725 3.0000 .339059 .339059 4.0000 -0.660433E-01 -0.660433E-01 5.0000 -.327579 -.327579 6.0000 -.276684 -.276684 7.0000 -0.468280E-02 -0.468291E-02 8.0000 .234636 .234637 9.0000 .245312 .245312 10.0000 0.434728E-01 0.434728E-01 11.0000 -.176785 -.176785 12.0000 -.223447 -.223448 13.0000 -0.703181E-01 -0.703180E-01 14.0000 .133375 .133375 15.0000 .205104 .205105 TEST07: BESJ evaluates the Bessel J function. NU X Exact F BESJ(NU)(X) 2 1.0000 .114903 .114903 2 2.0000 .352834 .352834 2 5.0000 0.465651E-01 0.465651E-01 2 10.0000 .254630 .254630 2 50.0000 -0.597128E-01 -0.597128E-01 5 1.0000 0.249758E-03 0.249758E-03 5 2.0000 0.703963E-02 0.703963E-02 5 5.0000 .261141 .261140 5 10.0000 -.234062 -.234061 5 50.0000 -0.814002E-01 -0.814001E-01 10 1.0000 0.263062E-09 0.263061E-09 10 2.0000 0.251539E-06 0.251538E-06 10 5.0000 0.146780E-02 0.146780E-02 10 10.0000 .207486 .207486 10 50.0000 -.113848 -.113848 20 1.0000 0.387350E-24 0.387349E-24 20 2.0000 0.391897E-18 0.391897E-18 20 5.0000 0.277033E-10 0.277033E-10 20 10.0000 0.115134E-04 0.115134E-04 20 50.0000 -.116704 -.116704 TEST08: BETA evaluates the Beta function. X Y Exact F BETA(X) 0.2000 1.0000 5.00000 5.00000 0.4000 1.0000 2.50000 2.50000 0.6000 1.0000 1.66667 1.66667 0.8000 1.0000 1.25000 1.25000 1.0000 0.2000 5.00000 5.00000 1.0000 0.4000 2.50000 2.50000 1.0000 1.0000 1.00000 1.00000 2.0000 2.0000 .166667 .166667 3.0000 3.0000 0.333333E-01 0.333333E-01 4.0000 4.0000 0.714286E-02 0.714286E-02 5.0000 5.0000 0.158730E-02 0.158730E-02 6.0000 2.0000 0.238095E-01 0.238095E-01 6.0000 3.0000 0.595238E-02 0.595238E-02 6.0000 4.0000 0.198413E-02 0.198413E-02 6.0000 5.0000 0.793651E-03 0.793651E-03 6.0000 6.0000 0.360750E-03 0.360750E-03 7.0000 7.0000 0.832501E-04 0.832501E-04 TEST09: BRATIO evaluates the normalized incomplete Beta function BETAI(A,B,X). A B X Exact F BETAI(A,B,X) 0.5000 0.5000 0.0100 0.637686E-01 0.637686E-01 0.5000 0.5000 0.1000 .204833 .204833 0.5000 0.5000 1.0000 1.00000 1.00000 1.0000 0.5000 0.0100 0.501260E-02 0.501256E-02 1.0000 0.5000 0.1000 0.513167E-01 0.513167E-01 1.0000 0.5000 1.0000 1.00000 1.00000 1.0000 1.0000 0.5000 .500000 .500000 5.0000 5.0000 0.5000 .500000 .500000 10.0000 0.5000 0.9000 .151641 .151641 10.0000 5.0000 0.5000 0.897827E-01 0.897827E-01 10.0000 5.0000 1.0000 1.00000 1.00000 10.0000 10.0000 0.5000 .500000 .500000 20.0000 5.0000 0.8000 .459877 .459877 20.0000 10.0000 0.6000 .214682 .214682 20.0000 10.0000 0.8000 .950737 .950736 20.0000 20.0000 0.5000 .500000 .500000 20.0000 20.0000 0.6000 .897941 .897942 30.0000 10.0000 0.7000 .224130 .224130 30.0000 10.0000 0.8000 .758641 .758640 40.0000 20.0000 0.7000 .700178 .700178 TEST10: CIN evaluates the cosine integral function. X Exact F CIN(X) 0.5000 0.618526E-01 0.618526E-01 0.6000 0.886608E-01 0.886607E-01 0.7000 .120026 .120026 0.8000 .155794 .155794 0.9000 .195787 .195787 1.0000 .239812 .239812 1.2000 .339078 .339078 1.4000 .451681 .451681 1.6000 .575487 .575487 1.8000 .708191 .708191 2.0000 .847382 .847382 2.5000 1.20764 1.20764 3.0000 1.55620 1.55620 3.5000 1.86211 1.86211 4.0000 2.10449 2.10449 4.5000 2.27478 2.27478 TEST11: DAWSON evaluates Dawson's integral. X Exact F DAWSON(X) 0.0000 .000000 .000000 0.1000 0.993360E-01 0.993360E-01 0.2000 .194751 .194751 0.3000 .282632 .282632 0.4000 .359943 .359943 0.5000 .424436 .424436 0.6000 .474763 .474763 0.7000 .510504 .510504 0.8000 .532102 .532102 0.9000 .540724 .540724 1.0000 .538080 .538080 1.1000 .526207 .526207 1.2000 .507273 .507273 1.3000 .483398 .483398 1.4000 .456507 .456507 1.5000 .428249 .428249 1.6000 .399940 .399940 1.7000 .372559 .372559 1.8000 .346773 .346773 1.9000 .322974 .322974 2.0000 .301340 .301340 TEST12: DILOGARITHM evaluates the DILOGARITHM function. X Exact F DILOGARITHM(X) 0.0000 1.64493 1.64493 0.0500 1.44063 1.44063 0.1000 1.29971 1.29971 0.1500 1.18058 1.18058 0.2000 1.07479 1.07479 0.2500 .978469 .978469 0.3000 .889378 .889378 0.3500 .806083 .806083 0.4000 .727586 .727586 0.4500 .653158 .653158 0.5000 .582241 .582241 TEST13: ELLPF evaluates the Jacobi elliptic function SN. A X Exact F SN(A,X) 0.0000 0.1000 0.998330E-01 0.998334E-01 0.0000 0.2000 .198670 .198669 0.0000 0.5000 .479430 .479426 0.0000 1.0000 .841470 .841471 0.0000 2.0000 .909300 .909297 0.5000 0.1000 0.997510E-01 0.997920E-01 0.5000 0.2000 .198020 .198345 0.5000 0.5000 .470750 .475083 0.5000 1.0000 .803000 .822636 0.5000 2.0000 .994660 .962898 1.0000 0.1000 0.996680E-01 0.996680E-01 1.0000 0.2000 .197380 .197375 1.0000 0.5000 .462120 .462117 1.0000 1.0000 .761590 .761594 1.0000 2.0000 .964030 .964028 1.0000 4.0000 .999330 .999329 1.0000 -0.2000 -.197380 -.197375 1.0000 -0.5000 -.462120 -.462117 1.0000 -1.0000 -.761590 -.761594 1.0000 -2.0000 -.964030 -.964028 TEST14: ERF evaluates the ERF function. X Exact F ERF(X) 0.1000 .112463 .112463 0.2000 .222703 .222703 0.3000 .328627 .328627 0.4000 .428392 .428392 0.5000 .520500 .520500 0.6000 .603856 .603856 0.7000 .677801 .677801 0.8000 .742101 .742101 0.9000 .796908 .796908 1.0000 .842701 .842701 1.1000 .880205 .880205 1.2000 .910314 .910314 1.3000 .934008 .934008 1.4000 .952285 .952285 1.5000 .966105 .966105 1.6000 .976348 .976348 1.7000 .983790 .983790 1.8000 .989091 .989091 1.9000 .992790 .992790 2.0000 .995322 .995322 TEST15: EXPLI evaluates the exponential integral function EI(X). X Exact F EI(X) 0.5000 .454220 .454220 0.6000 .769881 .769881 0.7000 1.06491 1.06491 0.8000 1.34740 1.34740 0.9000 1.62281 1.62281 1.0000 1.89512 1.89512 1.1000 2.16738 2.16738 1.2000 2.44209 2.44209 1.3000 2.72140 2.72140 1.4000 3.00721 3.00721 1.5000 3.30129 3.30129 1.6000 3.60532 3.60532 1.7000 3.92096 3.92096 1.8000 4.24987 4.24987 1.9000 4.59371 4.59371 2.0000 4.95423 4.95423 TEST16: EXPLI evaluates the exponential integral function E1(X). X Exact F E1(X) 0.5000 .559774 .559774 0.6000 .454379 .454379 0.7000 .373769 .373769 0.8000 .310597 .310597 0.9000 .260184 .260184 1.0000 .219384 .219384 1.1000 .185991 .185991 1.2000 .158408 .158408 1.3000 .135451 .135451 1.4000 .116219 .116219 1.5000 .100020 .100020 1.6000 0.863083E-01 0.863083E-01 1.7000 0.746546E-01 0.746546E-01 1.8000 0.647131E-01 0.647131E-01 1.9000 0.562044E-01 0.562044E-01 2.0000 0.489005E-01 0.489005E-01 TEST17: FRESNEL evaluates the Fresnel cosine integral function. X Exact F C(X) 0.0000 .000000 .000000 0.2000 .199921 .199921 0.4000 .397481 .397481 0.6000 .581095 .581095 0.8000 .722844 .722844 1.0000 .779893 .779893 1.2000 .715438 .715438 1.4000 .543096 .543096 1.6000 .365462 .365462 1.8000 .333633 .333633 2.0000 .488253 .488253 2.2000 .636286 .636286 2.4000 .554961 .554961 2.6000 .388938 .388938 2.8000 .467492 .467492 3.0000 .605721 .605721 FRESNEL evaluates the Fresnel sine integral function. X Exact F S(X) 0.0000 .000000 .000000 0.2000 0.418760E-02 0.418761E-02 0.4000 0.333594E-01 0.333594E-01 0.6000 .110540 .110540 0.8000 .249341 .249341 1.0000 .438259 .438259 1.2000 .623401 .623401 1.4000 .713525 .713525 1.6000 .638888 .638888 1.8000 .450939 .450939 2.0000 .343416 .343416 2.2000 .455705 .455705 2.4000 .619690 .619690 2.6000 .549989 .549989 2.8000 .391528 .391528 3.0000 .496313 .496313 TEST18: GAMMA evaluates the Gamma function. X Exact F GAMMA(X) 0.2000 4.59085 4.59084 0.4000 2.21816 2.21816 0.6000 1.48919 1.48919 0.8000 1.16423 1.16423 1.0000 1.00000 1.00000 1.1000 .951351 .951351 1.2000 .918169 .918169 1.3000 .897471 .897471 1.4000 .887264 .887264 1.5000 .886227 .886227 1.6000 .893515 .893515 1.7000 .908639 .908639 1.8000 .931384 .931384 1.9000 .961766 .961766 2.0000 1.00000 1.00000 10.0000 362880. 362880. 20.0000 0.121645E+18 0.121645E+18 30.0000 0.884176E+31 0.884176E+31 TEST19: GAMMP evaluates the normalized incomplete Gamma function P(A,X). A X Exact F GAMMP(A,X) 0.1000 0.0316 .742026 .742026 0.1000 0.3162 .911975 .911975 0.1000 1.5811 .989896 .989897 0.5000 0.0707 .293128 .293128 0.5000 0.7071 .765642 .765642 0.5000 3.5355 .992166 .992166 1.0000 0.1000 0.951626E-01 0.951626E-01 1.0000 1.0000 .632121 .632121 1.0000 5.0000 .993262 .993262 1.1000 0.1049 0.757471E-01 0.757471E-01 1.1000 1.0488 .607646 .607644 1.1000 5.2440 .993343 .993342 2.0000 0.1414 0.910540E-02 0.910536E-02 2.0000 1.4142 .413064 .413062 2.0000 7.0711 .993145 .993145 6.0000 2.4495 0.387318E-01 0.387315E-01 6.0000 12.2474 .982594 .982594 11.0000 16.5831 .940427 .940428 26.0000 25.4951 .486387 .486381 41.0000 44.8219 .735971 .735968 TEST20 MEXP computes the matrix exponential. The matrix A: 1 2 1 2.00000 2.00000 2 .000000 2.00000 The matrix exponential exp(A): 1 2 1 7.38906 14.7781 2 .000000 7.38906 TEST21: PSI evaluates the PSI function. X Exact F PSI(X) 1.0000 -.577216 -.577216 1.1000 -.423755 -.423755 1.2000 -.289040 -.289040 1.3000 -.169191 -.169191 1.4000 -0.613845E-01 -0.613846E-01 1.5000 -0.364900E-01 0.364900E-01 1.6000 .126047 .126047 1.7000 .208548 .208548 1.8000 .284991 .284991 1.9000 .356184 .356184 2.0000 .422784 .422784 TEST22: SI evaluates the sine integral function. X Exact F SI(X) 0.5000 .493107 .493107 0.6000 .588129 .588129 0.7000 .681222 .681222 0.8000 .772096 .772096 0.9000 .860471 .860471 1.0000 .946083 .946083 1.2000 1.10805 1.10805 1.4000 1.25623 1.25623 1.6000 1.38918 1.38918 1.8000 1.50582 1.50582 2.0000 1.60541 1.60541 2.5000 1.77852 1.77852 3.0000 1.84865 1.84865 3.5000 1.83313 1.83313 4.0000 1.75820 1.75820 4.5000 1.65414 1.65414 NSWC_PRB Normal end of execution. May 17 2007 2:48:12.378 PM