(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.3' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 41828, 1180] NotebookOptionsPosition[ 39676, 1116] NotebookOutlinePosition[ 40033, 1132] CellTagsIndexPosition[ 39990, 1129] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ "This", " ", "notebook", " ", "contains", " ", "the", " ", "modal", " ", "solution", " ", "to", " ", RowBox[{"Lilley", "'"}], "s", " ", "equation", " ", "for", " ", "a", " ", "point", " ", "source", " ", RowBox[{ RowBox[{"(", "1", ")"}], ".", " ", "Also"}], " ", "included", " ", "here", " ", "are", " ", "the", " ", "symmetric", " ", "ring", " ", "source", " ", RowBox[{"(", "2", ")"}], "and", " ", "weighted", " ", "ring", " ", "source", " ", RowBox[{ RowBox[{"(", "3", ")"}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", "This"}], " ", "work", " ", "is", " ", "used", " ", "in", " ", "the", " ", "paper", " ", "\"\\"", "\[IndentingNewLine]", "\[IndentingNewLine]", "In", " ", "all", " ", "cases", " ", "the", " ", "single", " ", "simple", " ", "turning"}], "-", RowBox[{"point", " ", "solution", " ", "is", " ", RowBox[{"used", ".", " ", "Extensions"}], " ", "can", " ", "be", " ", "made", " ", "using", " ", "the", " ", "methods", " ", "described", " ", "in", " ", "the", " ", "above", " ", RowBox[{"paper", "."}]}]}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", "*)"}]], "Input", CellChangeTimes->{{3.6895085630628357`*^9, 3.689508727797757*^9}, { 3.689509425490583*^9, 3.689509436136877*^9}, {3.689578663522025*^9, 3.689578667182661*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"-", " ", "1."}]}]}]}]}]}]}]}]}]}]}]}]}], " ", "Point", " ", "source", " ", RowBox[{ RowBox[{"calculation", " ", "--"}], "--"}]}], "-", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", "-"}]}]}]}]}]}]}]}]}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689509361666357*^9, 3.6895094052457457`*^9}, { 3.689509439747902*^9, 3.689509446317203*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "The", " ", "modal", " ", "solution", " ", "is", " ", "calculated", " ", "quickly", " ", "by", " ", "noting", " ", "the", " ", "separable", " ", "nature", " ", "of", " ", "the", " ", RowBox[{"sum", ".", " ", "Compute"}], " ", "the", " ", "azimuthal", " ", "components", " ", "first", " ", RowBox[{"using", "."}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689508979926484*^9, 3.689508996809946*^9}, { 3.689509073286088*^9, 3.6895090961145163`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"ExpFunc", "[", RowBox[{"\[CapitalDelta]\[Phi]_", ",", "N_"}], "]"}], ":=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Exp", "[", RowBox[{"\[ImaginaryI]", "*", "n", "*", "\[CapitalDelta]\[Phi]"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", RowBox[{"-", "N"}], ",", "N"}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Compute", " ", "the", " ", "radial", " ", "solutions"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689509112455064*^9, 3.689509118484796*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"myArgm", "[", "x_", "]"}], ":=", RowBox[{"Mod", "[", RowBox[{ RowBox[{"Arg", "[", "x", "]"}], ",", RowBox[{"2", " ", "\[Pi]"}], ",", RowBox[{"-", FractionBox["\[Pi]", "2"]}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"q", "[", RowBox[{ RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"r_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}]}], "]"}], ":=", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"1", "-", RowBox[{"MJ", " ", SuperscriptBox[ RowBox[{"Sech", "[", RowBox[{"2", " ", "r"}], "]"}], "2"], " ", RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}], "2"], "-", SuperscriptBox[ RowBox[{"Cos", "[", "\[Theta]", "]"}], "2"]}]]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Q", "[", RowBox[{ RowBox[{"\[Nu]_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"r_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}]}], "]"}], ":=", FractionBox[ RowBox[{ SqrtBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ SuperscriptBox["r", "2"], " ", SuperscriptBox[ RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "r", ",", "MJ"}], "]"}], "2"]}], "-", SuperscriptBox["\[Nu]", "2"]}], "]"}]], " ", RowBox[{"Exp", "[", RowBox[{ FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ RowBox[{"myArgm", "[", RowBox[{ RowBox[{"r", " ", RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "r", ",", "MJ"}], "]"}]}], "-", "\[Nu]"}], "]"}], "+", RowBox[{"myArgm", "[", RowBox[{ RowBox[{"r", " ", RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "r", ",", "MJ"}], "]"}]}], "+", "\[Nu]"}], "]"}]}], ")"}]}], "]"}]}], "r"]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"r\[Delta]", "[", RowBox[{ RowBox[{"\[Nu]_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"M_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"p", "/.", "\[VeryThinSpace]", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{ SuperscriptBox["p", "2"], " ", SuperscriptBox[ RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "p", ",", "M"}], "]"}], "2"]}], "-", SuperscriptBox["\[Nu]", "2"]}], "]"}], " ", RowBox[{"Exp", "[", RowBox[{"\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ RowBox[{"myArgm", "[", RowBox[{ RowBox[{"p", " ", RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "p", ",", "M"}], "]"}]}], "-", "\[Nu]"}], "]"}], "+", RowBox[{"myArgm", "[", RowBox[{ RowBox[{"p", " ", RowBox[{"q", "[", RowBox[{"\[Theta]", ",", "p", ",", "M"}], "]"}]}], "+", "\[Nu]"}], "]"}]}], ")"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"p", ",", "1"}], "}"}]}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"r\[Delta]3", "[", RowBox[{"\[Nu]_", ",", "\[Theta]_", ",", "M_"}], "]"}], ":=", RowBox[{"Abs", "[", RowBox[{"r\[Delta]", "[", RowBox[{"\[Nu]", ",", "\[Theta]", ",", "M"}], "]"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"DiffP", "[", RowBox[{ RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "x"], RowBox[{"(", RowBox[{ SuperscriptBox["x", "2"], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"1", "-", RowBox[{"MJ", " ", SuperscriptBox[ RowBox[{"Sech", "[", RowBox[{"2", " ", "x"}], "]"}], "2"], " ", RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}], "2"], "-", SuperscriptBox[ RowBox[{"Cos", "[", "\[Theta]", "]"}], "2"]}], ")"}]}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"x", "\[Rule]", "rs"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"P", "[", RowBox[{ RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}]}], "]"}], ":=", FractionBox["1", RowBox[{"rs", " ", RowBox[{"DiffP", "[", RowBox[{"\[Theta]", ",", "rs", ",", "MJ"}], "]"}]}]]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"ZetaN", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{"Q", "[", RowBox[{"n", ",", "\[Theta]", ",", "r", ",", "MJ"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"r", ",", RowBox[{"r\[Delta]3", "[", RowBox[{"n", ",", "\[Theta]", ",", "MJ"}], "]"}], ",", "rs"}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"EtaN", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"sign", "=", RowBox[{"Sign", "[", RowBox[{"rs", "-", RowBox[{"r\[Delta]", "[", RowBox[{"n", ",", "\[Theta]", ",", "MJ"}], "]"}]}], "]"}]}], "}"}], ",", RowBox[{ RowBox[{"-", "sign"}], " ", SuperscriptBox[ RowBox[{"(", RowBox[{ FractionBox["3", "2"], " ", "sign", " ", "k", " ", RowBox[{"ZetaN", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ"}], "]"}]}], ")"}], RowBox[{"2", "/", "3"}]]}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"AmplitudeFunc", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"sign", "=", RowBox[{"Sign", "[", RowBox[{"rs", "-", RowBox[{"r\[Delta]", "[", RowBox[{"n", ",", "\[Theta]", ",", "MJ"}], "]"}]}], "]"}]}], "}"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"sign", "\[Equal]", "0"}], ",", RowBox[{ SqrtBox[ FractionBox["2", "k"]], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"k", " ", RowBox[{"P", "[", RowBox[{"\[Theta]", ",", "rs", ",", "MJ"}], "]"}]}], ")"}], RowBox[{"1", "/", "3"}]]}], ",", RowBox[{ SqrtBox[ FractionBox["2", "k"]], " ", SqrtBox[ FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", "sign"}], " ", RowBox[{"EtaN", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}]}]], RowBox[{"Abs", "[", RowBox[{"rs", " ", RowBox[{"Q", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ"}], "]"}]}], "]"}]]]}]}], "]"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"AiMode", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"AiryAi", "[", RowBox[{"EtaN", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"ReverseFunc", "[", RowBox[{"array1_", ",", " ", "array0_"}], "]"}], ":=", RowBox[{"(*", " ", RowBox[{ "array1", " ", "is", " ", "a", " ", "sequence", " ", "of", " ", "terms", " ", RowBox[{"a", "[", "1", "]"}], " ", "to", " ", RowBox[{"a", "[", "N", "]"}]}], "*)"}], RowBox[{"Flatten", "[", RowBox[{"{", RowBox[{ RowBox[{"Reverse", "[", "array1", "]"}], ",", "array0", ",", "array1"}], "}"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"NTermQuick", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{ RowBox[{"AmplitudeFunc", "[", RowBox[{ FractionBox["n", "k"], ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}], " ", RowBox[{"AiMode", "[", RowBox[{ FractionBox["n", "k"], ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}], " ", RowBox[{"Exp", "[", RowBox[{"\[ImaginaryI]", " ", "k", " ", RowBox[{"(", RowBox[{ RowBox[{"ZetaN", "[", RowBox[{ FractionBox["n", "k"], ",", "\[Theta]", ",", RowBox[{"100", " ", RowBox[{"Sin", "[", "\[Theta]", "]"}]}], ",", FractionBox["9", "10"]}], "]"}], "-", RowBox[{"100", " ", SuperscriptBox[ RowBox[{"Sin", "[", "\[Theta]", "]"}], "2"]}]}], ")"}]}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"NTermQuickArray", "[", RowBox[{ RowBox[{"N_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"Table", "[", RowBox[{ RowBox[{"NTermQuick", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", "N"}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"NTermsForAPoint", "[", RowBox[{ RowBox[{"N_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"ReverseFunc", "[", RowBox[{ RowBox[{"NTermQuickArray", "[", RowBox[{"N", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}], ",", " ", RowBox[{"NTermQuick", "[", RowBox[{"0", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", " ", "Example"}]}]}]}]}]}], " ", "computations", " ", "of", " ", "point", " ", RowBox[{ RowBox[{"source", " ", "--"}], "--"}]}], "-", RowBox[{"--", RowBox[{"--", RowBox[{"--", "-"}]}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.68950934022993*^9, 3.6895093566414623`*^9}, { 3.6895094149302*^9, 3.689509419726527*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"N", ".", "B", ".", " ", "Use"}], " ", "radians"}], ",", " ", RowBox[{"not", " ", RowBox[{"degrees", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689509621915798*^9, 3.689509630852298*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"k0", "=", "28.27"}], ";", " ", RowBox[{"rs", " ", "=", " ", "0.75"}], ";", " ", RowBox[{"Ns", " ", "=", " ", "30"}], ";", " ", RowBox[{"MJ", " ", "=", " ", "0.9"}], ";", " ", RowBox[{"\[CapitalDelta]\[Phi]", " ", "=", " ", "0"}], ";", RowBox[{"cInf", " ", "=", " ", "340"}], ";"}]], "Input", CellChangeTimes->{{3.689575992911243*^9, 3.689576026088739*^9}, { 3.689576101709401*^9, 3.6895761056568947`*^9}, {3.6895763973288803`*^9, 3.689576403664977*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"radialSol", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"(", " ", FractionBox["\[ImaginaryI]", RowBox[{"cInf", "*", "k0", "*", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "-", RowBox[{"MJ", "*", SuperscriptBox[ RowBox[{"Sech", "[", RowBox[{"2", "*", "rs"}], "]"}], "2"], "*", RowBox[{"Cos", "[", RowBox[{"i", "*", RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}]}], ")"}], "2"]}]], ")"}], "*", RowBox[{"NTermsForAPoint", "[", RowBox[{"Ns", ",", RowBox[{"i", "*", RowBox[{"\[Pi]", "/", "180"}]}], ",", "rs", ",", "MJ", ",", "k0"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", " ", "2", ",", "100", ",", " ", "1"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.689509481573625*^9, 3.689509482695612*^9}, { 3.689509545632063*^9, 3.689509620044639*^9}, {3.689509688920383*^9, 3.689509694559883*^9}, {3.6895759647125187`*^9, 3.6895759719647083`*^9}, { 3.689576031754681*^9, 3.689576042401538*^9}, {3.6895761619034452`*^9, 3.689576162059105*^9}, {3.689576355093238*^9, 3.689576409320458*^9}, 3.689576617582552*^9}], Cell[BoxData[ RowBox[{ RowBox[{"aziAngles", " ", "=", " ", RowBox[{"ExpFunc", "[", RowBox[{"\[CapitalDelta]\[Phi]", ",", "30"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.6895760864207067`*^9, 3.6895761116440773`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"totSol\[CapitalDelta]\[Phi]", " ", "=", " ", RowBox[{"radialSol", ".", "aziAngles"}]}], " ", ";"}]], "Input", CellChangeTimes->{{3.689576082706661*^9, 3.6895760840699387`*^9}, { 3.689576123468993*^9, 3.689576145912661*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListLinePlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"i", "*", RowBox[{"\[Pi]", "/", "180"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "100", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"Abs", "[", "totSol\[CapitalDelta]\[Phi]", "]"}]}], "}"}], "]"}], ",", " ", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{ "\"\<\[Theta]\>\"", ",", "\"\\""}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.6895761479398518`*^9, 3.689576155494587*^9}, { 3.689576424527082*^9, 3.6895764501661243`*^9}, {3.689577245142013*^9, 3.689577279993988*^9}}], Cell[BoxData[ GraphicsBox[{{}, {{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.016666666666666666`], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJwBMQPO/CFib1JlAgAAADIAAAACAAAAOZ1SokbfkT/25BQKA5VmPtbre/Pp zqo/PqlJ2+zieT6HROdKGFe2P12e3YzHIY0+JJMQnLtGvz9WcX+cgXCfPuDw nHYvG8Q/nwVcg3lRsD4umDEfAZPIP2FkESflVcA+fT/Gx9IKzT+HTJw1e5LP PmZzLThSwdA/DwV3y4Bz3T4Nx3cMO/3SP5JFe2MEf+o+tBrC4CM51T/BrBYp lPT2PltuDLUMddc//5dcvlEcAz8CwlaJ9bDZP12OrQFkfQ4/qRWhXd7s2z9k 3LGEfzgXP1Bp6zHHKN4/ZbltgrPNID983hoDWDLgP60MvkpR9iY/TwhAbUxQ 4T/eyxSRl14tPyMyZddAbuI/f/4mA89bMT/2W4pBNYzjP0WZfvVClzI/yoWv qymq5D+/wbQQV2wxP56v1BUeyOU/CIRW5DydKj9x2fl/EubmP6DCVcQYSRw/ RQMf6gYE6D/kWdiDSnQXPxgtRFT7Iek/tzb8/EbIJj/sVmm+7z/qP5TM1wyM 1yo/v4COKORd6z8T6FToLoAjP5Oqs5LYe+w/fOyE5La4GD9n1Nj8zJntP7Mq kHOH1yQ/Ov79ZsG37j+cnV53PXAjPw4oI9G11e8/JgbiPtAeHz/xKKQd1Xnw P6xBwDRhPCE/2r22Us8I8T+fjTIOlrwgP8RSyYfJl/E/FEL+mit5ID+u59u8 wybyP9O2n8ocOyA/mHzu8b218j/3dv6P//sfP4IRASe4RPM/Og1M47mEHz9r phNcstPzPzjRzYdPEB8/VTsmkaxi9D8/WM3Zn54ePz/QOMam8fQ/0z97zL4v Hj8pZUv7oID1P2d3pdB0wx0/EvpdMJsP9j8NPG3KKlkdP/yOcGWVnvY/mnLd cyDxHD/mI4Oajy33Px2NI9MZixw/0LiVz4m89z9Q2y9cTCccP7pNqASES/g/ 1wIlrmPFGz+j4ro5ftr4P3KX0wBQZRs/jXfNbnhp+T92o9zmTAcbP3cM4KNy +Pk/Is6VqdeqGj9hofLYbIf6P9nrlXstUBo/SzYFDmcW+z8HM8ZENPcZPzTL F0Nhpfs/KYwrk9+fGT8bOn4N "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"\[Theta]\"", TraditionalForm], FormBox[ "\"Abs[\\!\\(\\*SubscriptBox[\\(G\\), \ \\(\[Omega]\\)]\\)/\\!\\(\\*SubscriptBox[\\(g\\), \\(\[Omega]\\)]\\)]\"", TraditionalForm]}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{}, PlotRange->{{0, 1.7278759594743862`}, {0, 0.00028367409458453604`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.689576156242044*^9, 3.689576611932705*^9, {3.689577267823657*^9, 3.689577280357973*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"-", " ", "2."}]}]}]}]}]}]}]}]}]}]}]}]}], " ", "Ring", " ", "source", " ", RowBox[{ RowBox[{"calculation", " ", "--"}], "--"}]}], "-", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", "-"}]}]}]}]}]}]}]}]}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689509361666357*^9, 3.6895094052457457`*^9}, { 3.689509439747902*^9, 3.689509446317203*^9}, {3.689576332868601*^9, 3.689576334958552*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "The", " ", "ring", " ", "source", " ", "is", " ", "calculated", " ", "using", " ", "the", " ", "point", " ", "source", " ", "solution", " ", "of", " ", "the", " ", "previous", " ", RowBox[{"section", "."}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689576455426899*^9, 3.6895764998857737`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"AmplitudeFuncRS", "[", RowBox[{ RowBox[{"n_", "?", "NumericQ"}], ",", RowBox[{"\[Theta]_", "?", "NumericQ"}], ",", RowBox[{"rs_", "?", "NumericQ"}], ",", RowBox[{"MJ_", "?", "NumericQ"}], ",", RowBox[{"k_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"sign", "=", RowBox[{"Sign", "[", RowBox[{"rs", "-", RowBox[{"r\[Delta]", "[", RowBox[{"n", ",", "\[Theta]", ",", "MJ"}], "]"}]}], "]"}]}], "}"}], ",", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", "sign"}], " ", RowBox[{"EtaN", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k"}], "]"}]}]], RowBox[{"Abs", "[", RowBox[{"rs", " ", RowBox[{"Q", "[", RowBox[{"n", ",", "\[Theta]", ",", "rs", ",", "MJ"}], "]"}]}], "]"}]]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"D\[Omega]", "[", RowBox[{"\[Theta]_", ",", " ", "rs_", ",", " ", "MJ_", ",", " ", "k0_"}], "]"}], ":=", RowBox[{ FractionBox[ RowBox[{"4", "*", "\[Pi]"}], RowBox[{ SuperscriptBox["cInf", "2"], "*", SuperscriptBox["k0", "3"], "*", SuperscriptBox[ RowBox[{"(", RowBox[{"1", " ", "-", " ", RowBox[{"MJ", "*", SuperscriptBox[ RowBox[{"Sech", "[", RowBox[{"2", "*", "rs"}], "]"}], "2"], "*", RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ")"}], "4"]}]], "*", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"AmplitudeFuncRS", "[", RowBox[{ RowBox[{"n", "/", "k0"}], ",", "\[Theta]", ",", "rs", ",", "MJ", ",", "k0"}], "]"}], "*", SuperscriptBox[ RowBox[{"AiMode", "[", RowBox[{ RowBox[{"n", "/", "k0"}], ",", " ", "\[Theta]", ",", "rs", ",", "MJ", ",", "k0"}], "]"}], "2"]}], ",", RowBox[{"{", RowBox[{"n", ",", RowBox[{"-", "30"}], ",", " ", "30"}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.6895766442776546`*^9, 3.68957665898383*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"ringSource", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"D\[Omega]", "[", RowBox[{"#", ",", " ", "rs", ",", " ", "MJ", ",", " ", "k0"}], "]"}], "&"}], "/@", RowBox[{"Table", "[", RowBox[{ RowBox[{"i", "*", RowBox[{"\[Pi]", "/", "180"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", " ", "158", ",", "0.5"}], "}"}]}], "]"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.689576533934832*^9, 3.689576538461903*^9}, { 3.689576663957909*^9, 3.6895766702275267`*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.689508979926484*^9, 3.689508996809946*^9}, { 3.689509073286088*^9, 3.6895090961145163`*^9}, 3.689576339603603*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"-", " ", "3."}]}]}]}]}]}]}]}]}]}], " ", "Weighted", " ", "ring", " ", "source", " ", RowBox[{ RowBox[{"calculation", " ", "--"}], "--"}]}], "-", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", "-"}]}]}]}]}]}]}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689509361666357*^9, 3.6895094052457457`*^9}, { 3.689509439747902*^9, 3.689509446317203*^9}, {3.689576332868601*^9, 3.689576334958552*^9}, {3.689576847019689*^9, 3.689576860796979*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "We", " ", "need", " ", "the", " ", "radial", " ", "solutions", " ", "from", " ", "the", " ", "point", " ", "source", " ", RowBox[{"solution", ".", " ", "I", ".", "e", "."}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"needs", " ", "'"}], RowBox[{"radialSol", "'"}], " ", "from", " ", "section", " ", "1.", " ", "This", " ", "code", " ", "generates", " ", "one", " ", "cross"}], "-", RowBox[{"section", " ", "at", " ", "\[Phi]"}]}], "=", RowBox[{ RowBox[{"const", ".", " ", "of"}], " ", "the", " ", "weighted", " ", "ring", " ", "source", " ", RowBox[{"calculation", "."}]}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689576912862864*^9, 3.689576932750504*^9}, { 3.689577123795817*^9, 3.689577137942379*^9}, {3.689577169137891*^9, 3.689577191779233*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Example", " ", "parameters"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.689577042130295*^9, 3.6895770467720003`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Alpha]T", " ", "=", " ", "1"}], ";", " ", RowBox[{"\[Phi]", " ", "=", " ", RowBox[{ RowBox[{"\[Pi]", "/", "2"}], "+", RowBox[{"28", "*", RowBox[{"\[Pi]", "/", "180"}]}]}]}], ";"}]], "Input", CellChangeTimes->{{3.6895770241422253`*^9, 3.689577040082822*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"YNewFunc", "[", " ", RowBox[{"\[Alpha]T_", ",", "\[Phi]_"}], " ", "]"}], ":=", RowBox[{"(*", " ", RowBox[{"\[Phi]P", " ", "is", " ", "at", " ", RowBox[{"\[Pi]", "/", "2"}]}], " ", "*)"}], RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ FractionBox[ SqrtBox["\[Pi]"], RowBox[{"2", " ", SqrtBox["\[Alpha]T"]}]], "\[ImaginaryI]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"m", "-", "n"}], ")"}], " ", RowBox[{"(", RowBox[{"m", "-", "n", "-", RowBox[{"4", " ", "\[ImaginaryI]", " ", "*", "\[Alpha]T", "*", RowBox[{"(", RowBox[{"\[Phi]", " ", "-", " ", RowBox[{"\[Pi]", "/", "2"}]}], ")"}]}]}], ")"}]}], RowBox[{"4", " ", "*", "\[Alpha]T"}]]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"Erfi", "[", FractionBox[ RowBox[{ RowBox[{"-", "m"}], "+", "n", "+", RowBox[{ "2", " ", "\[ImaginaryI]", " ", "*", "\[Alpha]T", "*", " ", RowBox[{"(", RowBox[{"\[Pi]", "-", RowBox[{"\[Pi]", "/", "2"}]}], ")"}]}]}], RowBox[{"2", " ", SqrtBox["\[Alpha]T"]}]], "]"}]}], "+", RowBox[{"Erfi", "[", FractionBox[ RowBox[{ RowBox[{"-", "m"}], "+", "n", "-", RowBox[{"2", " ", "\[ImaginaryI]", " ", "*", "\[Alpha]T", "*", RowBox[{"(", RowBox[{"\[Pi]", "+", RowBox[{"\[Pi]", "/", "2"}]}], ")"}]}]}], RowBox[{"2", " ", SqrtBox["\[Alpha]T"]}]], "]"}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"n", ",", RowBox[{"-", "40"}], ",", "40"}], "}"}], ",", RowBox[{"{", RowBox[{"m", ",", RowBox[{"-", "40"}], ",", "40"}], "}"}]}], "]"}], "//", "N"}]}], ";"}]], "Input", CellChangeTimes->{{3.674196092866037*^9, 3.674196172849228*^9}, 3.674196217165133*^9, {3.6741962662609*^9, 3.6741962746192427`*^9}, { 3.674196463909836*^9, 3.674196464138926*^9}, {3.674196516208467*^9, 3.674196525611731*^9}, {3.674198430572459*^9, 3.6741984441590557`*^9}, { 3.6741986534733677`*^9, 3.674198676966761*^9}, {3.674213322494532*^9, 3.674213327196993*^9}, {3.6742135555262003`*^9, 3.6742135619146357`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"siftFunc", "[", RowBox[{"a_", ",", "n_"}], "]"}], ":=", RowBox[{ RowBox[{"a", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"D1", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", " ", RowBox[{ "radialSol", "\[LeftDoubleBracket]", "#", "\[RightDoubleBracket]"}], ",", RowBox[{"Conjugate", "[", RowBox[{ "radialSol", "\[LeftDoubleBracket]", "#", "\[RightDoubleBracket]"}], "]"}]}], "]"}], "&"}], "/@", RowBox[{"Table", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", RowBox[{"Length", "[", "radialSol", "]"}]}], "}"}]}], "]"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.6209964113597126`*^9, 3.620996416897266*^9}, { 3.6209964505066266`*^9, 3.6209964506666427`*^9}, {3.6209964837549305`*^9, 3.620996530302931*^9}, {3.620996660816931*^9, 3.6209966610089307`*^9}, { 3.620997204962243*^9, 3.620997214614208*^9}, {3.620997859482929*^9, 3.620997859649929*^9}, {3.6209980746459293`*^9, 3.620998075045929*^9}, { 3.6209985345479293`*^9, 3.620998544687929*^9}, {3.620998770701929*^9, 3.620998802880929*^9}, {3.6895769909367743`*^9, 3.689577001187439*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"D1", " ", "=", " ", RowBox[{"Flatten", "[", RowBox[{"D1", ",", "1"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.620996873083931*^9, 3.6209969094579306`*^9}, 3.6209972020809546`*^9, 3.620998381259929*^9}], Cell[BoxData[ RowBox[{ RowBox[{"Yarrays", " ", "=", " ", RowBox[{"YNewFunc", "[", RowBox[{"\[Alpha]T", ",", "\[Phi]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.689577050794903*^9, 3.689577060509108*^9}}], Cell[BoxData[ RowBox[{"d\[Phi]Example", " ", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"D1", "\[LeftDoubleBracket]", "#", "\[RightDoubleBracket]"}], "*", "Yarrays"}], "&"}], "/@", RowBox[{"Table", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "radialSol", "]"}]}], "}"}]}], "]"}]}], "}"}], ";;"}]}]], "Input", CellChangeTimes->{{3.689577103529965*^9, 3.6895771209854727`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"d\[Phi]Example2", "=", RowBox[{"Flatten", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"Total", "[", RowBox[{"Flatten", "[", RowBox[{"siftFunc", "[", RowBox[{"d\[Phi]Example", ",", "#"}], "]"}], "]"}], "]"}], "}"}], "&"}], "/@", RowBox[{"Table", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "radialSol", "]"}]}], "}"}]}], "]"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.6895771500908623`*^9, 3.689577165102434*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListLinePlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"i", "*", RowBox[{"\[Pi]", "/", "180"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "100", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"Abs", "[", "d\[Phi]Example2", "]"}]}], "}"}], "]"}], " ", ",", " ", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{ "\"\<\[Theta]\>\"", ",", "\"\\""}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.674195567639019*^9, 3.674195574945075*^9}, { 3.6741989151333303`*^9, 3.674198920293249*^9}, {3.674200582740075*^9, 3.674200583225944*^9}, {3.689577222579186*^9, 3.689577224417762*^9}, { 3.6895772896881847`*^9, 3.68957735007836*^9}}], Cell[BoxData[ GraphicsBox[{{}, {{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011111111111111112`], AbsoluteThickness[1.6], LineBox[CompressedData[" 1:eJwtyg0sFHAYx3E08zLm5U5i1rzccnTJnULn3P3uzrvDleqIsZSXRV7HbApF qtEWs6W8hchLzcsW85KR3nShvESNbiYptZhqq6VUu+e/PXv2eb5/m5iUQ7Fa Ghoa0n/zf6vfuui0t46UxekQ0AGzazNlhokjZF0Mu4jyR7dPkY1Rf1nYGaSv IjPRn7QacEv7A3kHVo0NSqtTNshWKCwb5sSYbZKtkTvGSX3C1/ZU2xaRyuZg hsqIzEJpTj1v4KE5eRc6mO+WCmpsyGxE3O3zHPzpQHaEDs+mS8/AhcxBt/71 5weVHmQnJPxws55z8iI7I9Q8xqR+TxCZi7rZ2zXntw6TeSjva2prMY8iu2D7 9FvRVHYseR+MX4XFLxQkk/fDlWHO9a7IJLvCPmR+ejQ1l+wGi1+aZ45HXSC7 Q/UySpq5eIV8AJ/SF0c+jpeT+eBIZv5kLFSRPVCj+ModXW8gC9D+vt6pJbWN 7Il4zs6JZ2NdZCFqV/gJPa29ZBFCS5bcLdOHyQCz17J6d80TtfOB/uVMOXd2 groY2vciKzOiZ6iLIW4XKNll89QlYP65c7R72xJ1CZQRbIcV3kfqUhxrZW81 1a5Tl+K+srjj6cx36l4Q+yT5z5VsUvdCcHGmqW+5llDdvdEsKxCcMNFVO98b qr2mz+t0jKj7YMNr3ObLCoO6D6q7B49onLKg7ovXgbZxMpY1dV9M2n1rzIhj UffDUGW484tGB+p+kPcoqhP1nKn7Y2BBAtaGC3V/tKW6qkrb3KkHIFxzlsVW CKkHYGzgEy+7REI9EIyonLWsYh/qgbjqOiS1EgRSl2HZbuXNepicugxZn0sU nXmh1IMQbiAqbLdUUA9CxJm4m4+cI6gHo8VNdq3icTT1YDzgr5muJsdQD8GI 7bazHH4s9RD0NrT+7hInUJdjwc7xxsm8RLUhR7T9OTPmajL9l6Nn07Cxip8m nCwrcuy5KMelIj2uX3qa8C91aCIL "]], LineBox[CompressedData[" 1:eJw11QlQU9ceBvAUqy0WF7QqWtSAoFgtUseq1Sd8UB8JFYSyBEhYEpKwZU9A XKBeBCyoVYQiiigR6oZY0KIipXIRqLbIYlGwVZ8BixQVwZWqrT7e85/MZDK/ Oec73zk3Nzd2URp/uQWHw+kffv/vc1Fpcp3VQCAyN1t+zNfrXDn/fwVBcOry rE292jdGENJm3iqcOFr9xkwQUFjgKLsc/8ZsEHwMGpukoFjKC2DvPOLevapo yguwIKW8dKQVmREgwnWCMMWLzApwcVm+cEVQDOWD8fftDU/T/6H1EIw7R8Jq a15RHxOM60eXnnlqqaJ8MD7M1nSUNtP+OCGoiPdWPXytoXwIPDe27NYsovMw IYgzfSY13qNxNgTTx1lf3bPNnA+Fl8B6yySDkvKh+FX87GB8TxzlQ2GVsSwx yI32y4bitK2r24rbMsoLcapky2/+TlLKC8Fzq927gx9FeSE2XjqUHOVPZoW4 0uLj/NKO5nNEYMc2x4nktB5E8LjWuSelTE55EcoDZ/WOTDdfPxH+5q1+uTaP zAnDgGTD06IGmo8wXCyPXFLjRusxw+OH+uZ1rTX3hyHn8Nhv1o+QUD4ce8s0 S1KEYsqHQ+XoMMJ2O5kJx7/ubjuz6nsyG44g5h//KSfInAhc222fL8s15yOQ VOL4dUSsOR+B0p6Rn1d+as5HYFKi71cfWZnzkeC/2mZRcT6S8pHo4L1WifRk JhKBNVu+c59NZiPR2FnjvexmBOXF6E3r2hS7i8wVoyysNdjjCzLEODZ3/7oN Y8hiMaoXl1j2/xJO64vR4J4jSMsiG8VYLEj6JWkFmRUjaHHTioeWZJMYmTmp a041hVG/BLu+drmRsJ3MlaDn8I3vHPzJkGB78zupZyeTxRI0VXJ2zL8hon4J ON/qZiQUk40SdGuas1LjyKwELXsfJgd8TDZJkNL/3uru50Lqj8Ik7yfrHOrJ 3Cj8p+VpwfSvyYhC1vXyofpgsjgKxXfPjp04i8xEIXehyHPMQCj1R6G34dVc 4w9kNgqDC6blXMokm6LQvfZSZq6AzJHC4ph/yh0HMleKA7szqi4+DqF+KdoT /7yxuJ4slmJq3u3U+blkRoofph+WFUvJRimaWlv99i0is1LM9qsYNXkU2SRF ZAQ/d9y1YOqXwXrXzg/TSslcGV4MHhsvTyFDBknNoGu1H1kswxjl2krGgczI cDNx7cXTfwmof3j8edvBL5rJrAzXk4aOC4vJJhmwta7wShKZI8e109l3z/iQ uXI4uA78ZeFAhhzH53bOq3wRRP1y+BX6XG+7TGbkqJySIfA7SjbK0XQrZ9zc VDIrx/RBuwWSULJJjpeGRuHQx2RONF60H+jpGU3mRqPXMmH6wj8CqT8a7xTs X3XzR7I4GvzjoRld+WRmeL67PHy5nmyMhv9thVOfN5mNhtRqZdL9OWRTNHiR ogLvEWRODK7Zq7tG3gqg/hh8cNdP80c1GTF4Ma3guWc+WRyDR/4X+3sMZCYG e/1OGNv9yMYY/KJfbzvOmczGIGiaTJf/HtkUg85qm1OKPn/qj0WU4m2LrRfI 3FjI6s+seXKQjFhUTwn55HA6WRwLr1szpMVSMhOLuJ975/R4kI2x+Cmn5bja nszGQvTod66rBdkUi8sDucX+3V9QfxzWCWasP3aezI1Dd/KzJl4JGXGwPb/l 59npZHEcnGI+O8iTk5k4HPlWkXPUk2yMw41TMxtXOZHZ4fXfddN/NJpsisOS aZJWn/t+1B+PothLE8payNx4lK9SG3gnyIjH/NR2p5nfkMXxSFnOVSxMIjPx eBD4roERko3x8M5UNP66y5f64+HUOP7sz3Y0boqHi9BaUT/S3K/A5L4n+/ra aT5XgcfJL+ZVF5KhQPURq/pncrJYAZdoh+D9C8iMAtPGuHDLn6+mfgU2O+7P m9NAZhW4HPDIedQOskmBui+LpV5CMkeJfbzHpQOOZK4SnyiubBh66EP9Sljt rPGOPEcWK7H7h4xcx61kRgmPfK9GfgjZqETaNlfBBUcyq4RtQt7Zosfe1K/E TYsQSVsdmaNCsjv/XHA2matCQXid26JIMlSoytwTIHMmi1Xw9h8ZfeefVdSv Quf15K6qZrJRBdHpYuff95FZFTQN8yv5arJJhVFMa4eFG5mjRtfOec9HjSdz 1ZAkn9sQ0PU59ash3dxRc+ckWazGN9ezbGrTyYwaK7Un//xdQDaqccTmQfrS uWRWjT7LjfeuvvSifjUcdO2OJ1vIHA08zl1xbDpA5mrw1Ldo2cxEMjQIKt1w ooxPFmvw8v64Dp0tmdFg0p6yJ8pBPvVrUHMhz7+wgcxq0Fkyadmr3WSTBttf Rz/OVpE5WnxoP+2Cvwd5vBYt7IwX7lPIXC2SN83tCrvPe2MXLcr/nPzTgToy tMg+G/1ybD7ZT4uEjU43ipVksRZNr5+zIg+yVgs79tmkpTZkRgtP19u+ix94 vnG2FokWZY8CG8hGLdKb5+XlFJArtBjBr+3p15JZLTKC28vlPHKbFq7tR7r/ mk42abHzWd3eQ0/+/caDWvSvsZ+lbiJzdGgMXzrVt5g8XoeZO5RRXuvIXB3G yKa4h/uRXXRIa3/15dY5ZOgQ4PPj+quvVtL5daivrvJY0UEW61D1VudQ43Gy Vofa4IIzhgwyo4PlOd8PvMPJ2TqMwposwSdkow41vGvlhWPIFTrkne63m3fn Mzq/DkUuu2ZPqCW36fDWzrL66N1kkw62O4RTV+rJgzoknn3wfrM3maPHk+8H ep2dyOP1SMpJTy95m8zVI2VTa0BitwedX49vfcdljThPhh7bllRYXSgh++nh fuH85JBMsliPiTbfFyZoyVo9rs7qt3onnMzoMfq3sKT5vuRsPXglNkPL+WSj HnEH1k5p9iJX6HFwoWg0G0Rm9Zjqxaj5anLb8P74b4f9nUc26aHl/sTebSUP Du+HV1h5eqb5/AZcNhatKcsyn9+AtLSEPkdrut5cA6o6SiJ+qyW7GHDy4e0J E/PM3/9w3u7H8Js76X7zM2BqKv/8wSPm+9+A1Yx15tNj9HvSGnBp6JD1H4Hm 368Bx77M+/TKUXK2AZ5X7y/P7ScbDVDGuH1wh0fPkwoDdoUnv1/baH6+GFA0 6CNKiKXnV5sB9iULVD0NK+j8Bnz13rqll7zo/2bQAK/1V/fHWMD1vyMcWu0= "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{}, PlotRange->{{0., 319.}, {0, 6.499295702449204*^-8}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.6741989207069473`*^9, {3.6741992150414953`*^9, 3.674199217263775*^9}, 3.6742004146462593`*^9, 3.674200451249134*^9, 3.674200583596402*^9, 3.6742008473663197`*^9}] }, Open ]] }, WindowSize->{808, 604}, WindowMargins->{{36, Automatic}, {Automatic, 1}}, FrontEndVersion->"10.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (December 10, \ 2015)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 1559, 30, 250, "Input"], Cell[2120, 52, 918, 31, 28, "Input"], Cell[3041, 85, 522, 10, 46, "Input"], Cell[3566, 97, 374, 11, 28, "Input"], Cell[3943, 110, 194, 4, 28, "Input"], Cell[4140, 116, 257, 8, 45, "Input"], Cell[4400, 126, 581, 18, 57, "Input"], Cell[4984, 146, 1156, 35, 80, "Input"], Cell[6143, 183, 1247, 36, 90, "Input"], Cell[7393, 221, 248, 6, 28, "Input"], Cell[7644, 229, 856, 26, 54, "Input"], Cell[8503, 257, 345, 10, 51, "Input"], Cell[8851, 269, 604, 18, 46, "Input"], Cell[9458, 289, 882, 26, 67, "Input"], Cell[10343, 317, 1464, 45, 118, "Input"], Cell[11810, 364, 444, 12, 46, "Input"], Cell[12257, 378, 496, 14, 46, "Input"], Cell[12756, 394, 1076, 31, 102, "Input"], Cell[13835, 427, 546, 15, 46, "Input"], Cell[14384, 444, 609, 16, 46, "Input"], Cell[14996, 462, 542, 17, 28, "Input"], Cell[15541, 481, 279, 7, 28, "Input"], Cell[15823, 490, 509, 10, 28, "Input"], Cell[16335, 502, 1322, 33, 109, "Input"], Cell[17660, 537, 236, 5, 28, "Input"], Cell[17899, 544, 265, 5, 28, "Input"], Cell[CellGroupData[{ Cell[18189, 553, 841, 23, 46, "Input"], Cell[19033, 578, 2283, 51, 241, "Output"] }, Open ]], Cell[21331, 632, 966, 32, 28, "Input"], Cell[22300, 666, 360, 7, 46, "Input"], Cell[22663, 675, 951, 29, 73, "Input"], Cell[23617, 706, 1185, 34, 82, "Input"], Cell[24805, 742, 569, 16, 28, "Input"], Cell[25377, 760, 166, 2, 28, "Input"], Cell[25546, 764, 876, 27, 28, "Input"], Cell[26425, 793, 903, 20, 63, "Input"], Cell[27331, 815, 167, 3, 28, "Input"], Cell[27501, 820, 317, 8, 28, "Input"], Cell[27821, 830, 2561, 65, 155, "Input"], Cell[30385, 897, 250, 6, 28, "Input"], Cell[30638, 905, 1181, 27, 46, "Input"], Cell[31822, 934, 257, 6, 28, "Input"], Cell[32082, 942, 223, 5, 28, "Input"], Cell[32308, 949, 520, 15, 28, "Input"], Cell[32831, 966, 610, 18, 63, "Input"], Cell[CellGroupData[{ Cell[33466, 988, 883, 23, 46, "Input"], Cell[34352, 1013, 5308, 100, 229, "Output"] }, Open ]] } ] *) (* End of internal cache information *)