24 else if( dx < 1.0e-5 )
26 exp1_v = ((-x*0.5) - 1.0)*x;
28 else if( dx < 1.0e-3 )
30 exp1_v = (((-x*0.1666666667) - 0.5)*x - 1.0)*x;
34 exp1_v = 1.0 - exp(x);
116 static double ab[
NFREQ],
127 static int first =
true;
128 static double zero = 0.;
129 static double one = 1.;
130 static double two = 2.;
131 static double four = 4.;
132 static double ten = 10.;
133 static double con1 = 0.8559;
164 fprintf( stderr,
" increase ncm in hypho to at least%5ld\n",
172 for( i=1; i < 7000; i++ )
175 alo[i-1] = log10(ai);
176 fal[i] = alo[i-1] + fal[i-1];
197 con2 = con1*
POW2(fn/ zed);
200 gn0 = 2.3052328943 - 2.302585093*fal[n+n-1] - fn*0.61370563888 +
201 alo[n-1]*(fn + one)*2.30258093;
202 lmax1 =
MIN2(lmax,n-1);
210 for( i=0; i < lmin; i++ )
217 gnt = 2.*(
POW2( (
double)n ) - gl);
221 for( i=0; i < ncell; i++ )
225 freq[i] = anu[i]*zed2;
226 for( j=0; j < 2; j++ )
235 freq[0] =
MAX2(freq[0],en);
239 for( il=1; il <= ilmax; il++ )
245 con3 = con2/(two*al + one);
248 for( ifp=0; ifp < ncell; ifp++ )
250 if( freq[ifp] < fth )
261 se = freq[ifp] - fth;
264 fprintf( stderr,
" %4ld%12.4e%12.4e\n", ifp,
276 p = 0.5*log10(
exp1(-ee));
285 a = two*(fn - atan(fn*e)/e);
292 ab[m] = (gn0 + a)/2.302585 - p - (fn + two)*
295 gn1e = x*gne/(fn + fn);
298 g2[1][m] = gne*fn*x*(fn + fn - one);
299 g2[0][m] = gn1e*(fn + fn - one)*(four +
341 sl =
POW2( (
double)ll );
344 g12 = (sn4 - sl4 + (two*sl - fll)*x)*g22 -
345 sn4*(sn - sl)*(one +
POW2(fll + one)*se)*
350 sm =
POW2( (
double)lm );
353 g11 = (sn4 - sm4 + (two*sm + flm)*x)*g21 -
354 sn4*(sn -
POW2(flm + one))*(one +
393 alfac = fal[n+l] - fal[n-l-1] + two*(al - fn)*
403 for( i=1; i <= llm; i++ )
406 p1 *= one + ai*ai*se;
420 for( i=llk; i <= lll; i++ )
423 p2 *= one + ai*ai*se;
441 if( fabs(sum) < 35. )
451 anl = fac*x*con3*(g11*al + g12*(al + 1.));
452 anl *= 2.*(2.*al + 1.);
454 bfnu[ifp] += (
realnum)(anl*1e-18);
464 for( i=0; i < ncell; i++ )