(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 103708, 2447] NotebookOptionsPosition[ 101194, 2360] NotebookOutlinePosition[ 101552, 2376] CellTagsIndexPosition[ 101509, 2373] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["FOR PROCESSING", "Subtitle", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ StyleBox[ RowBox[{ RowBox[{"nParticle", "=", StyleBox[ RowBox[{"0.1748", "+", RowBox[{"3.4618", " ", "I"}]}], FontSize->12]}], ";"}], FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"particles", " ", "complex", " ", "refractive", " ", "index"}], ",", " ", RowBox[{ RowBox[{"-", "3.4618"}], " ", "I"}]}], "*)"}], FontColor->RGBColor[1, 0, 0]]}], "\[IndentingNewLine]", RowBox[{ StyleBox[ RowBox[{"R", "=", StyleBox[ RowBox[{"3", "0"}]]}], FontColor->RGBColor[1, 0, 0]], StyleBox[";", FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{"(*", RowBox[{ "Particles", " ", "radius", " ", "R", " ", "should", " ", "have", " ", "same", " ", "units", " ", "as", " ", RowBox[{"\[Lambda]", "!!"}]}], "*)"}], FontColor->RGBColor[1, 0, 0]], "\[IndentingNewLine]", StyleBox[ RowBox[{"nmed", "=", "1.46"}], FontColor->RGBColor[1, 0, 0]], StyleBox[";", FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"mediums", "'"}], " ", "real", " ", "refractive", " ", "index"}], "*)"}], FontColor->RGBColor[1, 0, 0]], "\[IndentingNewLine]", StyleBox[ RowBox[{"\[Lambda]", "=", "635"}], FontColor->RGBColor[1, 0, 0]], StyleBox[";", FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{ "probing", " ", "beam", " ", "wavelength", " ", "in", " ", "vacuum"}], ",", " ", RowBox[{"used", " ", "for", " ", "scatteruing", " ", "calculations"}]}], "*)"}], FontColor->RGBColor[1, 0, 0]], "\[IndentingNewLine]", RowBox[{"Mn", "=", RowBox[{"nParticle", "/", "nmed"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k", "=", RowBox[{"nmed", "*", "2", RowBox[{"\[Pi]", "/", "\[Lambda]"}]}]}], ";", RowBox[{"(*", RowBox[{"wave", "-", RowBox[{"vector", " ", "in", " ", "the", " ", "medium"}]}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"nmax", "=", "50"}], ";", RowBox[{"(*", RowBox[{ RowBox[{ "maximum", " ", "multipole", " ", "order", " ", "for", " ", "non"}], "-", RowBox[{"multilayer", " ", "Mie"}], "-", RowBox[{ RowBox[{"coefficients", ".", " ", "Multilayer"}], " ", "coefs", " ", "use", " ", "computed", " ", "nmax"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"npmax", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[Alpha]0", "=", RowBox[{"4", "\[Pi]", " ", SuperscriptBox["R", "3"], RowBox[{"(", FractionBox[ RowBox[{ SuperscriptBox["Mn", "2"], "-", "1"}], RowBox[{ SuperscriptBox["Mn", "2"], "+", "2"}]], ")"}]}]}], ",", RowBox[{ SuperscriptBox["k", "3"], "\[Alpha]0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"\[Alpha]a1", "=", FractionBox[ RowBox[{"6", "\[Pi]", " ", RowBox[{"aMie", "[", RowBox[{"1", ",", RowBox[{"k", " ", "R"}], ",", "Mn"}], "]"}]}], RowBox[{"I", " ", SuperscriptBox["k", "3"]}]]}], ",", RowBox[{ SuperscriptBox["k", "3"], "\[Alpha]a1"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"\[Alpha]RR2", "=", FractionBox["\[Alpha]0", RowBox[{"1", "+", RowBox[{ FractionBox["1", RowBox[{"6", "\[Pi]"}]], "I", " ", SuperscriptBox["k", "3"], "\[Alpha]0"}]}]]}], ",", RowBox[{ SuperscriptBox["k", "3"], "\[Alpha]RR2"}]}], "}"}]}], "}"}]}], "Input", CellChangeTimes->{{3.601903700113559*^9, 3.601903775895965*^9}, { 3.601903862804308*^9, 3.6019038872914963`*^9}, {3.6019118211065817`*^9, 3.601911848756424*^9}, {3.601911918007708*^9, 3.601912022700201*^9}, { 3.6019120943288393`*^9, 3.60191213374644*^9}, {3.601912314268362*^9, 3.601912314502075*^9}, {3.6019126997375383`*^9, 3.6019128860461206`*^9}, { 3.601916462647184*^9, 3.6019165222387466`*^9}, {3.6019173049944677`*^9, 3.601917400282938*^9}, {3.60191865083187*^9, 3.601918654894009*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"614608.52787048`", "\[VeryThinSpace]", "+", RowBox[{"43327.352432756416`", " ", "\[ImaginaryI]"}]}], ",", RowBox[{"1.8530006790699598`", "\[VeryThinSpace]", "+", RowBox[{"0.13062886347896585`", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"791747.788183802`", "\[VeryThinSpace]", "-", RowBox[{"25701.255069766918`", " ", "\[ImaginaryI]"}]}], ",", RowBox[{"2.3870628581090827`", "\[VeryThinSpace]", "-", RowBox[{"0.07748744271778885`", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"617168.7321178211`", "\[VeryThinSpace]", "-", RowBox[{"17464.287674589214`", " ", "\[ImaginaryI]"}]}], ",", RowBox[{"1.8607195114547272`", "\[VeryThinSpace]", "-", RowBox[{"0.05265357614319766`", " ", "\[ImaginaryI]"}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.6019173879364243`*^9, 3.6019174007944202`*^9}, { 3.601918651342744*^9, 3.601918655360819*^9}, 3.601919701797957*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Bessel", " ", "related", " ", "functions", " ", "needed", " ", "for", " ", "Mie", " ", RowBox[{"coeff", "."}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"psi", "[", RowBox[{"n_", ",", "x_"}], "]"}], ":=", RowBox[{"x", " ", RowBox[{"SphericalBesselJ", "[", RowBox[{"n", ",", "x"}], "]"}]}]}], RowBox[{"(*", RowBox[{"Riccarti", " ", "\.06Bessel", " ", "\.06Functions"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xi", "[", RowBox[{"n_", ",", "x_"}], "]"}], ":=", RowBox[{"x", " ", RowBox[{"SphericalHankelH2", "[", RowBox[{"n", ",", "x"}], "]"}]}]}], "\n", RowBox[{ RowBox[{ RowBox[{"psidev", "[", RowBox[{"n_", ",", "x_"}], "]"}], ":=", " ", RowBox[{ RowBox[{ RowBox[{"Derivative", "[", RowBox[{"0", ",", "1"}], "]"}], "[", "psi", "]"}], "[", RowBox[{"n", ",", "x"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"xidev", "[", RowBox[{"n_", ",", "x_"}], "]"}], ":=", " ", RowBox[{ RowBox[{ RowBox[{"Derivative", "[", RowBox[{"0", ",", "1"}], "]"}], "[", "xi", "]"}], "[", RowBox[{"n", ",", "x"}], "]"}]}], ";"}], "\n", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Mie", " ", RowBox[{"Coeff", ".", " ", RowBox[{"{", RowBox[{ SubscriptBox["a", "n"], ",", SubscriptBox["b", "n"], ",", SubscriptBox["c", "n"], ",", SubscriptBox["d", "n"]}], "}"}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"aMie", "[", RowBox[{"n_", ",", "x_", ",", "m_"}], "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"m", " ", RowBox[{"psi", "[", RowBox[{"n", ",", RowBox[{"m", " ", "x"}]}], "]"}], RowBox[{"psidev", "[", RowBox[{"n", ",", "x"}], "]"}]}], "-", RowBox[{ RowBox[{"psi", "[", RowBox[{"n", ",", "x"}], "]"}], RowBox[{"psidev", "[", RowBox[{"n", ",", RowBox[{"m", " ", "x"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"m", " ", RowBox[{"psi", "[", RowBox[{"n", ",", RowBox[{"m", " ", "x"}]}], "]"}], RowBox[{"xidev", "[", RowBox[{"n", ",", "x"}], "]"}]}], "-", RowBox[{ RowBox[{"xi", "[", RowBox[{"n", ",", "x"}], "]"}], RowBox[{"psidev", "[", RowBox[{"n", ",", RowBox[{"m", " ", "x"}]}], "]"}]}]}], ")"}]}]}], ";"}], "\n"}]}]], "Input", InitializationCell->True, CellChangeTimes->{{3.497679293332241*^9, 3.497679311135788*^9}, 3.4976910585174065`*^9, {3.497691089835766*^9, 3.497691127150219*^9}, { 3.4976911798367424`*^9, 3.4976911799968987`*^9}, {3.4976914407937737`*^9, 3.4976914413533435`*^9}, {3.4976914896785393`*^9, 3.4976914900261955`*^9}, {3.49772132537926*^9, 3.497721339464965*^9}, { 3.497721570260353*^9, 3.497721598022975*^9}, {3.497721749402769*^9, 3.497721762310053*^9}, {3.4977217972727537`*^9, 3.4977218026468554`*^9}, 3.4977218349709177`*^9, {3.497721891211123*^9, 3.497721903768449*^9}, 3.49772196643891*^9, {3.4977221601291847`*^9, 3.497722185887884*^9}, { 3.497722242014049*^9, 3.497722243310285*^9}, {3.497722275862769*^9, 3.497722275970105*^9}, 3.4977224348748493`*^9, {3.4977229731364822`*^9, 3.497723000014002*^9}, {3.4977238550572367`*^9, 3.4977238554144363`*^9}, { 3.497723910366666*^9, 3.497723977430134*^9}, {3.497726965139494*^9, 3.497726986421715*^9}, {3.497776478969555*^9, 3.4977764811267815`*^9}, { 3.497777430895336*^9, 3.497777433495922*^9}, 3.4977775617029533`*^9, { 3.49779038878318*^9, 3.497790415014697*^9}, {3.497794082126857*^9, 3.4977940969297132`*^9}, {3.497801845201638*^9, 3.497801940001289*^9}, 3.497801975350624*^9, {3.497802014252404*^9, 3.497802020197878*^9}, { 3.497802052282159*^9, 3.497802052396743*^9}, 3.497802310840019*^9, { 3.497802653906659*^9, 3.497802682161769*^9}, {3.497808600540502*^9, 3.497808607650961*^9}, 3.497808668699174*^9, {3.497808706986237*^9, 3.497808721164303*^9}, {3.4978088072476482`*^9, 3.4978088141990757`*^9}, { 3.497809053490217*^9, 3.4978090605601673`*^9}, {3.4978103014510717`*^9, 3.4978103770193453`*^9}, {3.497812709771573*^9, 3.497812714534457*^9}, { 3.49786188195744*^9, 3.497861885312722*^9}, 3.4978671955675383`*^9, { 3.498283607300857*^9, 3.498283618488854*^9}, {3.498286561151466*^9, 3.498286602612596*^9}, {3.498305987823595*^9, 3.498305999290082*^9}, 3.4984587480803432`*^9, {3.4984654932136803`*^9, 3.4984654937634706`*^9}, { 3.4984655289584017`*^9, 3.4984655293138227`*^9}, 3.498471095479938*^9, { 3.4984772889311857`*^9, 3.498477314656807*^9}, 3.498907445117175*^9, { 3.498908229192465*^9, 3.498908229379212*^9}, {3.498909133123087*^9, 3.498909137412986*^9}, {3.498910465415127*^9, 3.498910472713295*^9}, { 3.498913997988604*^9, 3.4989139981235943`*^9}, {3.498914044248953*^9, 3.4989140444040003`*^9}, {3.4989242481077633`*^9, 3.498924258711973*^9}, { 3.498924495285656*^9, 3.498924503175386*^9}, {3.4989248362605257`*^9, 3.498924844196879*^9}, {3.498926513767701*^9, 3.498926516569848*^9}, { 3.498926648318211*^9, 3.498926654989155*^9}, {3.4989896674195127`*^9, 3.4989896684066477`*^9}, {3.498990506010316*^9, 3.498990513117537*^9}, { 3.498990961287756*^9, 3.49899096745578*^9}, {3.4990624977519073`*^9, 3.49906251192861*^9}, {3.499063058123124*^9, 3.4990630592149067`*^9}, 3.499063347322259*^9, {3.4990636296587954`*^9, 3.49906363253699*^9}, { 3.4990637398070602`*^9, 3.499063745369041*^9}, 3.4990642800514517`*^9, 3.4990801117228603`*^9, {3.499080305501318*^9, 3.499080308325286*^9}, { 3.499160548274925*^9, 3.4991605521064043`*^9}, {3.499161023996531*^9, 3.499161026555059*^9}, {3.499161440224216*^9, 3.499161444143832*^9}, { 3.499163325826466*^9, 3.4991633288822613`*^9}, {3.499172034004579*^9, 3.499172037384995*^9}, {3.499172834333021*^9, 3.499172838850675*^9}, { 3.499173043551001*^9, 3.499173046009667*^9}, {3.499173157549073*^9, 3.499173160513274*^9}, {3.49942808085717*^9, 3.499428087390359*^9}, { 3.500206447540786*^9, 3.500206456023302*^9}, {3.508589187074369*^9, 3.508589190231892*^9}, {3.508589476850644*^9, 3.508589477348749*^9}, { 3.509170498779584*^9, 3.509170505949778*^9}, {3.509170589295842*^9, 3.509170598956603*^9}, {3.5091706523625507`*^9, 3.509170658059001*^9}, { 3.515318592297471*^9, 3.515318595836187*^9}, {3.515387946050088*^9, 3.515387950131812*^9}, {3.5165588160855722`*^9, 3.5165588188258343`*^9}, { 3.51659636067799*^9, 3.516596361320985*^9}, {3.516596438546509*^9, 3.516596439519425*^9}, {3.516597085231395*^9, 3.516597086470118*^9}, { 3.51659787660032*^9, 3.5165978936755457`*^9}, {3.517039463911162*^9, 3.517039474548794*^9}, {3.5237776182622347`*^9, 3.523777684667384*^9}, { 3.528442316587049*^9, 3.528442318165063*^9}, {3.5284424733858147`*^9, 3.528442475209852*^9}, {3.528442685207932*^9, 3.5284426857944717`*^9}, { 3.528442953987035*^9, 3.5284429846481256`*^9}, {3.5284507666223927`*^9, 3.5284507667840652`*^9}, {3.52845097245695*^9, 3.5284509856091967`*^9}, { 3.528451312452523*^9, 3.528451312669868*^9}, {3.529644872179673*^9, 3.5296449006552877`*^9}, {3.529645362971634*^9, 3.52964539217627*^9}, { 3.5296464753971043`*^9, 3.529646475591346*^9}, {3.529647169063978*^9, 3.529647169396578*^9}, {3.529649905201055*^9, 3.529649906343062*^9}, { 3.529655641621482*^9, 3.529655662949066*^9}, {3.529656954436404*^9, 3.529656955149273*^9}, 3.529686467808062*^9, {3.52969186833315*^9, 3.529691869421298*^9}, {3.529734606679757*^9, 3.5297346265927267`*^9}, { 3.529762118479937*^9, 3.529762119219027*^9}, {3.529780246934883*^9, 3.529780247266218*^9}, {3.529815738091066*^9, 3.529815747640449*^9}, { 3.5298162916389847`*^9, 3.5298163049870253`*^9}, {3.52981641462272*^9, 3.529816431643127*^9}, {3.529817467083869*^9, 3.529817468394849*^9}, { 3.529828554087121*^9, 3.529828557632923*^9}, {3.5298377116949368`*^9, 3.52983774266153*^9}, {3.529837784012268*^9, 3.5298377873534822`*^9}, { 3.530514165820072*^9, 3.5305141660926123`*^9}, {3.530519836088729*^9, 3.5305198363613567`*^9}, {3.531025988709428*^9, 3.531026013224771*^9}, { 3.531026163231715*^9, 3.531026170141305*^9}, {3.531026321071066*^9, 3.5310263375185328`*^9}, {3.531473128740703*^9, 3.531473157027091*^9}, 3.5315889720604343`*^9, {3.531632664894644*^9, 3.531632665486025*^9}, { 3.532106877147643*^9, 3.5321068809947853`*^9}, {3.5321069981023417`*^9, 3.532107109692912*^9}, {3.53210716548494*^9, 3.532107250402657*^9}, { 3.532107325394965*^9, 3.532107344065818*^9}, 3.532107391450946*^9, { 3.532108064517661*^9, 3.5321080783195753`*^9}, 3.532118597138898*^9, { 3.532275337607031*^9, 3.532275338893052*^9}, {3.532276030704335*^9, 3.5322760746648617`*^9}, {3.532276121560622*^9, 3.5322761225254917`*^9}, { 3.532276278969397*^9, 3.532276288388688*^9}, {3.532277454290394*^9, 3.532277455257657*^9}, {3.532281652034698*^9, 3.532281655600101*^9}, 3.532368085348022*^9, 3.532368344027742*^9, 3.532368376204671*^9, { 3.5323684466232157`*^9, 3.532368494004251*^9}, {3.532368634586609*^9, 3.532368644638282*^9}, {3.5323687247489643`*^9, 3.532368727968088*^9}, { 3.5323688393619337`*^9, 3.532368954385887*^9}, {3.532427670113476*^9, 3.532427671015189*^9}, {3.532427826514886*^9, 3.532427827176198*^9}, { 3.534743131478382*^9, 3.534743131981379*^9}, {3.534743163085504*^9, 3.534743163467737*^9}, {3.5347437584336777`*^9, 3.534743774984112*^9}, { 3.534743970576194*^9, 3.534743998366474*^9}, {3.534745813770913*^9, 3.5347458214501953`*^9}, {3.535178367294001*^9, 3.5351784018195887`*^9}, { 3.535178560434997*^9, 3.5351785605938587`*^9}, {3.535178654971407*^9, 3.535178699302227*^9}, {3.535179478708027*^9, 3.535179488178892*^9}, { 3.536678694995782*^9, 3.536678718167942*^9}, {3.538457139791945*^9, 3.538457141192031*^9}, {3.538457431973028*^9, 3.538457447216082*^9}, { 3.5389985087457438`*^9, 3.53899852265202*^9}, {3.540577983837792*^9, 3.540577985381318*^9}, {3.540578066338348*^9, 3.540578066586337*^9}, { 3.542353416707424*^9, 3.5423534249436626`*^9}, {3.542353567332198*^9, 3.5423535778438873`*^9}, {3.542354860256979*^9, 3.5423548701430817`*^9}, { 3.5423552085091867`*^9, 3.542355215120421*^9}, {3.5423569751576023`*^9, 3.542356989967235*^9}, {3.543900544366967*^9, 3.543900562480075*^9}, { 3.543900808527608*^9, 3.543900821669887*^9}, {3.5478310177935743`*^9, 3.547831018327734*^9}, {3.547831412738585*^9, 3.547831417078899*^9}, { 3.547832115820421*^9, 3.547832159728569*^9}, {3.547834492495185*^9, 3.547834495966104*^9}, {3.547882838738327*^9, 3.54788284384301*^9}, 3.5478899420787277`*^9, {3.547893203971627*^9, 3.5478932167280207`*^9}, { 3.547893471548398*^9, 3.547893474056078*^9}, {3.547893622810422*^9, 3.547893638827704*^9}, {3.547915057194215*^9, 3.5479150911172256`*^9}, { 3.5479154008001833`*^9, 3.547915442347324*^9}, {3.547915986969235*^9, 3.5479160016286697`*^9}, {3.547918165005031*^9, 3.547918175075529*^9}, { 3.547918557011546*^9, 3.5479185610698223`*^9}, {3.547919912481065*^9, 3.547919950553218*^9}, {3.54792005825746*^9, 3.547920059115597*^9}, { 3.547920251150463*^9, 3.5479202544878817`*^9}, {3.5479209990464687`*^9, 3.5479209991657877`*^9}, {3.547921351940798*^9, 3.5479213583834057`*^9}, { 3.547921401615552*^9, 3.547921413543582*^9}, {3.554279521975708*^9, 3.554279528053816*^9}, {3.554279685109273*^9, 3.5542797268268538`*^9}, { 3.554282214311639*^9, 3.5542822401436768`*^9}, {3.5542829294297047`*^9, 3.55428293812042*^9}, {3.5542848313677053`*^9, 3.5542848323264017`*^9}, { 3.554303300350402*^9, 3.554303307136869*^9}, 3.554433616137577*^9, 3.5544336930858603`*^9, {3.554435225628139*^9, 3.554435238852293*^9}, 3.554436484438078*^9, 3.5544375725780153`*^9, 3.554439097742404*^9, 3.554439881037758*^9, 3.554440811267008*^9, 3.554449306103919*^9, 3.55444964496318*^9, 3.554450556844346*^9, 3.554452281522999*^9, 3.5544542247641773`*^9, 3.55445454856222*^9, 3.5544553137896013`*^9, 3.5544554161466503`*^9, {3.554455672060656*^9, 3.554455677652313*^9}, 3.554455909811925*^9, 3.554456063023418*^9, 3.554456228139711*^9, { 3.554462456142919*^9, 3.5544624745308247`*^9}, 3.554472477794661*^9, { 3.554505674614406*^9, 3.554505675320919*^9}, {3.55450581476824*^9, 3.554505815429551*^9}, 3.554524526045466*^9, {3.592657697874805*^9, 3.592657698974822*^9}, {3.592675619277109*^9, 3.5926756197966223`*^9}, { 3.601903638278137*^9, 3.601903677895439*^9}, {3.6019122567367907`*^9, 3.601912280930421*^9}}], Cell[CellGroupData[{ Cell[TextData[{ "Angular Integrals (analytical expressions / recursive relations, (", Cell[BoxData[ FormBox[ SubsuperscriptBox[ RowBox[{"I2", ")"}], RowBox[{"n", ",", RowBox[{"n", "'"}]}], "m"], TraditionalForm]]], " = Integrals\[Theta]f[\[Theta]f][2][n][n\[CloseCurlyQuote]][m+1], (", Cell[BoxData[ FormBox[ SubsuperscriptBox[ RowBox[{"I1", ")"}], RowBox[{"n", ",", RowBox[{"n", "'"}]}], "m"], TraditionalForm]]], " = Integrals\[Theta]f[\[Theta]f][1][n][n\[CloseCurlyQuote]][m+1])" }], "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.601712314484749*^9, 3.6017123405804234`*^9}, {3.60178702671509*^9, 3.60178704659487*^9}, {3.601787114175496*^9, 3.60178714927952*^9}, { 3.601787184118429*^9, 3.601787230579612*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "n", ",", "HmmtT", ",", "HnmT", ",", "\[Pi]nmT", ",", "\[Tau]nmT", ",", "\[Theta]max", ",", "Hnmt", ",", "I1N", ",", "I2N", ",", "I1T", ",", "I2T", ",", "ScatterOAE", ",", "Integrals\[Theta]f"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"I1N", "[", RowBox[{"n_", ",", "np_", ",", "m_", ",", "\[Theta]f_"}], "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["m", "2"], RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", "n", ",", "m"}], "]"}], RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", "np", ",", "m"}], "]"}]}], "+", RowBox[{ RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]", ",", "n", ",", "m"}], "]"}], RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]", ",", "np", ",", "m"}], "]"}]}]}], ")"}], RowBox[{"Sin", "[", "\[Theta]", "]"}]}], ",", RowBox[{"{", RowBox[{"\[Theta]", ",", "0", ",", "\[Theta]f"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"I2N", "[", RowBox[{"n_", ",", "np_", ",", "m_", ",", "\[Theta]f_"}], "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", "n", ",", "m"}], "]"}], RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]", ",", "np", ",", "m"}], "]"}]}], "+", RowBox[{ RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]", ",", "n", ",", "m"}], "]"}], RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", "np", ",", "m"}], "]"}]}]}], ")"}], RowBox[{"Sin", "[", "\[Theta]", "]"}]}], ",", RowBox[{"{", RowBox[{"\[Theta]", ",", "0", ",", "\[Theta]f"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]_", ",", "n_", ",", "m_"}], "]"}], ":=", RowBox[{ RowBox[{"LegendreP", "[", RowBox[{"n", ",", "m", ",", RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}], "/", RowBox[{"Sin", "[", "\[Theta]", "]"}]}]}], ";", RowBox[{ RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]_", ",", "n_", ",", "m_"}], "]"}], ":=", RowBox[{ RowBox[{"n", " ", RowBox[{"Cos", "[", "\[Theta]", "]"}], RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", "n", ",", "m"}], "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"n", "+", "m"}], ")"}], RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]", ",", RowBox[{"n", "-", "1"}], ",", "m"}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Pi]n\[Tau]n\[Theta]", "[", RowBox[{"\[Theta]\[Theta]_", ",", "NNmax_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[Theta]", "=", RowBox[{"N", "[", "\[Theta]\[Theta]", "]"}]}], ",", RowBox[{"Nmax", "=", "NNmax"}], ",", "\[Pi]n", ",", "\[Tau]n", ",", "n"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Pi]n", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", RowBox[{"Nmax", "+", "2"}]}], "}"}]}], "]"}]}], ";", RowBox[{"\[Tau]n", "=", "\[Pi]n"}], ";", RowBox[{ RowBox[{ "\[Pi]n", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "=", "0"}], ";", RowBox[{ RowBox[{ "\[Pi]n", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}], "=", "1"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "3"}], ",", RowBox[{"n", "\[LessEqual]", RowBox[{"Nmax", "+", "1"}]}], ",", RowBox[{"n", "++"}], ",", RowBox[{ RowBox[{ "\[Pi]n", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], "=", RowBox[{ FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{"(", RowBox[{"n", "-", "1"}], ")"}]}], "-", "1"}], ")"}], " ", RowBox[{"Cos", "[", "\[Theta]", "]"}], " ", RowBox[{"\[Pi]n", "\[LeftDoubleBracket]", RowBox[{"n", "-", "1"}], "\[RightDoubleBracket]"}]}], RowBox[{"n", "-", "2"}]], "-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"n", "-", "1"}], ")"}], " ", RowBox[{"\[Pi]n", "\[LeftDoubleBracket]", RowBox[{"n", "-", "2"}], "\[RightDoubleBracket]"}]}], RowBox[{"n", "-", "2"}]]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "2"}], ",", RowBox[{"n", "\[LessEqual]", RowBox[{"Nmax", "+", "1"}]}], ",", RowBox[{"n", "++"}], ",", RowBox[{ RowBox[{ "\[Tau]n", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"n", "-", "1"}], ")"}], " ", RowBox[{"Cos", "[", "\[Theta]", "]"}], " ", RowBox[{ "\[Pi]n", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}]}], "-", RowBox[{"n", " ", RowBox[{"\[Pi]n", "\[LeftDoubleBracket]", RowBox[{"n", "-", "1"}], "\[RightDoubleBracket]"}]}]}]}]}], "]"}], ";", RowBox[{"{", RowBox[{"\[Pi]n", ",", "\[Tau]n"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Integrals\[Theta]f", "[", "\[Theta]ff_", "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[Theta]f", "=", RowBox[{"0.99999", "\[Theta]ff"}]}], ",", "\[Pi]nmT", ",", "\[Tau]nmT", ",", "HmmtT", ",", "HnmT", ",", "I1Tt", ",", "I2Tt", ",", "n"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Pi]nmT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"\[Pi]nm", "[", RowBox[{"\[Theta]f", ",", "n", ",", "m"}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "0", ",", RowBox[{"nmax", "+", "3"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Tau]nmT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"\[Tau]nm", "[", RowBox[{"\[Theta]f", ",", "n", ",", "m"}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "0", ",", RowBox[{"nmax", "+", "3"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"HmmtT", "=", RowBox[{"RecurrenceTable", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Hmmt", "[", "n", "]"}], "\[Equal]", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], "-", "1"}], ")"}], "!!"}], ")"}], "2"], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Cos", "[", "\[Theta]f", "]"}], SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]f", "]"}], RowBox[{"2", "n"}]]}], "+", FractionBox[ RowBox[{"2", "n", " ", RowBox[{"Hmmt", "[", RowBox[{"n", "-", "1"}], "]"}]}], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], "-", "1"}], ")"}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], "-", "3"}], ")"}], "!!"}], ")"}], "2"]}]]}], ")"}]}]}], ",", RowBox[{ RowBox[{"Hmmt", "[", "0", "]"}], "\[Equal]", RowBox[{ RowBox[{"Cos", "[", "\[Theta]f", "]"}], "-", "1"}]}]}], "}"}], ",", "Hmmt", ",", RowBox[{"{", RowBox[{"n", ",", "0", ",", RowBox[{"nmax", "+", "2"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"HnmT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"m", ",", "0", ",", RowBox[{"nmax", "+", "3"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", RowBox[{"nmax", "+", "3"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"m", "=", "0"}], ",", RowBox[{"m", "\[LessEqual]", "nmax"}], ",", RowBox[{"m", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"HnmT", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"m", "+", "1"}], ",", RowBox[{"m", "+", "1"}]}], "\[RightDoubleBracket]"}], "=", RowBox[{"HmmtT", "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "m"}], ",", RowBox[{"n", "\[LessEqual]", "nmax"}], ",", RowBox[{"n", "++"}], ",", RowBox[{ RowBox[{"HnmT", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"n", "+", "2"}], ",", RowBox[{"m", "+", "1"}]}], "\[RightDoubleBracket]"}], "=", RowBox[{"(", RowBox[{ FractionBox["1", RowBox[{"n", "-", "m", "+", "1"}]], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"n", "+", "m", "+", "1"}], ")"}], RowBox[{"HnmT", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"n", "+", "1"}], ",", RowBox[{"m", "+", "1"}]}], "\[RightDoubleBracket]"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]f", "]"}], "2"], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Cos", "[", "\[Theta]f", "]"}], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"n", "+", "m", "+", "1"}], ")"}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], ")"}], "2"]}], "+", RowBox[{ RowBox[{"(", RowBox[{"n", "-", "m", "+", "1"}], ")"}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "2"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], ")"}], "2"]}]}], ")"}]}], "-", RowBox[{"2", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "2"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}]}], ")"}]}]}], ")"}]}], ")"}]}]}], "\[IndentingNewLine]", "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"I1Tt", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"np", "\[Equal]", "n"}], ",", RowBox[{ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"n", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}], RowBox[{ RowBox[{"2", "n"}], "+", "1"}]]}], RowBox[{"HnmT", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"n", "+", "1"}], ",", RowBox[{"m", "+", "1"}]}], "\[RightDoubleBracket]"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]f", "]"}], "2"], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], RowBox[{ RowBox[{"\[Tau]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}]}], ",", FractionBox[ RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]f", "]"}], "2"]}], RowBox[{"(", RowBox[{ RowBox[{"np", RowBox[{"(", RowBox[{"np", "+", "1"}], ")"}], RowBox[{ RowBox[{"\[Tau]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"np", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}], "-", RowBox[{"n", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], RowBox[{ RowBox[{"\[Tau]nmT", "\[LeftDoubleBracket]", RowBox[{"np", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}]}], ")"}]}], RowBox[{ RowBox[{"n", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}], "-", RowBox[{"np", RowBox[{"(", RowBox[{"np", "+", "1"}], ")"}]}]}]]}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", "0", ",", RowBox[{"Min", "[", RowBox[{"np", ",", "n"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"np", ",", "1", ",", "npmax"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "nmax"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"I2Tt", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]f", "]"}], "2"], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"n", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}], RowBox[{ RowBox[{"\[Pi]nmT", "\[LeftDoubleBracket]", RowBox[{"np", "+", "1"}], "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", RowBox[{"m", "+", "1"}], "\[RightDoubleBracket]"}]}], ",", RowBox[{"{", RowBox[{"m", ",", "0", ",", RowBox[{"Min", "[", RowBox[{"np", ",", "n"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"np", ",", "1", ",", "npmax"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "nmax"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "I1Tt", ",", "I2Tt", ",", "HmmtT", ",", "HnmT", ",", "\[Pi]nmT", ",", "\[Tau]nmT"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.5926427792540216`*^9, 3.592642794568759*^9}, { 3.592642830152213*^9, 3.592642862546418*^9}, {3.592643009379498*^9, 3.592643024011511*^9}, {3.592643056415991*^9, 3.592643199214016*^9}, { 3.5926432342311*^9, 3.592643294596505*^9}, {3.592643333593403*^9, 3.592643348119581*^9}, {3.5926434709944353`*^9, 3.592643645452895*^9}, { 3.592644094368346*^9, 3.5926441315715303`*^9}, {3.592644530922711*^9, 3.592644534744465*^9}, {3.592644633303133*^9, 3.5926446336823397`*^9}, { 3.592645198622354*^9, 3.592645198947795*^9}, {3.5926452794559317`*^9, 3.592645364474959*^9}, {3.592645446726019*^9, 3.592645496670294*^9}, { 3.592646305857114*^9, 3.592646364459711*^9}, {3.592646399839666*^9, 3.592646470780963*^9}, {3.592646521059557*^9, 3.592646745466783*^9}, { 3.592646783963552*^9, 3.592646879382978*^9}, {3.592646912950347*^9, 3.59264691819639*^9}, {3.592647078043831*^9, 3.592647121508973*^9}, { 3.5926471825845337`*^9, 3.592647182754076*^9}, {3.592647349944558*^9, 3.592647352952985*^9}, {3.592647545496751*^9, 3.592647584379525*^9}, { 3.592647630771235*^9, 3.592647661298295*^9}, {3.592647802873008*^9, 3.5926478633326817`*^9}, {3.5926479201106243`*^9, 3.59264795627558*^9}, { 3.592648019437244*^9, 3.5926481194789753`*^9}, {3.592648158831698*^9, 3.592648189370249*^9}, {3.5926482312409697`*^9, 3.5926483151896353`*^9}, { 3.5926483668105097`*^9, 3.592648397466782*^9}, {3.592650257710258*^9, 3.5926503630527983`*^9}, {3.592650400737734*^9, 3.592650512787529*^9}, { 3.592650747910375*^9, 3.592650748212998*^9}, {3.592650835436901*^9, 3.592650843028594*^9}, {3.592650874567988*^9, 3.5926509759466333`*^9}, { 3.5926510620116158`*^9, 3.592651063379702*^9}, {3.592651227842004*^9, 3.592651260313429*^9}, {3.5926514778133793`*^9, 3.592651497455518*^9}, { 3.592651543597643*^9, 3.592651546077271*^9}, {3.5926515928048553`*^9, 3.5926515929195423`*^9}, {3.5926516364521646`*^9, 3.5926516374428663`*^9}, {3.592651794409423*^9, 3.5926518214644547`*^9}, {3.5926519279105186`*^9, 3.592651962284721*^9}, {3.5926520823809032`*^9, 3.5926520839768763`*^9}, { 3.592652201494802*^9, 3.5926522443608913`*^9}, {3.592652310903236*^9, 3.592652402966281*^9}, {3.5926524840368567`*^9, 3.592652487317605*^9}, { 3.59265270324225*^9, 3.592652764953279*^9}, {3.59265281117514*^9, 3.592653080414287*^9}, {3.592653119697412*^9, 3.592653294573493*^9}, { 3.592653363242373*^9, 3.5926535235393047`*^9}, {3.592653559158606*^9, 3.592653617177002*^9}, {3.592653743207964*^9, 3.5926537833746634`*^9}, { 3.59265383460118*^9, 3.5926538560246773`*^9}, {3.5926539075288687`*^9, 3.592654024852335*^9}, {3.59265406107618*^9, 3.5926541464246197`*^9}, { 3.5926543252943287`*^9, 3.5926549817123003`*^9}, {3.592655197839485*^9, 3.5926551984100723`*^9}, {3.59265550011287*^9, 3.5926555035544443`*^9}, { 3.592655565151979*^9, 3.592655593062307*^9}, {3.592655656659622*^9, 3.5926556593724546`*^9}, {3.5926557424711123`*^9, 3.592655827187357*^9}, { 3.592655870541663*^9, 3.592655921296207*^9}, {3.5926559747943687`*^9, 3.592655996817231*^9}, {3.59265612143718*^9, 3.5926561301584377`*^9}, { 3.592656166837747*^9, 3.592656476093549*^9}, {3.592656520799947*^9, 3.592656548496765*^9}, {3.592656796893553*^9, 3.592656852412436*^9}, { 3.592656945682494*^9, 3.5926570376573257`*^9}, {3.592657087142173*^9, 3.592657108760231*^9}, {3.592657201365192*^9, 3.592657300394661*^9}, { 3.592657335079443*^9, 3.592657385155903*^9}, {3.592657571636003*^9, 3.592657589082115*^9}, {3.592657719724681*^9, 3.5926577410174923`*^9}, { 3.5926577893719473`*^9, 3.592657859516509*^9}, {3.592657894040359*^9, 3.592657894492758*^9}, {3.592658095970468*^9, 3.592658100634058*^9}, 3.592658174096962*^9, {3.5926583526315804`*^9, 3.592658364240327*^9}, { 3.592658657363935*^9, 3.592658755184544*^9}, {3.592658791968504*^9, 3.5926588950187387`*^9}, {3.5926589479830437`*^9, 3.592659015936792*^9}, { 3.5926591075273733`*^9, 3.592659113497985*^9}, {3.592659157955681*^9, 3.592659169513733*^9}, {3.592659230379725*^9, 3.592659257466988*^9}, { 3.592659310855897*^9, 3.592659321985558*^9}, {3.592659400547645*^9, 3.592659429311173*^9}, {3.5926594814286337`*^9, 3.5926594815840387`*^9}, { 3.5926595727868013`*^9, 3.592659573046012*^9}, {3.592733943106533*^9, 3.592733966880623*^9}, {3.59385786976322*^9, 3.593857881564617*^9}, { 3.593860760665586*^9, 3.5938607631291103`*^9}, {3.593860912682075*^9, 3.5938609139556627`*^9}, {3.597407011613701*^9, 3.597407022203287*^9}, { 3.597411063896636*^9, 3.597411086526745*^9}, {3.597411567361825*^9, 3.597411571641241*^9}, {3.597487268429966*^9, 3.5974872697733717`*^9}, { 3.598119811636578*^9, 3.5981198704784403`*^9}, 3.598157766969006*^9, 3.5981586229992647`*^9, 3.598158668478476*^9, {3.598166957659184*^9, 3.598166960768156*^9}, {3.5981670974678097`*^9, 3.5981671016472683`*^9}, { 3.598167148234441*^9, 3.598167151924514*^9}, {3.598167197607127*^9, 3.5981672211867743`*^9}, {3.5981672538738117`*^9, 3.59816727447323*^9}, { 3.598781599588106*^9, 3.5987816165691643`*^9}, 3.598781646882066*^9, { 3.601787015915666*^9, 3.601787016598316*^9}, {3.601807411178842*^9, 3.6018074115108232`*^9}, {3.601811262658208*^9, 3.6018112634595633`*^9}, { 3.601811347205225*^9, 3.601811348286314*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Beam Shape Coefficients (BSCs) for a Gaussian Beam\n", StyleBox["\[OpenCurlyDoubleQuote]Plane-wave spectrum of electromagnetic \ beams\[CloseCurlyDoubleQuote], A. Doicu and T. Wriedt, Opt. Commun. 136, \ 114-124, 1997\n\[OpenCurlyDoubleQuote]Computation of the beam shape \ coefficients in the Generalized Lorenz-Mie Theory by using the translational \ addition theorem for spherical vector wave functions\[CloseCurlyDoubleQuote], \ A. Doicu and T. Wriedt, 1997\nx0 = x/", FontSize->12], Cell[BoxData[ FormBox[ SubscriptBox["w", "0"], TraditionalForm]], "None", FontSize->12], StyleBox[", y0 = y/", FontSize->12], Cell[BoxData[ FormBox[ SubscriptBox["w", "0"], TraditionalForm]], "None", FontSize->12], StyleBox[", ", FontSize->12], Cell[BoxData[ FormBox[ SubscriptBox["z", "0"], TraditionalForm]], "None", FontSize->12], StyleBox[" = z/l, ", FontSize->12], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"l", "=", RowBox[{ RowBox[{"2", "zR"}], "=", RowBox[{ RowBox[{"k", " ", SuperscriptBox["\[Omega]0", "2"]}], "=", RowBox[{ SubscriptBox["n", "med"], "med", " ", "2", "\[Pi]", " ", RowBox[{ SuperscriptBox[ SubscriptBox["w", "0"], "2"], "/", "\[Lambda]"}]}]}]}]}], ","}], TraditionalForm]], FontSize->12], StyleBox[" s=", FontSize->12], Cell[BoxData[ FormBox[ SubscriptBox["w", "0"], TraditionalForm]], "None", FontSize->12], StyleBox["/l=1/(k ", FontSize->12], Cell[BoxData[ FormBox[ SubscriptBox["w", "0"], TraditionalForm]], "None", FontSize->12], StyleBox["), \[Gamma]=2 zp/", FontSize->12], Cell[BoxData[ FormBox[ RowBox[{ SubscriptBox["w", "0"], "=", RowBox[{"2", RowBox[{"z0", "/", "s"}]}]}], TraditionalForm]], "None", FontSize->12] }], "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.601712314484749*^9, 3.6017123405804234`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Qb", "[", RowBox[{"n_", ",", "s_", ",", "\[Gamma]_"}], "]"}], ":=", FractionBox["1", RowBox[{"\[ImaginaryI]", "-", RowBox[{"s", " ", "\[Gamma]"}]}]]}], ";", RowBox[{ RowBox[{"\[Rho]n", "[", RowBox[{"n_", ",", "s_"}], "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{"n", "+", "0.5"}], ")"}], "s"}]}], ";", RowBox[{ RowBox[{"\[Rho]0", "[", RowBox[{"x0_", ",", "y0_"}], "]"}], ":=", RowBox[{"Sqrt", "[", RowBox[{ SuperscriptBox["x0", "2"], "+", SuperscriptBox["y0", "2"]}], "]"}]}], ";", RowBox[{"eps", "=", "0.000000001"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Knm", "[", RowBox[{"n_", ",", "m_"}], "]"}], ":=", RowBox[{"If", "[", RowBox[{ RowBox[{"m", "\[Equal]", "0"}], ",", RowBox[{"n", RowBox[{ RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}], "/", RowBox[{"(", RowBox[{"n", "+", "0.5"}], ")"}]}]}], ",", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "I"}], ")"}], RowBox[{"Abs", "[", "m", "]"}]], RowBox[{"I", "/", SuperscriptBox[ RowBox[{"(", RowBox[{"n", "+", "0.5"}], ")"}], RowBox[{ RowBox[{"Abs", "[", "m", "]"}], "-", "1"}]]}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"gnmTMAddT", "[", RowBox[{ "n_", ",", "m_", ",", "s_", ",", "\[Gamma]_", ",", "z0_", ",", "x0_", ",", "y0_"}], "]"}], ":=", RowBox[{ RowBox[{"Knm", "[", RowBox[{"n", ",", "m"}], "]"}], RowBox[{"(", RowBox[{"I", " ", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "I"}], " ", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], "2"], "+", SuperscriptBox[ RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}], "2"]}], ")"}]}], "]"}]}], ")"}], FractionBox[ RowBox[{"Exp", "[", RowBox[{"I", " ", RowBox[{"\[Gamma]", "/", RowBox[{"(", RowBox[{"2", "s"}], ")"}]}]}], "]"}], "2"], RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Exp", "[", RowBox[{"I", RowBox[{"(", RowBox[{"m", "-", "1"}], ")"}], "\[Phi]"}], "]"}], RowBox[{"BesselJ", "[", RowBox[{ RowBox[{"m", "-", "1"}], ",", RowBox[{"2", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}]}]}], "]"}]}], "+", RowBox[{ RowBox[{"Exp", "[", RowBox[{"I", RowBox[{"(", RowBox[{"m", "+", "1"}], ")"}], "\[Phi]"}], "]"}], RowBox[{"BesselJ", "[", RowBox[{ RowBox[{"m", "+", "1"}], ",", RowBox[{"2", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}]}]}], "]"}]}]}], ")"}], "/.", RowBox[{"{", RowBox[{"\[Phi]", "\[Rule]", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"y0", "+", "eps"}], ",", RowBox[{"x0", "+", "eps"}]}], "]"}]}], "}"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"gnmTEAddT", "[", RowBox[{ "n_", ",", "m_", ",", "s_", ",", "\[Gamma]_", ",", "z0_", ",", "x0_", ",", "y0_"}], "]"}], ":=", RowBox[{ SuperscriptBox["I", RowBox[{"-", "1"}]], RowBox[{"Knm", "[", RowBox[{"n", ",", "m"}], "]"}], RowBox[{"(", RowBox[{"I", " ", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "I"}], " ", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], "2"], "+", SuperscriptBox[ RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}], "2"]}], ")"}]}], "]"}]}], ")"}], FractionBox[ RowBox[{"Exp", "[", RowBox[{"I", " ", RowBox[{"\[Gamma]", "/", RowBox[{"(", RowBox[{"2", "s"}], ")"}]}]}], "]"}], "2"], RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Exp", "[", RowBox[{"I", RowBox[{"(", RowBox[{"m", "-", "1"}], ")"}], "\[Phi]"}], "]"}], RowBox[{"BesselJ", "[", RowBox[{ RowBox[{"m", "-", "1"}], ",", RowBox[{"2", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}]}]}], "]"}]}], "-", RowBox[{ RowBox[{"Exp", "[", RowBox[{"I", RowBox[{"(", RowBox[{"m", "+", "1"}], ")"}], "\[Phi]"}], "]"}], RowBox[{"BesselJ", "[", RowBox[{ RowBox[{"m", "+", "1"}], ",", RowBox[{"2", RowBox[{"Qb", "[", RowBox[{"n", ",", "s", ",", "\[Gamma]"}], "]"}], RowBox[{"\[Rho]0", "[", RowBox[{"x0", ",", "y0"}], "]"}], RowBox[{"\[Rho]n", "[", RowBox[{"n", ",", "s"}], "]"}]}]}], "]"}]}]}], ")"}], "/.", RowBox[{"{", RowBox[{"\[Phi]", "\[Rule]", RowBox[{"ArcTan", "[", RowBox[{ RowBox[{"y0", "+", "eps"}], ",", RowBox[{"x0", "+", "eps"}]}], "]"}]}], "}"}]}], ")"}]}]}], ";"}]}], "Input", InitializationCell->True, CellChangeTimes->CompressedData[" 1:eJxTTMoPSmViYGAQA2IQ/eRUZoB+8CvHhBenkkD0v9R3mSDaZyN7Hoiu6ZS+ B6LTXWrug+iwU7mvQHTo2eOvwep2iH4B0XPMZn0F0XEVoUEGQNpc42wIiD6+ KyMaRMdccYwD0fPTT6SBaAPrhmYQPcmXqQNEl8rf6QfRf9jZpoPoupDVs0C0 9s/2hSA6uD1+HYhWedO2GUSvmXRhD4g20jt0EUTL3PlzCUQvsux6BqIrPNa8 BdFipyKPvD//yvHYjAIw3aBoeg5Eix7LBdPPrveqfwDSHvabNED05R/rnUG0 rkOlJ4i+cy00FUTfsJ0MpvW8tJl+AuloiXww/UPnMweIvsO/nAtET5dREgbR KdNSwDQAwNW26Q== "]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Functions for the Data-set generation\ \>", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.601807420307987*^9, 3.6018074311131268`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Nn", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], "+", "1"}], ")"}], "/", RowBox[{"(", RowBox[{"n", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}], ")"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "10000"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nmax", "=", "150"}], ";", RowBox[{"npmax", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Processing1", "[", RowBox[{ "I1I22_", ",", "sii_", ",", "\[Theta]ii_", ",", "x00_", ",", "z00_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"I1I2", "=", "I1I22"}], ",", RowBox[{"si", "=", "sii"}], ",", RowBox[{"\[Theta]i", "=", "\[Theta]ii"}], ",", RowBox[{"x0", "=", "x00"}], ",", RowBox[{"z0", "=", RowBox[{"-", "z00"}]}], ",", "Ax", ",", "gnmTET", ",", "gnmTMT", ",", "I1T", ",", "I2T"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"I1T", "=", RowBox[{ RowBox[{ RowBox[{ "I1I2", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}]}], ";", RowBox[{"I2T", "=", RowBox[{ RowBox[{ RowBox[{ "I1I2", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}]}], ";", " ", RowBox[{"(*", RowBox[{"precomputed", " ", "Angular", " ", "integrals"}], "*)"}], "\[IndentingNewLine]", RowBox[{"gnmTET", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"gnmTEAddT", "[", RowBox[{"n", ",", "m", ",", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], ",", RowBox[{"2", RowBox[{"z0", "/", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}]}]}], ",", "z0", ",", "x0", ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", RowBox[{"nmax", "+", "2"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"gnmTMT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"gnmTMAddT", "[", RowBox[{"n", ",", "m", ",", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], ",", RowBox[{"2", RowBox[{"z0", "/", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}]}]}], ",", "z0", ",", "x0", ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"m", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", RowBox[{"nmax", "+", "2"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Ax", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ "Nn", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], "/", "2"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"gnmTMT", "\[LeftDoubleBracket]", RowBox[{"n", ",", RowBox[{"2", "+", "m"}]}], "\[RightDoubleBracket]"}], RowBox[{"Conjugate", "[", RowBox[{"gnmTMT", "\[LeftDoubleBracket]", RowBox[{"1", ",", RowBox[{"2", "+", "m"}]}], "\[RightDoubleBracket]"}], "]"}], RowBox[{"I1T", "\[LeftDoubleBracket]", RowBox[{"n", ",", "1", ",", RowBox[{ RowBox[{"Abs", "[", "m", "]"}], "+", "1"}]}], "\[RightDoubleBracket]"}]}], "+", RowBox[{"I", " ", "m", " ", RowBox[{"gnmTET", "\[LeftDoubleBracket]", RowBox[{"n", ",", RowBox[{"2", "+", "m"}]}], "\[RightDoubleBracket]"}], RowBox[{"Conjugate", "[", RowBox[{"gnmTMT", "\[LeftDoubleBracket]", RowBox[{"1", ",", RowBox[{"2", "+", "m"}]}], "\[RightDoubleBracket]"}], "]"}], RowBox[{"I2T", "\[LeftDoubleBracket]", RowBox[{"n", ",", "1", ",", RowBox[{ RowBox[{"Abs", "[", "m", "]"}], "+", "1"}]}], "\[RightDoubleBracket]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "nmax"}], "}"}], ",", RowBox[{"{", RowBox[{"m", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Re", "[", "Ax", "]"}], ",", RowBox[{"-", RowBox[{"Im", "[", "Ax", "]"}]}]}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ProcessingInt", "[", RowBox[{"sTT_", ",", "\[Theta]Maxby\[Theta]divTT_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"sT", "=", "sTT"}], ",", RowBox[{ "\[Theta]Maxby\[Theta]divT", "=", "\[Theta]Maxby\[Theta]divTT"}], ",", "I1TT", ",", "I2TT", ",", "ki", ",", "si", ",", "\[Theta]i", ",", "\[Theta]Maxby\[Theta]divi", ",", "Integrals\[Theta]"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"I1TT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"\[Theta]i", ",", "1", ",", "Ntheta"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"si", ",", "1", ",", "Ns"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"I2TT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"\[Theta]i", ",", "1", ",", "Ntheta"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"si", ",", "1", ",", "Ns"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"\[Theta]i", "=", "1"}], ",", RowBox[{"\[Theta]i", "\[LessEqual]", "Ntheta"}], ",", RowBox[{"\[Theta]i", "++"}], ",", RowBox[{"For", "[", RowBox[{ RowBox[{"si", "=", "1"}], ",", RowBox[{"si", "\[LessEqual]", "Ns"}], ",", RowBox[{"si", "++"}], ",", RowBox[{ RowBox[{"Integrals\[Theta]", "=", RowBox[{"Integrals\[Theta]f", "[", RowBox[{ RowBox[{ "\[Theta]Maxby\[Theta]divT", "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}], "*", "2", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}]}], "]"}]}], ";", RowBox[{ RowBox[{ RowBox[{ "I1TT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}], "=", RowBox[{ "Integrals\[Theta]", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], ";", RowBox[{ RowBox[{ RowBox[{ "I2TT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}], "=", RowBox[{ "Integrals\[Theta]", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}]}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"I1TT", ",", "I2TT"}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Processing", "[", RowBox[{"AllI1I22_", ",", "x00_", ",", "z00_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"AllI1I2", "=", "AllI1I22"}], ",", RowBox[{"x0", "=", "x00"}], ",", RowBox[{"z0", "=", RowBox[{"-", "z00"}]}], ",", "\[Theta]i", ",", "si", ",", "T"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Processing1", "[", RowBox[{ "AllI1I2", ",", "si", ",", "\[Theta]i", ",", "x0", ",", "z0"}], "]"}], ",", RowBox[{"{", RowBox[{"\[Theta]i", ",", "1", ",", "Ntheta"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"si", ",", "1", ",", "Ns"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";", "\[IndentingNewLine]", "T"}]}], "\[IndentingNewLine]", "]"}]}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.601470225519223*^9, 3.6014702586793547`*^9}, { 3.601470725947863*^9, 3.6014707271387053`*^9}, {3.601476073722938*^9, 3.601476146168597*^9}, {3.601477331417449*^9, 3.601477332521552*^9}, { 3.60147960740593*^9, 3.601479614079176*^9}, {3.60147964440301*^9, 3.601479645275077*^9}, 3.6014800899903727`*^9, {3.6014802446916437`*^9, 3.601480273891508*^9}, {3.601494864709074*^9, 3.6014948713118362`*^9}, { 3.601494954348072*^9, 3.6014949603403387`*^9}, {3.601496656506143*^9, 3.6014966709773817`*^9}, {3.601496915062845*^9, 3.60149691721797*^9}, 3.601542063801547*^9, {3.601542634830572*^9, 3.601542634941766*^9}, { 3.601544388643242*^9, 3.6015444081138763`*^9}, {3.601544485466201*^9, 3.601544488161111*^9}, {3.601544586527172*^9, 3.601544738666931*^9}, { 3.6015448627620792`*^9, 3.601544877420809*^9}, {3.6015449474417963`*^9, 3.6015453757001877`*^9}, {3.601545413854149*^9, 3.601545482759251*^9}, { 3.601545543681039*^9, 3.601545548679089*^9}, {3.601545584987825*^9, 3.601545585523888*^9}, 3.601545633774952*^9, 3.6015457065333433`*^9, { 3.6015488433835*^9, 3.601548881783435*^9}, {3.6015491051271563`*^9, 3.601549352327194*^9}, {3.6015494099582987`*^9, 3.601549620362048*^9}, { 3.601549729553275*^9, 3.601549807656212*^9}, {3.601549851093926*^9, 3.601550016144586*^9}, {3.601550359235393*^9, 3.601550431162404*^9}, 3.6015504873314333`*^9, {3.601550563225851*^9, 3.6015506769359913`*^9}, { 3.6015507138384*^9, 3.601550719668439*^9}, {3.601550770828*^9, 3.601550784403097*^9}, 3.6015509009756117`*^9, {3.601550998348071*^9, 3.6015510081897907`*^9}, {3.60155107053459*^9, 3.6015511355962057`*^9}, { 3.601551210526455*^9, 3.6015512293724747`*^9}, {3.6015512735961933`*^9, 3.601551354642069*^9}, {3.601551548701129*^9, 3.601551561284626*^9}, { 3.601563144311461*^9, 3.601563153311468*^9}, {3.601563276275199*^9, 3.601563277264164*^9}, {3.6015633566471643`*^9, 3.601563356992647*^9}, { 3.601567154309621*^9, 3.601567207123556*^9}, {3.601712052808775*^9, 3.60171205537523*^9}, {3.601712358443159*^9, 3.601712370709833*^9}, { 3.601712405147663*^9, 3.601712413169304*^9}, {3.601721126633636*^9, 3.601721128787632*^9}, {3.6017212493601522`*^9, 3.601721262117486*^9}, 3.601725723980811*^9, {3.601731583102006*^9, 3.601731609479726*^9}, { 3.601731977877617*^9, 3.6017319801726427`*^9}, {3.6017865185772343`*^9, 3.601786524995343*^9}, {3.6017868571259193`*^9, 3.601786857357114*^9}, { 3.6018073643662443`*^9, 3.601807387434441*^9}, {3.601807449282332*^9, 3.601807454071175*^9}, 3.6018075894218073`*^9, {3.601901394766964*^9, 3.601901429381723*^9}, {3.6019197123414497`*^9, 3.601919712876449*^9}, 3.601919793079715*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Parameters of the generated Data\ \>", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122126282387`*^9, 3.601712220915333*^9}, {3.6018073991212873`*^9, 3.601807400193982*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"{", RowBox[{ RowBox[{"\[Theta]Maxby\[Theta]divStart", "=", "0.01"}], ",", RowBox[{"\[Theta]Maxby\[Theta]divMax", "=", "2"}], ",", RowBox[{"d\[Theta]Maxby\[Theta]div", "=", "0.2"}], ",", RowBox[{"\[Theta]Maxby\[Theta]divT", "=", RowBox[{"Table", "[", RowBox[{"\[Theta]", ",", RowBox[{"{", RowBox[{ "\[Theta]", ",", "\[Theta]Maxby\[Theta]divStart", ",", "\[Theta]Maxby\[Theta]divMax", ",", "d\[Theta]Maxby\[Theta]div"}], "}"}]}], "]"}]}], ",", RowBox[{"Ntheta", "=", RowBox[{"Length", "[", "\[Theta]Maxby\[Theta]divT", "]"}]}]}], "}"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"sStart", "=", "0.05"}], ",", RowBox[{"sMax", "=", "0.3"}], ",", RowBox[{"dds", "=", "0.025"}], ",", RowBox[{"sT", "=", RowBox[{"Table", "[", RowBox[{"s", ",", RowBox[{"{", RowBox[{"s", ",", "sStart", ",", "sMax", ",", "dds"}], "}"}]}], "]"}]}], ",", RowBox[{"Ns", "=", RowBox[{"Length", "[", "sT", "]"}]}]}], "}"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"XMax", "=", "4"}], ";", RowBox[{"dx", "=", "0.2"}], ";", RowBox[{"Nx", "=", RowBox[{"XMax", "/", "dx"}]}], ";", RowBox[{"ZMax", "=", "3"}], ";", RowBox[{"dz", "=", "0.1"}], ";", RowBox[{"Nz", "=", RowBox[{"ZMax", "/", "dz"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Nx", "+", "1"}], ",", RowBox[{"Nz", "+", "1"}]}], "}"}]}], "Input", CellChangeTimes->{ 3.601919720278648*^9, {3.6019588187610397`*^9, 3.601958821920426*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.01`", ",", "2", ",", "0.2`", ",", RowBox[{"{", RowBox[{ "0.01`", ",", "0.21000000000000002`", ",", "0.41000000000000003`", ",", "0.6100000000000001`", ",", "0.81`", ",", "1.01`", ",", "1.2100000000000002`", ",", "1.4100000000000001`", ",", "1.61`", ",", "1.81`"}], "}"}], ",", "10"}], "}"}]], "Output", CellChangeTimes->{ 3.601807402666354*^9, 3.601807436781579*^9, 3.601807467003783*^9, 3.6018075923368692`*^9, 3.6018097218357067`*^9, 3.601810971969933*^9, 3.601811367703339*^9, 3.6019014376196127`*^9, {3.6019197161990833`*^9, 3.601919725019939*^9}, 3.601919795741974*^9, 3.601958822380002*^9, { 3.601958887098686*^9, 3.601958892700366*^9}, 3.6019952356218224`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.05`", ",", "0.3`", ",", "0.025`", ",", RowBox[{"{", RowBox[{ "0.05`", ",", "0.07500000000000001`", ",", "0.1`", ",", "0.125`", ",", "0.15000000000000002`", ",", "0.175`", ",", "0.2`", ",", "0.22500000000000003`", ",", "0.25`", ",", "0.275`", ",", "0.3`"}], "}"}], ",", "11"}], "}"}]], "Output", CellChangeTimes->{ 3.601807402666354*^9, 3.601807436781579*^9, 3.601807467003783*^9, 3.6018075923368692`*^9, 3.6018097218357067`*^9, 3.601810971969933*^9, 3.601811367703339*^9, 3.6019014376196127`*^9, {3.6019197161990833`*^9, 3.601919725019939*^9}, 3.601919795741974*^9, 3.601958822380002*^9, { 3.601958887098686*^9, 3.601958892700366*^9}, 3.601995235623186*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"21.`", ",", "31.`"}], "}"}]], "Output", CellChangeTimes->{ 3.601807402666354*^9, 3.601807436781579*^9, 3.601807467003783*^9, 3.6018075923368692`*^9, 3.6018097218357067`*^9, 3.601810971969933*^9, 3.601811367703339*^9, 3.6019014376196127`*^9, {3.6019197161990833`*^9, 3.601919725019939*^9}, 3.601919795741974*^9, 3.601958822380002*^9, { 3.601958887098686*^9, 3.601958892700366*^9}, 3.6019952356240873`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Structure of the generated Data\ \>", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122126282387`*^9, 3.601712220915333*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ "\[Theta]Maxby\[Theta]divT", "\[LeftDoubleBracket]", "\[Theta]i", "\[RightDoubleBracket]"}], ",", RowBox[{ "sT", "\[LeftDoubleBracket]", "si", "\[RightDoubleBracket]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"\[Theta]i", ",", "1", ",", "Ntheta"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"si", ",", "1", ",", "Ns"}], "}"}]}], "]"}], ",", "2"}], "]"}]], "Input", CellChangeTimes->{{3.601559599100268*^9, 3.6015596044295483`*^9}, { 3.6015609283277903`*^9, 3.601560932465993*^9}, {3.601712465592617*^9, 3.601712473155444*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.01`", ",", "0.05`", ",", "0.21000000000000002`", ",", "0.05`", ",", "0.41000000000000003`", ",", "0.05`", ",", "0.6100000000000001`", ",", "0.05`", ",", "0.81`", ",", "0.05`", ",", "1.01`", ",", "0.05`", ",", "1.2100000000000002`", ",", "0.05`", ",", "1.4100000000000001`", ",", "0.05`", ",", "1.61`", ",", "0.05`", ",", "1.81`", ",", "0.05`", ",", "0.01`", ",", "0.07500000000000001`", ",", "0.21000000000000002`", ",", "0.07500000000000001`", ",", "0.41000000000000003`", ",", "0.07500000000000001`", ",", "0.6100000000000001`", ",", "0.07500000000000001`", ",", "0.81`", ",", "0.07500000000000001`", ",", "1.01`", ",", "0.07500000000000001`", ",", "1.2100000000000002`", ",", "0.07500000000000001`", ",", "1.4100000000000001`", ",", "0.07500000000000001`", ",", "1.61`", ",", "0.07500000000000001`", ",", "1.81`", ",", "0.07500000000000001`", ",", "0.01`", ",", "0.1`", ",", "0.21000000000000002`", ",", "0.1`", ",", "0.41000000000000003`", ",", "0.1`", ",", "0.6100000000000001`", ",", "0.1`", ",", "0.81`", ",", "0.1`", ",", "1.01`", ",", "0.1`", ",", "1.2100000000000002`", ",", "0.1`", ",", "1.4100000000000001`", ",", "0.1`", ",", "1.61`", ",", "0.1`", ",", "1.81`", ",", "0.1`", ",", "0.01`", ",", "0.125`", ",", "0.21000000000000002`", ",", "0.125`", ",", "0.41000000000000003`", ",", "0.125`", ",", "0.6100000000000001`", ",", "0.125`", ",", "0.81`", ",", "0.125`", ",", "1.01`", ",", "0.125`", ",", "1.2100000000000002`", ",", "0.125`", ",", "1.4100000000000001`", ",", "0.125`", ",", "1.61`", ",", "0.125`", ",", "1.81`", ",", "0.125`", ",", "0.01`", ",", "0.15000000000000002`", ",", "0.21000000000000002`", ",", "0.15000000000000002`", ",", "0.41000000000000003`", ",", "0.15000000000000002`", ",", "0.6100000000000001`", ",", "0.15000000000000002`", ",", "0.81`", ",", "0.15000000000000002`", ",", "1.01`", ",", "0.15000000000000002`", ",", "1.2100000000000002`", ",", "0.15000000000000002`", ",", "1.4100000000000001`", ",", "0.15000000000000002`", ",", "1.61`", ",", "0.15000000000000002`", ",", "1.81`", ",", "0.15000000000000002`", ",", "0.01`", ",", "0.175`", ",", "0.21000000000000002`", ",", "0.175`", ",", "0.41000000000000003`", ",", "0.175`", ",", "0.6100000000000001`", ",", "0.175`", ",", "0.81`", ",", "0.175`", ",", "1.01`", ",", "0.175`", ",", "1.2100000000000002`", ",", "0.175`", ",", "1.4100000000000001`", ",", "0.175`", ",", "1.61`", ",", "0.175`", ",", "1.81`", ",", "0.175`", ",", "0.01`", ",", "0.2`", ",", "0.21000000000000002`", ",", "0.2`", ",", "0.41000000000000003`", ",", "0.2`", ",", "0.6100000000000001`", ",", "0.2`", ",", "0.81`", ",", "0.2`", ",", "1.01`", ",", "0.2`", ",", "1.2100000000000002`", ",", "0.2`", ",", "1.4100000000000001`", ",", "0.2`", ",", "1.61`", ",", "0.2`", ",", "1.81`", ",", "0.2`", ",", "0.01`", ",", "0.22500000000000003`", ",", "0.21000000000000002`", ",", "0.22500000000000003`", ",", "0.41000000000000003`", ",", "0.22500000000000003`", ",", "0.6100000000000001`", ",", "0.22500000000000003`", ",", "0.81`", ",", "0.22500000000000003`", ",", "1.01`", ",", "0.22500000000000003`", ",", "1.2100000000000002`", ",", "0.22500000000000003`", ",", "1.4100000000000001`", ",", "0.22500000000000003`", ",", "1.61`", ",", "0.22500000000000003`", ",", "1.81`", ",", "0.22500000000000003`", ",", "0.01`", ",", "0.25`", ",", "0.21000000000000002`", ",", "0.25`", ",", "0.41000000000000003`", ",", "0.25`", ",", "0.6100000000000001`", ",", "0.25`", ",", "0.81`", ",", "0.25`", ",", "1.01`", ",", "0.25`", ",", "1.2100000000000002`", ",", "0.25`", ",", "1.4100000000000001`", ",", "0.25`", ",", "1.61`", ",", "0.25`", ",", "1.81`", ",", "0.25`", ",", "0.01`", ",", "0.275`", ",", "0.21000000000000002`", ",", "0.275`", ",", "0.41000000000000003`", ",", "0.275`", ",", "0.6100000000000001`", ",", "0.275`", ",", "0.81`", ",", "0.275`", ",", "1.01`", ",", "0.275`", ",", "1.2100000000000002`", ",", "0.275`", ",", "1.4100000000000001`", ",", "0.275`", ",", "1.61`", ",", "0.275`", ",", "1.81`", ",", "0.275`", ",", "0.01`", ",", "0.3`", ",", "0.21000000000000002`", ",", "0.3`", ",", "0.41000000000000003`", ",", "0.3`", ",", "0.6100000000000001`", ",", "0.3`", ",", "0.81`", ",", "0.3`", ",", "1.01`", ",", "0.3`", ",", "1.2100000000000002`", ",", "0.3`", ",", "1.4100000000000001`", ",", "0.3`", ",", "1.61`", ",", "0.3`", ",", "1.81`", ",", "0.3`"}], "}"}]], "Output", CellChangeTimes->{ 3.601559605064817*^9, 3.6015609329692793`*^9, 3.601567078865581*^9, { 3.6017124539367456`*^9, 3.601712473473742*^9}, 3.6017126388881283`*^9, 3.601725626782836*^9, 3.601725728745777*^9, 3.601731825300561*^9, 3.601732027989829*^9, 3.601787476622573*^9, 3.601807594578669*^9, 3.601809724789007*^9, 3.601810975189519*^9, 3.601901439789981*^9, 3.601919727171582*^9, 3.6019197979309196`*^9, 3.601958898115006*^9, 3.601995239041801*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Compute all angular Integrals first (takes about 1min)", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.601712224219471*^9, 3.6017122317962017`*^9}, {3.601725767759439*^9, 3.6017257727915773`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"AllI1I2", "=", RowBox[{"ProcessingInt", "[", RowBox[{"sT", ",", "\[Theta]Maxby\[Theta]divT"}], "]"}]}], ";"}]], "Input",\ CellChangeTimes->{{3.601551467342655*^9, 3.6015514773028*^9}, 3.601563158397647*^9, 3.6015632912496023`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Test", "=", RowBox[{"Timing", "[", RowBox[{"Processing", "[", RowBox[{"AllI1I2", ",", "1", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Test", "[", RowBox[{"[", "1", "]"}], "]"}], "*", "Nx", "*", RowBox[{"Nz", "/", RowBox[{"(", RowBox[{"60", "*", "60"}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{"Test", "[", RowBox[{"[", "2", "]"}], "]"}]}], "Input", CellChangeTimes->{{3.6018110290867653`*^9, 3.601811030110854*^9}, { 3.601901457378273*^9, 3.601901518211092*^9}}], Cell[BoxData["3.548021`"], "Output", CellChangeTimes->{ 3.601811162685341*^9, 3.601811375229395*^9, 3.601811516137822*^9, { 3.6019014976324778`*^9, 3.60190152563608*^9}, 3.6019197793013983`*^9, 3.601919876150498*^9, 3.601958985667324*^9, 3.6019953342586184`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.00002776435083996115`", ",", RowBox[{"-", "0.000023706595516817882`"}], ",", "0.012136853733925682`", ",", RowBox[{"-", "0.009486498959300944`"}], ",", "0.04416318067679872`", ",", RowBox[{"-", "0.026994435285479353`"}], ",", "0.0859525019483467`", ",", RowBox[{"-", "0.03471894900753904`"}], ",", "0.11864643391616587`", ",", RowBox[{"-", "0.025079035198804674`"}], ",", "0.13034607079026403`", ",", RowBox[{"-", "0.010369343174992753`"}], ",", "0.13033645352047465`", ",", RowBox[{"-", "0.005392386023988529`"}], ",", "0.133407445230775`", ",", RowBox[{"-", "0.006177860058997012`"}], ",", "0.1374436651852409`", ",", RowBox[{"-", "0.0025695056777874166`"}], ",", "0.135737218047719`", ",", "0.0015322926577817505`", ",", "0.00002767970185703243`", ",", RowBox[{"-", "0.00002361186048020501`"}], ",", "0.012097648995364706`", ",", RowBox[{"-", "0.009449249046788979`"}], ",", "0.04400724560783016`", ",", RowBox[{"-", "0.026897244630998365`"}], ",", "0.08563765049267591`", ",", RowBox[{"-", "0.03463610818159659`"}], ",", "0.1182568023656207`", ",", RowBox[{"-", "0.02511767465169514`"}], ",", "0.13005239304640925`", ",", RowBox[{"-", "0.010486887925560575`"}], ",", "0.13014154645090462`", ",", RowBox[{"-", "0.005394998595914553`"}], ",", "0.13310396487622767`", ",", RowBox[{"-", "0.006035536795374657`"}], ",", "0.1369823187412342`", ",", RowBox[{"-", "0.002496221458261576`"}], ",", "0.13530911586147734`", ",", "0.001481530269444809`", ",", "0.000027559991578729243`", ",", RowBox[{"-", "0.00002348071109014563`"}], ",", "0.01204284240779699`", ",", RowBox[{"-", "0.009397486489227967`"}], ",", "0.0437900802601694`", ",", RowBox[{"-", "0.026762025721071087`"}], ",", "0.08519914323731065`", ",", RowBox[{"-", "0.03452121295704365`"}], ",", "0.11771510051177651`", ",", RowBox[{"-", "0.02517080968250714`"}], ",", "0.1296427646173523`", ",", RowBox[{"-", "0.01064918775106198`"}], ",", "0.1298682567452575`", ",", RowBox[{"-", "0.005397478769787425`"}], ",", "0.1326792090215724`", ",", RowBox[{"-", "0.005839307564718184`"}], ",", "0.13634129839350426`", ",", RowBox[{"-", "0.0023950259277810454`"}], ",", "0.13471316368993075`", ",", "0.001411517265188077`", ",", "0.00002740667204873792`", ",", RowBox[{"-", "0.000023313424920426328`"}], ",", "0.01197269400256546`", ",", RowBox[{"-", "0.00933148832553098`"}], ",", "0.04351253077722427`", ",", RowBox[{"-", "0.02658980331482465`"}], ",", "0.08463963050690225`", ",", RowBox[{"-", "0.03437505043710269`"}], ",", "0.11702457567682822`", ",", RowBox[{"-", "0.025237837181933738`"}], ",", "0.12911924632589308`", ",", RowBox[{"-", "0.010853577686930247`"}], ",", "0.12951552523712087`", ",", RowBox[{"-", "0.005399212518444577`"}], ",", "0.13213407112977518`", ",", RowBox[{"-", "0.005591532353555713`"}], ",", "0.1355245650459254`", ",", RowBox[{"-", "0.0022679766681755397`"}], ",", "0.13395359361425085`", ",", "0.0013240534788368982`", ",", "0.00002722018918008919`", ",", RowBox[{"-", "0.00002311098651475332`"}], ",", "0.011887441122369257`", ",", RowBox[{"-", "0.009251662480008558`"}], ",", "0.04317583338825912`", ",", RowBox[{"-", "0.026381773304797275`"}], ",", "0.08396227476103442`", ",", RowBox[{"-", "0.03419875004664153`"}], ",", "0.11618959161359387`", ",", RowBox[{"-", "0.025317788233711638`"}], ",", "0.1284841289979257`", ",", RowBox[{"-", "0.011096928073697946`"}], ",", "0.12908238381482032`", ",", RowBox[{"-", "0.005399173261209543`"}], ",", "0.13146933089064436`", ",", RowBox[{"-", "0.005295509587104344`"}], ",", "0.1345375618988263`", ",", RowBox[{"-", "0.0021172743651457493`"}], ",", "0.13303525908595001`", ",", "0.0012209653147695597`", ",", "0.000027001088560791487`", ",", RowBox[{"-", "0.000022874581268310592`"}], ",", "0.011787373180899959`", ",", RowBox[{"-", "0.009158499457979318`"}], ",", "0.04278148664966544`", ",", RowBox[{"-", "0.02613936782388599`"}], ",", "0.08317088691158889`", ",", RowBox[{"-", "0.03399364818368514`"}], ",", "0.11521537464167121`", ",", RowBox[{"-", "0.02540947696764624`"}], ",", "0.12774013761960937`", ",", RowBox[{"-", "0.011375502887282118`"}], ",", "0.12856770536811044`", ",", RowBox[{"-", "0.0053961759300114095`"}], ",", "0.13068599137099254`", ",", RowBox[{"-", "0.004955170821152213`"}], ",", "0.13338679612777082`", ",", RowBox[{"-", "0.0019455076111182567`"}], ",", "0.1319639438800789`", ",", "0.0011043899699612761`", ",", "0.000026750015639977614`", ",", RowBox[{"-", "0.000022605584934736977`"}], ",", "0.01167283119531601`", ",", RowBox[{"-", "0.009052567637017717`"}], ",", "0.04233124393876439`", ",", RowBox[{"-", "0.025864238158773936`"}], ",", "0.08226989434337284`", ",", RowBox[{"-", "0.03376126319143601`"}], ",", "0.11410795172207197`", ",", RowBox[{"-", "0.025511500653960576`"}], ",", "0.12689038924901852`", ",", RowBox[{"-", "0.011685011519281481`"}], ",", "0.12797023647279987`", ",", RowBox[{"-", "0.005388915668743429`"}], ",", "0.12978530603911986`", ",", RowBox[{"-", "0.004575022911714009`"}], ",", "0.13207975532282365`", ",", RowBox[{"-", "0.0017555964008383604`"}], ",", "0.13074627884408643`", ",", "0.0009767225343884326`", ",", "0.000026467715792624963`", ",", RowBox[{"-", "0.00002230555146363148`"}], ",", "0.011544207178111514`", ",", RowBox[{"-", "0.008934507821040867`"}], ",", "0.041827104525596535`", ",", RowBox[{"-", "0.025558235111524726`"}], ",", "0.08126430367973178`", ",", RowBox[{"-", "0.03350326699761971`"}], ",", "0.11287407938676626`", ",", RowBox[{"-", "0.02562224074038666`"}], ",", "0.1259383457779331`", ",", RowBox[{"-", "0.012020668444096156`"}], ",", "0.12728863509428845`", ",", RowBox[{"-", "0.005376010872139206`"}], ",", "0.12876880805786636`", ",", RowBox[{"-", "0.004160081547631649`"}], ",", "0.13062481211214694`", ",", RowBox[{"-", "0.0015507289802496`"}], ",", "0.12938964841125108`", ",", "0.0008405570812278774`", ",", "0.000026155034196362074`", ",", RowBox[{"-", "0.000021976199284969163`"}], ",", "0.011401943363299116`", ",", RowBox[{"-", "0.008805027115002747`"}], ",", "0.04127130320419395`", ",", RowBox[{"-", "0.025223387074775644`"}], ",", "0.08015965857820832`", ",", RowBox[{"-", "0.03322145399004032`"}], ",", "0.11152116453327252`", ",", RowBox[{"-", "0.025739865233542117`"}], ",", "0.1248877626847767`", ",", RowBox[{"-", "0.012377260018798754`"}], ",", "0.1265215128331767`", ",", RowBox[{"-", "0.005356049338996321`"}], ",", "0.12763834066136337`", ",", RowBox[{"-", "0.0037157970669176475`"}], ",", "0.12903111855239482`", ",", RowBox[{"-", "0.0013342936710355041`"}], ",", "0.12790208903607003`", ",", "0.0006986233693567935`", ",", "0.000025812915448505798`", ",", RowBox[{"-", "0.000021619396174938495`"}], ",", "0.011246531240284815`", ",", RowBox[{"-", "0.008664892185473276`"}], ",", "0.04066629846886551`", ",", RowBox[{"-", "0.024861876117744177`"}], ",", "0.07896199288801654`", ",", RowBox[{"-", "0.03291770775480815`"}], ",", "0.11005717821473862`", ",", RowBox[{"-", "0.02586233284322552`"}], ",", "0.12374263501432135`", ",", RowBox[{"-", "0.01274921755507034`"}], ",", "0.12566748115322268`", ",", RowBox[{"-", "0.005327636216851564`"}], ",", "0.12639608734897512`", ",", RowBox[{"-", "0.003247973603114108`"}], ",", "0.12730849203373407`", ",", RowBox[{"-", "0.0011098074051714593`"}], ",", "0.12629218183308125`", ",", "0.0005537208972830961`", ",", "0.000025442402846310843`", ",", RowBox[{"-", "0.000021237142851601276`"}], ",", "0.011078510367994421`", ",", RowBox[{"-", "0.008514921979329461`"}], ",", "0.04001475922790768`", ",", RowBox[{"-", "0.024476012408249054`"}], ",", "0.07767777954854503`", ",", RowBox[{"-", "0.03259396634938285`"}], ",", "0.10849056366100049`", ",", RowBox[{"-", "0.025987399308731694`"}], ",", "0.12250714188321853`", ",", RowBox[{"-", "0.013130695701429034`"}], ",", "0.12472520094845249`", ",", RowBox[{"-", "0.005289442352263016`"}], ",", "0.1250446005862211`", ",", RowBox[{"-", "0.002762682724796792`"}], ",", "0.12546729458079303`", ",", RowBox[{"-", "0.0008808427657946789`"}], ",", "0.12456894175310942`", ",", "0.0004086521145190754`"}], "}"}]], "Output", CellChangeTimes->{ 3.601811162685341*^9, 3.601811375229395*^9, 3.601811516137822*^9, { 3.6019014976324778`*^9, 3.60190152563608*^9}, 3.6019197793013983`*^9, 3.601919876150498*^9, 3.601958985667324*^9, 3.60199533430824*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Generate Data Set (takes about 1day for 30 parameter combinations and 21x31 \ points, depending on your machine)\ \>", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122424687347`*^9, 3.6017122635465612`*^9}, {3.60189429865724*^9, 3.601894365413218*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"OAAxAll", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"z0", ",", "x0", ",", RowBox[{"Processing", "[", RowBox[{"AllI1I2", ",", "x0", ",", "z0"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x0", ",", "0", ",", "XMax", ",", "dx"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"z0", ",", RowBox[{"-", "ZMax"}], ",", "0", ",", "dz"}], "}"}]}], "]"}], ",", "2"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.601544793206133*^9, 3.6015448603837223`*^9}, { 3.60154568496122*^9, 3.601545692113284*^9}, {3.601545723097732*^9, 3.6015457286750507`*^9}, 3.601556222986932*^9, {3.601558853384388*^9, 3.601558869231247*^9}, {3.601559070170805*^9, 3.601559074050811*^9}, 3.601567310441114*^9, {3.601567538490617*^9, 3.601567538825157*^9}, { 3.601613714650173*^9, 3.601613714921245*^9}, 3.601784983111292*^9}], Cell[BoxData[ RowBox[{ RowBox[{"OAAxAl3", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"OAAxAll", "[", RowBox[{"[", RowBox[{ RowBox[{"3", "i"}], "+", "1"}], "]"}], "]"}], ",", RowBox[{"OAAxAll", "[", RowBox[{"[", RowBox[{ RowBox[{"3", "i"}], "+", "2"}], "]"}], "]"}], ",", RowBox[{"OAAxAll", "[", RowBox[{"[", RowBox[{ RowBox[{"3", "i"}], "+", "3"}], "]"}], "]"}]}], "}"}], ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{ RowBox[{ RowBox[{"Length", "[", "OAAxAll", "]"}], "/", "3"}], "-", "1"}], ",", "1"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.6016135885876293`*^9, 3.6016136845086737`*^9}, { 3.601613760719355*^9, 3.601613820891623*^9}, {3.601613890261571*^9, 3.601613929246291*^9}, {3.60161408112043*^9, 3.601614081279594*^9}, { 3.601614165113434*^9, 3.601614192624515*^9}, {3.601614249778594*^9, 3.601614253906803*^9}, {3.601614380445856*^9, 3.601614387553987*^9}, { 3.601614468322934*^9, 3.601614475435995*^9}, {3.601614601522099*^9, 3.6016146025115623`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Export Data for Processing", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122424687347`*^9, 3.6017122635465612`*^9}, {3.6017257486970587`*^9, 3.601725754425177*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Export", "[", RowBox[{ "\"\\"", ",", "OAAxAl3", ",", "\"\\""}], "]"}], ";"}]], "Input", CellChangeTimes->{{3.6015457432529373`*^9, 3.601545746033461*^9}, 3.60155907846342*^9, {3.6015670934537287`*^9, 3.6015670957257643`*^9}, 3.6016134385940237`*^9, 3.601613543174852*^9, 3.601613631968341*^9, { 3.6017126747166157`*^9, 3.601712675144404*^9}, {3.601725807942718*^9, 3.601725808222493*^9}, 3.601807627286932*^9, {3.601811550352386*^9, 3.601811550760222*^9}, 3.6019197502920647`*^9, {3.6019589114895773`*^9, 3.601958915177458*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Show Data", "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122424687347`*^9, 3.6017122635465612`*^9}, {3.601725756561282*^9, 3.6017257580732603`*^9}}], Cell[BoxData[ RowBox[{"FigsT", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"ListDensityPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"OAAxAl3", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"OAAxAl3", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"OAAxAl3", "[", RowBox[{"[", "i", "]"}], "]"}], "[", RowBox[{"[", "k", "]"}], "]"}]}], "}"}], ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "OAAxAl3", "]"}]}], "}"}]}], "]"}], ",", RowBox[{"ColorFunction", "\[Rule]", RowBox[{"(", RowBox[{ RowBox[{"Blend", "[", RowBox[{ RowBox[{"{", RowBox[{"Yellow", ",", "Orange", ",", "Red", ",", "Black"}], "}"}], ",", "#"}], "]"}], "&"}], ")"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", RowBox[{"1600", "/", "4000"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "200"}], ",", RowBox[{"LabelStyle", "\[Rule]", "Larger"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "3", ",", RowBox[{"Length", "[", RowBox[{"OAAxAl3", "[", RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.601561996098805*^9, 3.6015622347105637`*^9}, { 3.6015622979893847`*^9, 3.601562298147431*^9}, {3.601562823675922*^9, 3.601562844889489*^9}, {3.6015630414311132`*^9, 3.601563042111162*^9}, { 3.6016414222343903`*^9, 3.601641432738742*^9}, {3.6017864317649317`*^9, 3.6017864569441338`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Simplify Scatter Image for the near-forward direction? (", Cell[BoxData[ FormBox[ RowBox[{ SubsuperscriptBox[ RowBox[{"I1", ")"}], RowBox[{"n", ",", "1"}], "1"], "=", FormBox[ SubsuperscriptBox[ RowBox[{"(", "I2", ")"}], RowBox[{"n", ",", "1"}], "1"], TraditionalForm]}], TraditionalForm]]], ", ", Cell[BoxData[ FormBox[ RowBox[{ FormBox[ SubsuperscriptBox[ RowBox[{"(", "I2", ")"}], RowBox[{"n", ",", "1"}], "0"], TraditionalForm], "=", "1"}], TraditionalForm]]], ", ", " ", Cell[BoxData[ FormBox[ RowBox[{ FormBox[ SubsuperscriptBox[ RowBox[{"(", "I1", ")"}], RowBox[{"n", ",", "1"}], "0"], TraditionalForm], "=", "0"}], TraditionalForm]]] }], "Subsubsection", CellChangeTimes->{{3.601470607014696*^9, 3.60147061531664*^9}, { 3.6017122424687347`*^9, 3.6017122635465612`*^9}, {3.601725756561282*^9, 3.6017257580732603`*^9}, {3.6018076466481037`*^9, 3.601807779820788*^9}}], Cell[BoxData[ RowBox[{"II", "=", RowBox[{"Integrals\[Theta]f", "[", "0.001", "]"}]}]], "Input", CellChangeTimes->{{3.601786905393454*^9, 3.6017869176899443`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "20"}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "=", "1"}], ",", RowBox[{ RowBox[{"|", "m", "|"}], "=", "0"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "20"}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "=", "1"}], ",", RowBox[{ RowBox[{"|", "m", "|"}], "=", "1"}]}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "2", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "20"}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "=", "1"}], ",", RowBox[{ RowBox[{"|", "m", "|"}], "=", "0"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "2", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "20"}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "=", "1"}], ",", RowBox[{ RowBox[{"|", "m", "|"}], "=", "1"}]}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "-", RowBox[{ RowBox[{ RowBox[{ RowBox[{"II", "[", RowBox[{"[", "2", "]"}], "]"}], "[", RowBox[{"[", "n", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "20"}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "=", "1"}], ",", RowBox[{ RowBox[{"|", "m", "|"}], "=", "1"}]}], "*)"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.601786920924073*^9, 3.6017869683098707`*^9}, { 3.601787247814352*^9, 3.601787383041438*^9}, {3.601787950002325*^9, 3.601787952553419*^9}, {3.601788278952433*^9, 3.601788284547599*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "2.4989954491691144`*^-13", ",", "7.500115883975883`*^-13", ",", "1.4999083660817378`*^-12", ",", "2.4998447767301665`*^-12", ",", "3.749721875683818`*^-12", ",", "5.249676647180086`*^-12", ",", "6.999567560955798`*^-12", ",", "8.999408684179703`*^-12", ",", "1.1249314164408837`*^-11", ",", "1.3749115804069499`*^-11", ",", "1.649891943872798`*^-11", ",", "1.9498794403207634`*^-11", ",", "2.2748615654250393`*^-11", ",", "2.624832720644075`*^-11", ",", "2.999809192873197`*^-11", ",", "3.39977746736867`*^-11", ",", "3.824736093340346`*^-11", ",", "4.274696459206707`*^-11", ",", "4.749655191952416`*^-11", ",", "5.249597145230771`*^-11"}], "}"}]], "Output", CellChangeTimes->{{3.601786925435277*^9, 3.6017869694041843`*^9}, { 3.6017872717781*^9, 3.601787277801465*^9}, {3.601787308198646*^9, 3.6017873308955*^9}, 3.60178736176403*^9, 3.6017879120751257`*^9, 3.6017879529250603`*^9, 3.601788285047949*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "9.999796669168078`*^-7", ",", "2.9999375006085114`*^-6", ",", "5.99987050140152`*^-6", ",", "9.999774169418359`*^-6", ",", "0.000014999642504921272`", ",", "0.000020999468008252124`", ",", "0.00002799924167984177`", ",", "0.000035998953020221336`", ",", "0.000044998590030035306`", ",", "0.000054998139210056505`", ",", "0.00006599758556120308`", ",", "0.00007799691258455704`", ",", "0.00009099610228138512`", ",", "0.00010499513515316114`", ",", "0.00011999399020159032`", ",", "0.00013599264492863573`", ",", "0.00015299107533654627`", ",", "0.00017098925592788676`", ",", "0.00018998715970556946`", ",", "0.00020998475817288849`"}], "}"}]], "Output", CellChangeTimes->{{3.601786925435277*^9, 3.6017869694041843`*^9}, { 3.6017872717781*^9, 3.601787277801465*^9}, {3.601787308198646*^9, 3.6017873308955*^9}, 3.60178736176403*^9, 3.6017879120751257`*^9, 3.6017879529250603`*^9, 3.60178828505014*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.9999990000203329`", ",", "0.9999980000412909`", ",", "0.9999965000736655`", ",", "0.9999945001185815`", ",", "0.999992000177539`", ",", "0.9999890002524129`", ",", "0.9999855003454527`", ",", "0.9999815004592842`", ",", "0.9999770005969064`", ",", "0.9999720007616948`", ",", "0.9999665009573989`", ",", "0.9999605011881433`", ",", "0.9999540014584273`", ",", "0.9999470017731261`", ",", "0.9999395021374892`", ",", "0.9999315025571408`", ",", "0.9999230030380813`", ",", "0.9999140035866849`", ",", "0.9999045042097002`", ",", "0.9998945049142525`"}], "}"}]], "Output", CellChangeTimes->{{3.601786925435277*^9, 3.6017869694041843`*^9}, { 3.6017872717781*^9, 3.601787277801465*^9}, {3.601787308198646*^9, 3.6017873308955*^9}, 3.60178736176403*^9, 3.6017879120751257`*^9, 3.6017879529250603`*^9, 3.601788285051111*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "9.999796668491707`*^-7", ",", "2.9999375006081366`*^-6", ",", "5.999870501400019`*^-6", ",", "9.999774169414193`*^-6", ",", "0.000014999642504911899`", ",", "0.000020999468008233747`", ",", "0.000027999241679809108`", ",", "0.00003599895302016734`", ",", "0.000044998590029950935`", ",", "0.000054998139209930487`", ",", "0.00006599758556102159`", ",", "0.00007799691258430355`", ",", "0.00009099610228104009`", ",", "0.00010499513515270178`", ",", "0.00011999399020099038`", ",", "0.00013599264492786514`", ",", "0.00015299107533557098`", ",", "0.0001709892559266685`", ",", "0.00018998715970406542`", ",", "0.00020998475817105111`"}], "}"}]], "Output", CellChangeTimes->{{3.601786925435277*^9, 3.6017869694041843`*^9}, { 3.6017872717781*^9, 3.601787277801465*^9}, {3.601787308198646*^9, 3.6017873308955*^9}, 3.60178736176403*^9, 3.6017879120751257`*^9, 3.6017879529250603`*^9, 3.60178828505201*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ "6.763706314591358`*^-17", ",", "3.748120791600279`*^-19", ",", "1.5009423825346202`*^-18", ",", "4.165708034596649`*^-18", ",", "9.373266594316088`*^-18", ",", "1.83772268236293`*^-17", ",", "3.266159044612582`*^-17", ",", "5.399326818977812`*^-17", ",", "8.437125781010635`*^-17", ",", "1.2601817376070579`*^-16", ",", "1.8149544367407344`*^-16", ",", "2.5348646792711094`*^-16", ",", "3.450337888663557`*^-16", ",", "4.593629079549522`*^-16", ",", "5.999432721448539`*^-16", ",", "7.705966940940723`*^-16", ",", "9.752940642593355`*^-16", ",", "1.218263771113337`*^-15", ",", "1.5040323587262039`*^-15", ",", "1.8373703166568722`*^-15"}], "}"}]], "Output", CellChangeTimes->{{3.601786925435277*^9, 3.6017869694041843`*^9}, { 3.6017872717781*^9, 3.601787277801465*^9}, {3.601787308198646*^9, 3.6017873308955*^9}, 3.60178736176403*^9, 3.6017879120751257`*^9, 3.6017879529250603`*^9, 3.60178828505308*^9}] }, Open ]] }, Closed]] }, Open ]] }, WindowSize->{1345, 837}, WindowMargins->{{Automatic, 64}, {68, Automatic}}, FrontEndVersion->"9.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 20, \ 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 99, 1, 51, "Subtitle"], Cell[CellGroupData[{ Cell[703, 27, 4179, 123, 216, "Input"], Cell[4885, 152, 1106, 25, 28, "Output"] }, Open ]], Cell[6006, 180, 12956, 223, 165, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[18987, 407, 805, 20, 40, "Subsubsection"], Cell[19795, 429, 23917, 537, 778, "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[43749, 971, 1970, 66, 87, "Subsubsection"], Cell[45722, 1039, 6811, 204, 220, "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[52570, 1248, 186, 4, 35, "Subsubsection"], Cell[52759, 1254, 12732, 304, 386, "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[65528, 1563, 232, 5, 35, "Subsubsection"], Cell[CellGroupData[{ Cell[65785, 1572, 1600, 44, 97, "Input"], Cell[67388, 1618, 758, 14, 28, "Output"], Cell[68149, 1634, 747, 14, 28, "Output"], Cell[68899, 1650, 462, 8, 28, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[69410, 1664, 180, 4, 35, "Subsubsection"], Cell[CellGroupData[{ Cell[69615, 1672, 798, 22, 28, "Input"], Cell[70416, 1696, 5104, 74, 165, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[75569, 1776, 246, 3, 35, "Subsubsection"], Cell[75818, 1781, 280, 7, 28, "Input"], Cell[CellGroupData[{ Cell[76123, 1792, 579, 16, 63, "Input"], Cell[76705, 1810, 272, 4, 28, "Output"], Cell[76980, 1816, 8808, 129, 318, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[85837, 1951, 311, 6, 35, "Subsubsection"], Cell[86151, 1959, 1069, 25, 28, "Input"], Cell[87223, 1986, 1306, 34, 28, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[88566, 2025, 220, 3, 35, "Subsubsection"], Cell[88789, 2030, 680, 12, 28, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[89506, 2047, 203, 3, 35, "Subsubsection"], Cell[89712, 2052, 1984, 50, 63, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[91733, 2107, 1015, 35, 40, "Subsubsection"], Cell[92751, 2144, 167, 3, 28, "Input"], Cell[CellGroupData[{ Cell[92943, 2151, 3300, 110, 153, "Input"], Cell[96246, 2263, 1000, 17, 57, "Output"], Cell[97249, 2282, 994, 17, 56, "Output"], Cell[98246, 2301, 910, 16, 46, "Output"], Cell[99159, 2319, 994, 17, 56, "Output"], Cell[100156, 2338, 998, 17, 57, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (* End of internal cache information *)