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