13 static const double dicoef[2][NDIM] = {
14 {2.54e-3,6.15e-3,1.62e-3,0.}, {4.42e-2,5.88e-2,0.343,0.}
16 static const double dite[2][NDIM] = {
17 {1.57e5,1.41e5,8.19e4,0.}, {3.74e5,1.41e5,1.59e5,0.}
19 static const double ditcrt[NDIM] = {1.2e4,1.2e4,1.1e4,1e20};
20 static const double aa[NDIM] = {2.3196,0.,0.,0.};
21 static const double bb[NDIM] = {10.7328,0.,0.,0.};
22 static const double cc[NDIM] = {6.8830,0.,0.,0.};
23 static const double dd[NDIM] = {-0.1824,0.,0.,0.};
24 static const double ff[NDIM] = {0.4101,0.1,0.1,0.};
63 fprintf(
ioQQQ,
" Beryli returns; frac=" );
69 fprintf(
ioQQQ,
"\n" );