htswanalysis/MACS/lib/gsl/gsl-1.11/integration/qng.h

   1 /* integration/qng.h
   2  *
   3  * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Brian Gough
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18  */
  19
  20 /* Gauss-Kronrod-Patterson quadrature coefficients for use in
  21    quadpack routine qng. These coefficients were calculated with
  22    101 decimal digit arithmetic by L. W. Fullerton, Bell Labs, Nov
  23    1981. */
  24
  25 /* x1, abscissae common to the 10-, 21-, 43- and 87-point rule */
  26 static const double x1[5] = {
  27   0.973906528517171720077964012084452,
  28   0.865063366688984510732096688423493,
  29   0.679409568299024406234327365114874,
  30   0.433395394129247190799265943165784,
  31   0.148874338981631210884826001129720
  32 } ;
  33
  34 /* w10, weights of the 10-point formula */
  35 static const double w10[5] = {
  36   0.066671344308688137593568809893332,
  37   0.149451349150580593145776339657697,
  38   0.219086362515982043995534934228163,
  39   0.269266719309996355091226921569469,
  40   0.295524224714752870173892994651338
  41 } ;
  42
  43 /* x2, abscissae common to the 21-, 43- and 87-point rule */
  44 static const double x2[5] = {
  45   0.995657163025808080735527280689003,
  46   0.930157491355708226001207180059508,
  47   0.780817726586416897063717578345042,
  48   0.562757134668604683339000099272694,
  49   0.294392862701460198131126603103866
  50 } ;
  51
  52 /* w21a, weights of the 21-point formula for abscissae x1 */
  53 static const double w21a[5] = {
  54   0.032558162307964727478818972459390,
  55   0.075039674810919952767043140916190,
  56   0.109387158802297641899210590325805,
  57   0.134709217311473325928054001771707,
  58   0.147739104901338491374841515972068
  59 } ;
  60
  61 /* w21b, weights of the 21-point formula for abscissae x2 */
  62 static const double w21b[6] = {
  63   0.011694638867371874278064396062192,
  64   0.054755896574351996031381300244580,
  65   0.093125454583697605535065465083366,
  66   0.123491976262065851077958109831074,
  67   0.142775938577060080797094273138717,
  68   0.149445554002916905664936468389821
  69 } ;
  70
  71 /* x3, abscissae common to the 43- and 87-point rule */
  72 static const double x3[11] = {
  73   0.999333360901932081394099323919911,
  74   0.987433402908088869795961478381209,
  75   0.954807934814266299257919200290473,
  76   0.900148695748328293625099494069092,
  77   0.825198314983114150847066732588520,
  78   0.732148388989304982612354848755461,
  79   0.622847970537725238641159120344323,
  80   0.499479574071056499952214885499755,
  81   0.364901661346580768043989548502644,
  82   0.222254919776601296498260928066212,
  83   0.074650617461383322043914435796506
  84 } ;
  85
  86 /* w43a, weights of the 43-point formula for abscissae x1, x3 */
  87 static const double w43a[10] = {
  88   0.016296734289666564924281974617663,
  89   0.037522876120869501461613795898115,
  90   0.054694902058255442147212685465005,
  91   0.067355414609478086075553166302174,
  92   0.073870199632393953432140695251367,
  93   0.005768556059769796184184327908655,
  94   0.027371890593248842081276069289151,
  95   0.046560826910428830743339154433824,
  96   0.061744995201442564496240336030883,
  97   0.071387267268693397768559114425516
  98 } ;
  99
 100 /* w43b, weights of the 43-point formula for abscissae x3 */
 101 static const double w43b[12] = {
 102   0.001844477640212414100389106552965,
 103   0.010798689585891651740465406741293,
 104   0.021895363867795428102523123075149,
 105   0.032597463975345689443882222526137,
 106   0.042163137935191811847627924327955,
 107   0.050741939600184577780189020092084,
 108   0.058379395542619248375475369330206,
 109   0.064746404951445885544689259517511,
 110   0.069566197912356484528633315038405,
 111   0.072824441471833208150939535192842,
 112   0.074507751014175118273571813842889,
 113   0.074722147517403005594425168280423
 114 } ;
 115
 116 /* x4, abscissae of the 87-point rule */
 117 static const double x4[22] = {
 118   0.999902977262729234490529830591582,
 119   0.997989895986678745427496322365960,
 120   0.992175497860687222808523352251425,
 121   0.981358163572712773571916941623894,
 122   0.965057623858384619128284110607926,
 123   0.943167613133670596816416634507426,
 124   0.915806414685507209591826430720050,
 125   0.883221657771316501372117548744163,
 126   0.845710748462415666605902011504855,
 127   0.803557658035230982788739474980964,
 128   0.757005730685495558328942793432020,
 129   0.706273209787321819824094274740840,
 130   0.651589466501177922534422205016736,
 131   0.593223374057961088875273770349144,
 132   0.531493605970831932285268948562671,
 133   0.466763623042022844871966781659270,
 134   0.399424847859218804732101665817923,
 135   0.329874877106188288265053371824597,
 136   0.258503559202161551802280975429025,
 137   0.185695396568346652015917141167606,
 138   0.111842213179907468172398359241362,
 139   0.037352123394619870814998165437704
 140 } ;
 141
 142 /* w87a, weights of the 87-point formula for abscissae x1, x2, x3 */
 143 static const double w87a[21] = {
 144   0.008148377384149172900002878448190,
 145   0.018761438201562822243935059003794,
 146   0.027347451050052286161582829741283,
 147   0.033677707311637930046581056957588,
 148   0.036935099820427907614589586742499,
 149   0.002884872430211530501334156248695,
 150   0.013685946022712701888950035273128,
 151   0.023280413502888311123409291030404,
 152   0.030872497611713358675466394126442,
 153   0.035693633639418770719351355457044,
 154   0.000915283345202241360843392549948,
 155   0.005399280219300471367738743391053,
 156   0.010947679601118931134327826856808,
 157   0.016298731696787335262665703223280,
 158   0.021081568889203835112433060188190,
 159   0.025370969769253827243467999831710,
 160   0.029189697756475752501446154084920,
 161   0.032373202467202789685788194889595,
 162   0.034783098950365142750781997949596,
 163   0.036412220731351787562801163687577,
 164   0.037253875503047708539592001191226
 165 } ;
 166
 167 /* w87b, weights of the 87-point formula for abscissae x4    */
 168 static const double w87b[23] = {
 169   0.000274145563762072350016527092881,
 170   0.001807124155057942948341311753254,
 171   0.004096869282759164864458070683480,
 172   0.006758290051847378699816577897424,
 173   0.009549957672201646536053581325377,
 174   0.012329447652244853694626639963780,
 175   0.015010447346388952376697286041943,
 176   0.017548967986243191099665352925900,
 177   0.019938037786440888202278192730714,
 178   0.022194935961012286796332102959499,
 179   0.024339147126000805470360647041454,
 180   0.026374505414839207241503786552615,
 181   0.028286910788771200659968002987960,
 182   0.030052581128092695322521110347341,
 183   0.031646751371439929404586051078883,
 184   0.033050413419978503290785944862689,
 185   0.034255099704226061787082821046821,
 186   0.035262412660156681033782717998428,
 187   0.036076989622888701185500318003895,
 188   0.036698604498456094498018047441094,
 189   0.037120549269832576114119958413599,
 190   0.037334228751935040321235449094698,
 191   0.037361073762679023410321241766599
 192 } ;
 193