[BACK]Return to so3_nc.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / hgm / so3 / src

Annotation of OpenXM/src/hgm/so3/src/so3_nc.c, Revision 1.7

1.7     ! takayama    1: /* $OpenXM: OpenXM/src/hgm/so3/src/so3_nc.c,v 1.6 2013/03/07 05:23:31 takayama Exp $ */
1.1       takayama    2: #include <stdio.h>
                      3: #include <math.h>
                      4: #ifdef USE_GSL_LIB
                      5: #include <gsl/gsl_errno.h>
                      6: #include <gsl/gsl_odeiv.h>
                      7: #else
                      8: #include "t-gsl_errno.h"
                      9: #include "t-gsl_odeiv.h"
                     10: #endif
                     11:
1.7     ! takayama   12: #include "oxprint.h"
        !            13:
1.6       takayama   14: #ifndef STANDALONE
                     15: void mh_check_intr(int n);
                     16: #endif
                     17:
1.1       takayama   18: /* gcc evalnc.c -lgsl -lblas -lm */
                     19: /* gcc evalnc.c `pkg-config --cflags gsl` `pkg-config --libs gsl` */
                     20:
                     21: int so3_func();
                     22: void so3_nc(double a[3],double t0,double y[4]);
                     23: void so3_evalByS(int deg,double a,double b,double c,double t,double f[4]);
1.6       takayama   24: int so3_usage(void);
1.1       takayama   25:
                     26: #define MDEG 20
                     27:
1.2       takayama   28: #define SO3_QUIET_DEFAULT 0
                     29: #define SO3_Deg_DEFAULT 10
                     30: int SO3_Quiet = SO3_QUIET_DEFAULT;
                     31: int SO3_Deg = SO3_Deg_DEFAULT;
1.1       takayama   32: #ifdef STANDALONE
                     33: main(int argc,char *argv[]) {
                     34:   double a[3];
                     35:   double y[4];
                     36:   double t0;
                     37:   int i,j;
                     38:
                     39:   t0 = 0.0001;  /* Small enough number seems to be good (Hueristics) */
                     40:   j = 0;
                     41:   for (i=1; i<argc; i++) {
                     42:        if (strcmp(argv[i],"--quiet")==0) SO3_Quiet = 1;
                     43:        else if (strcmp(argv[i],"--t0")==0) {
                     44:          i++; sscanf(argv[i],"%lg",&t0);
                     45:        }else if (strcmp(argv[i],"--deg")==0) {
                     46:          i++; sscanf(argv[i],"%d",&SO3_Deg);
                     47:          if ((SO3_Deg > MDEG-2) || (SO3_Deg < 0)) {
1.7     ! takayama   48:            oxprintfe("Error: deg should be less than %d\n",MDEG-2); return(-1);
1.1       takayama   49:          }
                     50:     } else if (j<3) sscanf(argv[i],"%lg",&(a[j++]));
                     51:   }
                     52:   if (j != 3) {
                     53:     so3_usage();
                     54:     return(-1);
                     55:   }
                     56:   so3_nc(a,t0,y);
1.7     ! takayama   57: // oxprintf("%lg, %lg, %lg, %lg\n",y[0],y[1],y[2],y[3]);
        !            58:   oxprintf("%.15e, %.15e, %.15e, %.15e\n",y[0],y[1],y[2],y[3]);
1.1       takayama   59: }
                     60: #else
                     61: void so3_main(double *in1,double *in2,double *in3,double *t0p,int *quiet,int *deg,double *out) {
                     62:   double a[3];
                     63:   double y[4];
                     64:   double t0;
                     65:   int i,j;
1.2       takayama   66:   SO3_Quiet = SO3_QUIET_DEFAULT;
                     67:   SO3_Deg = SO3_Deg_DEFAULT;
1.1       takayama   68:   if (*quiet) SO3_Quiet = 1;
                     69:   if (*deg) {
                     70:     SO3_Deg = *deg;
1.7     ! takayama   71:     if (!SO3_Quiet) oxprintfe("deg is set to %d\n",SO3_Deg);
1.1       takayama   72:   }
                     73:   t0 = 0.0001;  /* Small enough number seems to be good (Hueristics) */
                     74:   if (*t0p > 0.0) {
                     75:     t0 = *t0p;
1.7     ! takayama   76:     if (!SO3_Quiet) oxprintfe("t0 is set to %lf\n",t0);
1.1       takayama   77:   }
                     78:   j = 0;
                     79:   if ((SO3_Deg > MDEG-2) || (SO3_Deg < 0)) {
1.7     ! takayama   80:     oxprintfe("Error: deg should be less than %d\n",MDEG-2);
1.1       takayama   81:     *out = 0.0; return;
                     82:   }
                     83:   a[0] = *in1; a[1] = *in2; a[2] = *in3;
1.7     ! takayama   84:   //  oxprintfe("DEBUG: %lf,%lf,%lf,%lf\n",t0,a[0],a[1],a[2]);
1.1       takayama   85:   so3_nc(a,t0,y);
1.7     ! takayama   86: // oxprintf("%lg, %lg, %lg, %lg\n",y[0],y[1],y[2],y[3]);
        !            87:   if (!SO3_Quiet) oxprintf("%.15e, %.15e, %.15e, %.15e\n",y[0],y[1],y[2],y[3]);
1.1       takayama   88:   *out = y[0];
                     89: }
                     90: #endif
                     91:
1.6       takayama   92: int so3_usage(void) {
1.7     ! takayama   93:   oxprintfe("Usage: so3_nc a b c returns nc(a,b,c) and its gradients\n");
        !            94:   oxprintfe("   where nc is the normalization constant\n" );
        !            95:   oxprintfe("   of the Fisher distribution on SO(3) for the diagonal matrix diag(a,b,c).\n");
        !            96:   oxprintfe("   See http://arxiv.org/abs/1110.0721\n");
        !            97:   oxprintfe("Options:  --quiet  --t0 T0 --deg DEG\n");
        !            98:   oxprintfe("   Series is evaluated at T0*(a,b,c) and the value is extended to (a,b,c) by diff. eq.\n");
1.6       takayama   99:   return(0);
1.1       takayama  100: }
                    101:
                    102: /* Evaluate normalization constant */
                    103: double SO3_A[3];
                    104: double SO3_R; /* rho */
                    105: int so3_func();
                    106: void so3_nc(double a[3],double t0,double y[4]) {
                    107:   int i;
                    108:   int deg;
                    109:   double r;
                    110:   double y0[4];
                    111:   double myerr,myerr2;
1.3       takayama  112:   double aa;
1.1       takayama  113:   deg = SO3_Deg;
1.3       takayama  114:   for (i=0; i<3; i++) SO3_A[i]=a[i];
                    115:
                    116:   /* When the argument is small, eval it only by series */
                    117:   aa = 0.0; for (i=0; i<3; i++) aa += a[i]*a[i];
                    118:   if (aa < 0.01) {
                    119:     so3_evalByS(deg,a[0],a[1],a[2], 1.0, y);
                    120:     return;
                    121:   }
                    122:
1.1       takayama  123:   SO3_R = 0.0;
                    124:   r = a[0]-a[1]-a[2]; if (r > SO3_R) SO3_R = r;
                    125:   r = -a[0]+a[1]-a[2]; if (r > SO3_R) SO3_R = r;
                    126:   r = -a[0]-a[1]+a[2]; if (r > SO3_R) SO3_R = r;
                    127:   r = a[0]+a[1]+a[2]; if (r > SO3_R) SO3_R = r;
1.7     ! takayama  128:   if (!SO3_Quiet) oxprintfe("SO3_R=%lg, exp(SO3_R t) is the asymptotics of nc.\n",SO3_R);
1.1       takayama  129:
                    130:   const gsl_odeiv_step_type *T = gsl_odeiv_step_rkf45;
                    131:   gsl_odeiv_step *s = gsl_odeiv_step_alloc(T, 4); /* rank4 */
                    132:
                    133:   /*
                    134:    * Absolute error 1e-6
                    135:    * Relative error 0.0.
                    136:    */
                    137:   gsl_odeiv_control *c = gsl_odeiv_control_y_new(1e-6, 0.0);
                    138:
                    139:   /* rank 4 system */
                    140:   gsl_odeiv_evolve *e = gsl_odeiv_evolve_alloc(4);
                    141:   gsl_odeiv_system sys = {so3_func, NULL, 4, NULL};
                    142:   /* t : start, t1 : goal */
                    143:   double t = t0, t1 = 1.0;
                    144:   double h = 1e-6;
                    145:   if (!SO3_Quiet) {
1.7     ! takayama  146:        oxprintfe("Set initial values at t0*(a,b,c) by evaluating series and find relevant t0.\n");
        !           147:        oxprintfe("t0=%lf, a=a[0]=%lf, b=a[1]=%lf, c=a[2]=%lf\n",t0,a[0],a[1],a[2]);
1.1       takayama  148:   }
                    149:   do {
1.7     ! takayama  150:     if (!SO3_Quiet) oxprintfe("t0=%lf\n",t0);
1.1       takayama  151:     so3_evalByS(deg,a[0],a[1],a[2], t0, y0);
1.7     ! takayama  152:     if (!SO3_Quiet) oxprintfe("[%2d]: %lg,%lg,%lg,%lg\n",deg,y0[0],y0[1],y0[2],y0[3]);
1.1       takayama  153:     so3_evalByS(deg+1,a[0],a[1],a[2], t0, y);
1.7     ! takayama  154:     if (!SO3_Quiet) oxprintfe("[%2d]:  %lg,%lg,%lg,%lg\n",deg+1,y[0],y[1],y[2],y[3]);
1.1       takayama  155:     myerr=0.0;
                    156:     for (i=0; i<4; i++) {
                    157:       myerr2 = (y0[i]-y[i])/y0[i]; if (myerr2 <0) myerr2 = -myerr2;
                    158:       if (myerr2 > myerr) myerr = myerr2;  /* sup norm */
                    159:     }
                    160:     if (myerr < 1e-6) break;  /* should take smaller value? */
                    161:     t0 = t0/2.0;
                    162:   } while (1);
                    163:   t=t0;
                    164:   for (i=0; i<4; i++) y[i]=y[i]*exp(-SO3_R*t0);
1.7     ! takayama  165:   if (!SO3_Quiet) oxprintfe("[%2d]*exp(-SO3_R*t0): %lg,%lg,%lg,%lg\n",deg+1,y[0],y[1],y[2],y[3]);
        !           166:   if (!SO3_Quiet) oxprintfe("Result by HGM (solving ODE) ------ \n");
1.1       takayama  167:
                    168:   while (t < t1) {
                    169:     int status = gsl_odeiv_evolve_apply(e, c, s, &sys, &t, t1, &h, y);
                    170:     if (status != GSL_SUCCESS) break;
                    171:   }
1.7     ! takayama  172:   if (!SO3_Quiet) oxprintfe("t and V   : t=%.5e %.5e %.5e, %.5e %.5e\n", t, y[0], y[1],y[2],y[3]);
1.1       takayama  173:   for (i=0; i<4; i++) y[i]=y[i]*exp(SO3_R*t);
1.7     ! takayama  174:   if (!SO3_Quiet) oxprintfe("V*exp(SO3_R*1): t= %.5e %.5e %.5e, %.5e %.5e\n", t, y[0], y[1],y[2],y[3]);
        !           175:   if (!SO3_Quiet) oxprintfe("Returned value is V=[so3_nc(a,b,c),  c_a, c_b, c_c]\n");
1.1       takayama  176:
                    177:   gsl_odeiv_evolve_free(e);
                    178:   gsl_odeiv_control_free(c);
                    179:   gsl_odeiv_step_free(s);
                    180: }
                    181:
                    182: /* From FB2/Prog/gls_ode_test2b.c */
                    183:
                    184: /*
                    185: d/dt y_i = f_i(t, y_1, ..., y_n),  f_i=f[i]
                    186:  y[0] = y, y[1]=y'.
                    187:
                    188: Ref: Note: See the Corollary 1 of the paper rotation.c on this function func.
                    189:
                    190: */
                    191: int so3_func(double t, const double y[], double f[], void *params)
                    192: {
                    193:   extern double SO3_A[3];
1.5       takayama  194: #ifndef STANDALONE
                    195:   mh_check_intr(100);
                    196: #endif
1.1       takayama  197:   f[0] = SO3_A[0]*y[1]+SO3_A[1]*y[2]+SO3_A[2]*y[3]               -SO3_R*y[0];
                    198:   f[1] = SO3_A[0]*y[0]+SO3_A[2]*y[2]+SO3_A[1]*y[3] - (2/t)*y[1]  -SO3_R*y[1];
                    199:   f[2] = SO3_A[1]*y[0]+SO3_A[2]*y[1]+SO3_A[0]*y[3] - (2/t)*y[2]  -SO3_R*y[2];
                    200:   f[3] = SO3_A[2]*y[0]+SO3_A[1]*y[1]+SO3_A[0]*y[2] - (2/t)*y[3]  -SO3_R*y[3];
                    201:   return GSL_SUCCESS;
                    202: }
                    203:
                    204: /* Evaluation of nc by a series */
                    205: void so3_evalByS(int deg,double a,double b,double c,double t,double f[4]) {
                    206:   double Tnc[MDEG][MDEG][MDEG];
                    207:   double ex[MDEG][MDEG][MDEG];
                    208:   int i,j,k;
                    209:   if (deg >= MDEG) {
1.7     ! takayama  210:     oxprintfe("Error: degree is too high\n");
1.1       takayama  211:   }
                    212:   for (i=0; i<MDEG; i++)for (j=0; j<MDEG; j++)for (k=0; k<MDEG; k++)Tnc[i][j][k]=0.0;
                    213:   /* load("sord.rr"); forC(20); */
                    214:   /* begin table */
                    215: Tnc[19][19][19]=1.8373e-54;
                    216: Tnc[19][19][17]=6.62551e-52;
                    217: Tnc[19][17][19]=6.62551e-52;
                    218: Tnc[17][19][19]=6.62551e-52;
                    219: Tnc[18][18][18]=1.36444e-50;
                    220: Tnc[19][19][15]=1.9059e-49;
                    221: Tnc[19][17][17]=2.39276e-49;
                    222: Tnc[17][19][17]=2.39276e-49;
                    223: Tnc[19][15][19]=1.9059e-49;
                    224: Tnc[17][17][19]=2.39276e-49;
                    225: Tnc[15][19][19]=1.9059e-49;
                    226: Tnc[18][18][16]=4.41547e-48;
                    227: Tnc[18][16][18]=4.41547e-48;
                    228: Tnc[16][18][18]=4.41547e-48;
                    229: Tnc[19][19][13]=4.24702e-47;
                    230: Tnc[19][17][15]=6.89438e-47;
                    231: Tnc[17][19][15]=6.89438e-47;
                    232: Tnc[19][15][17]=6.89438e-47;
                    233: Tnc[17][17][17]=8.65399e-47;
                    234: Tnc[15][19][17]=6.89438e-47;
                    235: Tnc[19][13][19]=4.24702e-47;
                    236: Tnc[17][15][19]=6.89438e-47;
                    237: Tnc[15][17][19]=6.89438e-47;
                    238: Tnc[13][19][19]=4.24702e-47;
                    239: Tnc[18][18][14]=1.12446e-45;
                    240: Tnc[18][16][16]=1.43125e-45;
                    241: Tnc[16][18][16]=1.43125e-45;
                    242: Tnc[18][14][18]=1.12446e-45;
                    243: Tnc[16][16][18]=1.43125e-45;
                    244: Tnc[14][18][18]=1.12446e-45;
                    245: Tnc[19][19][11]=7.05695e-45;
                    246: Tnc[19][17][13]=1.53917e-44;
                    247: Tnc[17][19][13]=1.53917e-44;
                    248: Tnc[19][15][15]=1.99021e-44;
                    249: Tnc[17][17][15]=2.49763e-44;
                    250: Tnc[15][19][15]=1.99021e-44;
                    251: Tnc[19][13][17]=1.53917e-44;
                    252: Tnc[17][15][17]=2.49763e-44;
                    253: Tnc[15][17][17]=2.49763e-44;
                    254: Tnc[13][19][17]=1.53917e-44;
                    255: Tnc[19][11][19]=7.05695e-45;
                    256: Tnc[17][13][19]=1.53917e-44;
                    257: Tnc[15][15][19]=1.99021e-44;
                    258: Tnc[13][17][19]=1.53917e-44;
                    259: Tnc[11][19][19]=7.05695e-45;
                    260: Tnc[18][18][12]=2.17975e-43;
                    261: Tnc[18][16][14]=3.65163e-43;
                    262: Tnc[16][18][14]=3.65163e-43;
                    263: Tnc[18][14][16]=3.65163e-43;
                    264: Tnc[16][16][16]=4.64694e-43;
                    265: Tnc[14][18][16]=3.65163e-43;
                    266: Tnc[18][12][18]=2.17975e-43;
                    267: Tnc[16][14][18]=3.65163e-43;
                    268: Tnc[14][16][18]=3.65163e-43;
                    269: Tnc[12][18][18]=2.17975e-43;
                    270: Tnc[19][19][9]=8.30395e-43;
                    271: Tnc[19][17][11]=2.56293e-42;
                    272: Tnc[17][19][11]=2.56293e-42;
                    273: Tnc[19][15][13]=4.45255e-42;
                    274: Tnc[17][17][13]=5.58632e-42;
                    275: Tnc[15][19][13]=4.45255e-42;
                    276: Tnc[19][13][15]=4.45255e-42;
                    277: Tnc[17][15][15]=7.22178e-42;
                    278: Tnc[15][17][15]=7.22178e-42;
                    279: Tnc[13][19][15]=4.45255e-42;
                    280: Tnc[19][11][17]=2.56293e-42;
                    281: Tnc[17][13][17]=5.58632e-42;
                    282: Tnc[15][15][17]=7.22178e-42;
                    283: Tnc[13][17][17]=5.58632e-42;
                    284: Tnc[11][19][17]=2.56293e-42;
                    285: Tnc[19][9][19]=8.30395e-43;
                    286: Tnc[17][11][19]=2.56293e-42;
                    287: Tnc[15][13][19]=4.45255e-42;
                    288: Tnc[13][15][19]=4.45255e-42;
                    289: Tnc[11][17][19]=2.56293e-42;
                    290: Tnc[9][19][19]=8.30395e-43;
                    291: Tnc[18][18][10]=3.07779e-41;
                    292: Tnc[18][16][12]=7.09362e-41;
                    293: Tnc[16][18][12]=7.09362e-41;
                    294: Tnc[18][14][14]=9.3363e-41;
                    295: Tnc[16][16][14]=1.1878e-40;
                    296: Tnc[14][18][14]=9.3363e-41;
                    297: Tnc[18][12][16]=7.09362e-41;
                    298: Tnc[16][14][16]=1.1878e-40;
                    299: Tnc[14][16][16]=1.1878e-40;
                    300: Tnc[12][18][16]=7.09362e-41;
                    301: Tnc[18][10][18]=3.07779e-41;
                    302: Tnc[16][12][18]=7.09362e-41;
                    303: Tnc[14][14][18]=9.3363e-41;
                    304: Tnc[12][16][18]=7.09362e-41;
                    305: Tnc[10][18][18]=3.07779e-41;
                    306: Tnc[19][19][7]=6.42744e-41;
                    307: Tnc[19][17][9]=3.02313e-40;
                    308: Tnc[17][19][9]=3.02313e-40;
                    309: Tnc[19][15][11]=7.43206e-40;
                    310: Tnc[17][17][11]=9.3216e-40;
                    311: Tnc[15][19][11]=7.43206e-40;
                    312: Tnc[19][13][13]=9.98548e-40;
                    313: Tnc[17][15][13]=1.61866e-39;
                    314: Tnc[15][17][13]=1.61866e-39;
                    315: Tnc[13][19][13]=9.98548e-40;
                    316: Tnc[19][11][15]=7.43206e-40;
                    317: Tnc[17][13][15]=1.61866e-39;
                    318: Tnc[15][15][15]=2.09201e-39;
                    319: Tnc[13][17][15]=1.61866e-39;
                    320: Tnc[11][19][15]=7.43206e-40;
                    321: Tnc[19][9][17]=3.02313e-40;
                    322: Tnc[17][11][17]=9.3216e-40;
                    323: Tnc[15][13][17]=1.61866e-39;
                    324: Tnc[13][15][17]=1.61866e-39;
                    325: Tnc[11][17][17]=9.3216e-40;
                    326: Tnc[9][19][17]=3.02313e-40;
                    327: Tnc[19][7][19]=6.42744e-41;
                    328: Tnc[17][9][19]=3.02313e-40;
                    329: Tnc[15][11][19]=7.43206e-40;
                    330: Tnc[13][13][19]=9.98548e-40;
                    331: Tnc[11][15][19]=7.43206e-40;
                    332: Tnc[9][17][19]=3.02313e-40;
                    333: Tnc[7][19][19]=6.42744e-41;
                    334: Tnc[18][18][8]=2.97769e-39;
                    335: Tnc[18][16][10]=1.00404e-38;
                    336: Tnc[16][18][10]=1.00404e-38;
                    337: Tnc[18][14][12]=1.81805e-38;
                    338: Tnc[16][16][12]=2.31226e-38;
                    339: Tnc[14][18][12]=1.81805e-38;
                    340: Tnc[18][12][14]=1.81805e-38;
                    341: Tnc[16][14][14]=3.04251e-38;
                    342: Tnc[14][16][14]=3.04251e-38;
                    343: Tnc[12][18][14]=1.81805e-38;
                    344: Tnc[18][10][16]=1.00404e-38;
                    345: Tnc[16][12][16]=2.31226e-38;
                    346: Tnc[14][14][16]=3.04251e-38;
                    347: Tnc[12][16][16]=2.31226e-38;
                    348: Tnc[10][18][16]=1.00404e-38;
                    349: Tnc[18][8][18]=2.97769e-39;
                    350: Tnc[16][10][18]=1.00404e-38;
                    351: Tnc[14][12][18]=1.81805e-38;
                    352: Tnc[12][14][18]=1.81805e-38;
                    353: Tnc[10][16][18]=1.00404e-38;
                    354: Tnc[8][18][18]=2.97769e-39;
                    355: Tnc[19][19][5]=2.91875e-39;
                    356: Tnc[19][17][7]=2.34655e-38;
                    357: Tnc[17][19][7]=2.34655e-38;
                    358: Tnc[19][15][9]=8.79116e-38;
                    359: Tnc[17][17][9]=1.1022e-37;
                    360: Tnc[15][19][9]=8.79116e-38;
                    361: Tnc[19][13][11]=1.67142e-37;
                    362: Tnc[17][15][11]=2.70751e-37;
                    363: Tnc[15][17][11]=2.70751e-37;
                    364: Tnc[13][19][11]=1.67142e-37;
                    365: Tnc[19][11][13]=1.67142e-37;
                    366: Tnc[17][13][13]=3.63678e-37;
                    367: Tnc[15][15][13]=4.69881e-37;
                    368: Tnc[13][17][13]=3.63678e-37;
                    369: Tnc[11][19][13]=1.67142e-37;
                    370: Tnc[19][9][15]=8.79116e-38;
                    371: Tnc[17][11][15]=2.70751e-37;
                    372: Tnc[15][13][15]=4.69881e-37;
                    373: Tnc[13][15][15]=4.69881e-37;
                    374: Tnc[11][17][15]=2.70751e-37;
                    375: Tnc[9][19][15]=8.79116e-38;
                    376: Tnc[19][7][17]=2.34655e-38;
                    377: Tnc[17][9][17]=1.1022e-37;
                    378: Tnc[15][11][17]=2.70751e-37;
                    379: Tnc[13][13][17]=3.63678e-37;
                    380: Tnc[11][15][17]=2.70751e-37;
                    381: Tnc[9][17][17]=1.1022e-37;
                    382: Tnc[7][19][17]=2.34655e-38;
                    383: Tnc[19][5][19]=2.91875e-39;
                    384: Tnc[17][7][19]=2.34655e-38;
                    385: Tnc[15][9][19]=8.79116e-38;
                    386: Tnc[13][11][19]=1.67142e-37;
                    387: Tnc[11][13][19]=1.67142e-37;
                    388: Tnc[9][15][19]=8.79116e-38;
                    389: Tnc[7][17][19]=2.34655e-38;
                    390: Tnc[5][19][19]=2.91875e-39;
                    391: Tnc[18][18][6]=1.80281e-37;
                    392: Tnc[18][16][8]=9.74107e-37;
                    393: Tnc[16][18][8]=9.74107e-37;
                    394: Tnc[18][14][10]=2.58048e-36;
                    395: Tnc[16][16][10]=3.2807e-36;
                    396: Tnc[14][18][10]=2.58048e-36;
                    397: Tnc[18][12][12]=3.55017e-36;
                    398: Tnc[16][14][12]=5.93708e-36;
                    399: Tnc[14][16][12]=5.93708e-36;
                    400: Tnc[12][18][12]=3.55017e-36;
                    401: Tnc[18][10][14]=2.58048e-36;
                    402: Tnc[16][12][14]=5.93708e-36;
                    403: Tnc[14][14][14]=7.80966e-36;
                    404: Tnc[12][16][14]=5.93708e-36;
                    405: Tnc[10][18][14]=2.58048e-36;
                    406: Tnc[18][8][16]=9.74107e-37;
                    407: Tnc[16][10][16]=3.2807e-36;
                    408: Tnc[14][12][16]=5.93708e-36;
                    409: Tnc[12][14][16]=5.93708e-36;
                    410: Tnc[10][16][16]=3.2807e-36;
                    411: Tnc[8][18][16]=9.74107e-37;
                    412: Tnc[18][6][18]=1.80281e-37;
                    413: Tnc[16][8][18]=9.74107e-37;
                    414: Tnc[14][10][18]=2.58048e-36;
                    415: Tnc[12][12][18]=3.55017e-36;
                    416: Tnc[10][14][18]=2.58048e-36;
                    417: Tnc[8][16][18]=9.74107e-37;
                    418: Tnc[6][18][18]=1.80281e-37;
                    419: Tnc[19][19][3]=6.35428e-38;
                    420: Tnc[19][17][5]=1.0691e-36;
                    421: Tnc[17][19][5]=1.0691e-36;
                    422: Tnc[19][15][7]=6.84611e-36;
                    423: Tnc[17][17][7]=8.57921e-36;
                    424: Tnc[15][19][7]=6.84611e-36;
                    425: Tnc[19][13][9]=1.98356e-35;
                    426: Tnc[17][15][9]=3.21035e-35;
                    427: Tnc[15][17][9]=3.21035e-35;
                    428: Tnc[13][19][9]=1.98356e-35;
                    429: Tnc[19][11][11]=2.80689e-35;
                    430: Tnc[17][13][11]=6.10019e-35;
                    431: Tnc[15][15][11]=7.87855e-35;
                    432: Tnc[13][17][11]=6.10019e-35;
                    433: Tnc[11][19][11]=2.80689e-35;
                    434: Tnc[19][9][13]=1.98356e-35;
                    435: Tnc[17][11][13]=6.10019e-35;
                    436: Tnc[15][13][13]=1.05793e-34;
                    437: Tnc[13][15][13]=1.05793e-34;
                    438: Tnc[11][17][13]=6.10019e-35;
                    439: Tnc[9][19][13]=1.98356e-35;
                    440: Tnc[19][7][15]=6.84611e-36;
                    441: Tnc[17][9][15]=3.21035e-35;
                    442: Tnc[15][11][15]=7.87855e-35;
                    443: Tnc[13][13][15]=1.05793e-34;
                    444: Tnc[11][15][15]=7.87855e-35;
                    445: Tnc[9][17][15]=3.21035e-35;
                    446: Tnc[7][19][15]=6.84611e-36;
                    447: Tnc[19][5][17]=1.0691e-36;
                    448: Tnc[17][7][17]=8.57921e-36;
                    449: Tnc[15][9][17]=3.21035e-35;
                    450: Tnc[13][11][17]=6.10019e-35;
                    451: Tnc[11][13][17]=6.10019e-35;
                    452: Tnc[9][15][17]=3.21035e-35;
                    453: Tnc[7][17][17]=8.57921e-36;
                    454: Tnc[5][19][17]=1.0691e-36;
                    455: Tnc[19][3][19]=6.35428e-38;
                    456: Tnc[17][5][19]=1.0691e-36;
                    457: Tnc[15][7][19]=6.84611e-36;
                    458: Tnc[13][9][19]=1.98356e-35;
                    459: Tnc[11][11][19]=2.80689e-35;
                    460: Tnc[9][13][19]=1.98356e-35;
                    461: Tnc[7][15][19]=6.84611e-36;
                    462: Tnc[5][17][19]=1.0691e-36;
                    463: Tnc[3][19][19]=6.35428e-38;
                    464: Tnc[18][18][4]=5.88673e-36;
                    465: Tnc[18][16][6]=5.91702e-35;
                    466: Tnc[16][18][6]=5.91702e-35;
                    467: Tnc[18][14][8]=2.51178e-34;
                    468: Tnc[16][16][8]=3.1918e-34;
                    469: Tnc[14][18][8]=2.51178e-34;
                    470: Tnc[18][12][10]=5.05553e-34;
                    471: Tnc[16][14][10]=8.44721e-34;
                    472: Tnc[14][16][10]=8.44721e-34;
                    473: Tnc[12][18][10]=5.05553e-34;
                    474: Tnc[18][10][12]=5.05553e-34;
                    475: Tnc[16][12][12]=1.16178e-33;
                    476: Tnc[14][14][12]=1.52762e-33;
                    477: Tnc[12][16][12]=1.16178e-33;
                    478: Tnc[10][18][12]=5.05553e-34;
                    479: Tnc[18][8][14]=2.51178e-34;
                    480: Tnc[16][10][14]=8.44721e-34;
                    481: Tnc[14][12][14]=1.52762e-33;
                    482: Tnc[12][14][14]=1.52762e-33;
                    483: Tnc[10][16][14]=8.44721e-34;
                    484: Tnc[8][18][14]=2.51178e-34;
                    485: Tnc[18][6][16]=5.91702e-35;
                    486: Tnc[16][8][16]=3.1918e-34;
                    487: Tnc[14][10][16]=8.44721e-34;
                    488: Tnc[12][12][16]=1.16178e-33;
                    489: Tnc[10][14][16]=8.44721e-34;
                    490: Tnc[8][16][16]=3.1918e-34;
                    491: Tnc[6][18][16]=5.91702e-35;
                    492: Tnc[18][4][18]=5.88673e-36;
                    493: Tnc[16][6][18]=5.91702e-35;
                    494: Tnc[14][8][18]=2.51178e-34;
                    495: Tnc[12][10][18]=5.05553e-34;
                    496: Tnc[10][12][18]=5.05553e-34;
                    497: Tnc[8][14][18]=2.51178e-34;
                    498: Tnc[6][16][18]=5.91702e-35;
                    499: Tnc[4][18][18]=5.88673e-36;
                    500: Tnc[19][19][1]=4.18362e-37;
                    501: Tnc[19][17][3]=2.33662e-35;
                    502: Tnc[17][19][3]=2.33662e-35;
                    503: Tnc[19][15][5]=3.13132e-34;
                    504: Tnc[17][17][5]=3.9216e-34;
                    505: Tnc[15][19][5]=3.13132e-34;
                    506: Tnc[19][13][7]=1.55073e-33;
                    507: Tnc[17][15][7]=2.50705e-33;
                    508: Tnc[15][17][7]=2.50705e-33;
                    509: Tnc[13][19][7]=1.55073e-33;
                    510: Tnc[19][11][9]=3.34407e-33;
                    511: Tnc[17][13][9]=7.25682e-33;
                    512: Tnc[15][15][9]=9.36783e-33;
                    513: Tnc[13][17][9]=7.25682e-33;
                    514: Tnc[11][19][9]=3.34407e-33;
                    515: Tnc[19][9][11]=3.34407e-33;
                    516: Tnc[17][11][11]=1.02657e-32;
                    517: Tnc[15][13][11]=1.77879e-32;
                    518: Tnc[13][15][11]=1.77879e-32;
                    519: Tnc[11][17][11]=1.02657e-32;
                    520: Tnc[9][19][11]=3.34407e-33;
                    521: Tnc[19][7][13]=1.55073e-33;
                    522: Tnc[17][9][13]=7.25682e-33;
                    523: Tnc[15][11][13]=1.77879e-32;
                    524: Tnc[13][13][13]=2.38763e-32;
                    525: Tnc[11][15][13]=1.77879e-32;
                    526: Tnc[9][17][13]=7.25682e-33;
                    527: Tnc[7][19][13]=1.55073e-33;
                    528: Tnc[19][5][15]=3.13132e-34;
                    529: Tnc[17][7][15]=2.50705e-33;
                    530: Tnc[15][9][15]=9.36783e-33;
                    531: Tnc[13][11][15]=1.77879e-32;
                    532: Tnc[11][13][15]=1.77879e-32;
                    533: Tnc[9][15][15]=9.36783e-33;
                    534: Tnc[7][17][15]=2.50705e-33;
                    535: Tnc[5][19][15]=3.13132e-34;
                    536: Tnc[19][3][17]=2.33662e-35;
                    537: Tnc[17][5][17]=3.9216e-34;
                    538: Tnc[15][7][17]=2.50705e-33;
                    539: Tnc[13][9][17]=7.25682e-33;
                    540: Tnc[11][11][17]=1.02657e-32;
                    541: Tnc[9][13][17]=7.25682e-33;
                    542: Tnc[7][15][17]=2.50705e-33;
                    543: Tnc[5][17][17]=3.9216e-34;
                    544: Tnc[3][19][17]=2.33662e-35;
                    545: Tnc[19][1][19]=4.18362e-37;
                    546: Tnc[17][3][19]=2.33662e-35;
                    547: Tnc[15][5][19]=3.13132e-34;
                    548: Tnc[13][7][19]=1.55073e-33;
                    549: Tnc[11][9][19]=3.34407e-33;
                    550: Tnc[9][11][19]=3.34407e-33;
                    551: Tnc[7][13][19]=1.55073e-33;
                    552: Tnc[5][15][19]=3.13132e-34;
                    553: Tnc[3][17][19]=2.33662e-35;
                    554: Tnc[1][19][19]=4.18362e-37;
                    555: Tnc[18][18][2]=7.7507e-35;
                    556: Tnc[18][16][4]=1.93965e-33;
                    557: Tnc[16][18][4]=1.93965e-33;
                    558: Tnc[18][14][6]=1.53168e-32;
                    559: Tnc[16][16][6]=1.94516e-32;
                    560: Tnc[14][18][6]=1.53168e-32;
                    561: Tnc[18][12][8]=4.94009e-32;
                    562: Tnc[16][14][8]=8.24521e-32;
                    563: Tnc[14][16][8]=8.24521e-32;
                    564: Tnc[12][18][8]=4.94009e-32;
                    565: Tnc[18][10][10]=7.22723e-32;
                    566: Tnc[16][12][10]=1.65837e-31;
                    567: Tnc[14][14][10]=2.17953e-31;
                    568: Tnc[12][16][10]=1.65837e-31;
                    569: Tnc[10][18][10]=7.22723e-32;
                    570: Tnc[18][8][12]=4.94009e-32;
                    571: Tnc[16][10][12]=1.65837e-31;
                    572: Tnc[14][12][12]=2.99644e-31;
                    573: Tnc[12][14][12]=2.99644e-31;
                    574: Tnc[10][16][12]=1.65837e-31;
                    575: Tnc[8][18][12]=4.94009e-32;
                    576: Tnc[18][6][14]=1.53168e-32;
                    577: Tnc[16][8][14]=8.24521e-32;
                    578: Tnc[14][10][14]=2.17953e-31;
                    579: Tnc[12][12][14]=2.99644e-31;
                    580: Tnc[10][14][14]=2.17953e-31;
                    581: Tnc[8][16][14]=8.24521e-32;
                    582: Tnc[6][18][14]=1.53168e-32;
                    583: Tnc[18][4][16]=1.93965e-33;
                    584: Tnc[16][6][16]=1.94516e-32;
                    585: Tnc[14][8][16]=8.24521e-32;
                    586: Tnc[12][10][16]=1.65837e-31;
                    587: Tnc[10][12][16]=1.65837e-31;
                    588: Tnc[8][14][16]=8.24521e-32;
                    589: Tnc[6][16][16]=1.94516e-32;
                    590: Tnc[4][18][16]=1.93965e-33;
                    591: Tnc[18][2][18]=7.7507e-35;
                    592: Tnc[16][4][18]=1.93965e-33;
                    593: Tnc[14][6][18]=1.53168e-32;
                    594: Tnc[12][8][18]=4.94009e-32;
                    595: Tnc[10][10][18]=7.22723e-32;
                    596: Tnc[8][12][18]=4.94009e-32;
                    597: Tnc[6][14][18]=1.53168e-32;
                    598: Tnc[4][16][18]=1.93965e-33;
                    599: Tnc[2][18][18]=7.7507e-35;
                    600: Tnc[19][17][1]=1.54574e-34;
                    601: Tnc[17][19][1]=1.54574e-34;
                    602: Tnc[19][15][3]=6.87619e-33;
                    603: Tnc[17][17][3]=8.6046e-33;
                    604: Tnc[15][19][3]=6.87619e-33;
                    605: Tnc[19][13][5]=7.12626e-32;
                    606: Tnc[17][15][5]=1.15045e-31;
                    607: Tnc[15][17][5]=1.15045e-31;
                    608: Tnc[13][19][5]=7.12626e-32;
                    609: Tnc[19][11][7]=2.62665e-31;
                    610: Tnc[17][13][7]=5.68909e-31;
                    611: Tnc[15][15][7]=7.3395e-31;
                    612: Tnc[13][17][7]=5.68909e-31;
                    613: Tnc[11][19][7]=2.62665e-31;
                    614: Tnc[19][9][9]=4.00278e-31;
                    615: Tnc[17][11][9]=1.22596e-30;
                    616: Tnc[15][13][9]=2.12194e-30;
                    617: Tnc[13][15][9]=2.12194e-30;
                    618: Tnc[11][17][9]=1.22596e-30;
                    619: Tnc[9][19][9]=4.00278e-31;
                    620: Tnc[19][7][11]=2.62665e-31;
                    621: Tnc[17][9][11]=1.22596e-30;
                    622: Tnc[15][11][11]=3.00058e-30;
                    623: Tnc[13][13][11]=4.02567e-30;
                    624: Tnc[11][15][11]=3.00058e-30;
                    625: Tnc[9][17][11]=1.22596e-30;
                    626: Tnc[7][19][11]=2.62665e-31;
                    627: Tnc[19][5][13]=7.12626e-32;
                    628: Tnc[17][7][13]=5.68909e-31;
                    629: Tnc[15][9][13]=2.12194e-30;
                    630: Tnc[13][11][13]=4.02567e-30;
                    631: Tnc[11][13][13]=4.02567e-30;
                    632: Tnc[9][15][13]=2.12194e-30;
                    633: Tnc[7][17][13]=5.68909e-31;
                    634: Tnc[5][19][13]=7.12626e-32;
                    635: Tnc[19][3][15]=6.87619e-33;
                    636: Tnc[17][5][15]=1.15045e-31;
                    637: Tnc[15][7][15]=7.3395e-31;
                    638: Tnc[13][9][15]=2.12194e-30;
                    639: Tnc[11][11][15]=3.00058e-30;
                    640: Tnc[9][13][15]=2.12194e-30;
                    641: Tnc[7][15][15]=7.3395e-31;
                    642: Tnc[5][17][15]=1.15045e-31;
                    643: Tnc[3][19][15]=6.87619e-33;
                    644: Tnc[19][1][17]=1.54574e-34;
                    645: Tnc[17][3][17]=8.6046e-33;
                    646: Tnc[15][5][17]=1.15045e-31;
                    647: Tnc[13][7][17]=5.68909e-31;
                    648: Tnc[11][9][17]=1.22596e-30;
                    649: Tnc[9][11][17]=1.22596e-30;
                    650: Tnc[7][13][17]=5.68909e-31;
                    651: Tnc[5][15][17]=1.15045e-31;
                    652: Tnc[3][17][17]=8.6046e-33;
                    653: Tnc[1][19][17]=1.54574e-34;
                    654: Tnc[17][1][19]=1.54574e-34;
                    655: Tnc[15][3][19]=6.87619e-33;
                    656: Tnc[13][5][19]=7.12626e-32;
                    657: Tnc[11][7][19]=2.62665e-31;
                    658: Tnc[9][9][19]=4.00278e-31;
                    659: Tnc[7][11][19]=2.62665e-31;
                    660: Tnc[5][13][19]=7.12626e-32;
                    661: Tnc[3][15][19]=6.87619e-33;
                    662: Tnc[1][17][19]=1.54574e-34;
                    663: Tnc[18][18][0]=1.71749e-34;
                    664: Tnc[18][16][2]=2.56592e-32;
                    665: Tnc[16][18][2]=2.56592e-32;
                    666: Tnc[18][14][4]=5.04466e-31;
                    667: Tnc[16][16][4]=6.40127e-31;
                    668: Tnc[14][18][4]=5.04466e-31;
                    669: Tnc[18][12][6]=3.02663e-30;
                    670: Tnc[16][14][6]=5.04436e-30;
                    671: Tnc[14][16][6]=5.04436e-30;
                    672: Tnc[12][18][6]=3.02663e-30;
                    673: Tnc[18][10][8]=7.09536e-30;
                    674: Tnc[16][12][8]=1.625e-29;
                    675: Tnc[14][14][8]=2.13435e-29;
                    676: Tnc[12][16][8]=1.625e-29;
                    677: Tnc[10][18][8]=7.09536e-30;
                    678: Tnc[18][8][10]=7.09536e-30;
                    679: Tnc[16][10][10]=2.37641e-29;
                    680: Tnc[14][12][10]=4.28908e-29;
                    681: Tnc[12][14][10]=4.28908e-29;
                    682: Tnc[10][16][10]=2.37641e-29;
                    683: Tnc[8][18][10]=7.09536e-30;
                    684: Tnc[18][6][12]=3.02663e-30;
                    685: Tnc[16][8][12]=1.625e-29;
                    686: Tnc[14][10][12]=4.28908e-29;
                    687: Tnc[12][12][12]=5.8938e-29;
                    688: Tnc[10][14][12]=4.28908e-29;
                    689: Tnc[8][16][12]=1.625e-29;
                    690: Tnc[6][18][12]=3.02663e-30;
                    691: Tnc[18][4][14]=5.04466e-31;
                    692: Tnc[16][6][14]=5.04436e-30;
                    693: Tnc[14][8][14]=2.13435e-29;
                    694: Tnc[12][10][14]=4.28908e-29;
                    695: Tnc[10][12][14]=4.28908e-29;
                    696: Tnc[8][14][14]=2.13435e-29;
                    697: Tnc[6][16][14]=5.04436e-30;
                    698: Tnc[4][18][14]=5.04466e-31;
                    699: Tnc[18][2][16]=2.56592e-32;
                    700: Tnc[16][4][16]=6.40127e-31;
                    701: Tnc[14][6][16]=5.04436e-30;
                    702: Tnc[12][8][16]=1.625e-29;
                    703: Tnc[10][10][16]=2.37641e-29;
                    704: Tnc[8][12][16]=1.625e-29;
                    705: Tnc[6][14][16]=5.04436e-30;
                    706: Tnc[4][16][16]=6.40127e-31;
                    707: Tnc[2][18][16]=2.56592e-32;
                    708: Tnc[18][0][18]=1.71749e-34;
                    709: Tnc[16][2][18]=2.56592e-32;
                    710: Tnc[14][4][18]=5.04466e-31;
                    711: Tnc[12][6][18]=3.02663e-30;
                    712: Tnc[10][8][18]=7.09536e-30;
                    713: Tnc[8][10][18]=7.09536e-30;
                    714: Tnc[6][12][18]=3.02663e-30;
                    715: Tnc[4][14][18]=5.04466e-31;
                    716: Tnc[2][16][18]=2.56592e-32;
                    717: Tnc[0][18][18]=1.71749e-34;
                    718: Tnc[19][15][1]=4.57538e-32;
                    719: Tnc[17][17][1]=5.71923e-32;
                    720: Tnc[15][19][1]=4.57538e-32;
                    721: Tnc[19][13][3]=1.57398e-30;
                    722: Tnc[17][15][3]=2.5362e-30;
                    723: Tnc[15][17][3]=2.5362e-30;
                    724: Tnc[13][19][3]=1.57398e-30;
                    725: Tnc[19][11][5]=1.21405e-29;
                    726: Tnc[17][13][5]=2.62293e-29;
                    727: Tnc[15][15][5]=3.38111e-29;
                    728: Tnc[13][17][5]=2.62293e-29;
                    729: Tnc[11][19][5]=1.21405e-29;
                    730: Tnc[19][9][7]=3.16223e-29;
                    731: Tnc[17][11][7]=9.65617e-29;
                    732: Tnc[15][13][7]=1.66894e-28;
                    733: Tnc[13][15][7]=1.66894e-28;
                    734: Tnc[11][17][7]=9.65617e-29;
                    735: Tnc[9][19][7]=3.16223e-29;
                    736: Tnc[19][7][9]=3.16223e-29;
                    737: Tnc[17][9][9]=1.47094e-28;
                    738: Tnc[15][11][9]=3.59329e-28;
                    739: Tnc[13][13][9]=4.81788e-28;
                    740: Tnc[11][15][9]=3.59329e-28;
                    741: Tnc[9][17][9]=1.47094e-28;
                    742: Tnc[7][19][9]=3.16223e-29;
                    743: Tnc[19][5][11]=1.21405e-29;
                    744: Tnc[17][7][11]=9.65617e-29;
                    745: Tnc[15][9][11]=3.59329e-28;
                    746: Tnc[13][11][11]=6.80953e-28;
                    747: Tnc[11][13][11]=6.80953e-28;
                    748: Tnc[9][15][11]=3.59329e-28;
                    749: Tnc[7][17][11]=9.65617e-29;
                    750: Tnc[5][19][11]=1.21405e-29;
                    751: Tnc[19][3][13]=1.57398e-30;
                    752: Tnc[17][5][13]=2.62293e-29;
                    753: Tnc[15][7][13]=1.66894e-28;
                    754: Tnc[13][9][13]=4.81788e-28;
                    755: Tnc[11][11][13]=6.80953e-28;
                    756: Tnc[9][13][13]=4.81788e-28;
                    757: Tnc[7][15][13]=1.66894e-28;
                    758: Tnc[5][17][13]=2.62293e-29;
                    759: Tnc[3][19][13]=1.57398e-30;
                    760: Tnc[19][1][15]=4.57538e-32;
                    761: Tnc[17][3][15]=2.5362e-30;
                    762: Tnc[15][5][15]=3.38111e-29;
                    763: Tnc[13][7][15]=1.66894e-28;
                    764: Tnc[11][9][15]=3.59329e-28;
                    765: Tnc[9][11][15]=3.59329e-28;
                    766: Tnc[7][13][15]=1.66894e-28;
                    767: Tnc[5][15][15]=3.38111e-29;
                    768: Tnc[3][17][15]=2.5362e-30;
                    769: Tnc[1][19][15]=4.57538e-32;
                    770: Tnc[17][1][17]=5.71923e-32;
                    771: Tnc[15][3][17]=2.5362e-30;
                    772: Tnc[13][5][17]=2.62293e-29;
                    773: Tnc[11][7][17]=9.65617e-29;
                    774: Tnc[9][9][17]=1.47094e-28;
                    775: Tnc[7][11][17]=9.65617e-29;
                    776: Tnc[5][13][17]=2.62293e-29;
                    777: Tnc[3][15][17]=2.5362e-30;
                    778: Tnc[1][17][17]=5.71923e-32;
                    779: Tnc[15][1][19]=4.57538e-32;
                    780: Tnc[13][3][19]=1.57398e-30;
                    781: Tnc[11][5][19]=1.21405e-29;
                    782: Tnc[9][7][19]=3.16223e-29;
                    783: Tnc[7][9][19]=3.16223e-29;
                    784: Tnc[5][11][19]=1.21405e-29;
                    785: Tnc[3][13][19]=1.57398e-30;
                    786: Tnc[1][15][19]=4.57538e-32;
                    787: Tnc[18][16][0]=5.71923e-32;
                    788: Tnc[16][18][0]=5.71923e-32;
                    789: Tnc[18][14][2]=6.71235e-30;
                    790: Tnc[16][16][2]=8.50819e-30;
                    791: Tnc[14][18][2]=6.71235e-30;
                    792: Tnc[18][12][4]=1.00261e-28;
                    793: Tnc[16][14][4]=1.66785e-28;
                    794: Tnc[14][16][4]=1.66785e-28;
                    795: Tnc[12][18][4]=1.00261e-28;
                    796: Tnc[18][10][6]=4.37222e-28;
                    797: Tnc[16][12][6]=9.98824e-28;
                    798: Tnc[14][14][6]=1.31084e-27;
                    799: Tnc[12][16][6]=9.98824e-28;
                    800: Tnc[10][18][6]=4.37222e-28;
                    801: Tnc[18][8][8]=7.00612e-28;
                    802: Tnc[16][10][8]=2.33948e-27;
                    803: Tnc[14][12][8]=4.21641e-27;
                    804: Tnc[12][14][8]=4.21641e-27;
                    805: Tnc[10][16][8]=2.33948e-27;
                    806: Tnc[8][18][8]=7.00612e-28;
                    807: Tnc[18][6][10]=4.37222e-28;
                    808: Tnc[16][8][10]=2.33948e-27;
                    809: Tnc[14][10][10]=6.16308e-27;
                    810: Tnc[12][12][10]=8.46369e-27;
                    811: Tnc[10][14][10]=6.16308e-27;
                    812: Tnc[8][16][10]=2.33948e-27;
                    813: Tnc[6][18][10]=4.37222e-28;
                    814: Tnc[18][4][12]=1.00261e-28;
                    815: Tnc[16][6][12]=9.98824e-28;
                    816: Tnc[14][8][12]=4.21641e-27;
                    817: Tnc[12][10][12]=8.46369e-27;
                    818: Tnc[10][12][12]=8.46369e-27;
                    819: Tnc[8][14][12]=4.21641e-27;
                    820: Tnc[6][16][12]=9.98824e-28;
                    821: Tnc[4][18][12]=1.00261e-28;
                    822: Tnc[18][2][14]=6.71235e-30;
                    823: Tnc[16][4][14]=1.66785e-28;
                    824: Tnc[14][6][14]=1.31084e-27;
                    825: Tnc[12][8][14]=4.21641e-27;
                    826: Tnc[10][10][14]=6.16308e-27;
                    827: Tnc[8][12][14]=4.21641e-27;
                    828: Tnc[6][14][14]=1.31084e-27;
                    829: Tnc[4][16][14]=1.66785e-28;
                    830: Tnc[2][18][14]=6.71235e-30;
                    831: Tnc[18][0][16]=5.71923e-32;
                    832: Tnc[16][2][16]=8.50819e-30;
                    833: Tnc[14][4][16]=1.66785e-28;
                    834: Tnc[12][6][16]=9.98824e-28;
                    835: Tnc[10][8][16]=2.33948e-27;
                    836: Tnc[8][10][16]=2.33948e-27;
                    837: Tnc[6][12][16]=9.98824e-28;
                    838: Tnc[4][14][16]=1.66785e-28;
                    839: Tnc[2][16][16]=8.50819e-30;
                    840: Tnc[0][18][16]=5.71923e-32;
                    841: Tnc[16][0][18]=5.71923e-32;
                    842: Tnc[14][2][18]=6.71235e-30;
                    843: Tnc[12][4][18]=1.00261e-28;
                    844: Tnc[10][6][18]=4.37222e-28;
                    845: Tnc[8][8][18]=7.00612e-28;
                    846: Tnc[6][10][18]=4.37222e-28;
                    847: Tnc[4][12][18]=1.00261e-28;
                    848: Tnc[2][14][18]=6.71235e-30;
                    849: Tnc[0][16][18]=5.71923e-32;
                    850: Tnc[19][13][1]=1.05503e-29;
                    851: Tnc[17][15][1]=1.69558e-29;
                    852: Tnc[15][17][1]=1.69558e-29;
                    853: Tnc[13][19][1]=1.05503e-29;
                    854: Tnc[19][11][3]=2.70111e-28;
                    855: Tnc[17][13][3]=5.81585e-28;
                    856: Tnc[15][15][3]=7.48882e-28;
                    857: Tnc[13][17][3]=5.81585e-28;
                    858: Tnc[11][19][3]=2.70111e-28;
                    859: Tnc[19][9][5]=1.47227e-27;
                    860: Tnc[17][11][5]=4.47766e-27;
                    861: Tnc[15][13][5]=7.72439e-27;
                    862: Tnc[13][15][5]=7.72439e-27;
                    863: Tnc[11][17][5]=4.47766e-27;
                    864: Tnc[9][19][5]=1.47227e-27;
                    865: Tnc[19][7][7]=2.5164e-27;
                    866: Tnc[17][9][7]=1.16526e-26;
                    867: Tnc[15][11][7]=2.8394e-26;
                    868: Tnc[13][13][7]=3.80398e-26;
                    869: Tnc[11][15][7]=2.8394e-26;
                    870: Tnc[9][17][7]=1.16526e-26;
                    871: Tnc[7][19][7]=2.5164e-27;
                    872: Tnc[19][5][9]=1.47227e-27;
                    873: Tnc[17][7][9]=1.16526e-26;
                    874: Tnc[15][9][9]=4.32318e-26;
                    875: Tnc[13][11][9]=8.181e-26;
                    876: Tnc[11][13][9]=8.181e-26;
                    877: Tnc[9][15][9]=4.32318e-26;
                    878: Tnc[7][17][9]=1.16526e-26;
                    879: Tnc[5][19][9]=1.47227e-27;
                    880: Tnc[19][3][11]=2.70111e-28;
                    881: Tnc[17][5][11]=4.47766e-27;
                    882: Tnc[15][7][11]=2.8394e-26;
                    883: Tnc[13][9][11]=8.181e-26;
                    884: Tnc[11][11][11]=1.15557e-25;
                    885: Tnc[9][13][11]=8.181e-26;
                    886: Tnc[7][15][11]=2.8394e-26;
                    887: Tnc[5][17][11]=4.47766e-27;
                    888: Tnc[3][19][11]=2.70111e-28;
                    889: Tnc[19][1][13]=1.05503e-29;
                    890: Tnc[17][3][13]=5.81585e-28;
                    891: Tnc[15][5][13]=7.72439e-27;
                    892: Tnc[13][7][13]=3.80398e-26;
                    893: Tnc[11][9][13]=8.181e-26;
                    894: Tnc[9][11][13]=8.181e-26;
                    895: Tnc[7][13][13]=3.80398e-26;
                    896: Tnc[5][15][13]=7.72439e-27;
                    897: Tnc[3][17][13]=5.81585e-28;
                    898: Tnc[1][19][13]=1.05503e-29;
                    899: Tnc[17][1][15]=1.69558e-29;
                    900: Tnc[15][3][15]=7.48882e-28;
                    901: Tnc[13][5][15]=7.72439e-27;
                    902: Tnc[11][7][15]=2.8394e-26;
                    903: Tnc[9][9][15]=4.32318e-26;
                    904: Tnc[7][11][15]=2.8394e-26;
                    905: Tnc[5][13][15]=7.72439e-27;
                    906: Tnc[3][15][15]=7.48882e-28;
                    907: Tnc[1][17][15]=1.69558e-29;
                    908: Tnc[15][1][17]=1.69558e-29;
                    909: Tnc[13][3][17]=5.81585e-28;
                    910: Tnc[11][5][17]=4.47766e-27;
                    911: Tnc[9][7][17]=1.16526e-26;
                    912: Tnc[7][9][17]=1.16526e-26;
                    913: Tnc[5][11][17]=4.47766e-27;
                    914: Tnc[3][13][17]=5.81585e-28;
                    915: Tnc[1][15][17]=1.69558e-29;
                    916: Tnc[13][1][19]=1.05503e-29;
                    917: Tnc[11][3][19]=2.70111e-28;
                    918: Tnc[9][5][19]=1.47227e-27;
                    919: Tnc[7][7][19]=2.5164e-27;
                    920: Tnc[5][9][19]=1.47227e-27;
                    921: Tnc[3][11][19]=2.70111e-28;
                    922: Tnc[1][13][19]=1.05503e-29;
                    923: Tnc[18][14][0]=1.50718e-29;
                    924: Tnc[16][16][0]=1.90753e-29;
                    925: Tnc[14][18][0]=1.50718e-29;
                    926: Tnc[18][12][2]=1.34384e-27;
                    927: Tnc[16][14][2]=2.22969e-27;
                    928: Tnc[14][16][2]=2.22969e-27;
                    929: Tnc[12][18][2]=1.34384e-27;
                    930: Tnc[18][10][4]=1.45893e-26;
                    931: Tnc[16][12][4]=3.32156e-26;
                    932: Tnc[14][14][4]=4.35443e-26;
                    933: Tnc[12][16][4]=3.32156e-26;
                    934: Tnc[10][18][4]=1.45893e-26;
                    935: Tnc[18][8][6]=4.34867e-26;
                    936: Tnc[16][10][6]=1.44628e-25;
                    937: Tnc[14][12][6]=2.60166e-25;
                    938: Tnc[12][14][6]=2.60166e-25;
                    939: Tnc[10][16][6]=1.44628e-25;
                    940: Tnc[8][18][6]=4.34867e-26;
                    941: Tnc[18][6][8]=4.34867e-26;
                    942: Tnc[16][8][8]=2.3164e-25;
                    943: Tnc[14][10][8]=6.08691e-25;
                    944: Tnc[12][12][8]=8.35231e-25;
                    945: Tnc[10][14][8]=6.08691e-25;
                    946: Tnc[8][16][8]=2.3164e-25;
                    947: Tnc[6][18][8]=4.34867e-26;
                    948: Tnc[18][4][10]=1.45893e-26;
                    949: Tnc[16][6][10]=1.44628e-25;
                    950: Tnc[14][8][10]=6.08691e-25;
                    951: Tnc[12][10][10]=1.22009e-24;
                    952: Tnc[10][12][10]=1.22009e-24;
                    953: Tnc[8][14][10]=6.08691e-25;
                    954: Tnc[6][16][10]=1.44628e-25;
                    955: Tnc[4][18][10]=1.45893e-26;
                    956: Tnc[18][2][12]=1.34384e-27;
                    957: Tnc[16][4][12]=3.32156e-26;
                    958: Tnc[14][6][12]=2.60166e-25;
                    959: Tnc[12][8][12]=8.35231e-25;
                    960: Tnc[10][10][12]=1.22009e-24;
                    961: Tnc[8][12][12]=8.35231e-25;
                    962: Tnc[6][14][12]=2.60166e-25;
                    963: Tnc[4][16][12]=3.32156e-26;
                    964: Tnc[2][18][12]=1.34384e-27;
                    965: Tnc[18][0][14]=1.50718e-29;
                    966: Tnc[16][2][14]=2.22969e-27;
                    967: Tnc[14][4][14]=4.35443e-26;
                    968: Tnc[12][6][14]=2.60166e-25;
                    969: Tnc[10][8][14]=6.08691e-25;
                    970: Tnc[8][10][14]=6.08691e-25;
                    971: Tnc[6][12][14]=2.60166e-25;
                    972: Tnc[4][14][14]=4.35443e-26;
                    973: Tnc[2][16][14]=2.22969e-27;
                    974: Tnc[0][18][14]=1.50718e-29;
                    975: Tnc[16][0][16]=1.90753e-29;
                    976: Tnc[14][2][16]=2.22969e-27;
                    977: Tnc[12][4][16]=3.32156e-26;
                    978: Tnc[10][6][16]=1.44628e-25;
                    979: Tnc[8][8][16]=2.3164e-25;
                    980: Tnc[6][10][16]=1.44628e-25;
                    981: Tnc[4][12][16]=3.32156e-26;
                    982: Tnc[2][14][16]=2.22969e-27;
                    983: Tnc[0][16][16]=1.90753e-29;
                    984: Tnc[14][0][18]=1.50718e-29;
                    985: Tnc[12][2][18]=1.34384e-27;
                    986: Tnc[10][4][18]=1.45893e-26;
                    987: Tnc[8][6][18]=4.34867e-26;
                    988: Tnc[6][8][18]=4.34867e-26;
                    989: Tnc[4][10][18]=1.45893e-26;
                    990: Tnc[2][12][18]=1.34384e-27;
                    991: Tnc[0][14][18]=1.50718e-29;
                    992: Tnc[19][11][1]=1.82784e-27;
                    993: Tnc[17][13][1]=3.9168e-27;
                    994: Tnc[15][15][1]=5.03588e-27;
                    995: Tnc[13][17][1]=3.9168e-27;
                    996: Tnc[11][19][1]=1.82784e-27;
                    997: Tnc[19][9][3]=3.30672e-26;
                    998: Tnc[17][11][3]=1.00009e-25;
                    999: Tnc[15][13][3]=1.72078e-25;
                   1000: Tnc[13][15][3]=1.72078e-25;
                   1001: Tnc[11][17][3]=1.00009e-25;
                   1002: Tnc[9][19][3]=3.30672e-26;
                   1003: Tnc[19][7][5]=1.18268e-25;
                   1004: Tnc[17][9][5]=5.44273e-25;
                   1005: Tnc[15][11][5]=1.32174e-24;
                   1006: Tnc[13][13][5]=1.76882e-24;
                   1007: Tnc[11][15][5]=1.32174e-24;
                   1008: Tnc[9][17][5]=5.44273e-25;
                   1009: Tnc[7][19][5]=1.18268e-25;
                   1010: Tnc[19][5][7]=1.18268e-25;
                   1011: Tnc[17][7][7]=9.29809e-25;
                   1012: Tnc[15][9][7]=3.43579e-24;
                   1013: Tnc[13][11][7]=6.48942e-24;
                   1014: Tnc[11][13][7]=6.48942e-24;
                   1015: Tnc[9][15][7]=3.43579e-24;
                   1016: Tnc[7][17][7]=9.29809e-25;
                   1017: Tnc[5][19][7]=1.18268e-25;
                   1018: Tnc[19][3][9]=3.30672e-26;
                   1019: Tnc[17][5][9]=5.44273e-25;
                   1020: Tnc[15][7][9]=3.43579e-24;
                   1021: Tnc[13][9][9]=9.8744e-24;
                   1022: Tnc[11][11][9]=1.39364e-23;
                   1023: Tnc[9][13][9]=9.8744e-24;
                   1024: Tnc[7][15][9]=3.43579e-24;
                   1025: Tnc[5][17][9]=5.44273e-25;
                   1026: Tnc[3][19][9]=3.30672e-26;
                   1027: Tnc[19][1][11]=1.82784e-27;
                   1028: Tnc[17][3][11]=1.00009e-25;
                   1029: Tnc[15][5][11]=1.32174e-24;
                   1030: Tnc[13][7][11]=6.48942e-24;
                   1031: Tnc[11][9][11]=1.39364e-23;
                   1032: Tnc[9][11][11]=1.39364e-23;
                   1033: Tnc[7][13][11]=6.48942e-24;
                   1034: Tnc[5][15][11]=1.32174e-24;
                   1035: Tnc[3][17][11]=1.00009e-25;
                   1036: Tnc[1][19][11]=1.82784e-27;
                   1037: Tnc[17][1][13]=3.9168e-27;
                   1038: Tnc[15][3][13]=1.72078e-25;
                   1039: Tnc[13][5][13]=1.76882e-24;
                   1040: Tnc[11][7][13]=6.48942e-24;
                   1041: Tnc[9][9][13]=9.8744e-24;
                   1042: Tnc[7][11][13]=6.48942e-24;
                   1043: Tnc[5][13][13]=1.76882e-24;
                   1044: Tnc[3][15][13]=1.72078e-25;
                   1045: Tnc[1][17][13]=3.9168e-27;
                   1046: Tnc[15][1][15]=5.03588e-27;
                   1047: Tnc[13][3][15]=1.72078e-25;
                   1048: Tnc[11][5][15]=1.32174e-24;
                   1049: Tnc[9][7][15]=3.43579e-24;
                   1050: Tnc[7][9][15]=3.43579e-24;
                   1051: Tnc[5][11][15]=1.32174e-24;
                   1052: Tnc[3][13][15]=1.72078e-25;
                   1053: Tnc[1][15][15]=5.03588e-27;
                   1054: Tnc[13][1][17]=3.9168e-27;
                   1055: Tnc[11][3][17]=1.00009e-25;
                   1056: Tnc[9][5][17]=5.44273e-25;
                   1057: Tnc[7][7][17]=9.29809e-25;
                   1058: Tnc[5][9][17]=5.44273e-25;
                   1059: Tnc[3][11][17]=1.00009e-25;
                   1060: Tnc[1][13][17]=3.9168e-27;
                   1061: Tnc[11][1][19]=1.82784e-27;
                   1062: Tnc[9][3][19]=3.30672e-26;
                   1063: Tnc[7][5][19]=1.18268e-25;
                   1064: Tnc[5][7][19]=1.18268e-25;
                   1065: Tnc[3][9][19]=3.30672e-26;
                   1066: Tnc[1][11][19]=1.82784e-27;
                   1067: Tnc[18][12][0]=3.0464e-27;
                   1068: Tnc[16][14][0]=5.03588e-27;
                   1069: Tnc[14][16][0]=5.03588e-27;
                   1070: Tnc[12][18][0]=3.0464e-27;
                   1071: Tnc[18][10][2]=1.97406e-25;
                   1072: Tnc[16][12][2]=4.47298e-25;
                   1073: Tnc[14][14][2]=5.85505e-25;
                   1074: Tnc[12][16][2]=4.47298e-25;
                   1075: Tnc[10][18][2]=1.97406e-25;
                   1076: Tnc[18][8][4]=1.46481e-24;
                   1077: Tnc[16][10][4]=4.84455e-24;
                   1078: Tnc[14][12][4]=8.69215e-24;
                   1079: Tnc[12][14][4]=8.69215e-24;
                   1080: Tnc[10][16][4]=4.84455e-24;
                   1081: Tnc[8][18][4]=1.46481e-24;
                   1082: Tnc[18][6][6]=2.72473e-24;
                   1083: Tnc[16][8][6]=1.4424e-23;
                   1084: Tnc[14][10][6]=3.77738e-23;
                   1085: Tnc[12][12][6]=5.17761e-23;
                   1086: Tnc[10][14][6]=3.77738e-23;
                   1087: Tnc[8][16][6]=1.4424e-23;
                   1088: Tnc[6][18][6]=2.72473e-24;
                   1089: Tnc[18][4][8]=1.46481e-24;
                   1090: Tnc[16][6][8]=1.4424e-23;
                   1091: Tnc[14][8][8]=6.04614e-23;
                   1092: Tnc[12][10][8]=1.20962e-22;
                   1093: Tnc[10][12][8]=1.20962e-22;
                   1094: Tnc[8][14][8]=6.04614e-23;
                   1095: Tnc[6][16][8]=1.4424e-23;
                   1096: Tnc[4][18][8]=1.46481e-24;
                   1097: Tnc[18][2][10]=1.97406e-25;
                   1098: Tnc[16][4][10]=4.84455e-24;
                   1099: Tnc[14][6][10]=3.77738e-23;
                   1100: Tnc[12][8][10]=1.20962e-22;
                   1101: Tnc[10][10][10]=1.76555e-22;
                   1102: Tnc[8][12][10]=1.20962e-22;
                   1103: Tnc[6][14][10]=3.77738e-23;
                   1104: Tnc[4][16][10]=4.84455e-24;
                   1105: Tnc[2][18][10]=1.97406e-25;
                   1106: Tnc[18][0][12]=3.0464e-27;
                   1107: Tnc[16][2][12]=4.47298e-25;
                   1108: Tnc[14][4][12]=8.69215e-24;
                   1109: Tnc[12][6][12]=5.17761e-23;
                   1110: Tnc[10][8][12]=1.20962e-22;
                   1111: Tnc[8][10][12]=1.20962e-22;
                   1112: Tnc[6][12][12]=5.17761e-23;
                   1113: Tnc[4][14][12]=8.69215e-24;
                   1114: Tnc[2][16][12]=4.47298e-25;
                   1115: Tnc[0][18][12]=3.0464e-27;
                   1116: Tnc[16][0][14]=5.03588e-27;
                   1117: Tnc[14][2][14]=5.85505e-25;
                   1118: Tnc[12][4][14]=8.69215e-24;
                   1119: Tnc[10][6][14]=3.77738e-23;
                   1120: Tnc[8][8][14]=6.04614e-23;
                   1121: Tnc[6][10][14]=3.77738e-23;
                   1122: Tnc[4][12][14]=8.69215e-24;
                   1123: Tnc[2][14][14]=5.85505e-25;
                   1124: Tnc[0][16][14]=5.03588e-27;
                   1125: Tnc[14][0][16]=5.03588e-27;
                   1126: Tnc[12][2][16]=4.47298e-25;
                   1127: Tnc[10][4][16]=4.84455e-24;
                   1128: Tnc[8][6][16]=1.4424e-23;
                   1129: Tnc[6][8][16]=1.4424e-23;
                   1130: Tnc[4][10][16]=4.84455e-24;
                   1131: Tnc[2][12][16]=4.47298e-25;
                   1132: Tnc[0][14][16]=5.03588e-27;
                   1133: Tnc[12][0][18]=3.0464e-27;
                   1134: Tnc[10][2][18]=1.97406e-25;
                   1135: Tnc[8][4][18]=1.46481e-24;
                   1136: Tnc[6][6][18]=2.72473e-24;
                   1137: Tnc[4][8][18]=1.46481e-24;
                   1138: Tnc[2][10][18]=1.97406e-25;
                   1139: Tnc[0][12][18]=3.0464e-27;
                   1140: Tnc[19][9][1]=2.26652e-25;
                   1141: Tnc[17][11][1]=6.79956e-25;
                   1142: Tnc[15][13][1]=1.16564e-24;
                   1143: Tnc[13][15][1]=1.16564e-24;
                   1144: Tnc[11][17][1]=6.79956e-25;
                   1145: Tnc[9][19][1]=2.26652e-25;
                   1146: Tnc[19][7][3]=2.69039e-24;
                   1147: Tnc[17][9][3]=1.22716e-23;
                   1148: Tnc[15][11][3]=2.9659e-23;
                   1149: Tnc[13][13][3]=3.96317e-23;
                   1150: Tnc[11][15][3]=2.9659e-23;
                   1151: Tnc[9][17][3]=1.22716e-23;
                   1152: Tnc[7][19][3]=2.69039e-24;
                   1153: Tnc[19][5][5]=5.62966e-24;
                   1154: Tnc[17][7][5]=4.38402e-23;
                   1155: Tnc[15][9][5]=1.61095e-22;
                   1156: Tnc[13][11][5]=3.03484e-22;
                   1157: Tnc[11][13][5]=3.03484e-22;
                   1158: Tnc[9][15][5]=1.61095e-22;
                   1159: Tnc[7][17][5]=4.38402e-23;
                   1160: Tnc[5][19][5]=5.62966e-24;
                   1161: Tnc[19][3][7]=2.69039e-24;
                   1162: Tnc[17][5][7]=4.38402e-23;
                   1163: Tnc[15][7][7]=2.75032e-22;
                   1164: Tnc[13][9][7]=7.87749e-22;
                   1165: Tnc[11][11][7]=1.11059e-21;
                   1166: Tnc[9][13][7]=7.87749e-22;
                   1167: Tnc[7][15][7]=2.75032e-22;
                   1168: Tnc[5][17][7]=4.38402e-23;
                   1169: Tnc[3][19][7]=2.69039e-24;
                   1170: Tnc[19][1][9]=2.26652e-25;
                   1171: Tnc[17][3][9]=1.22716e-23;
                   1172: Tnc[15][5][9]=1.61095e-22;
                   1173: Tnc[13][7][9]=7.87749e-22;
                   1174: Tnc[11][9][9]=1.68852e-21;
                   1175: Tnc[9][11][9]=1.68852e-21;
                   1176: Tnc[7][13][9]=7.87749e-22;
                   1177: Tnc[5][15][9]=1.61095e-22;
                   1178: Tnc[3][17][9]=1.22716e-23;
                   1179: Tnc[1][19][9]=2.26652e-25;
                   1180: Tnc[17][1][11]=6.79956e-25;
                   1181: Tnc[15][3][11]=2.9659e-23;
                   1182: Tnc[13][5][11]=3.03484e-22;
                   1183: Tnc[11][7][11]=1.11059e-21;
                   1184: Tnc[9][9][11]=1.68852e-21;
                   1185: Tnc[7][11][11]=1.11059e-21;
                   1186: Tnc[5][13][11]=3.03484e-22;
                   1187: Tnc[3][15][11]=2.9659e-23;
                   1188: Tnc[1][17][11]=6.79956e-25;
                   1189: Tnc[15][1][13]=1.16564e-24;
                   1190: Tnc[13][3][13]=3.96317e-23;
                   1191: Tnc[11][5][13]=3.03484e-22;
                   1192: Tnc[9][7][13]=7.87749e-22;
                   1193: Tnc[7][9][13]=7.87749e-22;
                   1194: Tnc[5][11][13]=3.03484e-22;
                   1195: Tnc[3][13][13]=3.96317e-23;
                   1196: Tnc[1][15][13]=1.16564e-24;
                   1197: Tnc[13][1][15]=1.16564e-24;
                   1198: Tnc[11][3][15]=2.9659e-23;
                   1199: Tnc[9][5][15]=1.61095e-22;
                   1200: Tnc[7][7][15]=2.75032e-22;
                   1201: Tnc[5][9][15]=1.61095e-22;
                   1202: Tnc[3][11][15]=2.9659e-23;
                   1203: Tnc[1][13][15]=1.16564e-24;
                   1204: Tnc[11][1][17]=6.79956e-25;
                   1205: Tnc[9][3][17]=1.22716e-23;
                   1206: Tnc[7][5][17]=4.38402e-23;
                   1207: Tnc[5][7][17]=4.38402e-23;
                   1208: Tnc[3][9][17]=1.22716e-23;
                   1209: Tnc[1][11][17]=6.79956e-25;
                   1210: Tnc[9][1][19]=2.26652e-25;
                   1211: Tnc[7][3][19]=2.69039e-24;
                   1212: Tnc[5][5][19]=5.62966e-24;
                   1213: Tnc[3][7][19]=2.69039e-24;
                   1214: Tnc[1][9][19]=2.26652e-25;
                   1215: Tnc[18][10][0]=4.53304e-25;
                   1216: Tnc[16][12][0]=1.01993e-24;
                   1217: Tnc[14][14][0]=1.33216e-24;
                   1218: Tnc[12][16][0]=1.01993e-24;
                   1219: Tnc[10][18][0]=4.53304e-25;
                   1220: Tnc[18][8][2]=2.00749e-23;
                   1221: Tnc[16][10][2]=6.58586e-23;
                   1222: Tnc[14][12][2]=1.17729e-22;
                   1223: Tnc[12][14][2]=1.17729e-22;
                   1224: Tnc[10][16][2]=6.58586e-23;
                   1225: Tnc[8][18][2]=2.00749e-23;
                   1226: Tnc[18][6][4]=9.29544e-23;
                   1227: Tnc[16][8][4]=4.87719e-22;
                   1228: Tnc[14][10][4]=1.27117e-21;
                   1229: Tnc[12][12][4]=1.73976e-21;
                   1230: Tnc[10][14][4]=1.27117e-21;
                   1231: Tnc[8][16][4]=4.87719e-22;
                   1232: Tnc[6][18][4]=9.29544e-23;
                   1233: Tnc[18][4][6]=9.29544e-23;
                   1234: Tnc[16][6][6]=9.06636e-22;
                   1235: Tnc[14][8][6]=3.77921e-21;
                   1236: Tnc[12][10][6]=7.54126e-21;
                   1237: Tnc[10][12][6]=7.54126e-21;
                   1238: Tnc[8][14][6]=3.77921e-21;
                   1239: Tnc[6][16][6]=9.06636e-22;
                   1240: Tnc[4][18][6]=9.29544e-23;
                   1241: Tnc[18][2][8]=2.00749e-23;
                   1242: Tnc[16][4][8]=4.87719e-22;
                   1243: Tnc[14][6][8]=3.77921e-21;
                   1244: Tnc[12][8][8]=1.20608e-20;
                   1245: Tnc[10][10][8]=1.75846e-20;
                   1246: Tnc[8][12][8]=1.20608e-20;
                   1247: Tnc[6][14][8]=3.77921e-21;
                   1248: Tnc[4][16][8]=4.87719e-22;
                   1249: Tnc[2][18][8]=2.00749e-23;
                   1250: Tnc[18][0][10]=4.53304e-25;
                   1251: Tnc[16][2][10]=6.58586e-23;
                   1252: Tnc[14][4][10]=1.27117e-21;
                   1253: Tnc[12][6][10]=7.54126e-21;
                   1254: Tnc[10][8][10]=1.75846e-20;
                   1255: Tnc[8][10][10]=1.75846e-20;
                   1256: Tnc[6][12][10]=7.54126e-21;
                   1257: Tnc[4][14][10]=1.27117e-21;
                   1258: Tnc[2][16][10]=6.58586e-23;
                   1259: Tnc[0][18][10]=4.53304e-25;
                   1260: Tnc[16][0][12]=1.01993e-24;
                   1261: Tnc[14][2][12]=1.17729e-22;
                   1262: Tnc[12][4][12]=1.73976e-21;
                   1263: Tnc[10][6][12]=7.54126e-21;
                   1264: Tnc[8][8][12]=1.20608e-20;
                   1265: Tnc[6][10][12]=7.54126e-21;
                   1266: Tnc[4][12][12]=1.73976e-21;
                   1267: Tnc[2][14][12]=1.17729e-22;
                   1268: Tnc[0][16][12]=1.01993e-24;
                   1269: Tnc[14][0][14]=1.33216e-24;
                   1270: Tnc[12][2][14]=1.17729e-22;
                   1271: Tnc[10][4][14]=1.27117e-21;
                   1272: Tnc[8][6][14]=3.77921e-21;
                   1273: Tnc[6][8][14]=3.77921e-21;
                   1274: Tnc[4][10][14]=1.27117e-21;
                   1275: Tnc[2][12][14]=1.17729e-22;
                   1276: Tnc[0][14][14]=1.33216e-24;
                   1277: Tnc[12][0][16]=1.01993e-24;
                   1278: Tnc[10][2][16]=6.58586e-23;
                   1279: Tnc[8][4][16]=4.87719e-22;
                   1280: Tnc[6][6][16]=9.06636e-22;
                   1281: Tnc[4][8][16]=4.87719e-22;
                   1282: Tnc[2][10][16]=6.58586e-23;
                   1283: Tnc[0][12][16]=1.01993e-24;
                   1284: Tnc[10][0][18]=4.53304e-25;
                   1285: Tnc[8][2][18]=2.00749e-23;
                   1286: Tnc[6][4][18]=9.29544e-23;
                   1287: Tnc[4][6][18]=9.29544e-23;
                   1288: Tnc[2][8][18]=2.00749e-23;
                   1289: Tnc[0][10][18]=4.53304e-25;
                   1290: Tnc[19][7][1]=1.87797e-23;
                   1291: Tnc[17][9][1]=8.45088e-23;
                   1292: Tnc[15][11][1]=2.02821e-22;
                   1293: Tnc[13][13][1]=2.70428e-22;
                   1294: Tnc[11][15][1]=2.02821e-22;
                   1295: Tnc[9][17][1]=8.45088e-23;
                   1296: Tnc[7][19][1]=1.87797e-23;
                   1297: Tnc[19][5][3]=1.30407e-22;
                   1298: Tnc[17][7][3]=1.0011e-21;
                   1299: Tnc[15][9][3]=3.64905e-21;
                   1300: Tnc[13][11][3]=6.84911e-21;
                   1301: Tnc[11][13][3]=6.84911e-21;
                   1302: Tnc[9][15][3]=3.64905e-21;
                   1303: Tnc[7][17][3]=1.0011e-21;
                   1304: Tnc[5][19][3]=1.30407e-22;
                   1305: Tnc[19][3][5]=1.30407e-22;
                   1306: Tnc[17][5][5]=2.09345e-21;
                   1307: Tnc[15][7][5]=1.3017e-20;
                   1308: Tnc[13][9][5]=3.71059e-20;
                   1309: Tnc[11][11][5]=5.22342e-20;
                   1310: Tnc[9][13][5]=3.71059e-20;
                   1311: Tnc[7][15][5]=1.3017e-20;
                   1312: Tnc[5][17][5]=2.09345e-21;
                   1313: Tnc[3][19][5]=1.30407e-22;
                   1314: Tnc[19][1][7]=1.87797e-23;
                   1315: Tnc[17][3][7]=1.0011e-21;
                   1316: Tnc[15][5][7]=1.3017e-20;
                   1317: Tnc[13][7][7]=6.32972e-20;
                   1318: Tnc[11][9][7]=1.35323e-19;
                   1319: Tnc[9][11][7]=1.35323e-19;
                   1320: Tnc[7][13][7]=6.32972e-20;
                   1321: Tnc[5][15][7]=1.3017e-20;
                   1322: Tnc[3][17][7]=1.0011e-21;
                   1323: Tnc[1][19][7]=1.87797e-23;
                   1324: Tnc[17][1][9]=8.45088e-23;
                   1325: Tnc[15][3][9]=3.64905e-21;
                   1326: Tnc[13][5][9]=3.71059e-20;
                   1327: Tnc[11][7][9]=1.35323e-19;
                   1328: Tnc[9][9][9]=2.05517e-19;
                   1329: Tnc[7][11][9]=1.35323e-19;
                   1330: Tnc[5][13][9]=3.71059e-20;
                   1331: Tnc[3][15][9]=3.64905e-21;
                   1332: Tnc[1][17][9]=8.45088e-23;
                   1333: Tnc[15][1][11]=2.02821e-22;
                   1334: Tnc[13][3][11]=6.84911e-21;
                   1335: Tnc[11][5][11]=5.22342e-20;
                   1336: Tnc[9][7][11]=1.35323e-19;
                   1337: Tnc[7][9][11]=1.35323e-19;
                   1338: Tnc[5][11][11]=5.22342e-20;
                   1339: Tnc[3][13][11]=6.84911e-21;
                   1340: Tnc[1][15][11]=2.02821e-22;
                   1341: Tnc[13][1][13]=2.70428e-22;
                   1342: Tnc[11][3][13]=6.84911e-21;
                   1343: Tnc[9][5][13]=3.71059e-20;
                   1344: Tnc[7][7][13]=6.32972e-20;
                   1345: Tnc[5][9][13]=3.71059e-20;
                   1346: Tnc[3][11][13]=6.84911e-21;
                   1347: Tnc[1][13][13]=2.70428e-22;
                   1348: Tnc[11][1][15]=2.02821e-22;
                   1349: Tnc[9][3][15]=3.64905e-21;
                   1350: Tnc[7][5][15]=1.3017e-20;
                   1351: Tnc[5][7][15]=1.3017e-20;
                   1352: Tnc[3][9][15]=3.64905e-21;
                   1353: Tnc[1][11][15]=2.02821e-22;
                   1354: Tnc[9][1][17]=8.45088e-23;
                   1355: Tnc[7][3][17]=1.0011e-21;
                   1356: Tnc[5][5][17]=2.09345e-21;
                   1357: Tnc[3][7][17]=1.0011e-21;
                   1358: Tnc[1][9][17]=8.45088e-23;
                   1359: Tnc[7][1][19]=1.87797e-23;
                   1360: Tnc[5][3][19]=1.30407e-22;
                   1361: Tnc[3][5][19]=1.30407e-22;
                   1362: Tnc[1][7][19]=1.87797e-23;
                   1363: Tnc[18][8][0]=4.69493e-23;
                   1364: Tnc[16][10][0]=1.52116e-22;
                   1365: Tnc[14][12][0]=2.70428e-22;
                   1366: Tnc[12][14][0]=2.70428e-22;
                   1367: Tnc[10][16][0]=1.52116e-22;
                   1368: Tnc[8][18][0]=4.69493e-23;
                   1369: Tnc[18][6][2]=1.29725e-21;
                   1370: Tnc[16][8][2]=6.7152e-21;
                   1371: Tnc[14][10][2]=1.73802e-20;
                   1372: Tnc[12][12][2]=2.37353e-20;
                   1373: Tnc[10][14][2]=1.73802e-20;
                   1374: Tnc[8][16][2]=6.7152e-21;
                   1375: Tnc[6][18][2]=1.29725e-21;
                   1376: Tnc[18][4][4]=3.22906e-21;
                   1377: Tnc[16][6][4]=3.10472e-20;
                   1378: Tnc[14][8][4]=1.28375e-19;
                   1379: Tnc[12][10][4]=2.55227e-19;
                   1380: Tnc[10][12][4]=2.55227e-19;
                   1381: Tnc[8][14][4]=1.28375e-19;
                   1382: Tnc[6][16][4]=3.10472e-20;
                   1383: Tnc[4][18][4]=3.22906e-21;
                   1384: Tnc[18][2][6]=1.29725e-21;
                   1385: Tnc[16][4][6]=3.10472e-20;
                   1386: Tnc[14][6][6]=2.38442e-19;
                   1387: Tnc[12][8][6]=7.57325e-19;
                   1388: Tnc[10][10][6]=1.10252e-18;
                   1389: Tnc[8][12][6]=7.57325e-19;
                   1390: Tnc[6][14][6]=2.38442e-19;
                   1391: Tnc[4][16][6]=3.10472e-20;
                   1392: Tnc[2][18][6]=1.29725e-21;
                   1393: Tnc[18][0][8]=4.69493e-23;
                   1394: Tnc[16][2][8]=6.7152e-21;
                   1395: Tnc[14][4][8]=1.28375e-19;
                   1396: Tnc[12][6][8]=7.57325e-19;
                   1397: Tnc[10][8][8]=1.76133e-18;
                   1398: Tnc[8][10][8]=1.76133e-18;
                   1399: Tnc[6][12][8]=7.57325e-19;
                   1400: Tnc[4][14][8]=1.28375e-19;
                   1401: Tnc[2][16][8]=6.7152e-21;
                   1402: Tnc[0][18][8]=4.69493e-23;
                   1403: Tnc[16][0][10]=1.52116e-22;
                   1404: Tnc[14][2][10]=1.73802e-20;
                   1405: Tnc[12][4][10]=2.55227e-19;
                   1406: Tnc[10][6][10]=1.10252e-18;
                   1407: Tnc[8][8][10]=1.76133e-18;
                   1408: Tnc[6][10][10]=1.10252e-18;
                   1409: Tnc[4][12][10]=2.55227e-19;
                   1410: Tnc[2][14][10]=1.73802e-20;
                   1411: Tnc[0][16][10]=1.52116e-22;
                   1412: Tnc[14][0][12]=2.70428e-22;
                   1413: Tnc[12][2][12]=2.37353e-20;
                   1414: Tnc[10][4][12]=2.55227e-19;
                   1415: Tnc[8][6][12]=7.57325e-19;
                   1416: Tnc[6][8][12]=7.57325e-19;
                   1417: Tnc[4][10][12]=2.55227e-19;
                   1418: Tnc[2][12][12]=2.37353e-20;
                   1419: Tnc[0][14][12]=2.70428e-22;
                   1420: Tnc[12][0][14]=2.70428e-22;
                   1421: Tnc[10][2][14]=1.73802e-20;
                   1422: Tnc[8][4][14]=1.28375e-19;
                   1423: Tnc[6][6][14]=2.38442e-19;
                   1424: Tnc[4][8][14]=1.28375e-19;
                   1425: Tnc[2][10][14]=1.73802e-20;
                   1426: Tnc[0][12][14]=2.70428e-22;
                   1427: Tnc[10][0][16]=1.52116e-22;
                   1428: Tnc[8][2][16]=6.7152e-21;
                   1429: Tnc[6][4][16]=3.10472e-20;
                   1430: Tnc[4][6][16]=3.10472e-20;
                   1431: Tnc[2][8][16]=6.7152e-21;
                   1432: Tnc[0][10][16]=1.52116e-22;
                   1433: Tnc[8][0][18]=4.69493e-23;
                   1434: Tnc[6][2][18]=1.29725e-21;
                   1435: Tnc[4][4][18]=3.22906e-21;
                   1436: Tnc[2][6][18]=1.29725e-21;
                   1437: Tnc[0][8][18]=4.69493e-23;
                   1438: Tnc[19][5][1]=9.36097e-22;
                   1439: Tnc[17][7][1]=7.02073e-21;
                   1440: Tnc[15][9][1]=2.52746e-20;
                   1441: Tnc[13][11][1]=4.71793e-20;
                   1442: Tnc[11][13][1]=4.71793e-20;
                   1443: Tnc[9][15][1]=2.52746e-20;
                   1444: Tnc[7][17][1]=7.02073e-21;
                   1445: Tnc[5][19][1]=9.36097e-22;
                   1446: Tnc[19][3][3]=3.10614e-21;
                   1447: Tnc[17][5][3]=4.86771e-20;
                   1448: Tnc[15][7][3]=2.98615e-19;
                   1449: Tnc[13][9][3]=8.45296e-19;
                   1450: Tnc[11][11][3]=1.18735e-18;
                   1451: Tnc[9][13][3]=8.45296e-19;
                   1452: Tnc[7][15][3]=2.98615e-19;
                   1453: Tnc[5][17][3]=4.86771e-20;
                   1454: Tnc[3][19][3]=3.10614e-21;
                   1455: Tnc[19][1][5]=9.36097e-22;
                   1456: Tnc[17][3][5]=4.86771e-20;
                   1457: Tnc[15][5][5]=6.23917e-19;
                   1458: Tnc[13][7][5]=3.00947e-18;
                   1459: Tnc[11][9][5]=6.41031e-18;
                   1460: Tnc[9][11][5]=6.41031e-18;
                   1461: Tnc[7][13][5]=3.00947e-18;
                   1462: Tnc[5][15][5]=6.23917e-19;
                   1463: Tnc[3][17][5]=4.86771e-20;
                   1464: Tnc[1][19][5]=9.36097e-22;
                   1465: Tnc[17][1][7]=7.02073e-21;
                   1466: Tnc[15][3][7]=2.98615e-19;
                   1467: Tnc[13][5][7]=3.00947e-18;
                   1468: Tnc[11][7][7]=1.0923e-17;
                   1469: Tnc[9][9][7]=1.65638e-17;
                   1470: Tnc[7][11][7]=1.0923e-17;
                   1471: Tnc[5][13][7]=3.00947e-18;
                   1472: Tnc[3][15][7]=2.98615e-19;
                   1473: Tnc[1][17][7]=7.02073e-21;
                   1474: Tnc[15][1][9]=2.52746e-20;
                   1475: Tnc[13][3][9]=8.45296e-19;
                   1476: Tnc[11][5][9]=6.41031e-18;
                   1477: Tnc[9][7][9]=1.65638e-17;
                   1478: Tnc[7][9][9]=1.65638e-17;
                   1479: Tnc[5][11][9]=6.41031e-18;
                   1480: Tnc[3][13][9]=8.45296e-19;
                   1481: Tnc[1][15][9]=2.52746e-20;
                   1482: Tnc[13][1][11]=4.71793e-20;
                   1483: Tnc[11][3][11]=1.18735e-18;
                   1484: Tnc[9][5][11]=6.41031e-18;
                   1485: Tnc[7][7][11]=1.0923e-17;
                   1486: Tnc[5][9][11]=6.41031e-18;
                   1487: Tnc[3][11][11]=1.18735e-18;
                   1488: Tnc[1][13][11]=4.71793e-20;
                   1489: Tnc[11][1][13]=4.71793e-20;
                   1490: Tnc[9][3][13]=8.45296e-19;
                   1491: Tnc[7][5][13]=3.00947e-18;
                   1492: Tnc[5][7][13]=3.00947e-18;
                   1493: Tnc[3][9][13]=8.45296e-19;
                   1494: Tnc[1][11][13]=4.71793e-20;
                   1495: Tnc[9][1][15]=2.52746e-20;
                   1496: Tnc[7][3][15]=2.98615e-19;
                   1497: Tnc[5][5][15]=6.23917e-19;
                   1498: Tnc[3][7][15]=2.98615e-19;
                   1499: Tnc[1][9][15]=2.52746e-20;
                   1500: Tnc[7][1][17]=7.02073e-21;
                   1501: Tnc[5][3][17]=4.86771e-20;
                   1502: Tnc[3][5][17]=4.86771e-20;
                   1503: Tnc[1][7][17]=7.02073e-21;
                   1504: Tnc[5][1][19]=9.36097e-22;
                   1505: Tnc[3][3][19]=3.10614e-21;
                   1506: Tnc[1][5][19]=9.36097e-22;
                   1507: Tnc[18][6][0]=3.12032e-21;
                   1508: Tnc[16][8][0]=1.57966e-20;
                   1509: Tnc[14][10][0]=4.04394e-20;
                   1510: Tnc[12][12][0]=5.50425e-20;
                   1511: Tnc[10][14][0]=4.04394e-20;
                   1512: Tnc[8][16][0]=1.57966e-20;
                   1513: Tnc[6][18][0]=3.12032e-21;
                   1514: Tnc[18][4][2]=4.63368e-20;
                   1515: Tnc[16][6][2]=4.35285e-19;
                   1516: Tnc[14][8][2]=1.77765e-18;
                   1517: Tnc[12][10][2]=3.51486e-18;
                   1518: Tnc[10][12][2]=3.51486e-18;
                   1519: Tnc[8][14][2]=1.77765e-18;
                   1520: Tnc[6][16][2]=4.35285e-19;
                   1521: Tnc[4][18][2]=4.63368e-20;
                   1522: Tnc[18][2][4]=4.63368e-20;
                   1523: Tnc[16][4][4]=1.08255e-18;
                   1524: Tnc[14][6][4]=8.20251e-18;
                   1525: Tnc[12][8][4]=2.58709e-17;
                   1526: Tnc[10][10][4]=3.75812e-17;
                   1527: Tnc[8][12][4]=2.58709e-17;
                   1528: Tnc[6][14][4]=8.20251e-18;
                   1529: Tnc[4][16][4]=1.08255e-18;
                   1530: Tnc[2][18][4]=4.63368e-20;
                   1531: Tnc[18][0][6]=3.12032e-21;
                   1532: Tnc[16][2][6]=4.35285e-19;
                   1533: Tnc[14][4][6]=8.20251e-18;
                   1534: Tnc[12][6][6]=4.79984e-17;
                   1535: Tnc[10][8][6]=1.1122e-16;
                   1536: Tnc[8][10][6]=1.1122e-16;
                   1537: Tnc[6][12][6]=4.79984e-17;
                   1538: Tnc[4][14][6]=8.20251e-18;
                   1539: Tnc[2][16][6]=4.35285e-19;
                   1540: Tnc[0][18][6]=3.12032e-21;
                   1541: Tnc[16][0][8]=1.57966e-20;
                   1542: Tnc[14][2][8]=1.77765e-18;
                   1543: Tnc[12][4][8]=2.58709e-17;
                   1544: Tnc[10][6][8]=1.1122e-16;
                   1545: Tnc[8][8][8]=1.7741e-16;
                   1546: Tnc[6][10][8]=1.1122e-16;
                   1547: Tnc[4][12][8]=2.58709e-17;
                   1548: Tnc[2][14][8]=1.77765e-18;
                   1549: Tnc[0][16][8]=1.57966e-20;
                   1550: Tnc[14][0][10]=4.04394e-20;
                   1551: Tnc[12][2][10]=3.51486e-18;
                   1552: Tnc[10][4][10]=3.75812e-17;
                   1553: Tnc[8][6][10]=1.1122e-16;
                   1554: Tnc[6][8][10]=1.1122e-16;
                   1555: Tnc[4][10][10]=3.75812e-17;
                   1556: Tnc[2][12][10]=3.51486e-18;
                   1557: Tnc[0][14][10]=4.04394e-20;
                   1558: Tnc[12][0][12]=5.50425e-20;
                   1559: Tnc[10][2][12]=3.51486e-18;
                   1560: Tnc[8][4][12]=2.58709e-17;
                   1561: Tnc[6][6][12]=4.79984e-17;
                   1562: Tnc[4][8][12]=2.58709e-17;
                   1563: Tnc[2][10][12]=3.51486e-18;
                   1564: Tnc[0][12][12]=5.50425e-20;
                   1565: Tnc[10][0][14]=4.04394e-20;
                   1566: Tnc[8][2][14]=1.77765e-18;
                   1567: Tnc[6][4][14]=8.20251e-18;
                   1568: Tnc[4][6][14]=8.20251e-18;
                   1569: Tnc[2][8][14]=1.77765e-18;
                   1570: Tnc[0][10][14]=4.04394e-20;
                   1571: Tnc[8][0][16]=1.57966e-20;
                   1572: Tnc[6][2][16]=4.35285e-19;
                   1573: Tnc[4][4][16]=1.08255e-18;
                   1574: Tnc[2][6][16]=4.35285e-19;
                   1575: Tnc[0][8][16]=1.57966e-20;
                   1576: Tnc[6][0][18]=3.12032e-21;
                   1577: Tnc[4][2][18]=4.63368e-20;
                   1578: Tnc[2][4][18]=4.63368e-20;
                   1579: Tnc[0][6][18]=3.12032e-21;
                   1580: Tnc[19][3][1]=2.34024e-20;
                   1581: Tnc[17][5][1]=3.51036e-19;
                   1582: Tnc[15][7][1]=2.10622e-18;
                   1583: Tnc[13][9][1]=5.89741e-18;
                   1584: Tnc[11][11][1]=8.25638e-18;
                   1585: Tnc[9][13][1]=5.89741e-18;
                   1586: Tnc[7][15][1]=2.10622e-18;
                   1587: Tnc[5][17][1]=3.51036e-19;
                   1588: Tnc[3][19][1]=2.34024e-20;
                   1589: Tnc[19][1][3]=2.34024e-20;
                   1590: Tnc[17][3][3]=1.16374e-18;
                   1591: Tnc[15][5][3]=1.45733e-17;
                   1592: Tnc[13][7][3]=6.94286e-17;
                   1593: Tnc[11][9][3]=1.47078e-16;
                   1594: Tnc[9][11][3]=1.47078e-16;
                   1595: Tnc[7][13][3]=6.94286e-17;
                   1596: Tnc[5][15][3]=1.45733e-17;
                   1597: Tnc[3][17][3]=1.16374e-18;
                   1598: Tnc[1][19][3]=2.34024e-20;
                   1599: Tnc[17][1][5]=3.51036e-19;
                   1600: Tnc[15][3][5]=1.45733e-17;
                   1601: Tnc[13][5][5]=1.44898e-16;
                   1602: Tnc[11][7][5]=5.22243e-16;
                   1603: Tnc[9][9][5]=7.90217e-16;
                   1604: Tnc[7][11][5]=5.22243e-16;
                   1605: Tnc[5][13][5]=1.44898e-16;
                   1606: Tnc[3][15][5]=1.45733e-17;
                   1607: Tnc[1][17][5]=3.51036e-19;
                   1608: Tnc[15][1][7]=2.10622e-18;
                   1609: Tnc[13][3][7]=6.94286e-17;
                   1610: Tnc[11][5][7]=5.22243e-16;
                   1611: Tnc[9][7][7]=1.34445e-15;
                   1612: Tnc[7][9][7]=1.34445e-15;
                   1613: Tnc[5][11][7]=5.22243e-16;
                   1614: Tnc[3][13][7]=6.94286e-17;
                   1615: Tnc[1][15][7]=2.10622e-18;
                   1616: Tnc[13][1][9]=5.89741e-18;
                   1617: Tnc[11][3][9]=1.47078e-16;
                   1618: Tnc[9][5][9]=7.90217e-16;
                   1619: Tnc[7][7][9]=1.34445e-15;
                   1620: Tnc[5][9][9]=7.90217e-16;
                   1621: Tnc[3][11][9]=1.47078e-16;
                   1622: Tnc[1][13][9]=5.89741e-18;
                   1623: Tnc[11][1][11]=8.25638e-18;
                   1624: Tnc[9][3][11]=1.47078e-16;
                   1625: Tnc[7][5][11]=5.22243e-16;
                   1626: Tnc[5][7][11]=5.22243e-16;
                   1627: Tnc[3][9][11]=1.47078e-16;
                   1628: Tnc[1][11][11]=8.25638e-18;
                   1629: Tnc[9][1][13]=5.89741e-18;
                   1630: Tnc[7][3][13]=6.94286e-17;
                   1631: Tnc[5][5][13]=1.44898e-16;
                   1632: Tnc[3][7][13]=6.94286e-17;
                   1633: Tnc[1][9][13]=5.89741e-18;
                   1634: Tnc[7][1][15]=2.10622e-18;
                   1635: Tnc[5][3][15]=1.45733e-17;
                   1636: Tnc[3][5][15]=1.45733e-17;
                   1637: Tnc[1][7][15]=2.10622e-18;
                   1638: Tnc[5][1][17]=3.51036e-19;
                   1639: Tnc[3][3][17]=1.16374e-18;
                   1640: Tnc[1][5][17]=3.51036e-19;
                   1641: Tnc[3][1][19]=2.34024e-20;
                   1642: Tnc[1][3][19]=2.34024e-20;
                   1643: Tnc[18][4][0]=1.17012e-19;
                   1644: Tnc[16][6][0]=1.05311e-18;
                   1645: Tnc[14][8][0]=4.21244e-18;
                   1646: Tnc[12][10][0]=8.25638e-18;
                   1647: Tnc[10][12][0]=8.25638e-18;
                   1648: Tnc[8][14][0]=4.21244e-18;
                   1649: Tnc[6][16][0]=1.05311e-18;
                   1650: Tnc[4][18][0]=1.17012e-19;
                   1651: Tnc[18][2][2]=6.97818e-19;
                   1652: Tnc[16][4][2]=1.56052e-17;
                   1653: Tnc[14][6][2]=1.15651e-16;
                   1654: Tnc[12][8][2]=3.60814e-16;
                   1655: Tnc[10][10][2]=5.22403e-16;
                   1656: Tnc[8][12][2]=3.60814e-16;
                   1657: Tnc[6][14][2]=1.15651e-16;
                   1658: Tnc[4][16][2]=1.56052e-17;
                   1659: Tnc[2][18][2]=6.97818e-19;
                   1660: Tnc[18][0][4]=1.17012e-19;
                   1661: Tnc[16][2][4]=1.56052e-17;
                   1662: Tnc[14][4][4]=2.87288e-16;
                   1663: Tnc[12][6][4]=1.66039e-15;
                   1664: Tnc[10][8][4]=3.82635e-15;
                   1665: Tnc[8][10][4]=3.82635e-15;
                   1666: Tnc[6][12][4]=1.66039e-15;
                   1667: Tnc[4][14][4]=2.87288e-16;
                   1668: Tnc[2][16][4]=1.56052e-17;
                   1669: Tnc[0][18][4]=1.17012e-19;
                   1670: Tnc[16][0][6]=1.05311e-18;
                   1671: Tnc[14][2][6]=1.15651e-16;
                   1672: Tnc[12][4][6]=1.66039e-15;
                   1673: Tnc[10][6][6]=7.08812e-15;
                   1674: Tnc[8][8][6]=1.12818e-14;
                   1675: Tnc[6][10][6]=7.08812e-15;
                   1676: Tnc[4][12][6]=1.66039e-15;
                   1677: Tnc[2][14][6]=1.15651e-16;
                   1678: Tnc[0][16][6]=1.05311e-18;
                   1679: Tnc[14][0][8]=4.21244e-18;
                   1680: Tnc[12][2][8]=3.60814e-16;
                   1681: Tnc[10][4][8]=3.82635e-15;
                   1682: Tnc[8][6][8]=1.12818e-14;
                   1683: Tnc[6][8][8]=1.12818e-14;
                   1684: Tnc[4][10][8]=3.82635e-15;
                   1685: Tnc[2][12][8]=3.60814e-16;
                   1686: Tnc[0][14][8]=4.21244e-18;
                   1687: Tnc[12][0][10]=8.25638e-18;
                   1688: Tnc[10][2][10]=5.22403e-16;
                   1689: Tnc[8][4][10]=3.82635e-15;
                   1690: Tnc[6][6][10]=7.08812e-15;
                   1691: Tnc[4][8][10]=3.82635e-15;
                   1692: Tnc[2][10][10]=5.22403e-16;
                   1693: Tnc[0][12][10]=8.25638e-18;
                   1694: Tnc[10][0][12]=8.25638e-18;
                   1695: Tnc[8][2][12]=3.60814e-16;
                   1696: Tnc[6][4][12]=1.66039e-15;
                   1697: Tnc[4][6][12]=1.66039e-15;
                   1698: Tnc[2][8][12]=3.60814e-16;
                   1699: Tnc[0][10][12]=8.25638e-18;
                   1700: Tnc[8][0][14]=4.21244e-18;
                   1701: Tnc[6][2][14]=1.15651e-16;
                   1702: Tnc[4][4][14]=2.87288e-16;
                   1703: Tnc[2][6][14]=1.15651e-16;
                   1704: Tnc[0][8][14]=4.21244e-18;
                   1705: Tnc[6][0][16]=1.05311e-18;
                   1706: Tnc[4][2][16]=1.56052e-17;
                   1707: Tnc[2][4][16]=1.56052e-17;
                   1708: Tnc[0][6][16]=1.05311e-18;
                   1709: Tnc[4][0][18]=1.17012e-19;
                   1710: Tnc[2][2][18]=6.97818e-19;
                   1711: Tnc[0][4][18]=1.17012e-19;
                   1712: Tnc[19][1][1]=1.95729e-19;
                   1713: Tnc[17][3][1]=8.80782e-18;
                   1714: Tnc[15][5][1]=1.05694e-16;
                   1715: Tnc[13][7][1]=4.93238e-16;
                   1716: Tnc[11][9][1]=1.0358e-15;
                   1717: Tnc[9][11][1]=1.0358e-15;
                   1718: Tnc[7][13][1]=4.93238e-16;
                   1719: Tnc[5][15][1]=1.05694e-16;
                   1720: Tnc[3][17][1]=8.80782e-18;
                   1721: Tnc[1][19][1]=1.95729e-19;
                   1722: Tnc[17][1][3]=8.80782e-18;
                   1723: Tnc[15][3][3]=3.49964e-16;
                   1724: Tnc[13][5][3]=3.40334e-15;
                   1725: Tnc[11][7][3]=1.21337e-14;
                   1726: Tnc[9][9][3]=1.82991e-14;
                   1727: Tnc[7][11][3]=1.21337e-14;
                   1728: Tnc[5][13][3]=3.40334e-15;
                   1729: Tnc[3][15][3]=3.49964e-16;
                   1730: Tnc[1][17][3]=8.80782e-18;
                   1731: Tnc[15][1][5]=1.05694e-16;
                   1732: Tnc[13][3][5]=3.40334e-15;
                   1733: Tnc[11][5][5]=2.52834e-14;
                   1734: Tnc[9][7][5]=6.47326e-14;
                   1735: Tnc[7][9][5]=6.47326e-14;
                   1736: Tnc[5][11][5]=2.52834e-14;
                   1737: Tnc[3][13][5]=3.40334e-15;
                   1738: Tnc[1][15][5]=1.05694e-16;
                   1739: Tnc[13][1][7]=4.93238e-16;
                   1740: Tnc[11][3][7]=1.21337e-14;
                   1741: Tnc[9][5][7]=6.47326e-14;
                   1742: Tnc[7][7][7]=1.09893e-13;
                   1743: Tnc[5][9][7]=6.47326e-14;
                   1744: Tnc[3][11][7]=1.21337e-14;
                   1745: Tnc[1][13][7]=4.93238e-16;
                   1746: Tnc[11][1][9]=1.0358e-15;
                   1747: Tnc[9][3][9]=1.82991e-14;
                   1748: Tnc[7][5][9]=6.47326e-14;
                   1749: Tnc[5][7][9]=6.47326e-14;
                   1750: Tnc[3][9][9]=1.82991e-14;
                   1751: Tnc[1][11][9]=1.0358e-15;
                   1752: Tnc[9][1][11]=1.0358e-15;
                   1753: Tnc[7][3][11]=1.21337e-14;
                   1754: Tnc[5][5][11]=2.52834e-14;
                   1755: Tnc[3][7][11]=1.21337e-14;
                   1756: Tnc[1][9][11]=1.0358e-15;
                   1757: Tnc[7][1][13]=4.93238e-16;
                   1758: Tnc[5][3][13]=3.40334e-15;
                   1759: Tnc[3][5][13]=3.40334e-15;
                   1760: Tnc[1][7][13]=4.93238e-16;
                   1761: Tnc[5][1][15]=1.05694e-16;
                   1762: Tnc[3][3][15]=3.49964e-16;
                   1763: Tnc[1][5][15]=1.05694e-16;
                   1764: Tnc[3][1][17]=8.80782e-18;
                   1765: Tnc[1][3][17]=8.80782e-18;
                   1766: Tnc[1][1][19]=1.95729e-19;
                   1767: Tnc[18][2][0]=1.95729e-18;
                   1768: Tnc[16][4][0]=3.96352e-17;
                   1769: Tnc[14][6][0]=2.8185e-16;
                   1770: Tnc[12][8][0]=8.63167e-16;
                   1771: Tnc[10][10][0]=1.24296e-15;
                   1772: Tnc[8][12][0]=8.63167e-16;
                   1773: Tnc[6][14][0]=2.8185e-16;
                   1774: Tnc[4][16][0]=3.96352e-17;
                   1775: Tnc[2][18][0]=1.95729e-18;
                   1776: Tnc[18][0][2]=1.95729e-18;
                   1777: Tnc[16][2][2]=2.3605e-16;
                   1778: Tnc[14][4][2]=4.16434e-15;
                   1779: Tnc[12][6][2]=2.35768e-14;
                   1780: Tnc[10][8][2]=5.38616e-14;
                   1781: Tnc[8][10][2]=5.38616e-14;
                   1782: Tnc[6][12][2]=2.35768e-14;
                   1783: Tnc[4][14][2]=4.16434e-15;
                   1784: Tnc[2][16][2]=2.3605e-16;
                   1785: Tnc[0][18][2]=1.95729e-18;
                   1786: Tnc[16][0][4]=3.96352e-17;
                   1787: Tnc[14][2][4]=4.16434e-15;
                   1788: Tnc[12][4][4]=5.84734e-14;
                   1789: Tnc[10][6][4]=2.46915e-13;
                   1790: Tnc[8][8][4]=3.91705e-13;
                   1791: Tnc[6][10][4]=2.46915e-13;
                   1792: Tnc[4][12][4]=5.84734e-14;
                   1793: Tnc[2][14][4]=4.16434e-15;
                   1794: Tnc[0][16][4]=3.96352e-17;
                   1795: Tnc[14][0][6]=2.8185e-16;
                   1796: Tnc[12][2][6]=2.35768e-14;
                   1797: Tnc[10][4][6]=2.46915e-13;
                   1798: Tnc[8][6][6]=7.24011e-13;
                   1799: Tnc[6][8][6]=7.24011e-13;
                   1800: Tnc[4][10][6]=2.46915e-13;
                   1801: Tnc[2][12][6]=2.35768e-14;
                   1802: Tnc[0][14][6]=2.8185e-16;
                   1803: Tnc[12][0][8]=8.63167e-16;
                   1804: Tnc[10][2][8]=5.38616e-14;
                   1805: Tnc[8][4][8]=3.91705e-13;
                   1806: Tnc[6][6][8]=7.24011e-13;
                   1807: Tnc[4][8][8]=3.91705e-13;
                   1808: Tnc[2][10][8]=5.38616e-14;
                   1809: Tnc[0][12][8]=8.63167e-16;
                   1810: Tnc[10][0][10]=1.24296e-15;
                   1811: Tnc[8][2][10]=5.38616e-14;
                   1812: Tnc[6][4][10]=2.46915e-13;
                   1813: Tnc[4][6][10]=2.46915e-13;
                   1814: Tnc[2][8][10]=5.38616e-14;
                   1815: Tnc[0][10][10]=1.24296e-15;
                   1816: Tnc[8][0][12]=8.63167e-16;
                   1817: Tnc[6][2][12]=2.35768e-14;
                   1818: Tnc[4][4][12]=5.84734e-14;
                   1819: Tnc[2][6][12]=2.35768e-14;
                   1820: Tnc[0][8][12]=8.63167e-16;
                   1821: Tnc[6][0][14]=2.8185e-16;
                   1822: Tnc[4][2][14]=4.16434e-15;
                   1823: Tnc[2][4][14]=4.16434e-15;
                   1824: Tnc[0][6][14]=2.8185e-16;
                   1825: Tnc[4][0][16]=3.96352e-17;
                   1826: Tnc[2][2][16]=2.3605e-16;
                   1827: Tnc[0][4][16]=3.96352e-17;
                   1828: Tnc[2][0][18]=1.95729e-18;
                   1829: Tnc[0][2][18]=1.95729e-18;
                   1830: Tnc[17][1][1]=7.39857e-17;
                   1831: Tnc[15][3][1]=2.66349e-15;
                   1832: Tnc[13][5][1]=2.48592e-14;
                   1833: Tnc[11][7][1]=8.70072e-14;
                   1834: Tnc[9][9][1]=1.30511e-13;
                   1835: Tnc[7][11][1]=8.70072e-14;
                   1836: Tnc[5][13][1]=2.48592e-14;
                   1837: Tnc[3][15][1]=2.66349e-15;
                   1838: Tnc[1][17][1]=7.39857e-17;
                   1839: Tnc[15][1][3]=2.66349e-15;
                   1840: Tnc[13][3][3]=8.21735e-14;
                   1841: Tnc[11][5][3]=5.98002e-13;
                   1842: Tnc[9][7][3]=1.51779e-12;
                   1843: Tnc[7][9][3]=1.51779e-12;
                   1844: Tnc[5][11][3]=5.98002e-13;
                   1845: Tnc[3][13][3]=8.21735e-14;
                   1846: Tnc[1][15][3]=2.66349e-15;
                   1847: Tnc[13][1][5]=2.48592e-14;
                   1848: Tnc[11][3][5]=5.98002e-13;
                   1849: Tnc[9][5][5]=3.1556e-12;
                   1850: Tnc[7][7][5]=5.33934e-12;
                   1851: Tnc[5][9][5]=3.1556e-12;
                   1852: Tnc[3][11][5]=5.98002e-13;
                   1853: Tnc[1][13][5]=2.48592e-14;
                   1854: Tnc[11][1][7]=8.70072e-14;
                   1855: Tnc[9][3][7]=1.51779e-12;
                   1856: Tnc[7][5][7]=5.33934e-12;
                   1857: Tnc[5][7][7]=5.33934e-12;
                   1858: Tnc[3][9][7]=1.51779e-12;
                   1859: Tnc[1][11][7]=8.70072e-14;
                   1860: Tnc[9][1][9]=1.30511e-13;
                   1861: Tnc[7][3][9]=1.51779e-12;
                   1862: Tnc[5][5][9]=3.1556e-12;
                   1863: Tnc[3][7][9]=1.51779e-12;
                   1864: Tnc[1][9][9]=1.30511e-13;
                   1865: Tnc[7][1][11]=8.70072e-14;
                   1866: Tnc[5][3][11]=5.98002e-13;
                   1867: Tnc[3][5][11]=5.98002e-13;
                   1868: Tnc[1][7][11]=8.70072e-14;
                   1869: Tnc[5][1][13]=2.48592e-14;
                   1870: Tnc[3][3][13]=8.21735e-14;
                   1871: Tnc[1][5][13]=2.48592e-14;
                   1872: Tnc[3][1][15]=2.66349e-15;
                   1873: Tnc[1][3][15]=2.66349e-15;
                   1874: Tnc[1][1][17]=7.39857e-17;
                   1875: Tnc[18][0][0]=8.22064e-18;
                   1876: Tnc[16][2][0]=6.65871e-16;
                   1877: Tnc[14][4][0]=1.06539e-14;
                   1878: Tnc[12][6][0]=5.80048e-14;
                   1879: Tnc[10][8][0]=1.30511e-13;
                   1880: Tnc[8][10][0]=1.30511e-13;
                   1881: Tnc[6][12][0]=5.80048e-14;
                   1882: Tnc[4][14][0]=1.06539e-14;
                   1883: Tnc[2][16][0]=6.65871e-16;
                   1884: Tnc[0][18][0]=8.22064e-18;
                   1885: Tnc[16][0][2]=6.65871e-16;
                   1886: Tnc[14][2][2]=6.33318e-14;
                   1887: Tnc[12][4][2]=8.53499e-13;
                   1888: Tnc[10][6][2]=3.53829e-12;
                   1889: Tnc[8][8][2]=5.58296e-12;
                   1890: Tnc[6][10][2]=3.53829e-12;
                   1891: Tnc[4][12][2]=8.53499e-13;
                   1892: Tnc[2][14][2]=6.33318e-14;
                   1893: Tnc[0][16][2]=6.65871e-16;
                   1894: Tnc[14][0][4]=1.06539e-14;
                   1895: Tnc[12][2][4]=8.53499e-13;
                   1896: Tnc[10][4][4]=8.75527e-12;
                   1897: Tnc[8][6][4]=2.54496e-11;
                   1898: Tnc[6][8][4]=2.54496e-11;
                   1899: Tnc[4][10][4]=8.75527e-12;
                   1900: Tnc[2][12][4]=8.53499e-13;
                   1901: Tnc[0][14][4]=1.06539e-14;
                   1902: Tnc[12][0][6]=5.80048e-14;
                   1903: Tnc[10][2][6]=3.53829e-12;
                   1904: Tnc[8][4][6]=2.54496e-11;
                   1905: Tnc[6][6][6]=4.68831e-11;
                   1906: Tnc[4][8][6]=2.54496e-11;
                   1907: Tnc[2][10][6]=3.53829e-12;
                   1908: Tnc[0][12][6]=5.80048e-14;
                   1909: Tnc[10][0][8]=1.30511e-13;
                   1910: Tnc[8][2][8]=5.58296e-12;
                   1911: Tnc[6][4][8]=2.54496e-11;
                   1912: Tnc[4][6][8]=2.54496e-11;
                   1913: Tnc[2][8][8]=5.58296e-12;
                   1914: Tnc[0][10][8]=1.30511e-13;
                   1915: Tnc[8][0][10]=1.30511e-13;
                   1916: Tnc[6][2][10]=3.53829e-12;
                   1917: Tnc[4][4][10]=8.75527e-12;
                   1918: Tnc[2][6][10]=3.53829e-12;
                   1919: Tnc[0][8][10]=1.30511e-13;
                   1920: Tnc[6][0][12]=5.80048e-14;
                   1921: Tnc[4][2][12]=8.53499e-13;
                   1922: Tnc[2][4][12]=8.53499e-13;
                   1923: Tnc[0][6][12]=5.80048e-14;
                   1924: Tnc[4][0][14]=1.06539e-14;
                   1925: Tnc[2][2][14]=6.33318e-14;
                   1926: Tnc[0][4][14]=1.06539e-14;
                   1927: Tnc[2][0][16]=6.65871e-16;
                   1928: Tnc[0][2][16]=6.65871e-16;
                   1929: Tnc[0][0][18]=8.22064e-18;
                   1930: Tnc[15][1][1]=2.24917e-14;
                   1931: Tnc[13][3][1]=6.29766e-13;
                   1932: Tnc[11][5][1]=4.40836e-12;
                   1933: Tnc[9][7][1]=1.10209e-11;
                   1934: Tnc[7][9][1]=1.10209e-11;
                   1935: Tnc[5][11][1]=4.40836e-12;
                   1936: Tnc[3][13][1]=6.29766e-13;
                   1937: Tnc[1][15][1]=2.24917e-14;
                   1938: Tnc[13][1][3]=6.29766e-13;
                   1939: Tnc[11][3][3]=1.45371e-11;
                   1940: Tnc[9][5][3]=7.53096e-11;
                   1941: Tnc[7][7][3]=1.2674e-10;
                   1942: Tnc[5][9][3]=7.53096e-11;
                   1943: Tnc[3][11][3]=1.45371e-11;
                   1944: Tnc[1][13][3]=6.29766e-13;
                   1945: Tnc[11][1][5]=4.40836e-12;
                   1946: Tnc[9][3][5]=7.53096e-11;
                   1947: Tnc[7][5][5]=2.62613e-10;
                   1948: Tnc[5][7][5]=2.62613e-10;
                   1949: Tnc[3][9][5]=7.53096e-11;
                   1950: Tnc[1][11][5]=4.40836e-12;
                   1951: Tnc[9][1][7]=1.10209e-11;
                   1952: Tnc[7][3][7]=1.2674e-10;
                   1953: Tnc[5][5][7]=2.62613e-10;
                   1954: Tnc[3][7][7]=1.2674e-10;
                   1955: Tnc[1][9][7]=1.10209e-11;
                   1956: Tnc[7][1][9]=1.10209e-11;
                   1957: Tnc[5][3][9]=7.53096e-11;
                   1958: Tnc[3][5][9]=7.53096e-11;
                   1959: Tnc[1][7][9]=1.10209e-11;
                   1960: Tnc[5][1][11]=4.40836e-12;
                   1961: Tnc[3][3][11]=1.45371e-11;
                   1962: Tnc[1][5][11]=4.40836e-12;
                   1963: Tnc[3][1][13]=6.29766e-13;
                   1964: Tnc[1][3][13]=6.29766e-13;
                   1965: Tnc[1][1][15]=2.24917e-14;
                   1966: Tnc[16][0][0]=2.81146e-15;
                   1967: Tnc[14][2][0]=1.79933e-13;
                   1968: Tnc[12][4][0]=2.20418e-12;
                   1969: Tnc[10][6][0]=8.81673e-12;
                   1970: Tnc[8][8][0]=1.37761e-11;
                   1971: Tnc[6][10][0]=8.81673e-12;
                   1972: Tnc[4][12][0]=2.20418e-12;
                   1973: Tnc[2][14][0]=1.79933e-13;
                   1974: Tnc[0][16][0]=2.81146e-15;
                   1975: Tnc[14][0][2]=1.79933e-13;
                   1976: Tnc[12][2][2]=1.30677e-11;
                   1977: Tnc[10][4][2]=1.28945e-10;
                   1978: Tnc[8][6][2]=3.69201e-10;
                   1979: Tnc[6][8][2]=3.69201e-10;
                   1980: Tnc[4][10][2]=1.28945e-10;
                   1981: Tnc[2][12][2]=1.30677e-11;
                   1982: Tnc[0][14][2]=1.79933e-13;
                   1983: Tnc[12][0][4]=2.20418e-12;
                   1984: Tnc[10][2][4]=1.28945e-10;
                   1985: Tnc[8][4][4]=9.10406e-10;
                   1986: Tnc[6][6][4]=1.66809e-09;
                   1987: Tnc[4][8][4]=9.10406e-10;
                   1988: Tnc[2][10][4]=1.28945e-10;
                   1989: Tnc[0][12][4]=2.20418e-12;
                   1990: Tnc[10][0][6]=8.81673e-12;
                   1991: Tnc[8][2][6]=3.69201e-10;
                   1992: Tnc[6][4][6]=1.66809e-09;
                   1993: Tnc[4][6][6]=1.66809e-09;
                   1994: Tnc[2][8][6]=3.69201e-10;
                   1995: Tnc[0][10][6]=8.81673e-12;
                   1996: Tnc[8][0][8]=1.37761e-11;
                   1997: Tnc[6][2][8]=3.69201e-10;
                   1998: Tnc[4][4][8]=9.10406e-10;
                   1999: Tnc[2][6][8]=3.69201e-10;
                   2000: Tnc[0][8][8]=1.37761e-11;
                   2001: Tnc[6][0][10]=8.81673e-12;
                   2002: Tnc[4][2][10]=1.28945e-10;
                   2003: Tnc[2][4][10]=1.28945e-10;
                   2004: Tnc[0][6][10]=8.81673e-12;
                   2005: Tnc[4][0][12]=2.20418e-12;
                   2006: Tnc[2][2][12]=1.30677e-11;
                   2007: Tnc[0][4][12]=2.20418e-12;
                   2008: Tnc[2][0][14]=1.79933e-13;
                   2009: Tnc[0][2][14]=1.79933e-13;
                   2010: Tnc[0][0][16]=2.81146e-15;
                   2011: Tnc[13][1][1]=5.35301e-12;
                   2012: Tnc[11][3][1]=1.12413e-10;
                   2013: Tnc[9][5][1]=5.62067e-10;
                   2014: Tnc[7][7][1]=9.36778e-10;
                   2015: Tnc[5][9][1]=5.62067e-10;
                   2016: Tnc[3][11][1]=1.12413e-10;
                   2017: Tnc[1][13][1]=5.35301e-12;
                   2018: Tnc[11][1][3]=1.12413e-10;
                   2019: Tnc[9][3][3]=1.84679e-09;
                   2020: Tnc[7][5][3]=6.34332e-09;
                   2021: Tnc[5][7][3]=6.34332e-09;
                   2022: Tnc[3][9][3]=1.84679e-09;
                   2023: Tnc[1][11][3]=1.12413e-10;
                   2024: Tnc[9][1][5]=5.62067e-10;
                   2025: Tnc[7][3][5]=6.34332e-09;
                   2026: Tnc[5][5][5]=1.30721e-08;
                   2027: Tnc[3][7][5]=6.34332e-09;
                   2028: Tnc[1][9][5]=5.62067e-10;
                   2029: Tnc[7][1][7]=9.36778e-10;
                   2030: Tnc[5][3][7]=6.34332e-09;
                   2031: Tnc[3][5][7]=6.34332e-09;
                   2032: Tnc[1][7][7]=9.36778e-10;
                   2033: Tnc[5][1][9]=5.62067e-10;
                   2034: Tnc[3][3][9]=1.84679e-09;
                   2035: Tnc[1][5][9]=5.62067e-10;
                   2036: Tnc[3][1][11]=1.12413e-10;
                   2037: Tnc[1][3][11]=1.12413e-10;
                   2038: Tnc[1][1][13]=5.35301e-12;
                   2039: Tnc[14][0][0]=7.64716e-13;
                   2040: Tnc[12][2][0]=3.74711e-11;
                   2041: Tnc[10][4][0]=3.3724e-10;
                   2042: Tnc[8][6][0]=9.36778e-10;
                   2043: Tnc[6][8][0]=9.36778e-10;
                   2044: Tnc[4][10][0]=3.3724e-10;
                   2045: Tnc[2][12][0]=3.74711e-11;
                   2046: Tnc[0][14][0]=7.64716e-13;
                   2047: Tnc[12][0][2]=3.74711e-11;
                   2048: Tnc[10][2][2]=1.99132e-09;
                   2049: Tnc[8][4][2]=1.35699e-08;
                   2050: Tnc[6][6][2]=2.46239e-08;
                   2051: Tnc[4][8][2]=1.35699e-08;
                   2052: Tnc[2][10][2]=1.99132e-09;
                   2053: Tnc[0][12][2]=3.74711e-11;
                   2054: Tnc[10][0][4]=3.3724e-10;
                   2055: Tnc[8][2][4]=1.35699e-08;
                   2056: Tnc[6][4][4]=6.0382e-08;
                   2057: Tnc[4][6][4]=6.0382e-08;
                   2058: Tnc[2][8][4]=1.35699e-08;
                   2059: Tnc[0][10][4]=3.3724e-10;
                   2060: Tnc[8][0][6]=9.36778e-10;
                   2061: Tnc[6][2][6]=2.46239e-08;
                   2062: Tnc[4][4][6]=6.0382e-08;
                   2063: Tnc[2][6][6]=2.46239e-08;
                   2064: Tnc[0][8][6]=9.36778e-10;
                   2065: Tnc[6][0][8]=9.36778e-10;
                   2066: Tnc[4][2][8]=1.35699e-08;
                   2067: Tnc[2][4][8]=1.35699e-08;
                   2068: Tnc[0][6][8]=9.36778e-10;
                   2069: Tnc[4][0][10]=3.3724e-10;
                   2070: Tnc[2][2][10]=1.99132e-09;
                   2071: Tnc[0][4][10]=3.3724e-10;
                   2072: Tnc[2][0][12]=3.74711e-11;
                   2073: Tnc[0][2][12]=3.74711e-11;
                   2074: Tnc[0][0][14]=7.64716e-13;
                   2075: Tnc[11][1][1]=9.63543e-10;
                   2076: Tnc[9][3][1]=1.44531e-08;
                   2077: Tnc[7][5][1]=4.81771e-08;
                   2078: Tnc[5][7][1]=4.81771e-08;
                   2079: Tnc[3][9][1]=1.44531e-08;
                   2080: Tnc[1][11][1]=9.63543e-10;
                   2081: Tnc[9][1][3]=1.44531e-08;
                   2082: Tnc[7][3][3]=1.57379e-07;
                   2083: Tnc[5][5][3]=3.21181e-07;
                   2084: Tnc[3][7][3]=1.57379e-07;
                   2085: Tnc[1][9][3]=1.44531e-08;
                   2086: Tnc[7][1][5]=4.81771e-08;
                   2087: Tnc[5][3][5]=3.21181e-07;
                   2088: Tnc[3][5][5]=3.21181e-07;
                   2089: Tnc[1][7][5]=4.81771e-08;
                   2090: Tnc[5][1][7]=4.81771e-08;
                   2091: Tnc[3][3][7]=1.57379e-07;
                   2092: Tnc[1][5][7]=4.81771e-08;
                   2093: Tnc[3][1][9]=1.44531e-08;
                   2094: Tnc[1][3][9]=1.44531e-08;
                   2095: Tnc[1][1][11]=9.63543e-10;
                   2096: Tnc[12][0][0]=1.6059e-10;
                   2097: Tnc[10][2][0]=5.78126e-09;
                   2098: Tnc[8][4][0]=3.61328e-08;
                   2099: Tnc[6][6][0]=6.42362e-08;
                   2100: Tnc[4][8][0]=3.61328e-08;
                   2101: Tnc[2][10][0]=5.78126e-09;
                   2102: Tnc[0][12][0]=1.6059e-10;
                   2103: Tnc[10][0][2]=5.78126e-09;
                   2104: Tnc[8][2][2]=2.11979e-07;
                   2105: Tnc[6][4][2]=9.15365e-07;
                   2106: Tnc[4][6][2]=9.15365e-07;
                   2107: Tnc[2][8][2]=2.11979e-07;
                   2108: Tnc[0][10][2]=5.78126e-09;
                   2109: Tnc[8][0][4]=3.61328e-08;
                   2110: Tnc[6][2][4]=9.15365e-07;
                   2111: Tnc[4][4][4]=2.22257e-06;
                   2112: Tnc[2][6][4]=9.15365e-07;
                   2113: Tnc[0][8][4]=3.61328e-08;
                   2114: Tnc[6][0][6]=6.42362e-08;
                   2115: Tnc[4][2][6]=9.15365e-07;
                   2116: Tnc[2][4][6]=9.15365e-07;
                   2117: Tnc[0][6][6]=6.42362e-08;
                   2118: Tnc[4][0][8]=3.61328e-08;
                   2119: Tnc[2][2][8]=2.11979e-07;
                   2120: Tnc[0][4][8]=3.61328e-08;
                   2121: Tnc[2][0][10]=5.78126e-09;
                   2122: Tnc[0][2][10]=5.78126e-09;
                   2123: Tnc[0][0][12]=1.6059e-10;
                   2124: Tnc[9][1][1]=1.25261e-07;
                   2125: Tnc[7][3][1]=1.25261e-06;
                   2126: Tnc[5][5][1]=2.50521e-06;
                   2127: Tnc[3][7][1]=1.25261e-06;
                   2128: Tnc[1][9][1]=1.25261e-07;
                   2129: Tnc[7][1][3]=1.25261e-06;
                   2130: Tnc[5][3][3]=8.10018e-06;
                   2131: Tnc[3][5][3]=8.10018e-06;
                   2132: Tnc[1][7][3]=1.25261e-06;
                   2133: Tnc[5][1][5]=2.50521e-06;
                   2134: Tnc[3][3][5]=8.10018e-06;
                   2135: Tnc[1][5][5]=2.50521e-06;
                   2136: Tnc[3][1][7]=1.25261e-06;
                   2137: Tnc[1][3][7]=1.25261e-06;
                   2138: Tnc[1][1][9]=1.25261e-07;
                   2139: Tnc[10][0][0]=2.50521e-08;
                   2140: Tnc[8][2][0]=6.26303e-07;
                   2141: Tnc[6][4][0]=2.50521e-06;
                   2142: Tnc[4][6][0]=2.50521e-06;
                   2143: Tnc[2][8][0]=6.26303e-07;
                   2144: Tnc[0][10][0]=2.50521e-08;
                   2145: Tnc[8][0][2]=6.26303e-07;
                   2146: Tnc[6][2][2]=1.45302e-05;
                   2147: Tnc[4][4][2]=3.45719e-05;
                   2148: Tnc[2][6][2]=1.45302e-05;
                   2149: Tnc[0][8][2]=6.26303e-07;
                   2150: Tnc[6][0][4]=2.50521e-06;
                   2151: Tnc[4][2][4]=3.45719e-05;
                   2152: Tnc[2][4][4]=3.45719e-05;
                   2153: Tnc[0][6][4]=2.50521e-06;
                   2154: Tnc[4][0][6]=2.50521e-06;
                   2155: Tnc[2][2][6]=1.45302e-05;
                   2156: Tnc[0][4][6]=2.50521e-06;
                   2157: Tnc[2][0][8]=6.26303e-07;
                   2158: Tnc[0][2][8]=6.26303e-07;
                   2159: Tnc[0][0][10]=2.50521e-08;
                   2160: Tnc[7][1][1]=1.10229e-05;
                   2161: Tnc[5][3][1]=6.61376e-05;
                   2162: Tnc[3][5][1]=6.61376e-05;
                   2163: Tnc[1][7][1]=1.10229e-05;
                   2164: Tnc[5][1][3]=6.61376e-05;
                   2165: Tnc[3][3][3]=0.000209436;
                   2166: Tnc[1][5][3]=6.61376e-05;
                   2167: Tnc[3][1][5]=6.61376e-05;
                   2168: Tnc[1][3][5]=6.61376e-05;
                   2169: Tnc[1][1][7]=1.10229e-05;
                   2170: Tnc[8][0][0]=2.75573e-06;
                   2171: Tnc[6][2][0]=4.40917e-05;
                   2172: Tnc[4][4][0]=9.92063e-05;
                   2173: Tnc[2][6][0]=4.40917e-05;
                   2174: Tnc[0][8][0]=2.75573e-06;
                   2175: Tnc[6][0][2]=4.40917e-05;
                   2176: Tnc[4][2][2]=0.000562169;
                   2177: Tnc[2][4][2]=0.000562169;
                   2178: Tnc[0][6][2]=4.40917e-05;
                   2179: Tnc[4][0][4]=9.92063e-05;
                   2180: Tnc[2][2][4]=0.000562169;
                   2181: Tnc[0][4][4]=9.92063e-05;
                   2182: Tnc[2][0][6]=4.40917e-05;
                   2183: Tnc[0][2][6]=4.40917e-05;
                   2184: Tnc[0][0][8]=2.75573e-06;
                   2185: Tnc[5][1][1]=0.000595238;
                   2186: Tnc[3][3][1]=0.00178571;
                   2187: Tnc[1][5][1]=0.000595238;
                   2188: Tnc[3][1][3]=0.00178571;
                   2189: Tnc[1][3][3]=0.00178571;
                   2190: Tnc[1][1][5]=0.000595238;
                   2191: Tnc[6][0][0]=0.000198413;
                   2192: Tnc[4][2][0]=0.00178571;
                   2193: Tnc[2][4][0]=0.00178571;
                   2194: Tnc[0][6][0]=0.000198413;
                   2195: Tnc[4][0][2]=0.00178571;
                   2196: Tnc[2][2][2]=0.00952381;
                   2197: Tnc[0][4][2]=0.00178571;
                   2198: Tnc[2][0][4]=0.00178571;
                   2199: Tnc[0][2][4]=0.00178571;
                   2200: Tnc[0][0][6]=0.000198413;
                   2201: Tnc[3][1][1]=0.0166667;
                   2202: Tnc[1][3][1]=0.0166667;
                   2203: Tnc[1][1][3]=0.0166667;
                   2204: Tnc[4][0][0]=0.00833333;
                   2205: Tnc[2][2][0]=0.0333333;
                   2206: Tnc[0][4][0]=0.00833333;
                   2207: Tnc[2][0][2]=0.0333333;
                   2208: Tnc[0][2][2]=0.0333333;
                   2209: Tnc[0][0][4]=0.00833333;
                   2210: Tnc[1][1][1]=0.166667;
                   2211: Tnc[2][0][0]=0.166667;
                   2212: Tnc[0][2][0]=0.166667;
                   2213: Tnc[0][0][2]=0.166667;
                   2214: Tnc[0][0][0]=1;
                   2215: /* end of table */
                   2216:  f[0]=f[1]=f[2]=f[3]=0.0; ex[0][0][0] = 1.0;
                   2217:  for (i=0; i<deg; i++) {
                   2218:    for (j=0; j<deg; j++) {
                   2219:      for (k=0; k<deg; k++) {
                   2220:        if (i>0) ex[i][j][k] = ex[i-1][j][k]*(a*t);
                   2221:        else if (j>0) ex[i][j][k] = ex[i][j-1][k]*(b*t);
                   2222:        else if (k>0) ex[i][j][k] = ex[i][j][k-1]*(c*t);
                   2223:        f[0] += (ex[i][j][k])*(Tnc[i][j][k]);
                   2224:        if (i>0) f[1] += ((double) i)*(ex[i-1][j][k])*(Tnc[i][j][k]);
                   2225:        if (j>0) f[2] += ((double) j)*(ex[i][j-1][k])*(Tnc[i][j][k]);
                   2226:        if (k>0) f[3] += ((double) k)*(ex[i][j][k-1])*(Tnc[i][j][k]);
                   2227:      }
                   2228:    }
                   2229:  }
                   2230: }
                   2231:
                   2232:

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>