ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/multipole/QuadrupoleTensor.nb
Revision: 4433
Committed: Mon Jul 11 13:23:51 2016 UTC (9 years, 1 month ago) by gezelter
File size: 127171 byte(s)
Log Message:
Almost done

File Contents

# User Rev Content
1 gezelter 4432 (*CacheID: 234*)
2     (* Internal cache information:
3     NotebookFileLineBreakTest
4     NotebookFileLineBreakTest
5     NotebookDataPosition[ 0, 0]
6 gezelter 4433 NotebookDataLength[ 127170, 3497]
7     NotebookOptionsPosition[ 125250, 3424]
8     NotebookOutlinePosition[ 125606, 3440]
9     CellTagsIndexPosition[ 125563, 3437]
10 gezelter 4432 WindowFrame->Normal*)
11    
12     (* Beginning of Notebook Content *)
13     Notebook[{
14     Cell[BoxData[{
15     RowBox[{
16     RowBox[{"rc", "=."}], ";"}], "\[IndentingNewLine]",
17     RowBox[{
18     RowBox[{"\[Kappa]", "=."}], ";"}], "\[IndentingNewLine]",
19     RowBox[{
20     RowBox[{
21     RowBox[{"f", "[",
22     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
23     RowBox[{
24     RowBox[{"Erfc", "[",
25     RowBox[{"\[Kappa]", " ", "r"}], "]"}], "/", "r"}]}],
26     ";"}], "\[IndentingNewLine]",
27     RowBox[{
28     RowBox[{
29     RowBox[{"g", "[",
30     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
31     RowBox[{
32     RowBox[{
33     RowBox[{"(",
34     RowBox[{
35     RowBox[{"-", "2"}], "*", "\[Kappa]"}], ")"}], "/",
36     RowBox[{"(",
37     RowBox[{
38     RowBox[{"E", "^",
39     RowBox[{"(",
40     RowBox[{
41     RowBox[{"r", "^", "2"}], "*",
42     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
43     RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-",
44     RowBox[{
45     RowBox[{"Erfc", "[",
46     RowBox[{"r", "*", "\[Kappa]"}], "]"}], "/",
47     RowBox[{"r", "^", "2"}]}]}]}], " ", ";"}], "\[IndentingNewLine]",
48     RowBox[{
49     RowBox[{
50     RowBox[{"h", "[",
51     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
52     RowBox[{
53     RowBox[{
54     RowBox[{"(",
55     RowBox[{"4", "*", "\[Kappa]"}], ")"}], "/",
56     RowBox[{"(",
57     RowBox[{
58     RowBox[{"E", "^",
59     RowBox[{"(",
60     RowBox[{
61     RowBox[{"r", "^", "2"}], "*",
62     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
63     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
64     RowBox[{"r", "^", "2"}]}], ")"}]}], "+",
65     RowBox[{
66     RowBox[{"(",
67     RowBox[{"4", "*",
68     RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
69     RowBox[{"(",
70     RowBox[{
71     RowBox[{"E", "^",
72     RowBox[{"(",
73     RowBox[{
74     RowBox[{"r", "^", "2"}], "*",
75     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
76     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+",
77     RowBox[{
78     RowBox[{"(",
79     RowBox[{"2", "*",
80     RowBox[{"Erfc", "[",
81     RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
82     RowBox[{"r", "^", "3"}]}]}]}], ";"}], "\[IndentingNewLine]",
83     RowBox[{
84     RowBox[{
85     RowBox[{"s", "[",
86     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
87     RowBox[{
88     RowBox[{
89     RowBox[{"(",
90     RowBox[{
91     RowBox[{"-", "12"}], "*", "\[Kappa]"}], ")"}], "/",
92     RowBox[{"(",
93     RowBox[{
94     RowBox[{"E", "^",
95     RowBox[{"(",
96     RowBox[{
97     RowBox[{"r", "^", "2"}], "*",
98     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
99     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
100     RowBox[{"r", "^", "3"}]}], ")"}]}], "-",
101     RowBox[{
102     RowBox[{"(",
103     RowBox[{"8", "*",
104     RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
105     RowBox[{"(",
106     RowBox[{
107     RowBox[{"E", "^",
108     RowBox[{"(",
109     RowBox[{
110     RowBox[{"r", "^", "2"}], "*",
111     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
112     RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-",
113     RowBox[{
114     RowBox[{"(",
115     RowBox[{"8", "*", "r", "*",
116     RowBox[{"\[Kappa]", "^", "5"}]}], ")"}], "/",
117     RowBox[{"(",
118     RowBox[{
119     RowBox[{"E", "^",
120     RowBox[{"(",
121     RowBox[{
122     RowBox[{"r", "^", "2"}], "*",
123     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
124     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "-",
125     RowBox[{
126     RowBox[{"(",
127     RowBox[{"6", "*",
128     RowBox[{"Erfc", "[",
129     RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
130     RowBox[{"r", "^", "4"}]}]}]}], ";"}], "\[IndentingNewLine]",
131     RowBox[{
132     RowBox[{
133     RowBox[{"t", "[",
134     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
135     RowBox[{
136     RowBox[{
137     RowBox[{"(",
138     RowBox[{"48", "*", "\[Kappa]"}], ")"}], "/",
139     RowBox[{"(",
140     RowBox[{
141     RowBox[{"E", "^",
142     RowBox[{"(",
143     RowBox[{
144     RowBox[{"r", "^", "2"}], "*",
145     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
146     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
147     RowBox[{"r", "^", "4"}]}], ")"}]}], "+",
148     RowBox[{
149     RowBox[{"(",
150     RowBox[{"32", "*",
151     RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
152     RowBox[{"(",
153     RowBox[{
154     RowBox[{"E", "^",
155     RowBox[{"(",
156     RowBox[{
157     RowBox[{"r", "^", "2"}], "*",
158     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
159     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
160     RowBox[{"r", "^", "2"}]}], ")"}]}], "+",
161     RowBox[{
162     RowBox[{"(",
163     RowBox[{"8", "*",
164     RowBox[{"\[Kappa]", "^", "5"}]}], ")"}], "/",
165     RowBox[{"(",
166     RowBox[{
167     RowBox[{"E", "^",
168     RowBox[{"(",
169     RowBox[{
170     RowBox[{"r", "^", "2"}], "*",
171     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
172     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+",
173     RowBox[{
174     RowBox[{"(",
175     RowBox[{"16", "*",
176     RowBox[{"r", "^", "2"}], "*",
177     RowBox[{"\[Kappa]", "^", "7"}]}], ")"}], "/",
178     RowBox[{"(",
179     RowBox[{
180     RowBox[{"E", "^",
181     RowBox[{"(",
182     RowBox[{
183     RowBox[{"r", "^", "2"}], "*",
184     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
185     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+",
186     RowBox[{
187     RowBox[{"(",
188     RowBox[{"24", "*",
189     RowBox[{"Erfc", "[",
190     RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
191     RowBox[{"r", "^", "5"}]}]}]}], ";"}], "\[IndentingNewLine]",
192     RowBox[{
193     RowBox[{
194     RowBox[{"u", "[",
195     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
196     RowBox[{
197     RowBox[{
198     RowBox[{"(",
199     RowBox[{
200     RowBox[{"-", "240"}], "*", "\[Kappa]"}], ")"}], "/",
201     RowBox[{"(",
202     RowBox[{
203     RowBox[{"E", "^",
204     RowBox[{"(",
205     RowBox[{
206     RowBox[{"r", "^", "2"}], "*",
207     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
208     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
209     RowBox[{"r", "^", "5"}]}], ")"}]}], "-",
210     RowBox[{
211     RowBox[{"(",
212     RowBox[{"160", "*",
213     RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/",
214     RowBox[{"(",
215     RowBox[{
216     RowBox[{"E", "^",
217     RowBox[{"(",
218     RowBox[{
219     RowBox[{"r", "^", "2"}], "*",
220     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
221     RowBox[{"Sqrt", "[", "Pi", "]"}], "*",
222     RowBox[{"r", "^", "3"}]}], ")"}]}], "-",
223     RowBox[{
224     RowBox[{"(",
225     RowBox[{"64", "*",
226     RowBox[{"\[Kappa]", "^", "5"}]}], ")"}], "/",
227     RowBox[{"(",
228     RowBox[{
229     RowBox[{"E", "^",
230     RowBox[{"(",
231     RowBox[{
232     RowBox[{"r", "^", "2"}], "*",
233     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
234     RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "+",
235     RowBox[{
236     RowBox[{"(",
237     RowBox[{"16", "*", "r", "*",
238     RowBox[{"\[Kappa]", "^", "7"}]}], ")"}], "/",
239     RowBox[{"(",
240     RowBox[{
241     RowBox[{"E", "^",
242     RowBox[{"(",
243     RowBox[{
244     RowBox[{"r", "^", "2"}], "*",
245     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
246     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "-",
247     RowBox[{
248     RowBox[{"(",
249     RowBox[{"32", "*",
250     RowBox[{"r", "^", "3"}], "*",
251     RowBox[{"\[Kappa]", "^", "9"}]}], ")"}], "/",
252     RowBox[{"(",
253     RowBox[{
254     RowBox[{"E", "^",
255     RowBox[{"(",
256     RowBox[{
257     RowBox[{"r", "^", "2"}], "*",
258     RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*",
259     RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "-",
260     RowBox[{
261     RowBox[{"(",
262     RowBox[{"120", "*",
263     RowBox[{"Erfc", "[",
264     RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/",
265     RowBox[{"r", "^", "6"}]}]}]}], ";"}], "\[IndentingNewLine]",
266     RowBox[{
267     RowBox[{
268     RowBox[{"f4", "[",
269     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
270     RowBox[{
271     RowBox[{"f", "[",
272     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
273     RowBox[{"f", "[",
274     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
275     RowBox[{
276     RowBox[{"(",
277     RowBox[{"r", "-", "rc"}], ")"}],
278     RowBox[{"g", "[",
279     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
280     RowBox[{
281     RowBox[{
282     RowBox[{"(",
283     RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}], " ",
284     RowBox[{
285     RowBox[{"h", "[",
286     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}], " ", "-", " ",
287     RowBox[{
288     RowBox[{
289     RowBox[{"(",
290     RowBox[{"r", "-", "rc"}], ")"}], "^", "3"}], " ",
291     RowBox[{
292     RowBox[{"s", "[",
293     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "6"}]}], " ", "-",
294     RowBox[{
295     RowBox[{
296     RowBox[{"(",
297     RowBox[{"r", "-", "rc"}], ")"}], "^", "4"}], " ",
298     RowBox[{
299     RowBox[{"t", "[",
300     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "24"}]}], " ", "-", " ",
301     RowBox[{
302     RowBox[{
303     RowBox[{"(",
304     RowBox[{"r", "-", "rc"}], ")"}], "^", "5"}], " ",
305     RowBox[{
306     RowBox[{"u", "[",
307     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "120"}]}]}]}], " ",
308     ";"}], "\[IndentingNewLine]",
309     RowBox[{
310     RowBox[{
311     RowBox[{"g4", "[",
312     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ",
313     RowBox[{
314     RowBox[{"g", "[",
315     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ",
316     RowBox[{"g", "[",
317     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
318     RowBox[{
319     RowBox[{"(",
320     RowBox[{"r", "-", "rc"}], ")"}],
321     RowBox[{"h", "[",
322     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
323     RowBox[{
324     RowBox[{
325     RowBox[{"(",
326     RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}], " ",
327     RowBox[{
328     RowBox[{"s", "[",
329     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}], " ", "-",
330     RowBox[{
331     RowBox[{
332     RowBox[{"(",
333     RowBox[{"r", "-", "rc"}], ")"}], "^", "3"}], " ",
334     RowBox[{
335     RowBox[{"t", "[",
336     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "6"}]}], "-",
337     RowBox[{
338     RowBox[{
339     RowBox[{"(",
340     RowBox[{"r", "-", "rc"}], ")"}], "^", "4"}], " ",
341     RowBox[{
342     RowBox[{"u", "[",
343     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "24"}]}]}]}],
344     ";"}], "\[IndentingNewLine]",
345     RowBox[{
346     RowBox[{
347     RowBox[{"h4", "[",
348     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
349     RowBox[{
350     RowBox[{"h", "[",
351     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
352     RowBox[{"h", "[",
353     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], " ", "-", " ",
354     RowBox[{
355     RowBox[{"(",
356     RowBox[{"r", "-", "rc"}], ")"}],
357     RowBox[{"s", "[",
358     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
359     RowBox[{
360     RowBox[{
361     RowBox[{"(",
362     RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}],
363     RowBox[{
364     RowBox[{"t", "[",
365     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}], " ", "-", " ",
366     RowBox[{
367     RowBox[{
368     RowBox[{"(",
369     RowBox[{"r", "-", "rc"}], ")"}], "^", "3"}], " ",
370     RowBox[{
371     RowBox[{"u", "[",
372     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "6"}]}]}]}],
373     ";"}], "\[IndentingNewLine]",
374     RowBox[{
375     RowBox[{
376     RowBox[{"s4", "[",
377     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
378     RowBox[{
379     RowBox[{"s", "[",
380     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
381     RowBox[{"s", "[",
382     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], " ", "-", " ",
383     RowBox[{
384     RowBox[{"(",
385     RowBox[{"r", "-", "rc"}], ")"}],
386     RowBox[{"t", "[",
387     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-",
388     RowBox[{
389     RowBox[{
390     RowBox[{"(",
391     RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}],
392     RowBox[{
393     RowBox[{"u", "[",
394     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}]}]}], " ",
395     ";"}], "\[IndentingNewLine]",
396     RowBox[{
397     RowBox[{
398     RowBox[{"t4", "[",
399     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
400     RowBox[{
401     RowBox[{"t", "[",
402     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
403     RowBox[{"t", "[",
404     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], " ", "-", " ",
405     RowBox[{
406     RowBox[{"(",
407     RowBox[{"r", "-", "rc"}], ")"}],
408     RowBox[{"u", "[",
409     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], " ",
410     ";"}], "\[IndentingNewLine]",
411     RowBox[{
412     RowBox[{
413     RowBox[{"u4", "[",
414     RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=",
415     RowBox[{
416     RowBox[{"u", "[",
417     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
418     RowBox[{"u", "[",
419     RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], " ",
420     ";"}], "\[IndentingNewLine]",
421     RowBox[{
422     RowBox[{
423     RowBox[{"v41", "[",
424     RowBox[{"r_", ",", "method_"}], "]"}], ":=",
425     RowBox[{"Switch", "[",
426     RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", " ",
427     RowBox[{
428     RowBox[{
429     RowBox[{"h", "[",
430     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
431     RowBox[{"r", "^", "2"}]}], " ", "-", " ",
432     RowBox[{
433     RowBox[{"g", "[",
434     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
435     RowBox[{"r", "^", "3"}]}]}], ",", " ", "\[IndentingNewLine]", "TSF",
436     ",", " ",
437     RowBox[{
438     RowBox[{
439     RowBox[{"h4", "[",
440     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
441     RowBox[{"r", "^", "2"}]}], " ", "-", " ",
442     RowBox[{
443     RowBox[{"g4", "[",
444     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
445     RowBox[{"r", "^", "3"}]}]}], ",", "\[IndentingNewLine]", "SP", ",",
446     RowBox[{
447     RowBox[{
448     RowBox[{"h", "[",
449     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
450     RowBox[{"r", "^", "2"}]}], " ", "-", " ",
451     RowBox[{
452     RowBox[{"g", "[",
453     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
454     RowBox[{"r", "^", "3"}]}], " ", "-", " ",
455     RowBox[{
456     RowBox[{"h", "[",
457     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
458     RowBox[{"rc", "^", "2"}]}], " ", "+", " ",
459     RowBox[{
460     RowBox[{"g", "[",
461     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
462     RowBox[{"rc", "^", "3"}]}]}], ",", "\[IndentingNewLine]", "GSF", ",",
463     " ",
464     RowBox[{
465     RowBox[{
466     RowBox[{"h", "[",
467     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
468     RowBox[{"r", "^", "2"}]}], " ", "-", " ",
469     RowBox[{
470     RowBox[{"g", "[",
471     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
472     RowBox[{"r", "^", "3"}]}], " ", "-", " ",
473     RowBox[{
474     RowBox[{"h", "[",
475     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
476     RowBox[{"rc", "^", "2"}]}], " ", "+", " ",
477     RowBox[{
478     RowBox[{"g", "[",
479     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
480     RowBox[{"rc", "^", "3"}]}], " ", "-", " ",
481     RowBox[{
482     RowBox[{"(",
483     RowBox[{"r", "-", "rc"}], ")"}],
484     RowBox[{"(",
485     RowBox[{
486     RowBox[{"3",
487     RowBox[{
488     RowBox[{"g", "[",
489     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
490     RowBox[{"rc", "^", "4"}]}]}], " ", "-",
491     RowBox[{"3",
492     RowBox[{
493     RowBox[{"h", "[",
494     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
495     RowBox[{"rc", "^", "3"}]}]}], "+",
496     RowBox[{
497     RowBox[{"s", "[",
498     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
499     RowBox[{"rc", "^", "2"}]}]}], ")"}]}]}]}], "\[IndentingNewLine]",
500     "]"}]}], ";"}], "\[IndentingNewLine]",
501     RowBox[{
502     RowBox[{
503     RowBox[{"v42", "[",
504     RowBox[{"r_", ",", "method_"}], "]"}], ":=",
505     RowBox[{"Switch", "[",
506     RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", " ",
507     RowBox[{
508     RowBox[{"3", " ",
509     RowBox[{
510     RowBox[{"g", "[",
511     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
512     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
513     RowBox[{"3",
514     RowBox[{
515     RowBox[{"h", "[",
516     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
517     RowBox[{"r", "^", "2"}]}]}], "+", " ",
518     RowBox[{
519     RowBox[{"s", "[",
520     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], ",", " ",
521     "\[IndentingNewLine]", "TSF", ",", " ",
522     RowBox[{
523     RowBox[{"3", " ",
524     RowBox[{
525     RowBox[{"g4", "[",
526     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
527     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
528     RowBox[{"3",
529     RowBox[{
530     RowBox[{"h4", "[",
531     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
532     RowBox[{"r", "^", "2"}]}]}], "+", " ",
533     RowBox[{
534     RowBox[{"s4", "[",
535     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], ",",
536     "\[IndentingNewLine]", "SP", ",", " ",
537     RowBox[{
538     RowBox[{"3", " ",
539     RowBox[{
540     RowBox[{"g", "[",
541     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
542     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
543     RowBox[{"3",
544     RowBox[{
545     RowBox[{"h", "[",
546     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
547     RowBox[{"r", "^", "2"}]}]}], "+", " ",
548     RowBox[{
549     RowBox[{"s", "[",
550     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-", " ",
551     RowBox[{"(",
552     RowBox[{
553     RowBox[{"3", " ",
554     RowBox[{
555     RowBox[{"g", "[",
556     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
557     RowBox[{"rc", "^", "3"}]}]}], " ", "-", " ",
558     RowBox[{"3",
559     RowBox[{
560     RowBox[{"h", "[",
561     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
562     RowBox[{"rc", "^", "2"}]}]}], "+", " ",
563     RowBox[{
564     RowBox[{"s", "[",
565     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], ")"}]}], ",",
566     "\[IndentingNewLine]", "GSF", ",", " ",
567     RowBox[{
568     RowBox[{"3", " ",
569     RowBox[{
570     RowBox[{"g", "[",
571     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
572     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
573     RowBox[{"3",
574     RowBox[{
575     RowBox[{"h", "[",
576     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
577     RowBox[{"r", "^", "2"}]}]}], "+", " ",
578     RowBox[{
579     RowBox[{"s", "[",
580     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}], " ", "-", " ",
581     RowBox[{"(",
582     RowBox[{
583     RowBox[{"3", " ",
584     RowBox[{
585     RowBox[{"g", "[",
586     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
587     RowBox[{"rc", "^", "3"}]}]}], " ", "-", " ",
588     RowBox[{"3",
589     RowBox[{
590     RowBox[{"h", "[",
591     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
592     RowBox[{"rc", "^", "2"}]}]}], "+", " ",
593     RowBox[{
594     RowBox[{"s", "[",
595     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], ")"}], "-",
596     " ",
597     RowBox[{
598     RowBox[{"(",
599     RowBox[{"r", "-", "rc"}], ")"}],
600     RowBox[{"(",
601     RowBox[{
602     RowBox[{
603     RowBox[{"t", "[",
604     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}], " ", "-", " ",
605     RowBox[{"4", " ",
606     RowBox[{
607     RowBox[{"s", "[",
608     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
609     RowBox[{"rc", "^", "2"}]}]}], " ", "+", " ",
610     RowBox[{"9", " ",
611     RowBox[{
612     RowBox[{"h", "[",
613     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
614     RowBox[{"rc", "^", "3"}]}]}], " ", "-", " ",
615     RowBox[{"9",
616     RowBox[{
617     RowBox[{"g", "[",
618     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
619     RowBox[{"rc", "^", "4"}]}]}]}], ")"}]}]}]}], "\[IndentingNewLine]",
620     "]"}]}], ";"}], "\[IndentingNewLine]",
621     RowBox[{
622     RowBox[{
623     RowBox[{"v43", "[",
624     RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ",
625     RowBox[{"Switch", "[",
626     RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",",
627     RowBox[{
628     RowBox[{"t", "[",
629     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
630     RowBox[{"6",
631     RowBox[{
632     RowBox[{"s", "[",
633     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], "+",
634     RowBox[{"15", " ",
635     RowBox[{
636     RowBox[{"h", "[",
637     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
638     RowBox[{"r", "^", "2"}]}]}], " ", "-", " ",
639     RowBox[{"15", " ",
640     RowBox[{
641     RowBox[{"g", "[",
642     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
643     RowBox[{"r", "^", "3"}]}]}]}], ",", " ", "\[IndentingNewLine]", "TSF",
644     ",", " ",
645     RowBox[{
646     RowBox[{"t4", "[",
647     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
648     RowBox[{"6",
649     RowBox[{
650     RowBox[{"s4", "[",
651     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], "+",
652     RowBox[{"15", " ",
653     RowBox[{
654     RowBox[{"h4", "[",
655     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
656     RowBox[{"r", "^", "2"}]}]}], " ", "-", " ",
657     RowBox[{"15", " ",
658     RowBox[{
659     RowBox[{"g4", "[",
660     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
661     RowBox[{"r", "^", "3"}]}]}]}], ",", " ", "\[IndentingNewLine]", "SP",
662     ",", " ",
663     RowBox[{
664     RowBox[{"t", "[",
665     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
666     RowBox[{"6",
667     RowBox[{
668     RowBox[{"s", "[",
669     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], "+",
670     RowBox[{"15", " ",
671     RowBox[{
672     RowBox[{"h", "[",
673     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
674     RowBox[{"r", "^", "2"}]}]}], " ", "-", " ",
675     RowBox[{"15", " ",
676     RowBox[{
677     RowBox[{"g", "[",
678     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
679     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
680     RowBox[{"(",
681     RowBox[{
682     RowBox[{"t", "[",
683     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
684     RowBox[{"6",
685     RowBox[{
686     RowBox[{"s", "[",
687     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], "+",
688     RowBox[{"15", " ",
689     RowBox[{
690     RowBox[{"h", "[",
691     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
692     RowBox[{"rc", "^", "2"}]}]}], " ", "-", " ",
693     RowBox[{"15", " ",
694     RowBox[{
695     RowBox[{"g", "[",
696     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
697     RowBox[{"rc", "^", "3"}]}]}]}], ")"}]}], ",", " ",
698     "\[IndentingNewLine]", "GSF", ",", " ",
699     RowBox[{
700     RowBox[{"t", "[",
701     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-",
702     RowBox[{"6",
703     RowBox[{
704     RowBox[{"s", "[",
705     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/", "r"}]}], "+",
706     RowBox[{"15", " ",
707     RowBox[{
708     RowBox[{"h", "[",
709     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
710     RowBox[{"r", "^", "2"}]}]}], " ", "-", " ",
711     RowBox[{"15", " ",
712     RowBox[{
713     RowBox[{"g", "[",
714     RowBox[{"r", ",", "\[Kappa]"}], "]"}], "/",
715     RowBox[{"r", "^", "3"}]}]}], " ", "-", " ",
716     RowBox[{"(",
717     RowBox[{
718     RowBox[{"t", "[",
719     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-",
720     RowBox[{"6",
721     RowBox[{
722     RowBox[{"s", "[",
723     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], "+",
724     RowBox[{"15", " ",
725     RowBox[{
726     RowBox[{"h", "[",
727     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
728     RowBox[{"rc", "^", "2"}]}]}], " ", "-", " ",
729     RowBox[{"15", " ",
730     RowBox[{
731     RowBox[{"g", "[",
732     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
733     RowBox[{"rc", "^", "3"}]}]}]}], ")"}], " ", "-", " ",
734     RowBox[{
735     RowBox[{"(",
736     RowBox[{"r", "-", "rc"}], ")"}],
737     RowBox[{"(",
738     RowBox[{
739     RowBox[{"u", "[",
740     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-", " ",
741     RowBox[{"6", " ",
742     RowBox[{
743     RowBox[{"t", "[",
744     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "rc"}]}], " ", "+",
745     " ",
746     RowBox[{"21",
747     RowBox[{
748     RowBox[{"s", "[",
749     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
750     RowBox[{"rc", "^", "2"}]}]}], "-",
751     RowBox[{"45",
752     RowBox[{
753     RowBox[{"h", "[",
754     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
755     RowBox[{"rc", "^", "3"}]}]}], " ", "+", " ",
756     RowBox[{"45",
757     RowBox[{
758     RowBox[{"g", "[",
759     RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/",
760     RowBox[{"rc", "^", "4"}]}]}]}], ")"}]}]}]}], "\[IndentingNewLine]",
761     "]"}]}], ";"}], "\[IndentingNewLine]",
762     RowBox[{
763     RowBox[{
764     RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=",
765     RowBox[{"r", " ",
766     RowBox[{"Switch", "[",
767     RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",",
768     RowBox[{
769     RowBox[{"Sin", "[", "\[Theta]", "]"}],
770     RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
771     "2", ",",
772     RowBox[{
773     RowBox[{"Sin", "[", "\[Theta]", "]"}],
774     RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]",
775     "3", ",",
776     RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}],
777     ";"}], "\[IndentingNewLine]",
778     RowBox[{
779     RowBox[{
780     RowBox[{"T", "[",
781     RowBox[{
782     "\[Alpha]_", ",", "\[Beta]_", ",", "\[Gamma]_", ",", "\[Delta]_", ",",
783     " ", "method_"}], "]"}], " ", ":=", " ",
784     RowBox[{
785     RowBox[{
786     RowBox[{"(",
787     RowBox[{
788     RowBox[{
789     RowBox[{"KroneckerDelta", "[",
790     RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*",
791     RowBox[{"KroneckerDelta", "[",
792     RowBox[{"\[Gamma]", ",", "\[Delta]"}], "]"}]}], " ", "+", " ",
793     RowBox[{
794     RowBox[{"KroneckerDelta", "[",
795     RowBox[{"\[Alpha]", ",", "\[Gamma]"}], "]"}], "*",
796     RowBox[{"KroneckerDelta", "[",
797     RowBox[{"\[Beta]", ",", "\[Delta]"}], "]"}]}], "+", " ",
798     RowBox[{
799     RowBox[{"KroneckerDelta", "[",
800     RowBox[{"\[Beta]", ",", "\[Gamma]"}], "]"}], "*",
801     RowBox[{"KroneckerDelta", "[",
802     RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}]}]}], ")"}], " ",
803     RowBox[{"v41", "[",
804     RowBox[{"r", ",", " ", "method"}], "]"}]}], "+", " ",
805     RowBox[{
806     RowBox[{"(",
807     RowBox[{
808     RowBox[{
809     RowBox[{"KroneckerDelta", "[",
810     RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}],
811     RowBox[{"R", "[", "\[Gamma]", "]"}],
812     RowBox[{"R", "[", "\[Delta]", "]"}]}], "+", " ",
813     RowBox[{
814     RowBox[{"KroneckerDelta", "[",
815     RowBox[{"\[Alpha]", ",", "\[Gamma]"}], "]"}],
816     RowBox[{"R", "[", "\[Beta]", "]"}],
817     RowBox[{"R", "[", "\[Delta]", "]"}]}], "+",
818     RowBox[{
819     RowBox[{"KroneckerDelta", "[",
820     RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}],
821     RowBox[{"R", "[", "\[Beta]", "]"}],
822     RowBox[{"R", "[", "\[Gamma]", "]"}]}], "+",
823     RowBox[{
824     RowBox[{"KroneckerDelta", "[",
825     RowBox[{"\[Beta]", ",", "\[Gamma]"}], "]"}],
826     RowBox[{"R", "[", "\[Alpha]", "]"}],
827     RowBox[{"R", "[", "\[Delta]", "]"}]}], "+",
828     RowBox[{
829     RowBox[{"KroneckerDelta", "[",
830     RowBox[{"\[Beta]", ",", "\[Delta]"}], "]"}],
831     RowBox[{"R", "[", "\[Alpha]", "]"}],
832     RowBox[{"R", "[", "\[Gamma]", "]"}]}], "+",
833     RowBox[{
834     RowBox[{"KroneckerDelta", "[",
835     RowBox[{"\[Gamma]", ",", "\[Delta]"}], "]"}],
836     RowBox[{"R", "[", "\[Alpha]", "]"}],
837     RowBox[{"R", "[", "\[Beta]", "]"}]}]}], ")"}],
838     RowBox[{"(",
839     RowBox[{
840     RowBox[{"v42", "[",
841     RowBox[{"r", ",", "method"}], "]"}], "/",
842     RowBox[{"r", "^", "2"}]}], ")"}]}], " ", "+", " ",
843     RowBox[{
844     RowBox[{"R", "[", "\[Alpha]", "]"}],
845     RowBox[{"R", "[", "\[Beta]", "]"}],
846     RowBox[{"R", "[", "\[Gamma]", "]"}],
847     RowBox[{"R", "[", "\[Delta]", "]"}],
848     RowBox[{"(",
849     RowBox[{
850     RowBox[{"v43", "[",
851     RowBox[{"r", ",", "method"}], "]"}], "/",
852     RowBox[{"r", "^", "4"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
853     RowBox[{
854     RowBox[{"TBare", " ", "=", " ",
855     RowBox[{"MatrixForm", "[",
856     RowBox[{"Table", "[",
857     RowBox[{
858     RowBox[{"Table", "[",
859     RowBox[{
860     RowBox[{"Table", "[",
861     RowBox[{
862     RowBox[{"Table", "[",
863     RowBox[{
864     RowBox[{"Integrate", "[",
865     RowBox[{
866     RowBox[{
867     RowBox[{"r", "^", "2"}],
868     RowBox[{"Integrate", "[",
869     RowBox[{
870     RowBox[{"Integrate", "[",
871     RowBox[{
872     RowBox[{
873     RowBox[{"Sin", "[", "\[Theta]", "]"}],
874     RowBox[{"T", "[",
875     RowBox[{"i", ",", "j", ",", "k", ",", "l", ",", "BARE"}],
876     "]"}]}], ",",
877     RowBox[{"{",
878     RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
879     ",",
880     RowBox[{"{",
881     RowBox[{"\[Phi]", ",", "0", ",",
882     RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
883     RowBox[{"{",
884     RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
885     RowBox[{"{",
886     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
887     RowBox[{"{",
888     RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
889     RowBox[{"{",
890     RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
891     RowBox[{"{",
892     RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
893     ";"}], "\[IndentingNewLine]",
894     RowBox[{
895     RowBox[{"TSP", " ", "=", " ",
896     RowBox[{"MatrixForm", "[",
897     RowBox[{"Table", "[",
898     RowBox[{
899     RowBox[{"Table", "[",
900     RowBox[{
901     RowBox[{"Table", "[",
902     RowBox[{
903     RowBox[{"Table", "[",
904     RowBox[{
905     RowBox[{"Integrate", "[",
906     RowBox[{
907     RowBox[{
908     RowBox[{"r", "^", "2"}],
909     RowBox[{"Integrate", "[",
910     RowBox[{
911     RowBox[{"Integrate", "[",
912     RowBox[{
913     RowBox[{
914     RowBox[{"Sin", "[", "\[Theta]", "]"}],
915     RowBox[{"T", "[",
916     RowBox[{"i", ",", "j", ",", "k", ",", "l", ",", "SP"}],
917     "]"}]}], ",",
918     RowBox[{"{",
919     RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
920     ",",
921     RowBox[{"{",
922     RowBox[{"\[Phi]", ",", "0", ",",
923     RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
924     RowBox[{"{",
925     RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
926     RowBox[{"{",
927     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
928     RowBox[{"{",
929     RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
930     RowBox[{"{",
931     RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
932     RowBox[{"{",
933     RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
934     ";"}], "\[IndentingNewLine]",
935     RowBox[{
936     RowBox[{"TGSF", " ", "=", " ",
937     RowBox[{"MatrixForm", "[",
938     RowBox[{"Table", "[",
939     RowBox[{
940     RowBox[{"Table", "[",
941     RowBox[{
942     RowBox[{"Table", "[",
943     RowBox[{
944     RowBox[{"Table", "[",
945     RowBox[{
946     RowBox[{"Integrate", "[",
947     RowBox[{
948     RowBox[{
949     RowBox[{"r", "^", "2"}],
950     RowBox[{"Integrate", "[",
951     RowBox[{
952     RowBox[{"Integrate", "[",
953     RowBox[{
954     RowBox[{
955     RowBox[{"Sin", "[", "\[Theta]", "]"}],
956     RowBox[{"T", "[",
957     RowBox[{"i", ",", "j", ",", "k", ",", "l", ",", "GSF"}],
958     "]"}]}], ",",
959     RowBox[{"{",
960     RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
961     ",",
962     RowBox[{"{",
963     RowBox[{"\[Phi]", ",", "0", ",",
964     RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
965     RowBox[{"{",
966     RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
967     RowBox[{"{",
968     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
969     RowBox[{"{",
970     RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
971     RowBox[{"{",
972     RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
973     RowBox[{"{",
974     RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
975     ";"}], "\[IndentingNewLine]",
976     RowBox[{
977     RowBox[{"TTSF", " ", "=", " ",
978     RowBox[{"MatrixForm", "[",
979     RowBox[{"Table", "[",
980     RowBox[{
981     RowBox[{"Table", "[",
982     RowBox[{
983     RowBox[{"Table", "[",
984     RowBox[{
985     RowBox[{"Table", "[",
986     RowBox[{
987     RowBox[{"Integrate", "[",
988     RowBox[{
989     RowBox[{
990     RowBox[{"r", "^", "2"}],
991     RowBox[{"Integrate", "[",
992     RowBox[{
993     RowBox[{"Integrate", "[",
994     RowBox[{
995     RowBox[{
996     RowBox[{"Sin", "[", "\[Theta]", "]"}],
997     RowBox[{"T", "[",
998     RowBox[{"i", ",", "j", ",", "k", ",", "l", ",", "TSF"}],
999     "]"}]}], ",",
1000     RowBox[{"{",
1001     RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}],
1002     ",",
1003     RowBox[{"{",
1004     RowBox[{"\[Phi]", ",", "0", ",",
1005     RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",",
1006     RowBox[{"{",
1007     RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",",
1008     RowBox[{"{",
1009     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
1010     RowBox[{"{",
1011     RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
1012     RowBox[{"{",
1013     RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], ",",
1014     RowBox[{"{",
1015     RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}],
1016     ";"}], "\[IndentingNewLine]",
1017     RowBox[{
1018     RowBox[{"QBARE", " ", "=",
1019 gezelter 4433 RowBox[{
1020     RowBox[{"Part", "[",
1021     RowBox[{"TBare", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}],
1022     "]"}], "/",
1023     RowBox[{"(",
1024     RowBox[{"12", " ", "Pi"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
1025 gezelter 4432 RowBox[{
1026     RowBox[{"QSP", " ", "=", " ",
1027 gezelter 4433 RowBox[{
1028     RowBox[{"Part", "[",
1029     RowBox[{"TSP", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}], "]"}],
1030     "/",
1031     RowBox[{"(",
1032     RowBox[{"12", " ", "Pi"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
1033 gezelter 4432 RowBox[{
1034     RowBox[{"QGSF", " ", "=", " ",
1035 gezelter 4433 RowBox[{
1036     RowBox[{"Part", "[",
1037     RowBox[{"TGSF", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}],
1038     "]"}], "/",
1039     RowBox[{"(",
1040     RowBox[{"12", " ", "Pi"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
1041 gezelter 4432 RowBox[{
1042     RowBox[{"QTSF", " ", "=", " ",
1043 gezelter 4433 RowBox[{
1044     RowBox[{"Part", "[",
1045     RowBox[{"TTSF", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}],
1046     "]"}], "/",
1047     RowBox[{"(",
1048     RowBox[{"12", " ", "Pi"}], ")"}]}]}], ";"}]}], "Input",
1049 gezelter 4432 CellChangeTimes->{{3.613731558054296*^9, 3.613732200209785*^9},
1050     3.6137322781282454`*^9, {3.61373233169027*^9, 3.613732397354877*^9}, {
1051     3.613732488348338*^9, 3.6137327713906393`*^9}, {3.613732810835518*^9,
1052     3.613732821957967*^9}, {3.613732862763562*^9, 3.613732913141625*^9}, {
1053     3.61373295595298*^9, 3.613733106904414*^9}, {3.613733177027177*^9,
1054     3.613733199740512*^9}, {3.61373326352699*^9, 3.6137332727080936`*^9}, {
1055     3.613733425731758*^9, 3.6137335719815893`*^9}, {3.613733605778666*^9,
1056     3.613733648768111*^9}, {3.613733776819376*^9, 3.6137338150569887`*^9}, {
1057     3.6137338699918118`*^9, 3.61373393057332*^9}, {3.613733985426153*^9,
1058     3.613733985628289*^9}, {3.613734054134171*^9, 3.6137340588514023`*^9}, {
1059     3.6137371511302147`*^9, 3.613737188118826*^9}, {3.613815304031159*^9,
1060     3.613815306899302*^9}, {3.613815415274508*^9, 3.613815450160325*^9}, {
1061     3.613815535717202*^9, 3.613815604127789*^9}, {3.613815665097225*^9,
1062     3.613815673383039*^9}, {3.613816023942026*^9, 3.6138160527230864`*^9}, {
1063     3.613825536655637*^9, 3.613825555245915*^9}, {3.613825660601698*^9,
1064     3.6138256692518787`*^9}, {3.613825968874196*^9, 3.613826000721004*^9}, {
1065     3.6138263115515537`*^9, 3.613826312408416*^9}, {3.613826919678602*^9,
1066     3.613826920319913*^9}, {3.613827178556149*^9, 3.6138272197490597`*^9}, {
1067     3.614594522437893*^9, 3.61459456139437*^9}, 3.623084229697304*^9, {
1068     3.6230843379442043`*^9, 3.623084427365108*^9}, {3.62308446151657*^9,
1069     3.623084468250936*^9}, {3.623084499724618*^9, 3.623084561021247*^9}, {
1070     3.62308463451595*^9, 3.6230847861885157`*^9}, {3.6230849025921383`*^9,
1071     3.6230849294629087`*^9}, {3.623085006652482*^9, 3.623085043440975*^9}, {
1072     3.6230850760671787`*^9, 3.623086190170168*^9}, {3.62308622785679*^9,
1073     3.6230862554078197`*^9}, {3.6230862883261147`*^9, 3.623086432442491*^9}, {
1074     3.623086496576763*^9, 3.623086520303379*^9}, {3.623089619380378*^9,
1075 gezelter 4433 3.623089630706038*^9}, {3.623148156646064*^9, 3.623148193039057*^9}, {
1076     3.676984904725418*^9, 3.676984954321043*^9}}],
1077 gezelter 4432
1078     Cell[CellGroupData[{
1079    
1080     Cell[BoxData[
1081     RowBox[{"Expand", "[", "QBARE", "]"}]], "Input",
1082     CellChangeTimes->{{3.6138263631875963`*^9, 3.613826368246337*^9},
1083     3.6139024316559362`*^9}],
1084    
1085     Cell[BoxData[
1086     RowBox[{"-",
1087     FractionBox[
1088 gezelter 4433 RowBox[{"8", " ",
1089 gezelter 4432 SuperscriptBox["\[ExponentialE]",
1090     RowBox[{
1091     RowBox[{"-",
1092     SuperscriptBox["rc", "2"]}], " ",
1093     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1094     SuperscriptBox["rc", "3"], " ",
1095     SuperscriptBox["\[Kappa]", "5"]}],
1096 gezelter 4433 RowBox[{"15", " ",
1097 gezelter 4432 SqrtBox["\[Pi]"]}]]}]], "Output",
1098     CellChangeTimes->{3.613826419489653*^9, 3.613827358701158*^9,
1099     3.613902416523858*^9, 3.6139057770507097`*^9, 3.613988354332242*^9,
1100     3.614594677651219*^9, 3.623089508432043*^9, 3.623092199109468*^9,
1101 gezelter 4433 3.623150621773239*^9, 3.676988586611196*^9}]
1102 gezelter 4432 }, Open ]],
1103    
1104     Cell[CellGroupData[{
1105    
1106     Cell[BoxData[
1107     RowBox[{"Expand", "[", "QSP", "]"}]], "Input",
1108     CellChangeTimes->{{3.613826369100939*^9, 3.613826374086033*^9},
1109     3.623148196853565*^9}],
1110    
1111     Cell[BoxData[
1112     RowBox[{"-",
1113     FractionBox[
1114     RowBox[{"16", " ",
1115     SuperscriptBox["\[ExponentialE]",
1116     RowBox[{
1117     RowBox[{"-",
1118     SuperscriptBox["rc", "2"]}], " ",
1119     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1120     SuperscriptBox["rc", "5"], " ",
1121     SuperscriptBox["\[Kappa]", "7"]}],
1122 gezelter 4433 RowBox[{"45", " ",
1123 gezelter 4432 SqrtBox["\[Pi]"]}]]}]], "Output",
1124     CellChangeTimes->{3.613988356366749*^9, 3.6145946819255857`*^9,
1125 gezelter 4433 3.623089513390497*^9, 3.6230922026382523`*^9, 3.623150624262762*^9,
1126     3.676988589883754*^9}]
1127 gezelter 4432 }, Open ]],
1128    
1129     Cell[CellGroupData[{
1130    
1131     Cell[BoxData[
1132     RowBox[{"Expand", "[", "QGSF", "]"}]], "Input",
1133     CellChangeTimes->{{3.6138263748828497`*^9, 3.613826377787304*^9}, {
1134     3.613830939552072*^9, 3.613830971126011*^9}, {3.613905875167797*^9,
1135     3.6139058760662394`*^9}}],
1136    
1137     Cell[BoxData[
1138     RowBox[{
1139     FractionBox[
1140     RowBox[{"4", " ",
1141     SuperscriptBox["\[ExponentialE]",
1142     RowBox[{
1143     RowBox[{"-",
1144     SuperscriptBox["rc", "2"]}], " ",
1145     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1146     SuperscriptBox["rc", "5"], " ",
1147     SuperscriptBox["\[Kappa]", "7"]}],
1148 gezelter 4433 RowBox[{"45", " ",
1149 gezelter 4432 SqrtBox["\[Pi]"]}]], "-",
1150     FractionBox[
1151     RowBox[{"8", " ",
1152     SuperscriptBox["\[ExponentialE]",
1153     RowBox[{
1154     RowBox[{"-",
1155     SuperscriptBox["rc", "2"]}], " ",
1156     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1157     SuperscriptBox["rc", "7"], " ",
1158     SuperscriptBox["\[Kappa]", "9"]}],
1159 gezelter 4433 RowBox[{"45", " ",
1160 gezelter 4432 SqrtBox["\[Pi]"]}]]}]], "Output",
1161     CellChangeTimes->{
1162     3.613826419640369*^9, 3.6138270566935883`*^9, 3.6138273666891203`*^9, {
1163     3.613830943408469*^9, 3.6138309716792517`*^9}, 3.613905783149455*^9,
1164     3.613905876406279*^9, 3.613988358482745*^9, 3.614594685371746*^9,
1165 gezelter 4433 3.623089533253147*^9, 3.623092206193384*^9, 3.623150626825612*^9,
1166     3.6769886124964857`*^9}]
1167 gezelter 4432 }, Open ]],
1168    
1169     Cell[CellGroupData[{
1170    
1171     Cell[BoxData[
1172     RowBox[{"Expand", "[", "QTSF", "]"}]], "Input",
1173     CellChangeTimes->{{3.61382637871056*^9, 3.61382638209606*^9}}],
1174    
1175     Cell[BoxData[
1176     RowBox[{
1177     FractionBox[
1178 gezelter 4433 RowBox[{"4", " ",
1179 gezelter 4432 SuperscriptBox["\[ExponentialE]",
1180     RowBox[{
1181     RowBox[{"-",
1182     SuperscriptBox["rc", "2"]}], " ",
1183     SuperscriptBox["\[Kappa]", "2"]}]], " ", "\[Kappa]"}],
1184     RowBox[{
1185     SqrtBox["\[Pi]"], " ", "rc"}]], "+",
1186     FractionBox[
1187 gezelter 4433 RowBox[{"8", " ",
1188 gezelter 4432 SuperscriptBox["\[ExponentialE]",
1189     RowBox[{
1190     RowBox[{"-",
1191     SuperscriptBox["rc", "2"]}], " ",
1192     SuperscriptBox["\[Kappa]", "2"]}]], " ", "rc", " ",
1193     SuperscriptBox["\[Kappa]", "3"]}],
1194 gezelter 4433 RowBox[{"3", " ",
1195     SqrtBox["\[Pi]"]}]], "+",
1196 gezelter 4432 FractionBox[
1197 gezelter 4433 RowBox[{"16", " ",
1198 gezelter 4432 SuperscriptBox["\[ExponentialE]",
1199     RowBox[{
1200     RowBox[{"-",
1201     SuperscriptBox["rc", "2"]}], " ",
1202     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1203     SuperscriptBox["rc", "3"], " ",
1204     SuperscriptBox["\[Kappa]", "5"]}],
1205 gezelter 4433 RowBox[{"15", " ",
1206 gezelter 4432 SqrtBox["\[Pi]"]}]], "+",
1207     FractionBox[
1208 gezelter 4433 RowBox[{"4", " ",
1209 gezelter 4432 SuperscriptBox["\[ExponentialE]",
1210     RowBox[{
1211     RowBox[{"-",
1212     SuperscriptBox["rc", "2"]}], " ",
1213     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1214     SuperscriptBox["rc", "5"], " ",
1215     SuperscriptBox["\[Kappa]", "7"]}],
1216 gezelter 4433 RowBox[{"15", " ",
1217 gezelter 4432 SqrtBox["\[Pi]"]}]], "+",
1218     FractionBox[
1219     RowBox[{"8", " ",
1220     SuperscriptBox["\[ExponentialE]",
1221     RowBox[{
1222     RowBox[{"-",
1223     SuperscriptBox["rc", "2"]}], " ",
1224     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1225     SuperscriptBox["rc", "7"], " ",
1226     SuperscriptBox["\[Kappa]", "9"]}],
1227 gezelter 4433 RowBox[{"45", " ",
1228 gezelter 4432 SqrtBox["\[Pi]"]}]], "+",
1229     FractionBox[
1230 gezelter 4433 RowBox[{"2", " ",
1231 gezelter 4432 RowBox[{"Erfc", "[",
1232     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
1233     SuperscriptBox["rc", "2"]]}]], "Output",
1234     CellChangeTimes->{
1235     3.613826419723484*^9, 3.6138273774842777`*^9, 3.6139057853066607`*^9,
1236     3.613988360907867*^9, 3.614594688339691*^9, {3.6230922148593197`*^9,
1237 gezelter 4433 3.623092222429852*^9}, 3.623150629315596*^9, 3.6769886233095837`*^9}]
1238 gezelter 4432 }, Open ]],
1239    
1240     Cell[CellGroupData[{
1241    
1242     Cell[BoxData[
1243     RowBox[{
1244     RowBox[{"QGSF", "/.",
1245     RowBox[{"rc", "\[Rule]", "11.5"}]}], " ", "/.",
1246     RowBox[{"\[Kappa]", "\[Rule]", "0.2"}]}]], "Input",
1247     CellChangeTimes->{{3.613826607232987*^9, 3.6138266182268*^9}}],
1248    
1249     Cell[BoxData[
1250     RowBox[{"-", "0.05612571955354159`"}]], "Output",
1251     CellChangeTimes->{3.613827062575906*^9, 3.613827398725811*^9,
1252     3.613905789167375*^9, 3.6139075965611973`*^9, 3.613907713580661*^9,
1253     3.6139883631950703`*^9, 3.623092228852077*^9, 3.6231506340948057`*^9}]
1254     }, Open ]],
1255    
1256     Cell[CellGroupData[{
1257    
1258     Cell[BoxData[
1259     RowBox[{"Limit", "[",
1260     RowBox[{"QBARE", ",",
1261     RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
1262     CellChangeTimes->{{3.6145954324899607`*^9, 3.614595438979292*^9}}],
1263    
1264     Cell[BoxData["0"], "Output",
1265     CellChangeTimes->{3.614595439312615*^9, 3.6230922500382757`*^9,
1266     3.623150636706971*^9}]
1267     }, Open ]],
1268    
1269     Cell[CellGroupData[{
1270    
1271     Cell[BoxData[
1272     RowBox[{"Limit", "[",
1273     RowBox[{"QSP", ",",
1274     RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
1275     CellChangeTimes->{{3.614595446545167*^9, 3.61459545529193*^9},
1276     3.623148773969664*^9}],
1277    
1278     Cell[BoxData["0"], "Output",
1279     CellChangeTimes->{3.614595456000258*^9, 3.623150639226756*^9}]
1280     }, Open ]],
1281    
1282     Cell[CellGroupData[{
1283    
1284     Cell[BoxData[
1285     RowBox[{"Limit", "[",
1286     RowBox[{"QGSF", ",",
1287     RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
1288     CellChangeTimes->{{3.614595464309478*^9, 3.614595469692217*^9}}],
1289    
1290     Cell[BoxData["0"], "Output",
1291     CellChangeTimes->{3.6145954699957027`*^9, 3.623092256810233*^9,
1292     3.623150641465064*^9}]
1293     }, Open ]],
1294    
1295     Cell[CellGroupData[{
1296    
1297     Cell[BoxData[
1298     RowBox[{"Limit", "[",
1299     RowBox[{"QTSF", ",",
1300     RowBox[{"\[Kappa]", "\[Rule]", "0"}]}], "]"}]], "Input",
1301     CellChangeTimes->{{3.614595479813127*^9, 3.614595487321061*^9}}],
1302    
1303     Cell[BoxData[
1304     FractionBox["18",
1305     SuperscriptBox["rc", "2"]]], "Output",
1306     CellChangeTimes->{3.614595487675913*^9, 3.623092263603072*^9,
1307     3.623150643794711*^9}]
1308     }, Open ]],
1309    
1310     Cell[CellGroupData[{
1311    
1312     Cell[BoxData["TBare"], "Input",
1313     CellChangeTimes->{{3.614595726057316*^9, 3.6145957369270573`*^9}, {
1314     3.623092292825489*^9, 3.623092304621344*^9}}],
1315    
1316     Cell[BoxData[
1317     TagBox[
1318     RowBox[{"(", "\[NoBreak]", GridBox[{
1319     {
1320     RowBox[{"(", "\[NoBreak]", GridBox[{
1321     {
1322     RowBox[{
1323     RowBox[{"-",
1324     FractionBox["32", "5"]}], " ",
1325     SuperscriptBox["\[ExponentialE]",
1326     RowBox[{
1327     RowBox[{"-",
1328     SuperscriptBox["rc", "2"]}], " ",
1329     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1330     SqrtBox["\[Pi]"], " ",
1331     SuperscriptBox["rc", "3"], " ",
1332     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"},
1333     {"0",
1334     RowBox[{
1335     RowBox[{"-",
1336     FractionBox["32", "15"]}], " ",
1337     SuperscriptBox["\[ExponentialE]",
1338     RowBox[{
1339     RowBox[{"-",
1340     SuperscriptBox["rc", "2"]}], " ",
1341     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1342     SqrtBox["\[Pi]"], " ",
1343     SuperscriptBox["rc", "3"], " ",
1344     SuperscriptBox["\[Kappa]", "5"]}], "0"},
1345     {"0", "0",
1346     RowBox[{
1347     RowBox[{"-",
1348     FractionBox["32", "15"]}], " ",
1349     SuperscriptBox["\[ExponentialE]",
1350     RowBox[{
1351     RowBox[{"-",
1352     SuperscriptBox["rc", "2"]}], " ",
1353     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1354     SqrtBox["\[Pi]"], " ",
1355     SuperscriptBox["rc", "3"], " ",
1356     SuperscriptBox["\[Kappa]", "5"]}]}
1357     },
1358     GridBoxAlignment->{
1359     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1360     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1361     GridBoxSpacings->{"Columns" -> {
1362     Offset[0.27999999999999997`], {
1363     Offset[0.7]},
1364     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1365     "Rows" -> {
1366     Offset[0.2], {
1367     Offset[0.4]},
1368     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1369     RowBox[{"(", "\[NoBreak]", GridBox[{
1370     {"0",
1371     RowBox[{
1372     RowBox[{"-",
1373     FractionBox["32", "15"]}], " ",
1374     SuperscriptBox["\[ExponentialE]",
1375     RowBox[{
1376     RowBox[{"-",
1377     SuperscriptBox["rc", "2"]}], " ",
1378     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1379     SqrtBox["\[Pi]"], " ",
1380     SuperscriptBox["rc", "3"], " ",
1381     SuperscriptBox["\[Kappa]", "5"]}], "0"},
1382     {
1383     RowBox[{
1384     RowBox[{"-",
1385     FractionBox["32", "15"]}], " ",
1386     SuperscriptBox["\[ExponentialE]",
1387     RowBox[{
1388     RowBox[{"-",
1389     SuperscriptBox["rc", "2"]}], " ",
1390     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1391     SqrtBox["\[Pi]"], " ",
1392     SuperscriptBox["rc", "3"], " ",
1393     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"},
1394     {"0", "0", "0"}
1395     },
1396     GridBoxAlignment->{
1397     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1398     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1399     GridBoxSpacings->{"Columns" -> {
1400     Offset[0.27999999999999997`], {
1401     Offset[0.7]},
1402     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1403     "Rows" -> {
1404     Offset[0.2], {
1405     Offset[0.4]},
1406     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1407     RowBox[{"(", "\[NoBreak]", GridBox[{
1408     {"0", "0",
1409     RowBox[{
1410     RowBox[{"-",
1411     FractionBox["32", "15"]}], " ",
1412     SuperscriptBox["\[ExponentialE]",
1413     RowBox[{
1414     RowBox[{"-",
1415     SuperscriptBox["rc", "2"]}], " ",
1416     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1417     SqrtBox["\[Pi]"], " ",
1418     SuperscriptBox["rc", "3"], " ",
1419     SuperscriptBox["\[Kappa]", "5"]}]},
1420     {"0", "0", "0"},
1421     {
1422     RowBox[{
1423     RowBox[{"-",
1424     FractionBox["32", "15"]}], " ",
1425     SuperscriptBox["\[ExponentialE]",
1426     RowBox[{
1427     RowBox[{"-",
1428     SuperscriptBox["rc", "2"]}], " ",
1429     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1430     SqrtBox["\[Pi]"], " ",
1431     SuperscriptBox["rc", "3"], " ",
1432     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"}
1433     },
1434     GridBoxAlignment->{
1435     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1436     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1437     GridBoxSpacings->{"Columns" -> {
1438     Offset[0.27999999999999997`], {
1439     Offset[0.7]},
1440     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1441     "Rows" -> {
1442     Offset[0.2], {
1443     Offset[0.4]},
1444     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
1445     {
1446     RowBox[{"(", "\[NoBreak]", GridBox[{
1447     {"0",
1448     RowBox[{
1449     RowBox[{"-",
1450     FractionBox["32", "15"]}], " ",
1451     SuperscriptBox["\[ExponentialE]",
1452     RowBox[{
1453     RowBox[{"-",
1454     SuperscriptBox["rc", "2"]}], " ",
1455     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1456     SqrtBox["\[Pi]"], " ",
1457     SuperscriptBox["rc", "3"], " ",
1458     SuperscriptBox["\[Kappa]", "5"]}], "0"},
1459     {
1460     RowBox[{
1461     RowBox[{"-",
1462     FractionBox["32", "15"]}], " ",
1463     SuperscriptBox["\[ExponentialE]",
1464     RowBox[{
1465     RowBox[{"-",
1466     SuperscriptBox["rc", "2"]}], " ",
1467     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1468     SqrtBox["\[Pi]"], " ",
1469     SuperscriptBox["rc", "3"], " ",
1470     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"},
1471     {"0", "0", "0"}
1472     },
1473     GridBoxAlignment->{
1474     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1475     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1476     GridBoxSpacings->{"Columns" -> {
1477     Offset[0.27999999999999997`], {
1478     Offset[0.7]},
1479     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1480     "Rows" -> {
1481     Offset[0.2], {
1482     Offset[0.4]},
1483     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1484     RowBox[{"(", "\[NoBreak]", GridBox[{
1485     {
1486     RowBox[{
1487     RowBox[{"-",
1488     FractionBox["32", "15"]}], " ",
1489     SuperscriptBox["\[ExponentialE]",
1490     RowBox[{
1491     RowBox[{"-",
1492     SuperscriptBox["rc", "2"]}], " ",
1493     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1494     SqrtBox["\[Pi]"], " ",
1495     SuperscriptBox["rc", "3"], " ",
1496     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"},
1497     {"0",
1498     RowBox[{
1499     RowBox[{"-",
1500     FractionBox["32", "5"]}], " ",
1501     SuperscriptBox["\[ExponentialE]",
1502     RowBox[{
1503     RowBox[{"-",
1504     SuperscriptBox["rc", "2"]}], " ",
1505     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1506     SqrtBox["\[Pi]"], " ",
1507     SuperscriptBox["rc", "3"], " ",
1508     SuperscriptBox["\[Kappa]", "5"]}], "0"},
1509     {"0", "0",
1510     RowBox[{
1511     RowBox[{"-",
1512     FractionBox["32", "15"]}], " ",
1513     SuperscriptBox["\[ExponentialE]",
1514     RowBox[{
1515     RowBox[{"-",
1516     SuperscriptBox["rc", "2"]}], " ",
1517     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1518     SqrtBox["\[Pi]"], " ",
1519     SuperscriptBox["rc", "3"], " ",
1520     SuperscriptBox["\[Kappa]", "5"]}]}
1521     },
1522     GridBoxAlignment->{
1523     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1524     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1525     GridBoxSpacings->{"Columns" -> {
1526     Offset[0.27999999999999997`], {
1527     Offset[0.7]},
1528     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1529     "Rows" -> {
1530     Offset[0.2], {
1531     Offset[0.4]},
1532     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1533     RowBox[{"(", "\[NoBreak]", GridBox[{
1534     {"0", "0", "0"},
1535     {"0", "0",
1536     RowBox[{
1537     RowBox[{"-",
1538     FractionBox["32", "15"]}], " ",
1539     SuperscriptBox["\[ExponentialE]",
1540     RowBox[{
1541     RowBox[{"-",
1542     SuperscriptBox["rc", "2"]}], " ",
1543     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1544     SqrtBox["\[Pi]"], " ",
1545     SuperscriptBox["rc", "3"], " ",
1546     SuperscriptBox["\[Kappa]", "5"]}]},
1547     {"0",
1548     RowBox[{
1549     RowBox[{"-",
1550     FractionBox["32", "15"]}], " ",
1551     SuperscriptBox["\[ExponentialE]",
1552     RowBox[{
1553     RowBox[{"-",
1554     SuperscriptBox["rc", "2"]}], " ",
1555     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1556     SqrtBox["\[Pi]"], " ",
1557     SuperscriptBox["rc", "3"], " ",
1558     SuperscriptBox["\[Kappa]", "5"]}], "0"}
1559     },
1560     GridBoxAlignment->{
1561     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1562     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1563     GridBoxSpacings->{"Columns" -> {
1564     Offset[0.27999999999999997`], {
1565     Offset[0.7]},
1566     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1567     "Rows" -> {
1568     Offset[0.2], {
1569     Offset[0.4]},
1570     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
1571     {
1572     RowBox[{"(", "\[NoBreak]", GridBox[{
1573     {"0", "0",
1574     RowBox[{
1575     RowBox[{"-",
1576     FractionBox["32", "15"]}], " ",
1577     SuperscriptBox["\[ExponentialE]",
1578     RowBox[{
1579     RowBox[{"-",
1580     SuperscriptBox["rc", "2"]}], " ",
1581     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1582     SqrtBox["\[Pi]"], " ",
1583     SuperscriptBox["rc", "3"], " ",
1584     SuperscriptBox["\[Kappa]", "5"]}]},
1585     {"0", "0", "0"},
1586     {
1587     RowBox[{
1588     RowBox[{"-",
1589     FractionBox["32", "15"]}], " ",
1590     SuperscriptBox["\[ExponentialE]",
1591     RowBox[{
1592     RowBox[{"-",
1593     SuperscriptBox["rc", "2"]}], " ",
1594     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1595     SqrtBox["\[Pi]"], " ",
1596     SuperscriptBox["rc", "3"], " ",
1597     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"}
1598     },
1599     GridBoxAlignment->{
1600     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1601     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1602     GridBoxSpacings->{"Columns" -> {
1603     Offset[0.27999999999999997`], {
1604     Offset[0.7]},
1605     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1606     "Rows" -> {
1607     Offset[0.2], {
1608     Offset[0.4]},
1609     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1610     RowBox[{"(", "\[NoBreak]", GridBox[{
1611     {"0", "0", "0"},
1612     {"0", "0",
1613     RowBox[{
1614     RowBox[{"-",
1615     FractionBox["32", "15"]}], " ",
1616     SuperscriptBox["\[ExponentialE]",
1617     RowBox[{
1618     RowBox[{"-",
1619     SuperscriptBox["rc", "2"]}], " ",
1620     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1621     SqrtBox["\[Pi]"], " ",
1622     SuperscriptBox["rc", "3"], " ",
1623     SuperscriptBox["\[Kappa]", "5"]}]},
1624     {"0",
1625     RowBox[{
1626     RowBox[{"-",
1627     FractionBox["32", "15"]}], " ",
1628     SuperscriptBox["\[ExponentialE]",
1629     RowBox[{
1630     RowBox[{"-",
1631     SuperscriptBox["rc", "2"]}], " ",
1632     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1633     SqrtBox["\[Pi]"], " ",
1634     SuperscriptBox["rc", "3"], " ",
1635     SuperscriptBox["\[Kappa]", "5"]}], "0"}
1636     },
1637     GridBoxAlignment->{
1638     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1639     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1640     GridBoxSpacings->{"Columns" -> {
1641     Offset[0.27999999999999997`], {
1642     Offset[0.7]},
1643     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1644     "Rows" -> {
1645     Offset[0.2], {
1646     Offset[0.4]},
1647     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1648     RowBox[{"(", "\[NoBreak]", GridBox[{
1649     {
1650     RowBox[{
1651     RowBox[{"-",
1652     FractionBox["32", "15"]}], " ",
1653     SuperscriptBox["\[ExponentialE]",
1654     RowBox[{
1655     RowBox[{"-",
1656     SuperscriptBox["rc", "2"]}], " ",
1657     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1658     SqrtBox["\[Pi]"], " ",
1659     SuperscriptBox["rc", "3"], " ",
1660     SuperscriptBox["\[Kappa]", "5"]}], "0", "0"},
1661     {"0",
1662     RowBox[{
1663     RowBox[{"-",
1664     FractionBox["32", "15"]}], " ",
1665     SuperscriptBox["\[ExponentialE]",
1666     RowBox[{
1667     RowBox[{"-",
1668     SuperscriptBox["rc", "2"]}], " ",
1669     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1670     SqrtBox["\[Pi]"], " ",
1671     SuperscriptBox["rc", "3"], " ",
1672     SuperscriptBox["\[Kappa]", "5"]}], "0"},
1673     {"0", "0",
1674     RowBox[{
1675     RowBox[{"-",
1676     FractionBox["32", "5"]}], " ",
1677     SuperscriptBox["\[ExponentialE]",
1678     RowBox[{
1679     RowBox[{"-",
1680     SuperscriptBox["rc", "2"]}], " ",
1681     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1682     SqrtBox["\[Pi]"], " ",
1683     SuperscriptBox["rc", "3"], " ",
1684     SuperscriptBox["\[Kappa]", "5"]}]}
1685     },
1686     GridBoxAlignment->{
1687     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1688     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1689     GridBoxSpacings->{"Columns" -> {
1690     Offset[0.27999999999999997`], {
1691     Offset[0.7]},
1692     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1693     "Rows" -> {
1694     Offset[0.2], {
1695     Offset[0.4]},
1696     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]}
1697     },
1698     GridBoxAlignment->{
1699     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
1700     "RowsIndexed" -> {}},
1701     GridBoxSpacings->{"Columns" -> {
1702     Offset[0.27999999999999997`], {
1703     Offset[0.7]},
1704     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
1705     Offset[0.2], {
1706     Offset[0.4]},
1707     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1708     Function[BoxForm`e$,
1709     MatrixForm[BoxForm`e$]]]], "Output",
1710     CellChangeTimes->{{3.6145957290402822`*^9, 3.614595737188229*^9},
1711     3.614595821337948*^9, {3.6230922980508423`*^9, 3.623092304965407*^9},
1712     3.623150648866949*^9}]
1713     }, Open ]],
1714    
1715     Cell[CellGroupData[{
1716    
1717     Cell[BoxData["TSP"], "Input",
1718     CellChangeTimes->{{3.623150651900529*^9, 3.623150653967461*^9}}],
1719    
1720     Cell[BoxData[
1721     TagBox[
1722     RowBox[{"(", "\[NoBreak]", GridBox[{
1723     {
1724     RowBox[{"(", "\[NoBreak]", GridBox[{
1725     {
1726     RowBox[{
1727     RowBox[{"-",
1728     FractionBox["64", "15"]}], " ",
1729     SuperscriptBox["\[ExponentialE]",
1730     RowBox[{
1731     RowBox[{"-",
1732     SuperscriptBox["rc", "2"]}], " ",
1733     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1734     SqrtBox["\[Pi]"], " ",
1735     SuperscriptBox["rc", "5"], " ",
1736     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"},
1737     {"0",
1738     RowBox[{
1739     RowBox[{"-",
1740     FractionBox["64", "45"]}], " ",
1741     SuperscriptBox["\[ExponentialE]",
1742     RowBox[{
1743     RowBox[{"-",
1744     SuperscriptBox["rc", "2"]}], " ",
1745     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1746     SqrtBox["\[Pi]"], " ",
1747     SuperscriptBox["rc", "5"], " ",
1748     SuperscriptBox["\[Kappa]", "7"]}], "0"},
1749     {"0", "0",
1750     RowBox[{
1751     RowBox[{"-",
1752     FractionBox["64", "45"]}], " ",
1753     SuperscriptBox["\[ExponentialE]",
1754     RowBox[{
1755     RowBox[{"-",
1756     SuperscriptBox["rc", "2"]}], " ",
1757     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1758     SqrtBox["\[Pi]"], " ",
1759     SuperscriptBox["rc", "5"], " ",
1760     SuperscriptBox["\[Kappa]", "7"]}]}
1761     },
1762     GridBoxAlignment->{
1763     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1764     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1765     GridBoxSpacings->{"Columns" -> {
1766     Offset[0.27999999999999997`], {
1767     Offset[0.7]},
1768     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1769     "Rows" -> {
1770     Offset[0.2], {
1771     Offset[0.4]},
1772     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1773     RowBox[{"(", "\[NoBreak]", GridBox[{
1774     {"0",
1775     RowBox[{
1776     RowBox[{"-",
1777     FractionBox["64", "45"]}], " ",
1778     SuperscriptBox["\[ExponentialE]",
1779     RowBox[{
1780     RowBox[{"-",
1781     SuperscriptBox["rc", "2"]}], " ",
1782     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1783     SqrtBox["\[Pi]"], " ",
1784     SuperscriptBox["rc", "5"], " ",
1785     SuperscriptBox["\[Kappa]", "7"]}], "0"},
1786     {
1787     RowBox[{
1788     RowBox[{"-",
1789     FractionBox["64", "45"]}], " ",
1790     SuperscriptBox["\[ExponentialE]",
1791     RowBox[{
1792     RowBox[{"-",
1793     SuperscriptBox["rc", "2"]}], " ",
1794     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1795     SqrtBox["\[Pi]"], " ",
1796     SuperscriptBox["rc", "5"], " ",
1797     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"},
1798     {"0", "0", "0"}
1799     },
1800     GridBoxAlignment->{
1801     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1802     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1803     GridBoxSpacings->{"Columns" -> {
1804     Offset[0.27999999999999997`], {
1805     Offset[0.7]},
1806     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1807     "Rows" -> {
1808     Offset[0.2], {
1809     Offset[0.4]},
1810     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1811     RowBox[{"(", "\[NoBreak]", GridBox[{
1812     {"0", "0",
1813     RowBox[{
1814     RowBox[{"-",
1815     FractionBox["64", "45"]}], " ",
1816     SuperscriptBox["\[ExponentialE]",
1817     RowBox[{
1818     RowBox[{"-",
1819     SuperscriptBox["rc", "2"]}], " ",
1820     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1821     SqrtBox["\[Pi]"], " ",
1822     SuperscriptBox["rc", "5"], " ",
1823     SuperscriptBox["\[Kappa]", "7"]}]},
1824     {"0", "0", "0"},
1825     {
1826     RowBox[{
1827     RowBox[{"-",
1828     FractionBox["64", "45"]}], " ",
1829     SuperscriptBox["\[ExponentialE]",
1830     RowBox[{
1831     RowBox[{"-",
1832     SuperscriptBox["rc", "2"]}], " ",
1833     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1834     SqrtBox["\[Pi]"], " ",
1835     SuperscriptBox["rc", "5"], " ",
1836     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"}
1837     },
1838     GridBoxAlignment->{
1839     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1840     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1841     GridBoxSpacings->{"Columns" -> {
1842     Offset[0.27999999999999997`], {
1843     Offset[0.7]},
1844     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1845     "Rows" -> {
1846     Offset[0.2], {
1847     Offset[0.4]},
1848     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
1849     {
1850     RowBox[{"(", "\[NoBreak]", GridBox[{
1851     {"0",
1852     RowBox[{
1853     RowBox[{"-",
1854     FractionBox["64", "45"]}], " ",
1855     SuperscriptBox["\[ExponentialE]",
1856     RowBox[{
1857     RowBox[{"-",
1858     SuperscriptBox["rc", "2"]}], " ",
1859     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1860     SqrtBox["\[Pi]"], " ",
1861     SuperscriptBox["rc", "5"], " ",
1862     SuperscriptBox["\[Kappa]", "7"]}], "0"},
1863     {
1864     RowBox[{
1865     RowBox[{"-",
1866     FractionBox["64", "45"]}], " ",
1867     SuperscriptBox["\[ExponentialE]",
1868     RowBox[{
1869     RowBox[{"-",
1870     SuperscriptBox["rc", "2"]}], " ",
1871     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1872     SqrtBox["\[Pi]"], " ",
1873     SuperscriptBox["rc", "5"], " ",
1874     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"},
1875     {"0", "0", "0"}
1876     },
1877     GridBoxAlignment->{
1878     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1879     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1880     GridBoxSpacings->{"Columns" -> {
1881     Offset[0.27999999999999997`], {
1882     Offset[0.7]},
1883     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1884     "Rows" -> {
1885     Offset[0.2], {
1886     Offset[0.4]},
1887     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1888     RowBox[{"(", "\[NoBreak]", GridBox[{
1889     {
1890     RowBox[{
1891     RowBox[{"-",
1892     FractionBox["64", "45"]}], " ",
1893     SuperscriptBox["\[ExponentialE]",
1894     RowBox[{
1895     RowBox[{"-",
1896     SuperscriptBox["rc", "2"]}], " ",
1897     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1898     SqrtBox["\[Pi]"], " ",
1899     SuperscriptBox["rc", "5"], " ",
1900     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"},
1901     {"0",
1902     RowBox[{
1903     RowBox[{"-",
1904     FractionBox["64", "15"]}], " ",
1905     SuperscriptBox["\[ExponentialE]",
1906     RowBox[{
1907     RowBox[{"-",
1908     SuperscriptBox["rc", "2"]}], " ",
1909     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1910     SqrtBox["\[Pi]"], " ",
1911     SuperscriptBox["rc", "5"], " ",
1912     SuperscriptBox["\[Kappa]", "7"]}], "0"},
1913     {"0", "0",
1914     RowBox[{
1915     RowBox[{"-",
1916     FractionBox["64", "45"]}], " ",
1917     SuperscriptBox["\[ExponentialE]",
1918     RowBox[{
1919     RowBox[{"-",
1920     SuperscriptBox["rc", "2"]}], " ",
1921     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1922     SqrtBox["\[Pi]"], " ",
1923     SuperscriptBox["rc", "5"], " ",
1924     SuperscriptBox["\[Kappa]", "7"]}]}
1925     },
1926     GridBoxAlignment->{
1927     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1928     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1929     GridBoxSpacings->{"Columns" -> {
1930     Offset[0.27999999999999997`], {
1931     Offset[0.7]},
1932     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1933     "Rows" -> {
1934     Offset[0.2], {
1935     Offset[0.4]},
1936     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
1937     RowBox[{"(", "\[NoBreak]", GridBox[{
1938     {"0", "0", "0"},
1939     {"0", "0",
1940     RowBox[{
1941     RowBox[{"-",
1942     FractionBox["64", "45"]}], " ",
1943     SuperscriptBox["\[ExponentialE]",
1944     RowBox[{
1945     RowBox[{"-",
1946     SuperscriptBox["rc", "2"]}], " ",
1947     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1948     SqrtBox["\[Pi]"], " ",
1949     SuperscriptBox["rc", "5"], " ",
1950     SuperscriptBox["\[Kappa]", "7"]}]},
1951     {"0",
1952     RowBox[{
1953     RowBox[{"-",
1954     FractionBox["64", "45"]}], " ",
1955     SuperscriptBox["\[ExponentialE]",
1956     RowBox[{
1957     RowBox[{"-",
1958     SuperscriptBox["rc", "2"]}], " ",
1959     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1960     SqrtBox["\[Pi]"], " ",
1961     SuperscriptBox["rc", "5"], " ",
1962     SuperscriptBox["\[Kappa]", "7"]}], "0"}
1963     },
1964     GridBoxAlignment->{
1965     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
1966     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
1967     GridBoxSpacings->{"Columns" -> {
1968     Offset[0.27999999999999997`], {
1969     Offset[0.7]},
1970     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
1971     "Rows" -> {
1972     Offset[0.2], {
1973     Offset[0.4]},
1974     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
1975     {
1976     RowBox[{"(", "\[NoBreak]", GridBox[{
1977     {"0", "0",
1978     RowBox[{
1979     RowBox[{"-",
1980     FractionBox["64", "45"]}], " ",
1981     SuperscriptBox["\[ExponentialE]",
1982     RowBox[{
1983     RowBox[{"-",
1984     SuperscriptBox["rc", "2"]}], " ",
1985     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1986     SqrtBox["\[Pi]"], " ",
1987     SuperscriptBox["rc", "5"], " ",
1988     SuperscriptBox["\[Kappa]", "7"]}]},
1989     {"0", "0", "0"},
1990     {
1991     RowBox[{
1992     RowBox[{"-",
1993     FractionBox["64", "45"]}], " ",
1994     SuperscriptBox["\[ExponentialE]",
1995     RowBox[{
1996     RowBox[{"-",
1997     SuperscriptBox["rc", "2"]}], " ",
1998     SuperscriptBox["\[Kappa]", "2"]}]], " ",
1999     SqrtBox["\[Pi]"], " ",
2000     SuperscriptBox["rc", "5"], " ",
2001     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"}
2002     },
2003     GridBoxAlignment->{
2004     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2005     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2006     GridBoxSpacings->{"Columns" -> {
2007     Offset[0.27999999999999997`], {
2008     Offset[0.7]},
2009     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2010     "Rows" -> {
2011     Offset[0.2], {
2012     Offset[0.4]},
2013     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2014     RowBox[{"(", "\[NoBreak]", GridBox[{
2015     {"0", "0", "0"},
2016     {"0", "0",
2017     RowBox[{
2018     RowBox[{"-",
2019     FractionBox["64", "45"]}], " ",
2020     SuperscriptBox["\[ExponentialE]",
2021     RowBox[{
2022     RowBox[{"-",
2023     SuperscriptBox["rc", "2"]}], " ",
2024     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2025     SqrtBox["\[Pi]"], " ",
2026     SuperscriptBox["rc", "5"], " ",
2027     SuperscriptBox["\[Kappa]", "7"]}]},
2028     {"0",
2029     RowBox[{
2030     RowBox[{"-",
2031     FractionBox["64", "45"]}], " ",
2032     SuperscriptBox["\[ExponentialE]",
2033     RowBox[{
2034     RowBox[{"-",
2035     SuperscriptBox["rc", "2"]}], " ",
2036     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2037     SqrtBox["\[Pi]"], " ",
2038     SuperscriptBox["rc", "5"], " ",
2039     SuperscriptBox["\[Kappa]", "7"]}], "0"}
2040     },
2041     GridBoxAlignment->{
2042     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2043     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2044     GridBoxSpacings->{"Columns" -> {
2045     Offset[0.27999999999999997`], {
2046     Offset[0.7]},
2047     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2048     "Rows" -> {
2049     Offset[0.2], {
2050     Offset[0.4]},
2051     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2052     RowBox[{"(", "\[NoBreak]", GridBox[{
2053     {
2054     RowBox[{
2055     RowBox[{"-",
2056     FractionBox["64", "45"]}], " ",
2057     SuperscriptBox["\[ExponentialE]",
2058     RowBox[{
2059     RowBox[{"-",
2060     SuperscriptBox["rc", "2"]}], " ",
2061     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2062     SqrtBox["\[Pi]"], " ",
2063     SuperscriptBox["rc", "5"], " ",
2064     SuperscriptBox["\[Kappa]", "7"]}], "0", "0"},
2065     {"0",
2066     RowBox[{
2067     RowBox[{"-",
2068     FractionBox["64", "45"]}], " ",
2069     SuperscriptBox["\[ExponentialE]",
2070     RowBox[{
2071     RowBox[{"-",
2072     SuperscriptBox["rc", "2"]}], " ",
2073     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2074     SqrtBox["\[Pi]"], " ",
2075     SuperscriptBox["rc", "5"], " ",
2076     SuperscriptBox["\[Kappa]", "7"]}], "0"},
2077     {"0", "0",
2078     RowBox[{
2079     RowBox[{"-",
2080     FractionBox["64", "15"]}], " ",
2081     SuperscriptBox["\[ExponentialE]",
2082     RowBox[{
2083     RowBox[{"-",
2084     SuperscriptBox["rc", "2"]}], " ",
2085     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2086     SqrtBox["\[Pi]"], " ",
2087     SuperscriptBox["rc", "5"], " ",
2088     SuperscriptBox["\[Kappa]", "7"]}]}
2089     },
2090     GridBoxAlignment->{
2091     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2092     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2093     GridBoxSpacings->{"Columns" -> {
2094     Offset[0.27999999999999997`], {
2095     Offset[0.7]},
2096     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2097     "Rows" -> {
2098     Offset[0.2], {
2099     Offset[0.4]},
2100     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]}
2101     },
2102     GridBoxAlignment->{
2103     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
2104     "RowsIndexed" -> {}},
2105     GridBoxSpacings->{"Columns" -> {
2106     Offset[0.27999999999999997`], {
2107     Offset[0.7]},
2108     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
2109     Offset[0.2], {
2110     Offset[0.4]},
2111     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2112     Function[BoxForm`e$,
2113     MatrixForm[BoxForm`e$]]]], "Output",
2114     CellChangeTimes->{3.6231506551186953`*^9}]
2115     }, Open ]],
2116    
2117     Cell[CellGroupData[{
2118    
2119     Cell[BoxData["TGSF"], "Input",
2120     CellChangeTimes->{{3.623150658203018*^9, 3.6231506595127583`*^9}}],
2121    
2122     Cell[BoxData[
2123     TagBox[
2124     RowBox[{"(", "\[NoBreak]", GridBox[{
2125     {
2126     RowBox[{"(", "\[NoBreak]", GridBox[{
2127     {
2128     RowBox[{
2129     RowBox[{"-",
2130     FractionBox["16", "15"]}], " ",
2131     SuperscriptBox["\[ExponentialE]",
2132     RowBox[{
2133     RowBox[{"-",
2134     SuperscriptBox["rc", "2"]}], " ",
2135     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2136     SqrtBox["\[Pi]"], " ",
2137     SuperscriptBox["rc", "5"], " ",
2138     SuperscriptBox["\[Kappa]", "7"], " ",
2139     RowBox[{"(",
2140     RowBox[{
2141     RowBox[{"-", "1"}], "+",
2142     RowBox[{"2", " ",
2143     SuperscriptBox["rc", "2"], " ",
2144     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"},
2145     {"0",
2146     RowBox[{
2147     RowBox[{"-",
2148     FractionBox["16", "45"]}], " ",
2149     SuperscriptBox["\[ExponentialE]",
2150     RowBox[{
2151     RowBox[{"-",
2152     SuperscriptBox["rc", "2"]}], " ",
2153     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2154     SqrtBox["\[Pi]"], " ",
2155     SuperscriptBox["rc", "5"], " ",
2156     SuperscriptBox["\[Kappa]", "7"], " ",
2157     RowBox[{"(",
2158     RowBox[{
2159     RowBox[{"-", "1"}], "+",
2160     RowBox[{"2", " ",
2161     SuperscriptBox["rc", "2"], " ",
2162     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"},
2163     {"0", "0",
2164     RowBox[{
2165     RowBox[{"-",
2166     FractionBox["16", "45"]}], " ",
2167     SuperscriptBox["\[ExponentialE]",
2168     RowBox[{
2169     RowBox[{"-",
2170     SuperscriptBox["rc", "2"]}], " ",
2171     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2172     SqrtBox["\[Pi]"], " ",
2173     SuperscriptBox["rc", "5"], " ",
2174     SuperscriptBox["\[Kappa]", "7"], " ",
2175     RowBox[{"(",
2176     RowBox[{
2177     RowBox[{"-", "1"}], "+",
2178     RowBox[{"2", " ",
2179     SuperscriptBox["rc", "2"], " ",
2180     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]}
2181     },
2182     GridBoxAlignment->{
2183     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2184     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2185     GridBoxSpacings->{"Columns" -> {
2186     Offset[0.27999999999999997`], {
2187     Offset[0.7]},
2188     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2189     "Rows" -> {
2190     Offset[0.2], {
2191     Offset[0.4]},
2192     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2193     RowBox[{"(", "\[NoBreak]", GridBox[{
2194     {"0",
2195     RowBox[{
2196     RowBox[{"-",
2197     FractionBox["16", "45"]}], " ",
2198     SuperscriptBox["\[ExponentialE]",
2199     RowBox[{
2200     RowBox[{"-",
2201     SuperscriptBox["rc", "2"]}], " ",
2202     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2203     SqrtBox["\[Pi]"], " ",
2204     SuperscriptBox["rc", "5"], " ",
2205     SuperscriptBox["\[Kappa]", "7"], " ",
2206     RowBox[{"(",
2207     RowBox[{
2208     RowBox[{"-", "1"}], "+",
2209     RowBox[{"2", " ",
2210     SuperscriptBox["rc", "2"], " ",
2211     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"},
2212     {
2213     RowBox[{
2214     RowBox[{"-",
2215     FractionBox["16", "45"]}], " ",
2216     SuperscriptBox["\[ExponentialE]",
2217     RowBox[{
2218     RowBox[{"-",
2219     SuperscriptBox["rc", "2"]}], " ",
2220     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2221     SqrtBox["\[Pi]"], " ",
2222     SuperscriptBox["rc", "5"], " ",
2223     SuperscriptBox["\[Kappa]", "7"], " ",
2224     RowBox[{"(",
2225     RowBox[{
2226     RowBox[{"-", "1"}], "+",
2227     RowBox[{"2", " ",
2228     SuperscriptBox["rc", "2"], " ",
2229     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"},
2230     {"0", "0", "0"}
2231     },
2232     GridBoxAlignment->{
2233     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2234     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2235     GridBoxSpacings->{"Columns" -> {
2236     Offset[0.27999999999999997`], {
2237     Offset[0.7]},
2238     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2239     "Rows" -> {
2240     Offset[0.2], {
2241     Offset[0.4]},
2242     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2243     RowBox[{"(", "\[NoBreak]", GridBox[{
2244     {"0", "0",
2245     RowBox[{
2246     RowBox[{"-",
2247     FractionBox["16", "45"]}], " ",
2248     SuperscriptBox["\[ExponentialE]",
2249     RowBox[{
2250     RowBox[{"-",
2251     SuperscriptBox["rc", "2"]}], " ",
2252     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2253     SqrtBox["\[Pi]"], " ",
2254     SuperscriptBox["rc", "5"], " ",
2255     SuperscriptBox["\[Kappa]", "7"], " ",
2256     RowBox[{"(",
2257     RowBox[{
2258     RowBox[{"-", "1"}], "+",
2259     RowBox[{"2", " ",
2260     SuperscriptBox["rc", "2"], " ",
2261     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]},
2262     {"0", "0", "0"},
2263     {
2264     RowBox[{
2265     RowBox[{"-",
2266     FractionBox["16", "45"]}], " ",
2267     SuperscriptBox["\[ExponentialE]",
2268     RowBox[{
2269     RowBox[{"-",
2270     SuperscriptBox["rc", "2"]}], " ",
2271     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2272     SqrtBox["\[Pi]"], " ",
2273     SuperscriptBox["rc", "5"], " ",
2274     SuperscriptBox["\[Kappa]", "7"], " ",
2275     RowBox[{"(",
2276     RowBox[{
2277     RowBox[{"-", "1"}], "+",
2278     RowBox[{"2", " ",
2279     SuperscriptBox["rc", "2"], " ",
2280     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"}
2281     },
2282     GridBoxAlignment->{
2283     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2284     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2285     GridBoxSpacings->{"Columns" -> {
2286     Offset[0.27999999999999997`], {
2287     Offset[0.7]},
2288     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2289     "Rows" -> {
2290     Offset[0.2], {
2291     Offset[0.4]},
2292     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
2293     {
2294     RowBox[{"(", "\[NoBreak]", GridBox[{
2295     {"0",
2296     RowBox[{
2297     RowBox[{"-",
2298     FractionBox["16", "45"]}], " ",
2299     SuperscriptBox["\[ExponentialE]",
2300     RowBox[{
2301     RowBox[{"-",
2302     SuperscriptBox["rc", "2"]}], " ",
2303     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2304     SqrtBox["\[Pi]"], " ",
2305     SuperscriptBox["rc", "5"], " ",
2306     SuperscriptBox["\[Kappa]", "7"], " ",
2307     RowBox[{"(",
2308     RowBox[{
2309     RowBox[{"-", "1"}], "+",
2310     RowBox[{"2", " ",
2311     SuperscriptBox["rc", "2"], " ",
2312     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"},
2313     {
2314     RowBox[{
2315     RowBox[{"-",
2316     FractionBox["16", "45"]}], " ",
2317     SuperscriptBox["\[ExponentialE]",
2318     RowBox[{
2319     RowBox[{"-",
2320     SuperscriptBox["rc", "2"]}], " ",
2321     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2322     SqrtBox["\[Pi]"], " ",
2323     SuperscriptBox["rc", "5"], " ",
2324     SuperscriptBox["\[Kappa]", "7"], " ",
2325     RowBox[{"(",
2326     RowBox[{
2327     RowBox[{"-", "1"}], "+",
2328     RowBox[{"2", " ",
2329     SuperscriptBox["rc", "2"], " ",
2330     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"},
2331     {"0", "0", "0"}
2332     },
2333     GridBoxAlignment->{
2334     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2335     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2336     GridBoxSpacings->{"Columns" -> {
2337     Offset[0.27999999999999997`], {
2338     Offset[0.7]},
2339     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2340     "Rows" -> {
2341     Offset[0.2], {
2342     Offset[0.4]},
2343     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2344     RowBox[{"(", "\[NoBreak]", GridBox[{
2345     {
2346     RowBox[{
2347     RowBox[{"-",
2348     FractionBox["16", "45"]}], " ",
2349     SuperscriptBox["\[ExponentialE]",
2350     RowBox[{
2351     RowBox[{"-",
2352     SuperscriptBox["rc", "2"]}], " ",
2353     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2354     SqrtBox["\[Pi]"], " ",
2355     SuperscriptBox["rc", "5"], " ",
2356     SuperscriptBox["\[Kappa]", "7"], " ",
2357     RowBox[{"(",
2358     RowBox[{
2359     RowBox[{"-", "1"}], "+",
2360     RowBox[{"2", " ",
2361     SuperscriptBox["rc", "2"], " ",
2362     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"},
2363     {"0",
2364     RowBox[{
2365     RowBox[{"-",
2366     FractionBox["16", "15"]}], " ",
2367     SuperscriptBox["\[ExponentialE]",
2368     RowBox[{
2369     RowBox[{"-",
2370     SuperscriptBox["rc", "2"]}], " ",
2371     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2372     SqrtBox["\[Pi]"], " ",
2373     SuperscriptBox["rc", "5"], " ",
2374     SuperscriptBox["\[Kappa]", "7"], " ",
2375     RowBox[{"(",
2376     RowBox[{
2377     RowBox[{"-", "1"}], "+",
2378     RowBox[{"2", " ",
2379     SuperscriptBox["rc", "2"], " ",
2380     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"},
2381     {"0", "0",
2382     RowBox[{
2383     RowBox[{"-",
2384     FractionBox["16", "45"]}], " ",
2385     SuperscriptBox["\[ExponentialE]",
2386     RowBox[{
2387     RowBox[{"-",
2388     SuperscriptBox["rc", "2"]}], " ",
2389     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2390     SqrtBox["\[Pi]"], " ",
2391     SuperscriptBox["rc", "5"], " ",
2392     SuperscriptBox["\[Kappa]", "7"], " ",
2393     RowBox[{"(",
2394     RowBox[{
2395     RowBox[{"-", "1"}], "+",
2396     RowBox[{"2", " ",
2397     SuperscriptBox["rc", "2"], " ",
2398     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]}
2399     },
2400     GridBoxAlignment->{
2401     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2402     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2403     GridBoxSpacings->{"Columns" -> {
2404     Offset[0.27999999999999997`], {
2405     Offset[0.7]},
2406     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2407     "Rows" -> {
2408     Offset[0.2], {
2409     Offset[0.4]},
2410     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2411     RowBox[{"(", "\[NoBreak]", GridBox[{
2412     {"0", "0", "0"},
2413     {"0", "0",
2414     RowBox[{
2415     RowBox[{"-",
2416     FractionBox["16", "45"]}], " ",
2417     SuperscriptBox["\[ExponentialE]",
2418     RowBox[{
2419     RowBox[{"-",
2420     SuperscriptBox["rc", "2"]}], " ",
2421     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2422     SqrtBox["\[Pi]"], " ",
2423     SuperscriptBox["rc", "5"], " ",
2424     SuperscriptBox["\[Kappa]", "7"], " ",
2425     RowBox[{"(",
2426     RowBox[{
2427     RowBox[{"-", "1"}], "+",
2428     RowBox[{"2", " ",
2429     SuperscriptBox["rc", "2"], " ",
2430     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]},
2431     {"0",
2432     RowBox[{
2433     RowBox[{"-",
2434     FractionBox["16", "45"]}], " ",
2435     SuperscriptBox["\[ExponentialE]",
2436     RowBox[{
2437     RowBox[{"-",
2438     SuperscriptBox["rc", "2"]}], " ",
2439     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2440     SqrtBox["\[Pi]"], " ",
2441     SuperscriptBox["rc", "5"], " ",
2442     SuperscriptBox["\[Kappa]", "7"], " ",
2443     RowBox[{"(",
2444     RowBox[{
2445     RowBox[{"-", "1"}], "+",
2446     RowBox[{"2", " ",
2447     SuperscriptBox["rc", "2"], " ",
2448     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"}
2449     },
2450     GridBoxAlignment->{
2451     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2452     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2453     GridBoxSpacings->{"Columns" -> {
2454     Offset[0.27999999999999997`], {
2455     Offset[0.7]},
2456     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2457     "Rows" -> {
2458     Offset[0.2], {
2459     Offset[0.4]},
2460     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
2461     {
2462     RowBox[{"(", "\[NoBreak]", GridBox[{
2463     {"0", "0",
2464     RowBox[{
2465     RowBox[{"-",
2466     FractionBox["16", "45"]}], " ",
2467     SuperscriptBox["\[ExponentialE]",
2468     RowBox[{
2469     RowBox[{"-",
2470     SuperscriptBox["rc", "2"]}], " ",
2471     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2472     SqrtBox["\[Pi]"], " ",
2473     SuperscriptBox["rc", "5"], " ",
2474     SuperscriptBox["\[Kappa]", "7"], " ",
2475     RowBox[{"(",
2476     RowBox[{
2477     RowBox[{"-", "1"}], "+",
2478     RowBox[{"2", " ",
2479     SuperscriptBox["rc", "2"], " ",
2480     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]},
2481     {"0", "0", "0"},
2482     {
2483     RowBox[{
2484     RowBox[{"-",
2485     FractionBox["16", "45"]}], " ",
2486     SuperscriptBox["\[ExponentialE]",
2487     RowBox[{
2488     RowBox[{"-",
2489     SuperscriptBox["rc", "2"]}], " ",
2490     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2491     SqrtBox["\[Pi]"], " ",
2492     SuperscriptBox["rc", "5"], " ",
2493     SuperscriptBox["\[Kappa]", "7"], " ",
2494     RowBox[{"(",
2495     RowBox[{
2496     RowBox[{"-", "1"}], "+",
2497     RowBox[{"2", " ",
2498     SuperscriptBox["rc", "2"], " ",
2499     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"}
2500     },
2501     GridBoxAlignment->{
2502     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2503     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2504     GridBoxSpacings->{"Columns" -> {
2505     Offset[0.27999999999999997`], {
2506     Offset[0.7]},
2507     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2508     "Rows" -> {
2509     Offset[0.2], {
2510     Offset[0.4]},
2511     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2512     RowBox[{"(", "\[NoBreak]", GridBox[{
2513     {"0", "0", "0"},
2514     {"0", "0",
2515     RowBox[{
2516     RowBox[{"-",
2517     FractionBox["16", "45"]}], " ",
2518     SuperscriptBox["\[ExponentialE]",
2519     RowBox[{
2520     RowBox[{"-",
2521     SuperscriptBox["rc", "2"]}], " ",
2522     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2523     SqrtBox["\[Pi]"], " ",
2524     SuperscriptBox["rc", "5"], " ",
2525     SuperscriptBox["\[Kappa]", "7"], " ",
2526     RowBox[{"(",
2527     RowBox[{
2528     RowBox[{"-", "1"}], "+",
2529     RowBox[{"2", " ",
2530     SuperscriptBox["rc", "2"], " ",
2531     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]},
2532     {"0",
2533     RowBox[{
2534     RowBox[{"-",
2535     FractionBox["16", "45"]}], " ",
2536     SuperscriptBox["\[ExponentialE]",
2537     RowBox[{
2538     RowBox[{"-",
2539     SuperscriptBox["rc", "2"]}], " ",
2540     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2541     SqrtBox["\[Pi]"], " ",
2542     SuperscriptBox["rc", "5"], " ",
2543     SuperscriptBox["\[Kappa]", "7"], " ",
2544     RowBox[{"(",
2545     RowBox[{
2546     RowBox[{"-", "1"}], "+",
2547     RowBox[{"2", " ",
2548     SuperscriptBox["rc", "2"], " ",
2549     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"}
2550     },
2551     GridBoxAlignment->{
2552     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2553     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2554     GridBoxSpacings->{"Columns" -> {
2555     Offset[0.27999999999999997`], {
2556     Offset[0.7]},
2557     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2558     "Rows" -> {
2559     Offset[0.2], {
2560     Offset[0.4]},
2561     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2562     RowBox[{"(", "\[NoBreak]", GridBox[{
2563     {
2564     RowBox[{
2565     RowBox[{"-",
2566     FractionBox["16", "45"]}], " ",
2567     SuperscriptBox["\[ExponentialE]",
2568     RowBox[{
2569     RowBox[{"-",
2570     SuperscriptBox["rc", "2"]}], " ",
2571     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2572     SqrtBox["\[Pi]"], " ",
2573     SuperscriptBox["rc", "5"], " ",
2574     SuperscriptBox["\[Kappa]", "7"], " ",
2575     RowBox[{"(",
2576     RowBox[{
2577     RowBox[{"-", "1"}], "+",
2578     RowBox[{"2", " ",
2579     SuperscriptBox["rc", "2"], " ",
2580     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0", "0"},
2581     {"0",
2582     RowBox[{
2583     RowBox[{"-",
2584     FractionBox["16", "45"]}], " ",
2585     SuperscriptBox["\[ExponentialE]",
2586     RowBox[{
2587     RowBox[{"-",
2588     SuperscriptBox["rc", "2"]}], " ",
2589     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2590     SqrtBox["\[Pi]"], " ",
2591     SuperscriptBox["rc", "5"], " ",
2592     SuperscriptBox["\[Kappa]", "7"], " ",
2593     RowBox[{"(",
2594     RowBox[{
2595     RowBox[{"-", "1"}], "+",
2596     RowBox[{"2", " ",
2597     SuperscriptBox["rc", "2"], " ",
2598     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}], "0"},
2599     {"0", "0",
2600     RowBox[{
2601     RowBox[{"-",
2602     FractionBox["16", "15"]}], " ",
2603     SuperscriptBox["\[ExponentialE]",
2604     RowBox[{
2605     RowBox[{"-",
2606     SuperscriptBox["rc", "2"]}], " ",
2607     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2608     SqrtBox["\[Pi]"], " ",
2609     SuperscriptBox["rc", "5"], " ",
2610     SuperscriptBox["\[Kappa]", "7"], " ",
2611     RowBox[{"(",
2612     RowBox[{
2613     RowBox[{"-", "1"}], "+",
2614     RowBox[{"2", " ",
2615     SuperscriptBox["rc", "2"], " ",
2616     SuperscriptBox["\[Kappa]", "2"]}]}], ")"}]}]}
2617     },
2618     GridBoxAlignment->{
2619     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2620     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2621     GridBoxSpacings->{"Columns" -> {
2622     Offset[0.27999999999999997`], {
2623     Offset[0.7]},
2624     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2625     "Rows" -> {
2626     Offset[0.2], {
2627     Offset[0.4]},
2628     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]}
2629     },
2630     GridBoxAlignment->{
2631     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
2632     "RowsIndexed" -> {}},
2633     GridBoxSpacings->{"Columns" -> {
2634     Offset[0.27999999999999997`], {
2635     Offset[0.7]},
2636     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
2637     Offset[0.2], {
2638     Offset[0.4]},
2639     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2640     Function[BoxForm`e$,
2641     MatrixForm[BoxForm`e$]]]], "Output",
2642     CellChangeTimes->{3.623150659866877*^9}]
2643     }, Open ]],
2644    
2645     Cell[CellGroupData[{
2646    
2647     Cell[BoxData["TTSF"], "Input",
2648     CellChangeTimes->{{3.623150661672976*^9, 3.6231506628875427`*^9}}],
2649    
2650     Cell[BoxData[
2651     TagBox[
2652     RowBox[{"(", "\[NoBreak]", GridBox[{
2653     {
2654     RowBox[{"(", "\[NoBreak]", GridBox[{
2655     {
2656     RowBox[{
2657     FractionBox[
2658     RowBox[{"16", " ",
2659     SuperscriptBox["\[ExponentialE]",
2660     RowBox[{
2661     RowBox[{"-",
2662     SuperscriptBox["rc", "2"]}], " ",
2663     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2664     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2665     RowBox[{"(",
2666     RowBox[{"45", "+",
2667     RowBox[{"30", " ",
2668     SuperscriptBox["rc", "2"], " ",
2669     SuperscriptBox["\[Kappa]", "2"]}], "+",
2670     RowBox[{"12", " ",
2671     SuperscriptBox["rc", "4"], " ",
2672     SuperscriptBox["\[Kappa]", "4"]}], "+",
2673     RowBox[{"3", " ",
2674     SuperscriptBox["rc", "6"], " ",
2675     SuperscriptBox["\[Kappa]", "6"]}], "+",
2676     RowBox[{"2", " ",
2677     SuperscriptBox["rc", "8"], " ",
2678     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2679     RowBox[{"15", " ", "rc"}]], "+",
2680     FractionBox[
2681     RowBox[{"24", " ", "\[Pi]", " ",
2682     RowBox[{"Erfc", "[",
2683     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2684     SuperscriptBox["rc", "2"]]}], "0", "0"},
2685     {"0",
2686     RowBox[{
2687     FractionBox[
2688     RowBox[{"16", " ",
2689     SuperscriptBox["\[ExponentialE]",
2690     RowBox[{
2691     RowBox[{"-",
2692     SuperscriptBox["rc", "2"]}], " ",
2693     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2694     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2695     RowBox[{"(",
2696     RowBox[{"45", "+",
2697     RowBox[{"30", " ",
2698     SuperscriptBox["rc", "2"], " ",
2699     SuperscriptBox["\[Kappa]", "2"]}], "+",
2700     RowBox[{"12", " ",
2701     SuperscriptBox["rc", "4"], " ",
2702     SuperscriptBox["\[Kappa]", "4"]}], "+",
2703     RowBox[{"3", " ",
2704     SuperscriptBox["rc", "6"], " ",
2705     SuperscriptBox["\[Kappa]", "6"]}], "+",
2706     RowBox[{"2", " ",
2707     SuperscriptBox["rc", "8"], " ",
2708     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2709     RowBox[{"45", " ", "rc"}]], "+",
2710     FractionBox[
2711     RowBox[{"8", " ", "\[Pi]", " ",
2712     RowBox[{"Erfc", "[",
2713     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2714     SuperscriptBox["rc", "2"]]}], "0"},
2715     {"0", "0",
2716     RowBox[{
2717     FractionBox[
2718     RowBox[{"16", " ",
2719     SuperscriptBox["\[ExponentialE]",
2720     RowBox[{
2721     RowBox[{"-",
2722     SuperscriptBox["rc", "2"]}], " ",
2723     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2724     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2725     RowBox[{"(",
2726     RowBox[{"45", "+",
2727     RowBox[{"30", " ",
2728     SuperscriptBox["rc", "2"], " ",
2729     SuperscriptBox["\[Kappa]", "2"]}], "+",
2730     RowBox[{"12", " ",
2731     SuperscriptBox["rc", "4"], " ",
2732     SuperscriptBox["\[Kappa]", "4"]}], "+",
2733     RowBox[{"3", " ",
2734     SuperscriptBox["rc", "6"], " ",
2735     SuperscriptBox["\[Kappa]", "6"]}], "+",
2736     RowBox[{"2", " ",
2737     SuperscriptBox["rc", "8"], " ",
2738     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2739     RowBox[{"45", " ", "rc"}]], "+",
2740     FractionBox[
2741     RowBox[{"8", " ", "\[Pi]", " ",
2742     RowBox[{"Erfc", "[",
2743     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2744     SuperscriptBox["rc", "2"]]}]}
2745     },
2746     GridBoxAlignment->{
2747     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2748     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2749     GridBoxSpacings->{"Columns" -> {
2750     Offset[0.27999999999999997`], {
2751     Offset[0.7]},
2752     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2753     "Rows" -> {
2754     Offset[0.2], {
2755     Offset[0.4]},
2756     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2757     RowBox[{"(", "\[NoBreak]", GridBox[{
2758     {"0",
2759     RowBox[{
2760     FractionBox[
2761     RowBox[{"16", " ",
2762     SuperscriptBox["\[ExponentialE]",
2763     RowBox[{
2764     RowBox[{"-",
2765     SuperscriptBox["rc", "2"]}], " ",
2766     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2767     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2768     RowBox[{"(",
2769     RowBox[{"45", "+",
2770     RowBox[{"30", " ",
2771     SuperscriptBox["rc", "2"], " ",
2772     SuperscriptBox["\[Kappa]", "2"]}], "+",
2773     RowBox[{"12", " ",
2774     SuperscriptBox["rc", "4"], " ",
2775     SuperscriptBox["\[Kappa]", "4"]}], "+",
2776     RowBox[{"3", " ",
2777     SuperscriptBox["rc", "6"], " ",
2778     SuperscriptBox["\[Kappa]", "6"]}], "+",
2779     RowBox[{"2", " ",
2780     SuperscriptBox["rc", "8"], " ",
2781     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2782     RowBox[{"45", " ", "rc"}]], "+",
2783     FractionBox[
2784     RowBox[{"8", " ", "\[Pi]", " ",
2785     RowBox[{"Erfc", "[",
2786     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2787     SuperscriptBox["rc", "2"]]}], "0"},
2788     {
2789     RowBox[{
2790     FractionBox[
2791     RowBox[{"16", " ",
2792     SuperscriptBox["\[ExponentialE]",
2793     RowBox[{
2794     RowBox[{"-",
2795     SuperscriptBox["rc", "2"]}], " ",
2796     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2797     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2798     RowBox[{"(",
2799     RowBox[{"45", "+",
2800     RowBox[{"30", " ",
2801     SuperscriptBox["rc", "2"], " ",
2802     SuperscriptBox["\[Kappa]", "2"]}], "+",
2803     RowBox[{"12", " ",
2804     SuperscriptBox["rc", "4"], " ",
2805     SuperscriptBox["\[Kappa]", "4"]}], "+",
2806     RowBox[{"3", " ",
2807     SuperscriptBox["rc", "6"], " ",
2808     SuperscriptBox["\[Kappa]", "6"]}], "+",
2809     RowBox[{"2", " ",
2810     SuperscriptBox["rc", "8"], " ",
2811     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2812     RowBox[{"45", " ", "rc"}]], "+",
2813     FractionBox[
2814     RowBox[{"8", " ", "\[Pi]", " ",
2815     RowBox[{"Erfc", "[",
2816     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2817     SuperscriptBox["rc", "2"]]}], "0", "0"},
2818     {"0", "0", "0"}
2819     },
2820     GridBoxAlignment->{
2821     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2822     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2823     GridBoxSpacings->{"Columns" -> {
2824     Offset[0.27999999999999997`], {
2825     Offset[0.7]},
2826     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2827     "Rows" -> {
2828     Offset[0.2], {
2829     Offset[0.4]},
2830     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2831     RowBox[{"(", "\[NoBreak]", GridBox[{
2832     {"0", "0",
2833     RowBox[{
2834     FractionBox[
2835     RowBox[{"16", " ",
2836     SuperscriptBox["\[ExponentialE]",
2837     RowBox[{
2838     RowBox[{"-",
2839     SuperscriptBox["rc", "2"]}], " ",
2840     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2841     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2842     RowBox[{"(",
2843     RowBox[{"45", "+",
2844     RowBox[{"30", " ",
2845     SuperscriptBox["rc", "2"], " ",
2846     SuperscriptBox["\[Kappa]", "2"]}], "+",
2847     RowBox[{"12", " ",
2848     SuperscriptBox["rc", "4"], " ",
2849     SuperscriptBox["\[Kappa]", "4"]}], "+",
2850     RowBox[{"3", " ",
2851     SuperscriptBox["rc", "6"], " ",
2852     SuperscriptBox["\[Kappa]", "6"]}], "+",
2853     RowBox[{"2", " ",
2854     SuperscriptBox["rc", "8"], " ",
2855     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2856     RowBox[{"45", " ", "rc"}]], "+",
2857     FractionBox[
2858     RowBox[{"8", " ", "\[Pi]", " ",
2859     RowBox[{"Erfc", "[",
2860     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2861     SuperscriptBox["rc", "2"]]}]},
2862     {"0", "0", "0"},
2863     {
2864     RowBox[{
2865     FractionBox[
2866     RowBox[{"16", " ",
2867     SuperscriptBox["\[ExponentialE]",
2868     RowBox[{
2869     RowBox[{"-",
2870     SuperscriptBox["rc", "2"]}], " ",
2871     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2872     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2873     RowBox[{"(",
2874     RowBox[{"45", "+",
2875     RowBox[{"30", " ",
2876     SuperscriptBox["rc", "2"], " ",
2877     SuperscriptBox["\[Kappa]", "2"]}], "+",
2878     RowBox[{"12", " ",
2879     SuperscriptBox["rc", "4"], " ",
2880     SuperscriptBox["\[Kappa]", "4"]}], "+",
2881     RowBox[{"3", " ",
2882     SuperscriptBox["rc", "6"], " ",
2883     SuperscriptBox["\[Kappa]", "6"]}], "+",
2884     RowBox[{"2", " ",
2885     SuperscriptBox["rc", "8"], " ",
2886     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2887     RowBox[{"45", " ", "rc"}]], "+",
2888     FractionBox[
2889     RowBox[{"8", " ", "\[Pi]", " ",
2890     RowBox[{"Erfc", "[",
2891     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2892     SuperscriptBox["rc", "2"]]}], "0", "0"}
2893     },
2894     GridBoxAlignment->{
2895     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2896     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2897     GridBoxSpacings->{"Columns" -> {
2898     Offset[0.27999999999999997`], {
2899     Offset[0.7]},
2900     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2901     "Rows" -> {
2902     Offset[0.2], {
2903     Offset[0.4]},
2904     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
2905     {
2906     RowBox[{"(", "\[NoBreak]", GridBox[{
2907     {"0",
2908     RowBox[{
2909     FractionBox[
2910     RowBox[{"16", " ",
2911     SuperscriptBox["\[ExponentialE]",
2912     RowBox[{
2913     RowBox[{"-",
2914     SuperscriptBox["rc", "2"]}], " ",
2915     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2916     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2917     RowBox[{"(",
2918     RowBox[{"45", "+",
2919     RowBox[{"30", " ",
2920     SuperscriptBox["rc", "2"], " ",
2921     SuperscriptBox["\[Kappa]", "2"]}], "+",
2922     RowBox[{"12", " ",
2923     SuperscriptBox["rc", "4"], " ",
2924     SuperscriptBox["\[Kappa]", "4"]}], "+",
2925     RowBox[{"3", " ",
2926     SuperscriptBox["rc", "6"], " ",
2927     SuperscriptBox["\[Kappa]", "6"]}], "+",
2928     RowBox[{"2", " ",
2929     SuperscriptBox["rc", "8"], " ",
2930     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2931     RowBox[{"45", " ", "rc"}]], "+",
2932     FractionBox[
2933     RowBox[{"8", " ", "\[Pi]", " ",
2934     RowBox[{"Erfc", "[",
2935     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2936     SuperscriptBox["rc", "2"]]}], "0"},
2937     {
2938     RowBox[{
2939     FractionBox[
2940     RowBox[{"16", " ",
2941     SuperscriptBox["\[ExponentialE]",
2942     RowBox[{
2943     RowBox[{"-",
2944     SuperscriptBox["rc", "2"]}], " ",
2945     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2946     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2947     RowBox[{"(",
2948     RowBox[{"45", "+",
2949     RowBox[{"30", " ",
2950     SuperscriptBox["rc", "2"], " ",
2951     SuperscriptBox["\[Kappa]", "2"]}], "+",
2952     RowBox[{"12", " ",
2953     SuperscriptBox["rc", "4"], " ",
2954     SuperscriptBox["\[Kappa]", "4"]}], "+",
2955     RowBox[{"3", " ",
2956     SuperscriptBox["rc", "6"], " ",
2957     SuperscriptBox["\[Kappa]", "6"]}], "+",
2958     RowBox[{"2", " ",
2959     SuperscriptBox["rc", "8"], " ",
2960     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
2961     RowBox[{"45", " ", "rc"}]], "+",
2962     FractionBox[
2963     RowBox[{"8", " ", "\[Pi]", " ",
2964     RowBox[{"Erfc", "[",
2965     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
2966     SuperscriptBox["rc", "2"]]}], "0", "0"},
2967     {"0", "0", "0"}
2968     },
2969     GridBoxAlignment->{
2970     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
2971     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
2972     GridBoxSpacings->{"Columns" -> {
2973     Offset[0.27999999999999997`], {
2974     Offset[0.7]},
2975     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
2976     "Rows" -> {
2977     Offset[0.2], {
2978     Offset[0.4]},
2979     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
2980     RowBox[{"(", "\[NoBreak]", GridBox[{
2981     {
2982     RowBox[{
2983     FractionBox[
2984     RowBox[{"16", " ",
2985     SuperscriptBox["\[ExponentialE]",
2986     RowBox[{
2987     RowBox[{"-",
2988     SuperscriptBox["rc", "2"]}], " ",
2989     SuperscriptBox["\[Kappa]", "2"]}]], " ",
2990     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
2991     RowBox[{"(",
2992     RowBox[{"45", "+",
2993     RowBox[{"30", " ",
2994     SuperscriptBox["rc", "2"], " ",
2995     SuperscriptBox["\[Kappa]", "2"]}], "+",
2996     RowBox[{"12", " ",
2997     SuperscriptBox["rc", "4"], " ",
2998     SuperscriptBox["\[Kappa]", "4"]}], "+",
2999     RowBox[{"3", " ",
3000     SuperscriptBox["rc", "6"], " ",
3001     SuperscriptBox["\[Kappa]", "6"]}], "+",
3002     RowBox[{"2", " ",
3003     SuperscriptBox["rc", "8"], " ",
3004     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3005     RowBox[{"45", " ", "rc"}]], "+",
3006     FractionBox[
3007     RowBox[{"8", " ", "\[Pi]", " ",
3008     RowBox[{"Erfc", "[",
3009     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3010     SuperscriptBox["rc", "2"]]}], "0", "0"},
3011     {"0",
3012     RowBox[{
3013     FractionBox[
3014     RowBox[{"16", " ",
3015     SuperscriptBox["\[ExponentialE]",
3016     RowBox[{
3017     RowBox[{"-",
3018     SuperscriptBox["rc", "2"]}], " ",
3019     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3020     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3021     RowBox[{"(",
3022     RowBox[{"45", "+",
3023     RowBox[{"30", " ",
3024     SuperscriptBox["rc", "2"], " ",
3025     SuperscriptBox["\[Kappa]", "2"]}], "+",
3026     RowBox[{"12", " ",
3027     SuperscriptBox["rc", "4"], " ",
3028     SuperscriptBox["\[Kappa]", "4"]}], "+",
3029     RowBox[{"3", " ",
3030     SuperscriptBox["rc", "6"], " ",
3031     SuperscriptBox["\[Kappa]", "6"]}], "+",
3032     RowBox[{"2", " ",
3033     SuperscriptBox["rc", "8"], " ",
3034     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3035     RowBox[{"15", " ", "rc"}]], "+",
3036     FractionBox[
3037     RowBox[{"24", " ", "\[Pi]", " ",
3038     RowBox[{"Erfc", "[",
3039     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3040     SuperscriptBox["rc", "2"]]}], "0"},
3041     {"0", "0",
3042     RowBox[{
3043     FractionBox[
3044     RowBox[{"16", " ",
3045     SuperscriptBox["\[ExponentialE]",
3046     RowBox[{
3047     RowBox[{"-",
3048     SuperscriptBox["rc", "2"]}], " ",
3049     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3050     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3051     RowBox[{"(",
3052     RowBox[{"45", "+",
3053     RowBox[{"30", " ",
3054     SuperscriptBox["rc", "2"], " ",
3055     SuperscriptBox["\[Kappa]", "2"]}], "+",
3056     RowBox[{"12", " ",
3057     SuperscriptBox["rc", "4"], " ",
3058     SuperscriptBox["\[Kappa]", "4"]}], "+",
3059     RowBox[{"3", " ",
3060     SuperscriptBox["rc", "6"], " ",
3061     SuperscriptBox["\[Kappa]", "6"]}], "+",
3062     RowBox[{"2", " ",
3063     SuperscriptBox["rc", "8"], " ",
3064     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3065     RowBox[{"45", " ", "rc"}]], "+",
3066     FractionBox[
3067     RowBox[{"8", " ", "\[Pi]", " ",
3068     RowBox[{"Erfc", "[",
3069     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3070     SuperscriptBox["rc", "2"]]}]}
3071     },
3072     GridBoxAlignment->{
3073     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
3074     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
3075     GridBoxSpacings->{"Columns" -> {
3076     Offset[0.27999999999999997`], {
3077     Offset[0.7]},
3078     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
3079     "Rows" -> {
3080     Offset[0.2], {
3081     Offset[0.4]},
3082     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
3083     RowBox[{"(", "\[NoBreak]", GridBox[{
3084     {"0", "0", "0"},
3085     {"0", "0",
3086     RowBox[{
3087     FractionBox[
3088     RowBox[{"16", " ",
3089     SuperscriptBox["\[ExponentialE]",
3090     RowBox[{
3091     RowBox[{"-",
3092     SuperscriptBox["rc", "2"]}], " ",
3093     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3094     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3095     RowBox[{"(",
3096     RowBox[{"45", "+",
3097     RowBox[{"30", " ",
3098     SuperscriptBox["rc", "2"], " ",
3099     SuperscriptBox["\[Kappa]", "2"]}], "+",
3100     RowBox[{"12", " ",
3101     SuperscriptBox["rc", "4"], " ",
3102     SuperscriptBox["\[Kappa]", "4"]}], "+",
3103     RowBox[{"3", " ",
3104     SuperscriptBox["rc", "6"], " ",
3105     SuperscriptBox["\[Kappa]", "6"]}], "+",
3106     RowBox[{"2", " ",
3107     SuperscriptBox["rc", "8"], " ",
3108     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3109     RowBox[{"45", " ", "rc"}]], "+",
3110     FractionBox[
3111     RowBox[{"8", " ", "\[Pi]", " ",
3112     RowBox[{"Erfc", "[",
3113     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3114     SuperscriptBox["rc", "2"]]}]},
3115     {"0",
3116     RowBox[{
3117     FractionBox[
3118     RowBox[{"16", " ",
3119     SuperscriptBox["\[ExponentialE]",
3120     RowBox[{
3121     RowBox[{"-",
3122     SuperscriptBox["rc", "2"]}], " ",
3123     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3124     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3125     RowBox[{"(",
3126     RowBox[{"45", "+",
3127     RowBox[{"30", " ",
3128     SuperscriptBox["rc", "2"], " ",
3129     SuperscriptBox["\[Kappa]", "2"]}], "+",
3130     RowBox[{"12", " ",
3131     SuperscriptBox["rc", "4"], " ",
3132     SuperscriptBox["\[Kappa]", "4"]}], "+",
3133     RowBox[{"3", " ",
3134     SuperscriptBox["rc", "6"], " ",
3135     SuperscriptBox["\[Kappa]", "6"]}], "+",
3136     RowBox[{"2", " ",
3137     SuperscriptBox["rc", "8"], " ",
3138     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3139     RowBox[{"45", " ", "rc"}]], "+",
3140     FractionBox[
3141     RowBox[{"8", " ", "\[Pi]", " ",
3142     RowBox[{"Erfc", "[",
3143     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3144     SuperscriptBox["rc", "2"]]}], "0"}
3145     },
3146     GridBoxAlignment->{
3147     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
3148     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
3149     GridBoxSpacings->{"Columns" -> {
3150     Offset[0.27999999999999997`], {
3151     Offset[0.7]},
3152     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
3153     "Rows" -> {
3154     Offset[0.2], {
3155     Offset[0.4]},
3156     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]},
3157     {
3158     RowBox[{"(", "\[NoBreak]", GridBox[{
3159     {"0", "0",
3160     RowBox[{
3161     FractionBox[
3162     RowBox[{"16", " ",
3163     SuperscriptBox["\[ExponentialE]",
3164     RowBox[{
3165     RowBox[{"-",
3166     SuperscriptBox["rc", "2"]}], " ",
3167     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3168     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3169     RowBox[{"(",
3170     RowBox[{"45", "+",
3171     RowBox[{"30", " ",
3172     SuperscriptBox["rc", "2"], " ",
3173     SuperscriptBox["\[Kappa]", "2"]}], "+",
3174     RowBox[{"12", " ",
3175     SuperscriptBox["rc", "4"], " ",
3176     SuperscriptBox["\[Kappa]", "4"]}], "+",
3177     RowBox[{"3", " ",
3178     SuperscriptBox["rc", "6"], " ",
3179     SuperscriptBox["\[Kappa]", "6"]}], "+",
3180     RowBox[{"2", " ",
3181     SuperscriptBox["rc", "8"], " ",
3182     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3183     RowBox[{"45", " ", "rc"}]], "+",
3184     FractionBox[
3185     RowBox[{"8", " ", "\[Pi]", " ",
3186     RowBox[{"Erfc", "[",
3187     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3188     SuperscriptBox["rc", "2"]]}]},
3189     {"0", "0", "0"},
3190     {
3191     RowBox[{
3192     FractionBox[
3193     RowBox[{"16", " ",
3194     SuperscriptBox["\[ExponentialE]",
3195     RowBox[{
3196     RowBox[{"-",
3197     SuperscriptBox["rc", "2"]}], " ",
3198     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3199     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3200     RowBox[{"(",
3201     RowBox[{"45", "+",
3202     RowBox[{"30", " ",
3203     SuperscriptBox["rc", "2"], " ",
3204     SuperscriptBox["\[Kappa]", "2"]}], "+",
3205     RowBox[{"12", " ",
3206     SuperscriptBox["rc", "4"], " ",
3207     SuperscriptBox["\[Kappa]", "4"]}], "+",
3208     RowBox[{"3", " ",
3209     SuperscriptBox["rc", "6"], " ",
3210     SuperscriptBox["\[Kappa]", "6"]}], "+",
3211     RowBox[{"2", " ",
3212     SuperscriptBox["rc", "8"], " ",
3213     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3214     RowBox[{"45", " ", "rc"}]], "+",
3215     FractionBox[
3216     RowBox[{"8", " ", "\[Pi]", " ",
3217     RowBox[{"Erfc", "[",
3218     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3219     SuperscriptBox["rc", "2"]]}], "0", "0"}
3220     },
3221     GridBoxAlignment->{
3222     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
3223     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
3224     GridBoxSpacings->{"Columns" -> {
3225     Offset[0.27999999999999997`], {
3226     Offset[0.7]},
3227     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
3228     "Rows" -> {
3229     Offset[0.2], {
3230     Offset[0.4]},
3231     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
3232     RowBox[{"(", "\[NoBreak]", GridBox[{
3233     {"0", "0", "0"},
3234     {"0", "0",
3235     RowBox[{
3236     FractionBox[
3237     RowBox[{"16", " ",
3238     SuperscriptBox["\[ExponentialE]",
3239     RowBox[{
3240     RowBox[{"-",
3241     SuperscriptBox["rc", "2"]}], " ",
3242     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3243     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3244     RowBox[{"(",
3245     RowBox[{"45", "+",
3246     RowBox[{"30", " ",
3247     SuperscriptBox["rc", "2"], " ",
3248     SuperscriptBox["\[Kappa]", "2"]}], "+",
3249     RowBox[{"12", " ",
3250     SuperscriptBox["rc", "4"], " ",
3251     SuperscriptBox["\[Kappa]", "4"]}], "+",
3252     RowBox[{"3", " ",
3253     SuperscriptBox["rc", "6"], " ",
3254     SuperscriptBox["\[Kappa]", "6"]}], "+",
3255     RowBox[{"2", " ",
3256     SuperscriptBox["rc", "8"], " ",
3257     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3258     RowBox[{"45", " ", "rc"}]], "+",
3259     FractionBox[
3260     RowBox[{"8", " ", "\[Pi]", " ",
3261     RowBox[{"Erfc", "[",
3262     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3263     SuperscriptBox["rc", "2"]]}]},
3264     {"0",
3265     RowBox[{
3266     FractionBox[
3267     RowBox[{"16", " ",
3268     SuperscriptBox["\[ExponentialE]",
3269     RowBox[{
3270     RowBox[{"-",
3271     SuperscriptBox["rc", "2"]}], " ",
3272     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3273     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3274     RowBox[{"(",
3275     RowBox[{"45", "+",
3276     RowBox[{"30", " ",
3277     SuperscriptBox["rc", "2"], " ",
3278     SuperscriptBox["\[Kappa]", "2"]}], "+",
3279     RowBox[{"12", " ",
3280     SuperscriptBox["rc", "4"], " ",
3281     SuperscriptBox["\[Kappa]", "4"]}], "+",
3282     RowBox[{"3", " ",
3283     SuperscriptBox["rc", "6"], " ",
3284     SuperscriptBox["\[Kappa]", "6"]}], "+",
3285     RowBox[{"2", " ",
3286     SuperscriptBox["rc", "8"], " ",
3287     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3288     RowBox[{"45", " ", "rc"}]], "+",
3289     FractionBox[
3290     RowBox[{"8", " ", "\[Pi]", " ",
3291     RowBox[{"Erfc", "[",
3292     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3293     SuperscriptBox["rc", "2"]]}], "0"}
3294     },
3295     GridBoxAlignment->{
3296     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
3297     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
3298     GridBoxSpacings->{"Columns" -> {
3299     Offset[0.27999999999999997`], {
3300     Offset[0.7]},
3301     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
3302     "Rows" -> {
3303     Offset[0.2], {
3304     Offset[0.4]},
3305     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
3306     RowBox[{"(", "\[NoBreak]", GridBox[{
3307     {
3308     RowBox[{
3309     FractionBox[
3310     RowBox[{"16", " ",
3311     SuperscriptBox["\[ExponentialE]",
3312     RowBox[{
3313     RowBox[{"-",
3314     SuperscriptBox["rc", "2"]}], " ",
3315     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3316     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3317     RowBox[{"(",
3318     RowBox[{"45", "+",
3319     RowBox[{"30", " ",
3320     SuperscriptBox["rc", "2"], " ",
3321     SuperscriptBox["\[Kappa]", "2"]}], "+",
3322     RowBox[{"12", " ",
3323     SuperscriptBox["rc", "4"], " ",
3324     SuperscriptBox["\[Kappa]", "4"]}], "+",
3325     RowBox[{"3", " ",
3326     SuperscriptBox["rc", "6"], " ",
3327     SuperscriptBox["\[Kappa]", "6"]}], "+",
3328     RowBox[{"2", " ",
3329     SuperscriptBox["rc", "8"], " ",
3330     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3331     RowBox[{"45", " ", "rc"}]], "+",
3332     FractionBox[
3333     RowBox[{"8", " ", "\[Pi]", " ",
3334     RowBox[{"Erfc", "[",
3335     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3336     SuperscriptBox["rc", "2"]]}], "0", "0"},
3337     {"0",
3338     RowBox[{
3339     FractionBox[
3340     RowBox[{"16", " ",
3341     SuperscriptBox["\[ExponentialE]",
3342     RowBox[{
3343     RowBox[{"-",
3344     SuperscriptBox["rc", "2"]}], " ",
3345     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3346     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3347     RowBox[{"(",
3348     RowBox[{"45", "+",
3349     RowBox[{"30", " ",
3350     SuperscriptBox["rc", "2"], " ",
3351     SuperscriptBox["\[Kappa]", "2"]}], "+",
3352     RowBox[{"12", " ",
3353     SuperscriptBox["rc", "4"], " ",
3354     SuperscriptBox["\[Kappa]", "4"]}], "+",
3355     RowBox[{"3", " ",
3356     SuperscriptBox["rc", "6"], " ",
3357     SuperscriptBox["\[Kappa]", "6"]}], "+",
3358     RowBox[{"2", " ",
3359     SuperscriptBox["rc", "8"], " ",
3360     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3361     RowBox[{"45", " ", "rc"}]], "+",
3362     FractionBox[
3363     RowBox[{"8", " ", "\[Pi]", " ",
3364     RowBox[{"Erfc", "[",
3365     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3366     SuperscriptBox["rc", "2"]]}], "0"},
3367     {"0", "0",
3368     RowBox[{
3369     FractionBox[
3370     RowBox[{"16", " ",
3371     SuperscriptBox["\[ExponentialE]",
3372     RowBox[{
3373     RowBox[{"-",
3374     SuperscriptBox["rc", "2"]}], " ",
3375     SuperscriptBox["\[Kappa]", "2"]}]], " ",
3376     SqrtBox["\[Pi]"], " ", "\[Kappa]", " ",
3377     RowBox[{"(",
3378     RowBox[{"45", "+",
3379     RowBox[{"30", " ",
3380     SuperscriptBox["rc", "2"], " ",
3381     SuperscriptBox["\[Kappa]", "2"]}], "+",
3382     RowBox[{"12", " ",
3383     SuperscriptBox["rc", "4"], " ",
3384     SuperscriptBox["\[Kappa]", "4"]}], "+",
3385     RowBox[{"3", " ",
3386     SuperscriptBox["rc", "6"], " ",
3387     SuperscriptBox["\[Kappa]", "6"]}], "+",
3388     RowBox[{"2", " ",
3389     SuperscriptBox["rc", "8"], " ",
3390     SuperscriptBox["\[Kappa]", "8"]}]}], ")"}]}],
3391     RowBox[{"15", " ", "rc"}]], "+",
3392     FractionBox[
3393     RowBox[{"24", " ", "\[Pi]", " ",
3394     RowBox[{"Erfc", "[",
3395     RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}],
3396     SuperscriptBox["rc", "2"]]}]}
3397     },
3398     GridBoxAlignment->{
3399     "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
3400     "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
3401     GridBoxSpacings->{"Columns" -> {
3402     Offset[0.27999999999999997`], {
3403     Offset[0.7]},
3404     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {},
3405     "Rows" -> {
3406     Offset[0.2], {
3407     Offset[0.4]},
3408     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}]}
3409     },
3410     GridBoxAlignment->{
3411     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
3412     "RowsIndexed" -> {}},
3413     GridBoxSpacings->{"Columns" -> {
3414     Offset[0.27999999999999997`], {
3415     Offset[0.7]},
3416     Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
3417     Offset[0.2], {
3418     Offset[0.4]},
3419     Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
3420     Function[BoxForm`e$,
3421     MatrixForm[BoxForm`e$]]]], "Output",
3422     CellChangeTimes->{3.6231506637160788`*^9}]
3423     }, Open ]]
3424     },
3425     WindowSize->{884, 1337},
3426     WindowMargins->{{Automatic, 772}, {Automatic, 0}},
3427 gezelter 4433 FrontEndVersion->"10.4 for Mac OS X x86 (32-bit, 64-bit Kernel) (April 11, \
3428     2016)",
3429 gezelter 4432 StyleDefinitions->"Default.nb"
3430     ]
3431     (* End of Notebook Content *)
3432    
3433     (* Internal cache information *)
3434     (*CellTagsOutline
3435     CellTagsIndex->{}
3436     *)
3437     (*CellTagsIndex
3438     CellTagsIndex->{}
3439     *)
3440     (*NotebookFileOutline
3441     Notebook[{
3442 gezelter 4433 Cell[400, 13, 37088, 1062, 1525, "Input"],
3443 gezelter 4432 Cell[CellGroupData[{
3444 gezelter 4433 Cell[37513, 1079, 159, 3, 28, "Input"],
3445     Cell[37675, 1084, 617, 16, 58, "Output"]
3446 gezelter 4432 }, Open ]],
3447     Cell[CellGroupData[{
3448 gezelter 4433 Cell[38329, 1105, 153, 3, 28, "Input"],
3449     Cell[38485, 1110, 529, 15, 58, "Output"]
3450 gezelter 4432 }, Open ]],
3451     Cell[CellGroupData[{
3452 gezelter 4433 Cell[39051, 1130, 230, 4, 28, "Input"],
3453     Cell[39284, 1136, 1027, 29, 58, "Output"]
3454 gezelter 4432 }, Open ]],
3455     Cell[CellGroupData[{
3456 gezelter 4433 Cell[40348, 1170, 126, 2, 28, "Input"],
3457     Cell[40477, 1174, 1983, 62, 92, "Output"]
3458 gezelter 4432 }, Open ]],
3459     Cell[CellGroupData[{
3460 gezelter 4433 Cell[42497, 1241, 218, 5, 28, "Input"],
3461     Cell[42718, 1248, 272, 4, 28, "Output"]
3462 gezelter 4432 }, Open ]],
3463     Cell[CellGroupData[{
3464 gezelter 4433 Cell[43027, 1257, 190, 4, 28, "Input"],
3465     Cell[43220, 1263, 119, 2, 28, "Output"]
3466 gezelter 4432 }, Open ]],
3467     Cell[CellGroupData[{
3468 gezelter 4433 Cell[43376, 1270, 211, 5, 28, "Input"],
3469     Cell[43590, 1277, 92, 1, 28, "Output"]
3470 gezelter 4432 }, Open ]],
3471     Cell[CellGroupData[{
3472 gezelter 4433 Cell[43719, 1283, 187, 4, 28, "Input"],
3473     Cell[43909, 1289, 119, 2, 28, "Output"]
3474 gezelter 4432 }, Open ]],
3475     Cell[CellGroupData[{
3476 gezelter 4433 Cell[44065, 1296, 187, 4, 28, "Input"],
3477     Cell[44255, 1302, 163, 4, 49, "Output"]
3478 gezelter 4432 }, Open ]],
3479     Cell[CellGroupData[{
3480 gezelter 4433 Cell[44455, 1311, 148, 2, 28, "Input"],
3481     Cell[44606, 1315, 14732, 396, 267, "Output"]
3482 gezelter 4432 }, Open ]],
3483     Cell[CellGroupData[{
3484 gezelter 4433 Cell[59375, 1716, 95, 1, 28, "Input"],
3485     Cell[59473, 1719, 14611, 394, 267, "Output"]
3486 gezelter 4432 }, Open ]],
3487     Cell[CellGroupData[{
3488 gezelter 4433 Cell[74121, 2118, 98, 1, 28, "Input"],
3489     Cell[74222, 2121, 19397, 520, 267, "Output"]
3490 gezelter 4432 }, Open ]],
3491     Cell[CellGroupData[{
3492 gezelter 4433 Cell[93656, 2646, 98, 1, 28, "Input"],
3493     Cell[93757, 2649, 31477, 772, 436, "Output"]
3494 gezelter 4432 }, Open ]]
3495     }
3496     ]
3497     *)
3498