Blog: Update 4/20: charge_exchange_equilibrium.nb

File charge_exchange_equilibrium.nb, 22.6 KB (added by adebrech, 5 years ago)
Line 
1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 12.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 158, 7]
13NotebookDataLength[ 22935, 513]
14NotebookOptionsPosition[ 21245, 477]
15NotebookOutlinePosition[ 21580, 492]
16CellTagsIndexPosition[ 21537, 489]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[{
22 RowBox[{
23 RowBox[{"Rp", "=",
24 RowBox[{"1.07", "*",
25 SuperscriptBox["10", "10"]}]}], ";"}], "\[IndentingNewLine]",
26 RowBox[{
27 RowBox[{"\[Nu]0", "=",
28 RowBox[{"2.4661", "*",
29 SuperscriptBox["10", "15"]}]}], ";"}], "\[IndentingNewLine]",
30 RowBox[{
31 RowBox[{"c", "=",
32 RowBox[{"3", "*",
33 SuperscriptBox["10", "10"]}]}], ";"}], "\[IndentingNewLine]",
34 RowBox[{
35 RowBox[{"mH", "=",
36 RowBox[{"1.67", "*",
37 SuperscriptBox["10",
38 RowBox[{"-", "24"}]]}]}], ";"}], "\[IndentingNewLine]",
39 RowBox[{
40 RowBox[{"kB", "=",
41 RowBox[{"1.38", "*",
42 SuperscriptBox["10",
43 RowBox[{"-", "16"}]]}]}], ";"}]}], "Input",
44 CellChangeTimes->{{3.796138014685924*^9, 3.796138017232429*^9}, {
45 3.796138077065893*^9, 3.796138084092134*^9}, {3.796138132383556*^9,
46 3.796138138000414*^9}, {3.7961382547993307`*^9, 3.79613827156409*^9}, {
47 3.796138368361528*^9, 3.796138418969956*^9}, {3.79613893568709*^9,
48 3.79613893835375*^9}, {3.7963812099424686`*^9, 3.7963812135603952`*^9}},
49 CellLabel->"In[40]:=",ExpressionUUID->"253f762e-8af6-4b50-b652-0f3b4a07e0c9"],
50
51Cell[CellGroupData[{
52
53Cell[BoxData[{
54 RowBox[{
55 RowBox[{"n0c", "=",
56 RowBox[{"3.1", "*",
57 SuperscriptBox["10", "4"]}]}], ";"}], "\[IndentingNewLine]",
58 RowBox[{
59 RowBox[{"npc", "=",
60 RowBox[{"5.4", "*",
61 SuperscriptBox["10", "5"]}]}], ";"}], "\[IndentingNewLine]",
62 RowBox[{
63 RowBox[{"n0hs", "=",
64 RowBox[{"2.4", "*",
65 SuperscriptBox["10",
66 RowBox[{"-", "2"}]]}]}], ";"}], "\[IndentingNewLine]",
67 RowBox[{
68 RowBox[{"nphs", "=",
69 RowBox[{"3.8", "*",
70 SuperscriptBox["10", "4"]}]}], ";"}], "\[IndentingNewLine]",
71 RowBox[{
72 RowBox[{"n0hw", "=",
73 RowBox[{"7.2", "*",
74 SuperscriptBox["10",
75 RowBox[{"-", "6"}]]}]}], ";"}], "\[IndentingNewLine]",
76 RowBox[{
77 RowBox[{"nphw", "=",
78 RowBox[{"7.1", "*",
79 SuperscriptBox["10", "2"]}]}], ";"}], "\[IndentingNewLine]",
80 RowBox[{
81 RowBox[{"nph", "=", "nphw"}], ";"}], "\[IndentingNewLine]",
82 RowBox[{
83 RowBox[{"n0h", "=", "n0hw"}], ";"}], "\[IndentingNewLine]",
84 RowBox[{"\[CapitalDelta]", "=",
85 FractionBox[
86 RowBox[{
87 RowBox[{"n0c", " ", "nph"}], "-",
88 RowBox[{"n0h", " ", "npc"}]}],
89 RowBox[{"n0c", "+", "nph", "+", "n0h", "+", "npc"}]]}]}], "Input",
90 CellChangeTimes->{{3.795538356211307*^9, 3.7955384496723967`*^9}, {
91 3.7955385238874807`*^9, 3.795538616326481*^9}, {3.795538741529235*^9,
92 3.795538792052253*^9}, {3.795538891253086*^9, 3.795538913574039*^9}, {
93 3.795540411120379*^9, 3.795540429260208*^9}, {3.795541613466384*^9,
94 3.7955417243127117`*^9}, {3.7955418043978567`*^9, 3.7955418079894543`*^9}},
95 CellLabel->"In[45]:=",ExpressionUUID->"de1a34d1-f494-4416-a345-371b7a127382"],
96
97Cell[BoxData["38.49853266817584`"], "Output",
98 CellChangeTimes->{{3.795538774354725*^9, 3.7955388113443537`*^9}, {
99 3.795538895221807*^9, 3.7955389148559523`*^9}, 3.795540431172904*^9, {
100 3.7955417326877604`*^9, 3.795541756871379*^9}, 3.795541808596273*^9,
101 3.79638099598412*^9, 3.796381215778967*^9},
102 CellLabel->"Out[53]=",ExpressionUUID->"14e26464-572f-47dc-b947-58dc631736c6"]
103}, Open ]],
104
105Cell[CellGroupData[{
106
107Cell[BoxData[{
108 RowBox[{"Clear", "[",
109 RowBox[{"n0h", ",", "nph", ",", "n0c", ",", "npc", ",", "\[CapitalDelta]"}],
110 "]"}], "\[IndentingNewLine]",
111 RowBox[{"Solve", "[",
112 RowBox[{
113 RowBox[{
114 RowBox[{
115 RowBox[{"(",
116 RowBox[{"nph", "-", "\[CapitalDelta]"}], ")"}],
117 RowBox[{"(",
118 RowBox[{"n0c", "-", "\[CapitalDelta]"}], ")"}]}], "\[Equal]",
119 RowBox[{
120 RowBox[{"(",
121 RowBox[{"n0h", "+", "\[CapitalDelta]"}], ")"}],
122 RowBox[{"(",
123 RowBox[{"npc", "+", "\[CapitalDelta]"}], ")"}]}]}], ",",
124 "\[CapitalDelta]"}], "]"}]}], "Input",
125 CellChangeTimes->{{3.795540980427651*^9, 3.7955410148753357`*^9}, {
126 3.795541096827845*^9, 3.7955411804285192`*^9}, {3.795541213028809*^9,
127 3.7955412138899603`*^9}, {3.795541387900722*^9, 3.79554138943513*^9}, {
128 3.795541736746655*^9, 3.795541753228962*^9}},
129 CellLabel->"In[54]:=",ExpressionUUID->"76bb0e2c-dc81-413a-bed5-2921b26ee82b"],
130
131Cell[BoxData[
132 RowBox[{"{",
133 RowBox[{"{",
134 RowBox[{"\[CapitalDelta]", "\[Rule]",
135 FractionBox[
136 RowBox[{
137 RowBox[{
138 RowBox[{"-", "n0h"}], " ", "npc"}], "+",
139 RowBox[{"n0c", " ", "nph"}]}],
140 RowBox[{"n0c", "+", "n0h", "+", "npc", "+", "nph"}]]}], "}"}],
141 "}"}]], "Output",
142 CellChangeTimes->{
143 3.7955412148583937`*^9, 3.795541390287533*^9, {3.79554174761963*^9,
144 3.795541753965826*^9}, 3.796380998423135*^9, 3.796381217021764*^9},
145 CellLabel->"Out[55]=",ExpressionUUID->"0994e65f-46ee-4409-8255-4333b6653546"]
146}, Open ]],
147
148Cell[CellGroupData[{
149
150Cell[BoxData[{
151 RowBox[{
152 RowBox[{"\[Sigma]\[Nu]0", "=",
153 RowBox[{"1.102", "*",
154 SuperscriptBox["10",
155 RowBox[{"-", "2"}]]}]}], ";"}], "\[IndentingNewLine]",
156 RowBox[{
157 RowBox[{"n", "=",
158 RowBox[{"2", "*",
159 SuperscriptBox["10", "3"]}]}], ";"}], "\[IndentingNewLine]",
160 RowBox[{
161 RowBox[{"d", "=", "Rp"}], ";"}], "\[IndentingNewLine]",
162 RowBox[{
163 RowBox[{"T", "=",
164 SuperscriptBox["10", "6"]}], ";"}], "\[IndentingNewLine]",
165 RowBox[{"\[Sigma]n0h", "=",
166 RowBox[{"n", " ", "d"}]}], "\[IndentingNewLine]",
167 RowBox[{"\[Tau]", "=",
168 RowBox[{"\[Sigma]\[Nu]0", " ", "\[Sigma]n0h", " ",
169 FractionBox["c", "\[Nu]0"],
170 SqrtBox[
171 FractionBox["mH",
172 RowBox[{"2", "\[Pi]", " ", "kB", " ", "T"}]]],
173 SuperscriptBox["\[ExponentialE]",
174 RowBox[{"-",
175 FractionBox[
176 RowBox[{
177 SuperscriptBox[
178 RowBox[{"(",
179 RowBox[{"\[Nu]", "-", "\[Nu]0"}], ")"}], "2"],
180 SuperscriptBox["c", "2"], "mH"}],
181 RowBox[{"2", " ",
182 SuperscriptBox["\[Nu]0", "2"], "kB", " ", "T"}]]}]]}]}]}], "Input",
183 CellChangeTimes->{{3.796137915853367*^9, 3.79613800867592*^9}, {
184 3.7961380878255157`*^9, 3.796138129243558*^9}, {3.7961382814147673`*^9,
185 3.79613830360544*^9}, 3.7961383659446*^9, {3.7961384317656727`*^9,
186 3.796138470838286*^9}, {3.796139078341855*^9, 3.796139080545848*^9},
187 3.796139301833722*^9, {3.7961417213644543`*^9, 3.796141724077568*^9}},
188 CellLabel->"In[56]:=",ExpressionUUID->"d01c3e8e-c5ee-41e5-a4e3-5e218bfbf95a"],
189
190Cell[BoxData["2.14`*^13"], "Output",
191 CellChangeTimes->{
192 3.7961384723373013`*^9, 3.796138956291709*^9, {3.796139054763336*^9,
193 3.796139080924264*^9}, 3.7961393021630774`*^9, 3.796141726486257*^9,
194 3.7963809991647987`*^9, 3.796381218016006*^9},
195 CellLabel->"Out[60]=",ExpressionUUID->"8fd7ff38-9a2a-419e-a87c-b8cff9652dbd"],
196
197Cell[BoxData[
198 RowBox[{"0.12590259965726647`", " ",
199 SuperscriptBox["\[ExponentialE]",
200 RowBox[{
201 RowBox[{"-", "8.954235908507858`*^-25"}], " ",
202 SuperscriptBox[
203 RowBox[{"(",
204 RowBox[{
205 RowBox[{"-", "2.4661`*^15"}], "+", "\[Nu]"}], ")"}],
206 "2"]}]]}]], "Output",
207 CellChangeTimes->{
208 3.7961384723373013`*^9, 3.796138956291709*^9, {3.796139054763336*^9,
209 3.796139080924264*^9}, 3.7961393021630774`*^9, 3.796141726486257*^9,
210 3.7963809991647987`*^9, 3.796381218018855*^9},
211 CellLabel->"Out[61]=",ExpressionUUID->"e014cc80-ca01-41ee-943c-882968fdb6b8"]
212}, Open ]],
213
214Cell[CellGroupData[{
215
216Cell[BoxData[{
217 RowBox[{"Integrate", "[",
218 RowBox[{"\[Tau]", ",",
219 RowBox[{"{",
220 RowBox[{"\[Nu]", ",",
221 RowBox[{"2.4634", "*",
222 SuperscriptBox["10", "15"]}], ",",
223 RowBox[{"2.4684", "*",
224 SuperscriptBox["10", "15"]}]}], "}"}]}], "]"}], "\[IndentingNewLine]",
225 RowBox[{"Plot", "[",
226 RowBox[{"\[Tau]", ",",
227 RowBox[{"{",
228 RowBox[{"\[Nu]", ",",
229 RowBox[{"\[Nu]0", "-",
230 RowBox[{
231 FractionBox["1", "1000"], "\[Nu]0"}]}], ",",
232 RowBox[{"\[Nu]0", "+",
233 RowBox[{
234 FractionBox["1", "1000"], "\[Nu]0"}]}]}], "}"}], ",",
235 RowBox[{"PlotRange", "\[Rule]", "Full"}]}], "]"}]}], "Input",
236 CellChangeTimes->{{3.796138476473015*^9, 3.796138494635365*^9}, {
237 3.796138974046515*^9, 3.796139047342565*^9}, {3.7961391665994673`*^9,
238 3.796139197749282*^9}, {3.796139264668872*^9, 3.7961392889136467`*^9}, {
239 3.7961410623868303`*^9, 3.79614119546525*^9}, {3.796141420908752*^9,
240 3.7961414446145287`*^9}, {3.796141483938963*^9, 3.796141502862116*^9}, {
241 3.796141585200045*^9, 3.796141641148716*^9}, {3.796141736758945*^9,
242 3.796141748942387*^9}, {3.796141803327044*^9, 3.796141875725492*^9}},
243 CellLabel->"In[62]:=",ExpressionUUID->"7f83b4b0-1af8-44d0-9571-c30cdb9bb101"],
244
245Cell[BoxData["2.355465450328643`*^11"], "Output",
246 CellChangeTimes->{
247 3.796138495756485*^9, {3.796138960221661*^9, 3.7961389756057653`*^9},
248 3.796139008292447*^9, {3.796139048401775*^9, 3.7961390572076674`*^9}, {
249 3.7961391715663977`*^9, 3.796139198209818*^9}, {3.7961392717440434`*^9,
250 3.7961392893209352`*^9}, {3.7961410651450157`*^9,
251 3.7961411231684437`*^9}, {3.79614116237171*^9, 3.796141195850403*^9},
252 3.7961414290095377`*^9, {3.7961414929333878`*^9, 3.7961415031613092`*^9},
253 3.796141644337055*^9, {3.796141852490193*^9, 3.796141876389929*^9},
254 3.7963810011419907`*^9, 3.796381220217395*^9},
255 CellLabel->"Out[62]=",ExpressionUUID->"4020ff30-3750-435c-9995-c6228e6e1901"],
256
257Cell[BoxData[
258 GraphicsBox[{{{}, {},
259 TagBox[
260 {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
261 1.], LineBox[CompressedData["
2621:eJxF2Xk81F0XAHBr2bckpEhSSik8ttQpIiSPNSL7LoWsKftOZN9D2ZU9EWNm
263lC1LqaQsT5YoiSSRFL3n99f713y+nztz595zzz33d2f22Ljp29PR0NBU09LQ
264EK82TUtf9UP2qAaaf+PX6T8Nwj5qhQZoQa7w+bdfTsOYTJaxITpzx0utJBYV
265MK4+RTVCD0RV3E/RUAHtkttJJmg6hVcH33eqgFzqYVkL9N+j/F7tfarA6u58
2663Rn9015XrumvGnQfbj3sgi4kfRx8JKYOkV+4pgjrf3WrXNBUBxrHJi1XNMvr
267y3KMaeqwar5VyA09YfZ9/MPhszBxrpjshY5RbTa+56IBD/dP0gejR3WMOrR4
268zsGTb5ZdhE04L2ZVqJyDgcfvY0PQkwrbZE08z8HCuTHuMLTYgu7VvDfnQNz9
269rXAkeijl+vTUHW3IaOo/fgvdHNs1+OCEDvhrNntmouckyvrr3utCNI+ifBZa
270iZnHUvKPLqSPNq4TfjEvyTIooAf1VxpCctDZA5FFVUZ6MJ9UczsP3Tcp+zjj
271hR5YjJTcL0KfuuVe496nD6cvp3yoQacyN2UwzhuCVs2u7Fr0tUf6/JbsRmC4
272Uqpbh4aTlgufjhiBUxCptR6tzRTuF33NCBJTZjIeocvPSpuw/jGCiRZ5bRJa
273LecJ5zl+YwhmG3vYhf5SX5tl4HsRYvUcLncT81es7GDOvwip6d/2PEMX9Nd3
274rnZehDIRxts9aOrBBjWfHaYwIHvEuR8d6bskw9NiCiLmwbteo/+bKbddZbsE
2751EqxyPdobveRu/CfBShoCiqOo1P+uWtnw2oJ1dOc84QLFVtG2xQtIW/nb71J
276Yn1anRblMyzhZszLXdPo+Jhbv8IOWoGCfWDDZ3SLCJuo4lcrqBF6N/0D/VtZ
2774/aBBBs40NSfsULEJ1ZM+Gy5DeQbPNVaRYuws5tFt9tAQmxVzU90+/VBR8/f
278NuC6Fn5zncjfp9cGFl1s4cDgMV6a0D2qaz1fme/p2kF+XJwqK7oklyX/u7wD
279OEk/YGRDL89PclmbOsDR4b4uwpVVTu83bjoAVZxDmwMt/c7Y6OkTB5hqSzTk
280Rue0LsmH6TqC+Fqa/Q50T92ME723EzywK4jai96sL9nr8cEFvFjbNMXQBTHx
281zznYLsOJuknWfejXNSSzSdnL8PyvaKI4+mqT4+q+qMvwLbs4UwJd5XRIXOaI
282K8i8rCiXQluzrtyNj7wCzScaeo+jZ6leUTNH3SHQajJYGa0oWuggctYdVMLY
2835U6gffkyzbzN3aG3277gJNqLwXfmSqw7jBnweZ8mPj81wp074w4bzr67z6Lj
284fmQ00t71gFPpih766Hf3mhtzJD2B4bG9uAHaYR/d0ll1T+geTRolnLh8SZPb
285yhP09sypGaE7hna2/E7xBJsH2YIm6MuONyylNzwh/Mnvp+boyNsHa5TfeEHX
28611Y+Z/QZGoqTb64PjAltHyX8M+T2pWskH/iu5ZrvQrT7GJxJGfMBoVLBA67o
287Ig15r7O7fMHd0lfRDS0SdfznzQJfEBg4eskbXZrxX5NilR9cri28G4KmcegP
288nPnkD8Hj6/ah6E/w+dU42w1IZ9c/GIZmfOajyyB9A6jONPXh6OaVe7p9N28A
289r6h5exTaXKR37yXem9CazPcpHv1MqkMo61wAcHrFSGYT+XGoKeP6aBAkr51b
290IJzMr3glmy4YeAM4qnLQBpk1ngk7gkEwMkXqDvotY87S49PBIJZ5R7oAXc0m
291HSCSEQyKpFqFEvR9vQAZJfUQsGUYOVOHZuoSFNz2KBQaUw+a96C7/a8fXRyP
292ALm0Na5edBTbYeD7FQEP0zrbCWtlCxsZ8kRCTbqNZD86YjmiZ7taJFRkZv9+
293geYZnTxTfT8S8nJZst+gC3+LNZ4NiILwormhCfT+t+1FjtIxoPewQvcnevfa
294+0bBT7dAhFq4bQ2tmjAXx8ocD996c4cI28z94dl5KB4SPiRcWke/qfb1SneL
295hz4eT6cN9FLwigDNejyoexwPpg/D+rErxe0+721QPtJXzYm+Km7UPVmTCOJl
296CxwH0Dxetsmfj6TAktT3dsJdP8QvRqmkAKlx1V8CrfMifVrmQgrod/39eBC9
297c7nI9HFACgR95KYcRi8/+qjo058C78Tk3GTQKXOHLVrcUiH2buCLk+gdMrKG
298CZQ0+JrNmWiE5vNo9uSJzwQ+ia2SoWi96vAmhrY8uB7XllyHltNmaB9SK4Tn
29986akcXSxG8f7VutSsBBtrOEK36O6vfuG/R3NB6CzPEs5hXabNrvM2FwDx57+
300SrqC/lSV/PWxYj3869h/JA9NCi/ReLDUAJYswQm9aF4Z1vZBuyaoCFlp+oP+
301ZbQtLXqoGVYqWfw30A93VbPHjDfDqRHh45voKwmaKuGzzfBGWqv1L7q/9d8n
302LuvNQDOd10YXgf1x0Ov/3NUCRuoaPUzoIf8KMya7FqBhzRnlRZdz7md996UF
303tORrcrejvW938Yn+aIE02w5zPvSCZ9C6/Z8WOEhaHN+B/quVfGSUjQRGV85M
304C6JrcvKj/Q+ToOLFwrwI2jK+nTXqCranntqURPcr91mUzZIgV9m67jC6T2xM
305tWCRBFPTIQ5H0IIvuXmTV0lw9Z+nz6XQD57KWdsytEL0kFq+NLpOTGo3SbgV
306SPznTimg1+SvV542agWxO0ahqmhrG7cz/i2t4KLmI3cGva4031/c1gq1C+lz
307hN1Gzxb1dLXCyZPv9NXRBw7OidEMtoLJuOleTXTk2p9lgflWiNtj1a6DNtmi
308/YxhJxm+F7tsMUVLVUaqZXqSIfbio2eEFSK0My9cJ8Medrp4M7Rmsow9WxAZ
309dL2ytpmjZSWKTCxiyVCp0i1ihaZhvMCXU0AGp/F9xx3QobqLkZ69ZKBJ8dgk
3103B7BJVQ6QIYM9dY2R3Q9T8boyzdk6Kw01HBG25fMdrJNkGHvzTAjV/SL0Te9
311O36Q4T3/pNs1tPSt4aRAQQr49EnKeKJPyqcViwpTgD3Yb5XwqfWInNa9FFD+
312xBngjS6rY9/1TpICWQ0nY/3QwsaltPknKWCon1sUSIzf6fBTdysK9MYZD8cS
31362ng7cJ8hwIWLVW344h4jCkFni2gwPc5RvVb6GvR/PIBhRTg13pYF482Y62N
314f1FOATsm7luJxHwejUVyPqLAZnjPyXT05qxy2HofBZIe7lkh3CXS1Vr7ggJi
3150373M4j5SnX4WL2igJbqfv4s9CVLK7OitxTIoA1fykFvP3C5vmOSAkcDTxTd
316Rb9V3vfO9QcFnlammt5Dj5Vx0eWsUuDCf1+4Con++3NjqWsUCFDODihCO11J
317q1n6Q4Ge9dULpehc6mkuHzoqmB/UYS9DG9WOR7nQU+HbxeKnhJ/Y7i8zZaAC
31832ODoxVEfEsKxw9voYKtbw1zJRGPH/e2FzBTYbVkK5WwZOqzVg8WKsQOWfhU
319oXervWECVirU/MP+oRp9duNPeB8bFf4sO5PqiP3XukFfz0mF23vbrtWj1acK
320a625qCBqwC/xkGhX19lg5aaCZl1nWgP60dUwVwMeKqS773VrQge20+ZW8lIh
321Y3Vmz2P0h6wQP7ntVMi8WTZIuFxzSLAFnRV9WKkF7cxpUd3IR4Vsjm/zhHc4
322PyuQ2oFOrcsnoaP3+l8oROcWyDOS0Xa0p9kD+dHi642ETfZUsM+g7zwguVDQ
3231fmbA2cFqJDXdHqASrSHyubTCFIh/yRDWBt6uTKs/AK6oL3znyfo0o4s7zL0
3243YFzOU/RAhI79WEnFe5d4NBpR/Nzc9qEEx4boOlA79k1Kt2BLrRJqSe8SqPd
325TitEhaJZI4dOIl96j/EpEb7KL9CFbvpkeeAquvjHSC/h+t2V63fQJf53ArvR
326rGnLt54RprE69gyd0cc1/g1dGik6TZi+eObbtl1UKGObSe8h6pvphSfShJNL
327NXvRmVFaBufR5fwufwh37K0stEVX5ElW96Fjwl3rvQmLLVr3o/eZ+4eGoe9X
3281PI+R/8i93LFox8c9eoivFvgokUy4Udy/i/QTxV2X05BVyr/khwg5q8pKJ+I
329rnrSMk746ut/n0QT1ghMfknsJ3I7WwC6+vkptVfoz7qhfFfQNYb0a4QdJdPG
330TQiPdFS8Jva/PrvTKXStVbT5IFHPRDer96LrPmpxvSHiZxzYSEfYlf0pYV/2
331ssD/MD713194D6ETmjMZH6If+iUfeIvekPHQjiS8aThKeO2+kYEhuiF8R8I7
332dHfNFcHd6EeJucvD6Peunz4X4Xo28lmWjBD1hqv8uzW6KXfPxVEiPvwWJEH0
33347IS8hiaeX9UShDmS/MRZ4//iHxd9b57iPDDQ2LvifOMfN7jNeYfiVoTM07U
334g4evjPkIq3sqT6AZQra71WK+tvb9s0i4M9hVUwNNftdsOEWcP0dKNZwx39u+
335tQvPoH/8ULrhh/uFJ7RukvDrXwnbPuL+sttWUPiRiC9fcp0Omknuxv5ZNG3e
336f1u4t1FB58YxqS/o+b/bBExx/xawCS8RDk/V+pmP+3vpDlv9PDr9947JCQ4q
337pFI/yX1Fz1qrjBiwU2GEMe/kEvrOhK3vJtYTyYw4mu/EeCQpnULogAPXnxDu
3389iDvl2Oigsg5Q/UfaHeeLaKmWI8cE1l0fhLnmdIZLkesZz8Efc030CKHtM8f
339WqaA+gM74U0iH9j4Nnq+YT09oT9JWNNZYtHqKwWUrA7b00TiffKjeKHbZwqE
340FH9wpUdvXVGP4x+nAMdR3ZvM6CszLkHZ3RSwajtxkgW99MjLM7+DArX6h2hY
3410Vq+tS+zn+B54r0lnA2dcuzbWy8SBbKbSXGc6L9ld8weVVNgv5pEznZ0RWbT
342mEcqBbz0RXX50D23c08JJVHgieVOxh1osneDNzmeApf82d340fzB+hZzkRRI
343rl5S2Ynee0DaZ8KPAhsCzXMi6CmxOWUNMwq8/qqlJIkuqBTW2Ynno8gf1UXC
344fLsNYjq2U+Aq84miw+jIf+d+OnJTgElMikMK3SiyKyCNiQInL26bOoZ22+yO
345SVklQ/nT0Rh5dBEndWD3SzIEZboOq6Adu1SsbwWTYXdt4h1V9Gu1oMBX/mQg
346PXtofQat0pNMy+VNhrX1P5/V0MxMzjoezmS4Zh6/roE+PxaT2a9HBjvRqp3/
347oicGPeuXRMig+WDR3AzNsPVWEufjVthGvTbliT75itVh+SUJdkvRFXqhdavu
348j9n1kOBgXpKtN/o5W3/QiyckULlZM+2Dfkk1s0+sJ8E1ha8fr6MH2/ht21JJ
349MFDj/CUIHZC2PjVjRIL4u1Y/4tBCD/LzHF+3wNYwnS1F6LAB65gXpGZYOMO8
350MID2XpEPlAhpAo7RyDzaKKwXH5nKKiYaYI/7Bx4p9KYR58f3RvWQWrD9nlkU
3518XuMz4rucA2oFrP8iUMrN4k7igU8gOmy9/yP0b2hihYzi6VADfammUOXmQ7P
352Zl0rhATjLxU7ovF52/QDT69fPsZTzv0MuutpPP2DfVkwROcrrYb2+9ijYsCW
353BYvHG38QJgsMn2VbzgTRKvnrZ9EnRoOz66mZEJOsEKqFNrkbJ/L2UiZcMFNK
3541UXH9o54u2VmwNL8iaZLaMXhxfsDQukgzqVG442eLDOLiryYAsodiQ2EaU21
355PL+opYCe/5iLD3qsvX6LpXQK3Jz2fOOL/pj2jucaawq8arpX4Y+2C1hRP0tO
356hgBrGqNgNP1YpVKFeDIM1rWUx6FT8obJx+iSINRI2uAe+vyC1KRpZzzMbi/f
357UYi24H7nfrMwHv4dEh4jXNsTcrU1OB6EjNnti9Gln1773TkeDw0msz5l6Eta
358kRKDdbdg1iwvuwpt+HvHIEdpHOjYsEw1o8tOT7WfL4mBh6IhJS3otc+LBRAV
359Azs//HQhoRsmjm3oOMXArO3Mcit6w3JjpeNgDITYUxnb0FEfA8S21UbDQycf
360iS60K8OagtDTKBBwm3J/jbZ2FR3z/hMB8xHtxwbRHop9zH6TEUDNLflO2Msr
361uDy1MwKcely8htATL1Y2hZIioGnfsu8wEc/xtFNO+yPgwhh90Dixfm48jKMm
3624ZCquS/+C3pecJaF+jwUnKy26swT8ZBe37+E92Fl388cC+hOjeJqpfxQmC6q
363TPyKTmyXaOVzDwUZGrnUJfSTuzqMSzyh8OqRes5PNEcK82dh8xDgEnMqp4/B
3645/XVCwPLeF+feGRkzoD2mZfyWnofBNWaqtyM6IjcWRnRpiDQcdt9fQt6QVCR
365fc0lCOKbBzWY0V1b9JWvvQoEFr3TsxzojNIrSjFlAcAQKLhfEK2irFMR5X4D
366BrmYRgkn11n5Hte5AYWFKwk70TV9p8d5JG+AyrOBVSH0xKwkn+SsP4TwRnUK
367o3OkR04x2vjDZsWy/T70Kkk+5/il6/DrbV/RUfRZ93/Y6Vx8QdBnxvEY+mqD
368c5Wbpi8o8W4elEafWPQzppfwhZu6UrUy6Kz2+pO5sz6w2Z3UKoduvXLnVaKz
369D9A1XxhSRnM2Dz92uuYNrLkTWzXRzb/D2JbzPEFS6VcP4bdks1MfQj1B+x13
370ghYxvimOgBUHT0jgVeXVRssGnkwNkvIErvhikX/Rupy7Xos/uQa8gS6KRsT7
37199RFlyx4wC6rHy7WaLXi+XNpl9xByZ9/jw1a45vKuKi6O1xIVX5LuC9w9nmP
372lDskdIer2qHLlcHWmd4daI5tF3JEmzGMR7Yau8EUrVy/K/F9aatnDtFfhbJC
373Xyk/9NgJvkyRoMvQ3pozQ7g/I7nc1+QyTL6l5FxH27mE3ds4dhkE2ZiYbqB/
374/7eenjbjAnFemRMB6PFFKP6m4wJX1R4nhaGLroyGSB50BtlP68sJ6JZz1/ti
375aB2BbBxbcptY39yMLxyTDqDRLXAxEa1YMfupheoAl8oVyEnoR04jI7eCHSDs
376sk90KnpIwdBmC50DvFpcEspGbzREssix2YP72qxaMTEes26rAGVbWHfyWyNs
377cni4o17EFsKHt94vQdf6FWmzM9pCZrM4Vxl6O9nriMRzG6DctBupQJuGNejR
3782NgAB+3E1Rr0qTdK79kSreE+89uMFqL/ExJOFQxWEBmRdJGEZn1t+33zjSVY
3790WrvbEWHZDZFL92zBL5f1Dwymq5qklJ1whJCZyuK24j1ND2zP9jHAky6Auu7
3800BdEU2TlVi4BQ4T4i9forbQ7JZvETGGcZjJxEG1135Ht1a+L0HwzR/8NuvG3
381ld7BFxfBzYtraIjIN0WF4VT/izBiuz42jI6lauVpDppAjcrzuXH0jMBK97dk
382YzCn8dkyj+b1g2LF40Ygu3a8hTBzp4jfdwEjYPlG676A/qu71jO8ZghN4/HD
383X9HWncq8Mo2GwEMpfrCEThRoU2SRM4TOgDf6P9FUNrb0ncoGIPVbNp8uFs+X
384n7blpk56sGV53YAe/Wdvk9lNbT0Ym6MyMaC/nx/c2nlUD2JHtK8xogUXOGIG
385funC9GM7NSY0k4LLKCleFzL9Ur+woxl9H7gmXf8XaFeX5QXQ+w6crutz1YZu
386TiqtIPpXo7XxDwVtSJC41Uu4FNqFPRi1QdBczFII3RltKu+Ufw5k2w0jhdFT
387GQ0u8FYLnJIbBvehQ8ilkwcvaMLAEV8PaXRCnoTwcLw6pGuoHpdBtw9bcHbY
388qoO5DSejLJqktLFlTVEd5tJKM/9BZ5UkuRz7pAYMG28pCuhJm581H8+ogXyv
389AgegV/3qlenYz8DmNMM7wmd4E8uPfFKF9s2Bu6fQ54sH+nPaVEFX2vkfFbQE
39053qTnK8quGRlXVJDD2dK5djNqMDRentxdaK/EYnTt9pU4GffsW+Ecz94Fy7d
391UYFwmp4wDbRfhNgBL2MV0BJMP6+JXh6RT0qUVQFuWZsdWujP9Ruqn7lV4N35
392I5OEM50WpgMWT0Oe43rFudj//7/1P+TdG7w=
393 "]]},
394 Annotation[#, "Charting`Private`Tag$5566#1"]& ]}, {}},
395 AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
396 Axes->{True, True},
397 AxesLabel->{None, None},
398 AxesOrigin->{2.4636339*^15, 0.0005432730994269633},
399 DisplayFunction->Identity,
400 Frame->{{False, False}, {False, False}},
401 FrameLabel->{{None, None}, {None, None}},
402 FrameTicks->{{Automatic,
403 Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic,
404 Charting`ScaledFrameTicks[{Identity, Identity}]}},
405 GridLines->{None, None},
406 GridLinesStyle->Directive[
407 GrayLevel[0.5, 0.4]],
408 ImagePadding->All,
409 Method->{
410 "DefaultBoundaryStyle" -> Automatic,
411 "DefaultGraphicsInteraction" -> {
412 "Version" -> 1.2, "TrackMousePosition" -> {True, False},
413 "Effects" -> {
414 "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2},
415 "Droplines" -> {
416 "freeformCursorMode" -> True,
417 "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" ->
418 AbsolutePointSize[6], "ScalingFunctions" -> None,
419 "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
420 (Identity[#]& )[
421 Part[#, 1]],
422 (Identity[#]& )[
423 Part[#, 2]]}& ), "CopiedValueFunction" -> ({
424 (Identity[#]& )[
425 Part[#, 1]],
426 (Identity[#]& )[
427 Part[#, 2]]}& )}},
428 PlotRange->{{2.4636339*^15, 2.4685661*^15}, {0.0005432730994269633,
429 0.1259025834119197}},
430 PlotRangeClipping->True,
431 PlotRangePadding->{{
432 Scaled[0.02],
433 Scaled[0.02]}, {
434 Scaled[0.05],
435 Scaled[0.05]}},
436 Ticks->{Automatic, Automatic}]], "Output",
437 CellChangeTimes->{
438 3.796138495756485*^9, {3.796138960221661*^9, 3.7961389756057653`*^9},
439 3.796139008292447*^9, {3.796139048401775*^9, 3.7961390572076674`*^9}, {
440 3.7961391715663977`*^9, 3.796139198209818*^9}, {3.7961392717440434`*^9,
441 3.7961392893209352`*^9}, {3.7961410651450157`*^9,
442 3.7961411231684437`*^9}, {3.79614116237171*^9, 3.796141195850403*^9},
443 3.7961414290095377`*^9, {3.7961414929333878`*^9, 3.7961415031613092`*^9},
444 3.796141644337055*^9, {3.796141852490193*^9, 3.796141876389929*^9},
445 3.7963810011419907`*^9, 3.796381220260585*^9},
446 CellLabel->"Out[63]=",ExpressionUUID->"77f4fc56-1f2f-4e33-bd0b-0f8815fb5d4f"]
447}, Open ]],
448
449Cell[CellGroupData[{
450
451Cell[BoxData[{
452 RowBox[{
453 RowBox[{"\[Nu]", "=", "\[Nu]0"}],
454 ";"}], "\[IndentingNewLine]", "\[Tau]", "\[IndentingNewLine]",
455 RowBox[{
456 RowBox[{"\[Nu]", "=",
457 RowBox[{"2.467", "*",
458 SuperscriptBox["10", "15"]}]}],
459 ";"}], "\[IndentingNewLine]", "\[Tau]", "\[IndentingNewLine]",
460 RowBox[{"Clear", "[", "\[Nu]", "]"}]}], "Input",
461 CellChangeTimes->{{3.796380978091817*^9, 3.796381014063077*^9}, {
462 3.796381075279804*^9, 3.7963810973865337`*^9}, 3.79638124976886*^9},
463 CellLabel->"In[69]:=",ExpressionUUID->"c77b2858-169c-429d-8b5e-3cab4c0f20a3"],
464
465Cell[BoxData["0.12590259965726647`"], "Output",
466 CellChangeTimes->{{3.796381002870018*^9, 3.796381015329191*^9}, {
467 3.7963810902155867`*^9, 3.7963810989719477`*^9}, {3.7963812242970457`*^9,
468 3.7963812535133543`*^9}},
469 CellLabel->"Out[70]=",ExpressionUUID->"6fc9b0d8-dd45-4b9e-8d36-905f69c2fcbb"],
470
471Cell[BoxData["0.0609598518343331`"], "Output",
472 CellChangeTimes->{{3.796381002870018*^9, 3.796381015329191*^9}, {
473 3.7963810902155867`*^9, 3.7963810989719477`*^9}, {3.7963812242970457`*^9,
474 3.796381253515643*^9}},
475 CellLabel->"Out[72]=",ExpressionUUID->"8db0efbc-bbda-4025-b366-151fb55f695c"]
476}, Open ]]
477},
478WindowSize->{808, 804},
479WindowMargins->{{506, Automatic}, {Automatic, 30}},
480FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
481StyleDefinitions->"Default.nb"
482]
483(* End of Notebook Content *)
484
485(* Internal cache information *)
486(*CellTagsOutline
487CellTagsIndex->{}
488*)
489(*CellTagsIndex
490CellTagsIndex->{}
491*)
492(*NotebookFileOutline
493Notebook[{
494Cell[558, 20, 1101, 28, 130, "Input",ExpressionUUID->"253f762e-8af6-4b50-b652-0f3b4a07e0c9"],
495Cell[CellGroupData[{
496Cell[1684, 52, 1593, 42, 244, "Input",ExpressionUUID->"de1a34d1-f494-4416-a345-371b7a127382"],
497Cell[3280, 96, 387, 5, 35, "Output",ExpressionUUID->"14e26464-572f-47dc-b947-58dc631736c6"]
498}, Open ]],
499Cell[CellGroupData[{
500Cell[3704, 106, 928, 22, 55, "Input",ExpressionUUID->"76bb0e2c-dc81-413a-bed5-2921b26ee82b"],
501Cell[4635, 130, 549, 14, 57, "Output",ExpressionUUID->"0994e65f-46ee-4409-8255-4333b6653546"]
502}, Open ]],
503Cell[CellGroupData[{
504Cell[5221, 149, 1507, 38, 181, "Input",ExpressionUUID->"d01c3e8e-c5ee-41e5-a4e3-5e218bfbf95a"],
505Cell[6731, 189, 330, 5, 37, "Output",ExpressionUUID->"8fd7ff38-9a2a-419e-a87c-b8cff9652dbd"],
506Cell[7064, 196, 589, 14, 44, "Output",ExpressionUUID->"e014cc80-ca01-41ee-943c-882968fdb6b8"]
507}, Open ]],
508Cell[CellGroupData[{
509Cell[7690, 215, 1235, 27, 80, "Input",ExpressionUUID->"7f83b4b0-1af8-44d0-9571-c30cdb9bb101"],
510Cell[8928, 244, 700, 10, 37, "Output",ExpressionUUID->"4020ff30-3750-435c-9995-c6228e6e1901"],
511Cell[9631, 256, 10403, 189, 240, "Output",ExpressionUUID->"77f4fc56-1f2f-4e33-bd0b-0f8815fb5d4f"]
512}, Open ]],
513Cell[CellGroupData[{
514Cell[20071, 450, 559, 12, 125, "Input",ExpressionUUID->"c77b2858-169c-429d-8b5e-3cab4c0f20a3"],
515Cell[20633, 464, 298, 4, 35, "Output",ExpressionUUID->"6fc9b0d8-dd45-4b9e-8d36-905f69c2fcbb"],
516Cell[20934, 470, 295, 4, 35, "Output",ExpressionUUID->"8db0efbc-bbda-4025-b366-151fb55f695c"]
517}, Open ]]
518}
519]
520*)
521