C ============================================================================= C GAUSSIAN POINTS & WEIGHTS BLOCK DATA GAUSS04 IMPLICIT REAL*8 (A-Z) COMMON /GAUS04/ WW04(4),XX04(4) DATA WW04( 1),XX04( 1) / 3.47854853E-01 , -8.61136317E-01 / DATA WW04( 2),XX04( 2) / 6.52145147E-01 , -3.39981049E-01 / DATA WW04( 3),XX04( 3) / 6.52145147E-01 , 3.39981049E-01 / DATA WW04( 4),XX04( 4) / 3.47854853E-01 , 8.61136317E-01 / END BLOCK DATA GAUSS08 IMPLICIT REAL*8 (A-Z) COMMON /GAUS08/ WW08(8),XX08(8) DATA WW08( 1),XX08( 1) / 1.01228535E-01 , -9.60289836E-01 / DATA WW08( 2),XX08( 2) / 2.22381040E-01 , -7.96666503E-01 / DATA WW08( 3),XX08( 3) / 3.13706636E-01 , -5.25532424E-01 / DATA WW08( 4),XX08( 4) / 3.62683773E-01 , -1.83434635E-01 / DATA WW08( 5),XX08( 5) / 3.62683773E-01 , 1.83434635E-01 / DATA WW08( 6),XX08( 6) / 3.13706636E-01 , 5.25532424E-01 / DATA WW08( 7),XX08( 7) / 2.22381040E-01 , 7.96666503E-01 / DATA WW08( 8),XX08( 8) / 1.01228535E-01 , 9.60289836E-01 / END BLOCK DATA GAUSS16 IMPLICIT REAL*8 (A-Z) COMMON /GAUS16/ WW16(16),XX16(16) DATA WW16( 1),XX16( 1) / 2.71524601E-02 , -9.89400923E-01 / DATA WW16( 2),XX16( 2) / 6.22535236E-02 , -9.44575012E-01 / DATA WW16( 3),XX16( 3) / 9.51585099E-02 , -8.65631223E-01 / DATA WW16( 4),XX16( 4) / 1.24628969E-01 , -7.55404413E-01 / DATA WW16( 5),XX16( 5) / 1.49595991E-01 , -6.17876232E-01 / DATA WW16( 6),XX16( 6) / 1.69156522E-01 , -4.58016783E-01 / DATA WW16( 7),XX16( 7) / 1.82603419E-01 , -2.81603545E-01 / DATA WW16( 8),XX16( 8) / 1.89450607E-01 , -9.50125083E-02 / DATA WW16( 9),XX16( 9) / 1.89450607E-01 , 9.50125083E-02 / DATA WW16(10),XX16(10) / 1.82603419E-01 , 2.81603545E-01 / DATA WW16(11),XX16(11) / 1.69156522E-01 , 4.58016783E-01 / DATA WW16(12),XX16(12) / 1.49595991E-01 , 6.17876232E-01 / DATA WW16(13),XX16(13) / 1.24628969E-01 , 7.55404413E-01 / DATA WW16(14),XX16(14) / 9.51585099E-02 , 8.65631223E-01 / DATA WW16(15),XX16(15) / 6.22535236E-02 , 9.44575012E-01 / DATA WW16(16),XX16(16) / 2.71524601E-02 , 9.89400923E-01 / END BLOCK DATA GAUSS24 IMPLICIT REAL*8 (A-Z) COMMON /GAUS24/ WW24(24),XX24(24) DATA WW24( 1),XX24( 1) / 1.23412302E-02 , -9.95187223E-01 / DATA WW24( 2),XX24( 2) / 2.85313893E-02 , -9.74728584E-01 / DATA WW24( 3),XX24( 3) / 4.42774370E-02 , -9.38274562E-01 / DATA WW24( 4),XX24( 4) / 5.92985861E-02 , -8.86415541E-01 / DATA WW24( 5),XX24( 5) / 7.33464807E-02 , -8.20001960E-01 / DATA WW24( 6),XX24( 6) / 8.61901641E-02 , -7.40124166E-01 / DATA WW24( 7),XX24( 7) / 9.76186544E-02 , -6.48093641E-01 / DATA WW24( 8),XX24( 8) / 1.07444271E-01 , -5.45421481E-01 / DATA WW24( 9),XX24( 9) / 1.15505666E-01 , -4.33793515E-01 / DATA WW24(10),XX24(10) / 1.21670470E-01 , -3.15042675E-01 / DATA WW24(11),XX24(11) / 1.25837460E-01 , -1.91118866E-01 / DATA WW24(12),XX24(12) / 1.27938196E-01 , -6.40568957E-02 / DATA WW24(13),XX24(13) / 1.27938196E-01 , 6.40568957E-02 / DATA WW24(14),XX24(14) / 1.25837460E-01 , 1.91118866E-01 / DATA WW24(15),XX24(15) / 1.21670470E-01 , 3.15042675E-01 / DATA WW24(16),XX24(16) / 1.15505666E-01 , 4.33793515E-01 / DATA WW24(17),XX24(17) / 1.07444271E-01 , 5.45421481E-01 / DATA WW24(18),XX24(18) / 9.76186544E-02 , 6.48093641E-01 / DATA WW24(19),XX24(19) / 8.61901641E-02 , 7.40124166E-01 / DATA WW24(20),XX24(20) / 7.33464807E-02 , 8.20001960E-01 / DATA WW24(21),XX24(21) / 5.92985861E-02 , 8.86415541E-01 / DATA WW24(22),XX24(22) / 4.42774370E-02 , 9.38274562E-01 / DATA WW24(23),XX24(23) / 2.85313893E-02 , 9.74728584E-01 / DATA WW24(24),XX24(24) / 1.23412302E-02 , 9.95187223E-01 / END BLOCK DATA GAUSS32 IMPLICIT REAL*8 (A-Z) COMMON /GAUS32/ WW32(32),XX32(32) DATA WW32( 1),XX32( 1) / 7.01860990E-03 , -9.97263849E-01 / DATA WW32( 2),XX32( 2) / 1.62743945E-02 , -9.85611498E-01 / DATA WW32( 3),XX32( 3) / 2.53920648E-02 , -9.64762270E-01 / DATA WW32( 4),XX32( 4) / 3.42738628E-02 , -9.34906065E-01 / DATA WW32( 5),XX32( 5) / 4.28358987E-02 , -8.96321177E-01 / DATA WW32( 6),XX32( 6) / 5.09980582E-02 , -8.49367619E-01 / DATA WW32( 7),XX32( 7) / 5.86840920E-02 , -7.94483781E-01 / DATA WW32( 8),XX32( 8) / 6.58222213E-02 , -7.32182145E-01 / DATA WW32( 9),XX32( 9) / 7.23457932E-02 , -6.63044274E-01 / DATA WW32(10),XX32(10) / 7.81938955E-02 , -5.87715745E-01 / DATA WW32(11),XX32(11) / 8.33119228E-02 , -5.06899893E-01 / DATA WW32(12),XX32(12) / 8.76520947E-02 , -4.21351284E-01 / DATA WW32(13),XX32(13) / 9.11738798E-02 , -3.31868589E-01 / DATA WW32(14),XX32(14) / 9.38443989E-02 , -2.39287362E-01 / DATA WW32(15),XX32(15) / 9.56387222E-02 , -1.44471958E-01 / DATA WW32(16),XX32(16) / 9.65400860E-02 , -4.83076647E-02 / DATA WW32(17),XX32(17) / 9.65400860E-02 , 4.83076647E-02 / DATA WW32(18),XX32(18) / 9.56387222E-02 , 1.44471958E-01 / DATA WW32(19),XX32(19) / 9.38443989E-02 , 2.39287362E-01 / DATA WW32(20),XX32(20) / 9.11738798E-02 , 3.31868589E-01 / DATA WW32(21),XX32(21) / 8.76520947E-02 , 4.21351284E-01 / DATA WW32(22),XX32(22) / 8.33119228E-02 , 5.06899893E-01 / DATA WW32(23),XX32(23) / 7.81938955E-02 , 5.87715745E-01 / DATA WW32(24),XX32(24) / 7.23457932E-02 , 6.63044274E-01 / DATA WW32(25),XX32(25) / 6.58222213E-02 , 7.32182145E-01 / DATA WW32(26),XX32(26) / 5.86840920E-02 , 7.94483781E-01 / DATA WW32(27),XX32(27) / 5.09980582E-02 , 8.49367619E-01 / DATA WW32(28),XX32(28) / 4.28358987E-02 , 8.96321177E-01 / DATA WW32(29),XX32(29) / 3.42738628E-02 , 9.34906065E-01 / DATA WW32(30),XX32(30) / 2.53920648E-02 , 9.64762270E-01 / DATA WW32(31),XX32(31) / 1.62743945E-02 , 9.85611498E-01 / DATA WW32(32),XX32(32) / 7.01860990E-03 , 9.97263849E-01 / END BLOCK DATA GAUSS64 IMPLICIT REAL*8 (A-Z) COMMON /GAUS64/ WW64(64),XX64(64) DATA WW64( 1),XX64( 1) / 1.78328075E-03 , -9.99305069E-01 / DATA WW64( 2),XX64( 2) / 4.14703321E-03 , -9.96340096E-01 / DATA WW64( 3),XX64( 3) / 6.50445791E-03 , -9.91013348E-01 / DATA WW64( 4),XX64( 4) / 8.84675980E-03 , -9.83336270E-01 / DATA WW64( 5),XX64( 5) / 1.11681391E-02 , -9.73326802E-01 / DATA WW64( 6),XX64( 6) / 1.34630483E-02 , -9.61008787E-01 / DATA WW64( 7),XX64( 7) / 1.57260299E-02 , -9.46411371E-01 / DATA WW64( 8),XX64( 8) / 1.79517157E-02 , -9.29569185E-01 / DATA WW64( 9),XX64( 9) / 2.01348234E-02 , -9.10522163E-01 / DATA WW64(10),XX64(10) / 2.22701747E-02 , -8.89315426E-01 / DATA WW64(11),XX64(11) / 2.43527032E-02 , -8.65999401E-01 / DATA WW64(12),XX64(12) / 2.63774693E-02 , -8.40629280E-01 / DATA WW64(13),XX64(13) / 2.83396728E-02 , -8.13265324E-01 / DATA WW64(14),XX64(14) / 3.02346572E-02 , -7.83972383E-01 / DATA WW64(15),XX64(15) / 3.20579298E-02 , -7.52819896E-01 / DATA WW64(16),XX64(16) / 3.38051617E-02 , -7.19881833E-01 / DATA WW64(17),XX64(17) / 3.54722142E-02 , -6.85236335E-01 / DATA WW64(18),XX64(18) / 3.70551273E-02 , -6.48965478E-01 / DATA WW64(19),XX64(19) / 3.85501534E-02 , -6.11155331E-01 / DATA WW64(20),XX64(20) / 3.99537422E-02 , -5.71895659E-01 / DATA WW64(21),XX64(21) / 4.12625633E-02 , -5.31279445E-01 / DATA WW64(22),XX64(22) / 4.24735136E-02 , -4.89403158E-01 / DATA WW64(23),XX64(23) / 4.35837246E-02 , -4.46366012E-01 / DATA WW64(24),XX64(24) / 4.45905589E-02 , -4.02270168E-01 / DATA WW64(25),XX64(25) / 4.54916283E-02 , -3.57220173E-01 / DATA WW64(26),XX64(26) / 4.62847948E-02 , -3.11322868E-01 / DATA WW64(27),XX64(27) / 4.69681844E-02 , -2.64687151E-01 / DATA WW64(28),XX64(28) / 4.75401655E-02 , -2.17423648E-01 / DATA WW64(29),XX64(29) / 4.79993895E-02 , -1.69644415E-01 / DATA WW64(30),XX64(30) / 4.83447611E-02 , -1.21462822E-01 / DATA WW64(31),XX64(31) / 4.85754684E-02 , -7.29931220E-02 / DATA WW64(32),XX64(32) / 4.86909561E-02 , -2.43502930E-02 / DATA WW64(33),XX64(33) / 4.86909561E-02 , 2.43502930E-02 / DATA WW64(34),XX64(34) / 4.85754684E-02 , 7.29931220E-02 / DATA WW64(35),XX64(35) / 4.83447611E-02 , 1.21462822E-01 / DATA WW64(36),XX64(36) / 4.79993895E-02 , 1.69644415E-01 / DATA WW64(37),XX64(37) / 4.75401655E-02 , 2.17423648E-01 / DATA WW64(38),XX64(38) / 4.69681844E-02 , 2.64687151E-01 / DATA WW64(39),XX64(39) / 4.62847948E-02 , 3.11322868E-01 / DATA WW64(40),XX64(40) / 4.54916283E-02 , 3.57220173E-01 / DATA WW64(41),XX64(41) / 4.45905589E-02 , 4.02270168E-01 / DATA WW64(42),XX64(42) / 4.35837246E-02 , 4.46366012E-01 / DATA WW64(43),XX64(43) / 4.24735136E-02 , 4.89403158E-01 / DATA WW64(44),XX64(44) / 4.12625633E-02 , 5.31279445E-01 / DATA WW64(45),XX64(45) / 3.99537422E-02 , 5.71895659E-01 / DATA WW64(46),XX64(46) / 3.85501534E-02 , 6.11155331E-01 / DATA WW64(47),XX64(47) / 3.70551273E-02 , 6.48965478E-01 / DATA WW64(48),XX64(48) / 3.54722142E-02 , 6.85236335E-01 / DATA WW64(49),XX64(49) / 3.38051617E-02 , 7.19881833E-01 / DATA WW64(50),XX64(50) / 3.20579298E-02 , 7.52819896E-01 / DATA WW64(51),XX64(51) / 3.02346572E-02 , 7.83972383E-01 / DATA WW64(52),XX64(52) / 2.83396728E-02 , 8.13265324E-01 / DATA WW64(53),XX64(53) / 2.63774693E-02 , 8.40629280E-01 / DATA WW64(54),XX64(54) / 2.43527032E-02 , 8.65999401E-01 / DATA WW64(55),XX64(55) / 2.22701747E-02 , 8.89315426E-01 / DATA WW64(56),XX64(56) / 2.01348234E-02 , 9.10522163E-01 / DATA WW64(57),XX64(57) / 1.79517157E-02 , 9.29569185E-01 / DATA WW64(58),XX64(58) / 1.57260299E-02 , 9.46411371E-01 / DATA WW64(59),XX64(59) / 1.34630483E-02 , 9.61008787E-01 / DATA WW64(60),XX64(60) / 1.11681391E-02 , 9.73326802E-01 / DATA WW64(61),XX64(61) / 8.84675980E-03 , 9.83336270E-01 / DATA WW64(62),XX64(62) / 6.50445791E-03 , 9.91013348E-01 / DATA WW64(63),XX64(63) / 4.14703321E-03 , 9.96340096E-01 / DATA WW64(64),XX64(64) / 1.78328075E-03 , 9.99305069E-01 / END BLOCK DATA GAUSS96 IMPLICIT REAL*8 (A-Z) COMMON /GAUS96/ WW96(96),XX96(96) DATA WW96( 1),XX96( 1) / 7.96792039E-04 , -9.99689519E-01 / DATA WW96( 2),XX96( 2) / 1.85396080E-03 , -9.98364389E-01 / DATA WW96( 3),XX96( 3) / 2.91073183E-03 , -9.95981872E-01 / DATA WW96( 4),XX96( 4) / 3.96455452E-03 , -9.92543876E-01 / DATA WW96( 5),XX96( 5) / 5.01420256E-03 , -9.88054097E-01 / DATA WW96( 6),XX96( 6) / 6.05854532E-03 , -9.82517242E-01 / DATA WW96( 7),XX96( 7) / 7.09647080E-03 , -9.75939155E-01 / DATA WW96( 8),XX96( 8) / 8.12687725E-03 , -9.68326807E-01 / DATA WW96( 9),XX96( 9) / 9.14867129E-03 , -9.59688306E-01 / DATA WW96(10),XX96(10) / 1.01607703E-02 , -9.50032711E-01 / DATA WW96(11),XX96(11) / 1.11621022E-02 , -9.39370334E-01 / DATA WW96(12),XX96(12) / 1.21516045E-02 , -9.27712440E-01 / DATA WW96(13),XX96(13) / 1.31282294E-02 , -9.15071428E-01 / DATA WW96(14),XX96(14) / 1.40909422E-02 , -9.01460648E-01 / DATA WW96(15),XX96(15) / 1.50387213E-02 , -8.86894524E-01 / DATA WW96(16),XX96(16) / 1.59705635E-02 , -8.71388495E-01 / DATA WW96(17),XX96(17) / 1.68854799E-02 , -8.54959011E-01 / DATA WW96(18),XX96(18) / 1.77825019E-02 , -8.37623537E-01 / DATA WW96(19),XX96(19) / 1.86606795E-02 , -8.19400311E-01 / DATA WW96(20),XX96(20) / 1.95190813E-02 , -8.00308764E-01 / DATA WW96(21),XX96(21) / 2.03567967E-02 , -7.80369043E-01 / DATA WW96(22),XX96(22) / 2.11729407E-02 , -7.59602368E-01 / DATA WW96(23),XX96(23) / 2.19666436E-02 , -7.38030672E-01 / DATA WW96(24),XX96(24) / 2.27370691E-02 , -7.15676785E-01 / DATA WW96(25),XX96(25) / 2.34833993E-02 , -6.92564547E-01 / DATA WW96(26),XX96(26) / 2.42048409E-02 , -6.68718338E-01 / DATA WW96(27),XX96(27) / 2.49006338E-02 , -6.44163430E-01 / DATA WW96(28),XX96(28) / 2.55700368E-02 , -6.18925869E-01 / DATA WW96(29),XX96(29) / 2.62123402E-02 , -5.93032360E-01 / DATA WW96(30),XX96(30) / 2.68268660E-02 , -5.66510439E-01 / DATA WW96(31),XX96(31) / 2.74129622E-02 , -5.39388120E-01 / DATA WW96(32),XX96(32) / 2.79700067E-02 , -5.11694193E-01 / DATA WW96(33),XX96(33) / 2.84974109E-02 , -4.83457983E-01 / DATA WW96(34),XX96(34) / 2.89946143E-02 , -4.54709411E-01 / DATA WW96(35),XX96(35) / 2.94610895E-02 , -4.25478995E-01 / DATA WW96(36),XX96(36) / 2.98963450E-02 , -3.95797640E-01 / DATA WW96(37),XX96(37) / 3.02999150E-02 , -3.65696847E-01 / DATA WW96(38),XX96(38) / 3.06713767E-02 , -3.35208535E-01 / DATA WW96(39),XX96(39) / 3.10103334E-02 , -3.04364949E-01 / DATA WW96(40),XX96(40) / 3.13164257E-02 , -2.73198813E-01 / DATA WW96(41),XX96(41) / 3.15893292E-02 , -2.41743162E-01 / DATA WW96(42),XX96(42) / 3.18287574E-02 , -2.10031316E-01 / DATA WW96(43),XX96(43) / 3.20344567E-02 , -1.78096876E-01 / DATA WW96(44),XX96(44) / 3.22062038E-02 , -1.45973712E-01 / DATA WW96(45),XX96(45) / 3.23438235E-02 , -1.13695852E-01 / DATA WW96(46),XX96(46) / 3.24471630E-02 , -8.12974945E-02 / DATA WW96(47),XX96(47) / 3.25161181E-02 , -4.88129854E-02 / DATA WW96(48),XX96(48) / 3.25506143E-02 , -1.62767451E-02 / DATA WW96(49),XX96(49) / 3.25506143E-02 , 1.62767451E-02 / DATA WW96(50),XX96(50) / 3.25161181E-02 , 4.88129854E-02 / DATA WW96(51),XX96(51) / 3.24471630E-02 , 8.12974945E-02 / DATA WW96(52),XX96(52) / 3.23438235E-02 , 1.13695852E-01 / DATA WW96(53),XX96(53) / 3.22062038E-02 , 1.45973712E-01 / DATA WW96(54),XX96(54) / 3.20344567E-02 , 1.78096876E-01 / DATA WW96(55),XX96(55) / 3.18287574E-02 , 2.10031316E-01 / DATA WW96(56),XX96(56) / 3.15893292E-02 , 2.41743162E-01 / DATA WW96(57),XX96(57) / 3.13164257E-02 , 2.73198813E-01 / DATA WW96(58),XX96(58) / 3.10103334E-02 , 3.04364949E-01 / DATA WW96(59),XX96(59) / 3.06713767E-02 , 3.35208535E-01 / DATA WW96(60),XX96(60) / 3.02999150E-02 , 3.65696847E-01 / DATA WW96(61),XX96(61) / 2.98963450E-02 , 3.95797640E-01 / DATA WW96(62),XX96(62) / 2.94610895E-02 , 4.25478995E-01 / DATA WW96(63),XX96(63) / 2.89946143E-02 , 4.54709411E-01 / DATA WW96(64),XX96(64) / 2.84974109E-02 , 4.83457983E-01 / DATA WW96(65),XX96(65) / 2.79700067E-02 , 5.11694193E-01 / DATA WW96(66),XX96(66) / 2.74129622E-02 , 5.39388120E-01 / DATA WW96(67),XX96(67) / 2.68268660E-02 , 5.66510439E-01 / DATA WW96(68),XX96(68) / 2.62123402E-02 , 5.93032360E-01 / DATA WW96(69),XX96(69) / 2.55700368E-02 , 6.18925869E-01 / DATA WW96(70),XX96(70) / 2.49006338E-02 , 6.44163430E-01 / DATA WW96(71),XX96(71) / 2.42048409E-02 , 6.68718338E-01 / DATA WW96(72),XX96(72) / 2.34833993E-02 , 6.92564547E-01 / DATA WW96(73),XX96(73) / 2.27370691E-02 , 7.15676785E-01 / DATA WW96(74),XX96(74) / 2.19666436E-02 , 7.38030672E-01 / DATA WW96(75),XX96(75) / 2.11729407E-02 , 7.59602368E-01 / DATA WW96(76),XX96(76) / 2.03567967E-02 , 7.80369043E-01 / DATA WW96(77),XX96(77) / 1.95190813E-02 , 8.00308764E-01 / DATA WW96(78),XX96(78) / 1.86606795E-02 , 8.19400311E-01 / DATA WW96(79),XX96(79) / 1.77825019E-02 , 8.37623537E-01 / DATA WW96(80),XX96(80) / 1.68854799E-02 , 8.54959011E-01 / DATA WW96(81),XX96(81) / 1.59705635E-02 , 8.71388495E-01 / DATA WW96(82),XX96(82) / 1.50387213E-02 , 8.86894524E-01 / DATA WW96(83),XX96(83) / 1.40909422E-02 , 9.01460648E-01 / DATA WW96(84),XX96(84) / 1.31282294E-02 , 9.15071428E-01 / DATA WW96(85),XX96(85) / 1.21516045E-02 , 9.27712440E-01 / DATA WW96(86),XX96(86) / 1.11621022E-02 , 9.39370334E-01 / DATA WW96(87),XX96(87) / 1.01607703E-02 , 9.50032711E-01 / DATA WW96(88),XX96(88) / 9.14867129E-03 , 9.59688306E-01 / DATA WW96(89),XX96(89) / 8.12687725E-03 , 9.68326807E-01 / DATA WW96(90),XX96(90) / 7.09647080E-03 , 9.75939155E-01 / DATA WW96(91),XX96(91) / 6.05854532E-03 , 9.82517242E-01 / DATA WW96(92),XX96(92) / 5.01420256E-03 , 9.88054097E-01 / DATA WW96(93),XX96(93) / 3.96455452E-03 , 9.92543876E-01 / DATA WW96(94),XX96(94) / 2.91073183E-03 , 9.95981872E-01 / DATA WW96(95),XX96(95) / 1.85396080E-03 , 9.98364389E-01 / DATA WW96(96),XX96(96) / 7.96792039E-04 , 9.99689519E-01 / END C =============================================================================