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