(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 12.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 14656, 383] NotebookOptionsPosition[ 13472, 354] NotebookOutlinePosition[ 13808, 369] CellTagsIndexPosition[ 13765, 366] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Integrate", "[", RowBox[{ RowBox[{ RowBox[{"SphericalBesselJ", "[", RowBox[{"1", ",", RowBox[{"k", "*", "r"}]}], "]"}], "*", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "a"}], "*", "r"}], "]"}], "*", RowBox[{"r", "^", "2"}]}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "\[Infinity]"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.8315437592165213`*^9, 3.831543812261026*^9}, { 3.831547736280718*^9, 3.83154774639454*^9}}, CellLabel->"In[28]:=",ExpressionUUID->"dd75e20e-0826-4112-9d0a-a15ec853d80f"], Cell[BoxData[ RowBox[{"ConditionalExpression", "[", RowBox[{ FractionBox[ RowBox[{"2", " ", "k"}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["a", "2"], "+", SuperscriptBox["k", "2"]}], ")"}], "2"]], ",", RowBox[{ RowBox[{ RowBox[{"Im", "[", "k", "]"}], "<", RowBox[{"Re", "[", "a", "]"}]}], "&&", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Im", "[", "k", "]"}], ">", "0"}], "&&", RowBox[{ RowBox[{"Re", "[", "a", "]"}], ">", "0"}], "&&", RowBox[{ RowBox[{"Re", "[", "k", "]"}], "\[GreaterEqual]", "0"}]}], ")"}], "||", RowBox[{ RowBox[{"Re", "[", "k", "]"}], ">", "0"}]}], ")"}]}]}], "]"}]], "Output", CellChangeTimes->{3.831543802568611*^9, 3.83154384093207*^9, 3.831547740322646*^9, 3.831547773392627*^9}, CellLabel->"Out[28]=",ExpressionUUID->"d09c14aa-ccb8-428e-ba0b-71e134000b50"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Re", "[", FractionBox["r", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"I", "*", "t"}]}], ")"}], "2"]}], ")"}], "2"]], "]"}], ",", RowBox[{"Im", "[", FractionBox["r", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"I", "*", "t"}]}], ")"}], "2"]}], ")"}], "2"]], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.3"}], ",", "0.5"}], "}"}]}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.831546274584421*^9, 3.831546418750455*^9}, { 3.831547399158922*^9, 3.831547405632867*^9}, {3.831547749630069*^9, 3.831547763666521*^9}, {3.831547852165608*^9, 3.831547861431858*^9}}, CellLabel->"In[34]:=",ExpressionUUID->"7e4f9d43-d282-41d1-b71a-2bd5df2471f1"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`t$$ = 7.19211483001709, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`t$$], 0, 10}}, Typeset`size$$ = {360., {110., 115.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`t$1122070$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`t$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`t$$, $CellContext`t$1122070$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[{ Re[$CellContext`r/($CellContext`r^2 + (1 + I $CellContext`t$$)^2)^2], Im[$CellContext`r/($CellContext`r^2 + (1 + I $CellContext`t$$)^2)^2]}, {$CellContext`r, 0, 10}, PlotRange -> {{0, 10}, {-0.3, 0.5}}, Frame -> True], "Specifications" :> {{$CellContext`t$$, 0, 10}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{405., {169., 175.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.831546325870152*^9, 3.831546353048429*^9}, { 3.831546384493498*^9, 3.831546421381673*^9}, {3.8315474061652603`*^9, 3.83154740986261*^9}, 3.831547773437468*^9, 3.831547808748073*^9, { 3.831547853889669*^9, 3.831547864205819*^9}, 3.8334612426451693`*^9}, CellLabel->"Out[34]=",ExpressionUUID->"0680386d-6bec-4e12-9b49-bfa391341d65"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"r", "^", "2"}], "*", SuperscriptBox[ RowBox[{"Abs", "[", FractionBox["r", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["r", "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"I", "*", "t"}]}], ")"}], "2"]}], ")"}], "2"]], "]"}], "2"]}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.01"}], ",", "0.1"}], "}"}]}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.831546274584421*^9, 3.831546418750455*^9}, { 3.831546659744111*^9, 3.831546683243444*^9}, {3.8315473779484243`*^9, 3.831547378170991*^9}, {3.83154777755763*^9, 3.831547788205557*^9}, { 3.8315478278817997`*^9, 3.831547835781007*^9}}, CellLabel->"In[32]:=",ExpressionUUID->"61472b1b-1939-4e03-a284-0cab90acc208"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`t$$ = 7.507577896118164, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`t$$], 0, 10}}, Typeset`size$$ = {360., {112., 116.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`t$1118830$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`t$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`t$$, $CellContext`t$1118830$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[$CellContext`r^2 Abs[$CellContext`r/($CellContext`r^2 + (1 + I $CellContext`t$$)^2)^2]^2, {$CellContext`r, 0, 10}, PlotRange -> {{0, 10}, {-0.01, 0.1}}, Frame -> True], "Specifications" :> {{$CellContext`t$$, 0, 10}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{405., {170., 176.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.831546669483881*^9, 3.831546687213415*^9}, { 3.8315473786378727`*^9, 3.831547381377433*^9}, {3.831547789311469*^9, 3.831547794477837*^9}, {3.831547829400641*^9, 3.831547838521158*^9}, 3.833461275225813*^9}, CellLabel->"Out[32]=",ExpressionUUID->"b45acdcd-bc47-481a-b859-fcb91a39c8e6"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"rad", "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ct", "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", RowBox[{"y", "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], "]"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.831547520673568*^9, 3.831547548760683*^9}}, CellLabel-> "In[612]:=",ExpressionUUID->"85561817-1a68-4253-bcb7-e1d8652fe26c"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"DensityPlot", "[", RowBox[{ RowBox[{"Abs", "[", FractionBox[ RowBox[{ RowBox[{"rad", "[", RowBox[{"x", ",", "y"}], "]"}], "*", RowBox[{"ct", "[", RowBox[{"x", ",", "y"}], "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"rad", "[", RowBox[{"x", ",", "y"}], "]"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"I", "*", "t"}]}], ")"}], "2"]}], ")"}], "2"]], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "10"}], ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "10"}], ",", "10"}], "}"}], ",", RowBox[{"PlotPoints", "\[Rule]", "30"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.831546274584421*^9, 3.831546418750455*^9}, { 3.831546659744111*^9, 3.831546683243444*^9}, {3.8315473779484243`*^9, 3.831547378170991*^9}, {3.83154751540266*^9, 3.831547516285693*^9}, { 3.831547568000345*^9, 3.8315476326753883`*^9}, {3.8315476678784437`*^9, 3.831547700511856*^9}}, CellLabel-> "In[614]:=",ExpressionUUID->"ffc5e311-a7e6-4145-9d42-d48c7394883d"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`t$$ = 5.105257987976074, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`t$$], 0, 10}}, Typeset`size$$ = {360., {176., 180.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`t$194911$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`t$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`t$$, $CellContext`t$194911$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> DensityPlot[ Abs[$CellContext`rad[$CellContext`x, $CellContext`y] \ $CellContext`ct[$CellContext`x, \ $CellContext`y]/($CellContext`rad[$CellContext`x, $CellContext`y]^2 + (1 + I $CellContext`t$$)^2)^2], {$CellContext`x, -10, 10}, {$CellContext`y, -10, 10}, PlotPoints -> 30], "Specifications" :> {{$CellContext`t$$, 0, 10}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{405., {234., 240.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.8315476346468897`*^9, 3.83154763807407*^9}, { 3.831547671391802*^9, 3.831547705395103*^9}, {3.833461292263242*^9, 3.833461294574404*^9}}, CellLabel-> "Out[614]=",ExpressionUUID->"c1ab6fc5-9c56-47db-8f1e-f062fe9930ce"] }, Open ]] }, WindowSize->{1095, 911}, WindowMargins->{{4, Automatic}, {Automatic, 4}}, FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 586, 15, 44, "Input",ExpressionUUID->"dd75e20e-0826-4112-9d0a-a15ec853d80f"], Cell[1169, 39, 977, 29, 58, "Output",ExpressionUUID->"d09c14aa-ccb8-428e-ba0b-71e134000b50"] }, Open ]], Cell[CellGroupData[{ Cell[2183, 73, 1585, 45, 97, "Input",ExpressionUUID->"7e4f9d43-d282-41d1-b71a-2bd5df2471f1"], Cell[3771, 120, 2136, 40, 363, "Output",ExpressionUUID->"0680386d-6bec-4e12-9b49-bfa391341d65"] }, Open ]], Cell[CellGroupData[{ Cell[5944, 165, 1330, 36, 99, "Input",ExpressionUUID->"61472b1b-1939-4e03-a284-0cab90acc208"], Cell[7277, 203, 2038, 40, 365, "Output",ExpressionUUID->"b45acdcd-bc47-481a-b859-fcb91a39c8e6"] }, Open ]], Cell[9330, 246, 649, 20, 80, "Input",ExpressionUUID->"85561817-1a68-4253-bcb7-e1d8652fe26c"], Cell[CellGroupData[{ Cell[10004, 270, 1390, 38, 68, "Input",ExpressionUUID->"ffc5e311-a7e6-4145-9d42-d48c7394883d"], Cell[11397, 310, 2059, 41, 527, "Output",ExpressionUUID->"c1ab6fc5-9c56-47db-8f1e-f062fe9930ce"] }, Open ]] } ] *)