79 cc = 2. * (
sqrt(13.) - 1.) / 3. ;
91 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
97 cout <<
"!!!!! WARNING: Not yet prepared !!!!!" << endl ;
103 double hh = 1./double(nn) ;
108 double x1, x2, x3, x4, x5 ;
116 for (
int i=0; i<mm; i++) {
118 x1 = hh * double(4*i) ;
119 x2 = hh * double(4*i+1) ;
120 x3 = hh * double(4*i+2) ;
121 x4 = hh * double(4*i+3) ;
122 x5 = hh * double(4*i+4) ;
124 integ += (hh/45.) * (14.*ff(x1,cc) + 64.*ff(x2,cc)
125 + 24.*ff(x3,cc) + 64.*ff(x4,cc)
130 rah = 2. *
exp(integ) ;