11 July 2008 11:00:26.358 PM TOMS453_PRB FORTRAN80 version Test TOMS algorithm 453, which computes the complex abscissas and weights for Gaussian quadrature formulas for Bromwich's integral. TEST01 Determine abscissas and weights for a variety of values of S and N. N = 6 S = 0.00000 XR XI WR WI 1 3 3.04821 -7.65179 -9.31577 -13.1079 2 2 5.46926 -4.51993 37.1570 124.634 3 1 6.48253 -1.49931 -27.8412 -315.065 4 1 6.48253 1.49931 -27.8412 315.065 5 2 5.46926 4.51993 37.1570 -124.634 6 3 3.04821 7.65179 -9.31577 13.1079 WR total 0.404484E-09 N = 6 S = 0.100000 XR XI WR WI 1 3 3.14713 -7.72402 -7.74421 -12.1904 2 2 5.56941 -4.55940 30.2273 114.988 3 1 6.58346 -1.51198 -22.4305 -290.584 4 1 6.58346 1.51198 -22.4305 290.584 5 2 5.56941 4.55940 30.2273 -114.988 6 3 3.14713 7.72402 -7.74421 12.1904 WR total 0.105114 N = 6 S = 1.00000 XR XI WR WI 1 3 4.03885 -8.34560 -0.683871 -4.88323 2 2 6.47051 -4.90012 -0.490540 46.7814 3 1 7.49064 -1.62150 1.67441 -122.163 4 1 7.49064 1.62150 1.67441 122.163 5 2 6.47051 4.90012 -0.490540 -46.7814 6 3 4.03885 8.34560 -0.683871 4.88323 WR total 1.00000 N = 6 S = 4.00000 XR XI WR WI 1 3 7.02212 -10.1435 0.353115E-01 -0.271604E-01 2 2 9.47263 -5.89402 -0.357778 0.458555 3 1 10.5053 -1.94212 0.405800 -1.70332 4 1 10.5053 1.94212 0.405800 1.70332 5 2 9.47263 5.89402 -0.357778 -0.458555 6 3 7.02212 10.1435 0.353115E-01 0.271604E-01 WR total 0.166667 N = 9 S = 0.00000 XR XI WR WI 1 5 3.97536 -13.1066 7.79286 -48.4906 2 4 7.28383 -9.61625 -650.747 710.452 3 3 9.20449 -6.34804 4931.37 -2600.30 4 2 10.2467 -3.15936 -14158.4 3267.65 5 1 10.5793 0.00000 19740.0 -0.00000 6 2 10.2467 3.15936 -14158.4 -3267.65 7 3 9.20449 6.34804 4931.37 2600.30 8 4 7.28383 9.61625 -650.747 -710.452 9 5 3.97536 13.1066 7.79286 48.4906 WR total -0.997278E-07 N = 9 S = 0.100000 XR XI WR WI 1 5 4.07435 -13.1844 8.82134 -40.9010 2 4 7.38342 -9.66975 -587.140 590.179 3 3 9.30475 -6.38204 4360.66 -2145.19 4 2 10.3474 -3.17594 -12442.5 2686.32 5 1 10.6802 0.00000 17320.4 -0.00000 6 2 10.3474 3.17594 -12442.5 -2686.32 7 3 9.30475 6.38204 4360.66 2145.19 8 4 7.38342 9.66975 -587.140 -590.179 9 5 4.07435 13.1844 8.82134 40.9010 WR total 0.105116 N = 9 S = 1.00000 XR XI WR WI 1 5 4.96613 -13.8647 5.70963 -7.30219 2 4 8.28004 -10.1384 -187.912 88.0612 3 3 10.2069 -6.68014 1272.27 -283.871 4 2 11.2533 -3.32134 -3542.61 328.959 5 1 11.5874 0.00000 4906.08 -0.00000 6 2 11.2533 3.32134 -3542.61 -328.959 7 3 10.2069 6.68014 1272.27 283.871 8 4 8.28004 10.1384 -187.912 -88.0612 9 5 4.96613 13.8647 5.70963 7.30219 WR total 1.00000 N = 9 S = 4.00000 XR XI WR WI 1 5 7.94671 -15.9216 0.276524E-01 0.931263E-02 2 4 11.2719 -11.5626 -0.771534 -0.429792 3 3 13.2119 -7.58890 6.10477 2.73269 4 2 14.2673 -3.76533 -19.7412 -4.62697 5 1 14.6045 0.00000 28.9272 -0.00000 6 2 14.2673 3.76533 -19.7412 4.62697 7 3 13.2119 7.58890 6.10477 -2.73269 8 4 11.2719 11.5626 -0.771534 0.429792 9 5 7.94671 15.9216 0.276524E-01 -0.931263E-02 WR total 0.166667 N = 12 S = 0.00000 XR XI WR WI 1 6 4.70228 -18.6891 84.9833 -77.5368 2 5 8.67056 -14.9254 -4375.78 -283.201 3 4 11.2246 -11.4674 40404.5 24658.7 4 3 12.9259 -8.13362 -144785. -182670. 5 2 13.9836 -4.86081 234619. 584473. 6 1 14.4930 -1.61733 -125948. -0.100930E+07 7 1 14.4930 1.61733 -125948. 0.100930E+07 8 2 13.9836 4.86081 234619. -584473. 9 3 12.9259 8.13362 -144785. 182670. 10 4 11.2246 11.4674 40404.5 -24658.7 11 5 8.67056 14.9254 -4375.78 283.201 12 6 4.70228 18.6891 84.9833 77.5368 WR total 0.209991E-01 N = 12 S = 0.100000 XR XI WR WI 1 6 4.80135 -18.7702 73.2458 -59.9454 2 5 8.76993 -14.9867 -3628.70 -400.168 3 4 11.3245 -11.5128 33151.5 21898.5 4 3 13.0262 -8.16504 -118182. -157306. 5 2 14.0843 -4.87931 190985. 498186. 6 1 14.5938 -1.62344 -102398. -857076. 7 1 14.5938 1.62344 -102398. 857076. 8 2 14.0843 4.87931 190985. -498186. 9 3 13.0262 8.16504 -118182. 157306. 10 4 11.3245 11.5128 33151.5 -21898.5 11 5 8.76993 14.9867 -3628.70 400.168 12 6 4.80135 18.7702 73.2458 59.9454 WR total 0.895583E-01 N = 12 S = 1.00000 XR XI WR WI 1 6 5.69358 -19.4846 15.2259 -3.43448 2 5 9.66460 -15.5270 -584.577 -318.896 3 4 12.2232 -11.9134 4876.46 6045.66 4 3 13.9287 -8.44250 -16513.1 -36822.0 5 2 14.9895 -5.04267 25882.9 109815. 6 1 15.5004 -1.67741 -13676.4 -184762. 7 1 15.5004 1.67741 -13676.4 184762. 8 2 14.9895 5.04267 25882.9 -109815. 9 3 13.9287 8.44250 -16513.1 36822.0 10 4 12.2232 11.9134 4876.46 -6045.66 11 5 9.66460 15.5270 -584.577 318.896 12 6 5.69358 19.4846 15.2259 3.43448 WR total 0.999222 N = 12 S = 4.00000 XR XI WR WI 1 6 8.67366 -21.6952 0.988037E-02 0.199247E-01 2 5 12.6508 -17.2045 0.498979E-01 -1.24831 3 4 15.2200 -13.1600 -3.45609 17.2573 4 3 16.9351 -9.30721 21.7632 -101.905 5 2 18.0029 -5.55230 -47.3026 313.027 6 1 18.5175 -1.84584 29.0190 -540.588 7 1 18.5175 1.84584 29.0190 540.588 8 2 18.0029 5.55230 -47.3026 -313.027 9 3 16.9351 9.30721 21.7632 101.905 10 4 15.2200 13.1600 -3.45609 -17.2573 11 5 12.6508 17.2045 0.498979E-01 1.24831 12 6 8.67366 21.6952 0.988037E-02 -0.199247E-01 WR total 0.166667 TOMS453_PRB Normal end of execution. 11 July 2008 11:00:26.364 PM