u/adebrech/Matlab/IonizationFront: ionization front.nb

File ionization front.nb, 18.9 KB (added by adebrech, 7 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 11.1' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 158, 7]
13NotebookDataLength[ 19149, 431]
14NotebookOptionsPosition[ 18273, 403]
15NotebookOutlinePosition[ 18607, 418]
16CellTagsIndexPosition[ 18564, 415]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell[BoxData[{
24 RowBox[{
25 RowBox[{"\[Sigma]H", "=",
26 RowBox[{"6.3", "*",
27 SuperscriptBox["10",
28 RowBox[{"-", "18"}]]}]}], ";"}], "\[IndentingNewLine]",
29 RowBox[{
30 RowBox[{"T", "=",
31 SuperscriptBox["10", "4"]}], ";"}], "\[IndentingNewLine]",
32 RowBox[{
33 RowBox[{"\[Alpha]Rec", "=",
34 RowBox[{"2.59", "*",
35 SuperscriptBox["10",
36 RowBox[{"-", "13"}]],
37 SuperscriptBox[
38 RowBox[{"(",
39 FractionBox["T",
40 SuperscriptBox["10", "4"]], ")"}],
41 RowBox[{"-", "0.7"}]]}]}], ";"}], "\[IndentingNewLine]",
42 RowBox[{
43 RowBox[{"nH", "=",
44 RowBox[{"5", "*",
45 SuperscriptBox["10", "7"]}]}], ";"}], "\[IndentingNewLine]",
46 RowBox[{
47 RowBox[{"Rp", "=",
48 RowBox[{"1.5", "*",
49 SuperscriptBox["10", "10"]}]}], ";"}], "\[IndentingNewLine]",
50 RowBox[{
51 RowBox[{"cells", "=", "256"}], ";"}], "\[IndentingNewLine]",
52 RowBox[{
53 RowBox[{"boxSize", "=", "2"}], ";"}], "\[IndentingNewLine]",
54 RowBox[{
55 RowBox[{"dx", "=",
56 FractionBox[
57 RowBox[{"Rp", " ", "boxSize"}], "cells"]}],
58 ";"}], "\[IndentingNewLine]",
59 RowBox[{"\[Lambda]", "=",
60 FractionBox["1",
61 RowBox[{"\[Sigma]H", " ", "0.5", " ", "nH"}]]}]}], "Input",
62 CellChangeTimes->{{3.713258846511506*^9, 3.713258939851328*^9}, {
63 3.7132590823904333`*^9, 3.713259126478526*^9}, {3.713259175948367*^9,
64 3.713259192038424*^9}, {3.713259223154564*^9, 3.713259226710137*^9}, {
65 3.713259283745267*^9, 3.7132593106941233`*^9}, {3.713520754920622*^9,
66 3.713520757575594*^9}, {3.7135208004449463`*^9, 3.713520866144994*^9}, {
67 3.713520953852779*^9, 3.7135209655016537`*^9}, {3.713521867250255*^9,
68 3.713521868976025*^9}, {3.71352190023104*^9, 3.713521900973702*^9}, {
69 3.713522196478245*^9, 3.713522248221025*^9}, {3.7135226077012987`*^9,
70 3.7135226248555613`*^9}, {3.7135235310986347`*^9, 3.713523559665166*^9}, {
71 3.713615204688394*^9, 3.713615247162218*^9}, {3.713615325930789*^9,
72 3.7136153583084517`*^9}, {3.713615402955405*^9, 3.713615408176056*^9}, {
73 3.713615474619404*^9, 3.713615553893057*^9}, {3.7136324946782503`*^9,
74 3.7136325289198837`*^9}, {3.7136335496534*^9, 3.713633553984888*^9}, {
75 3.713633602255455*^9, 3.713633602374426*^9}, {3.713633788040331*^9,
76 3.713633788534278*^9}, {3.713693377788597*^9, 3.713693418352685*^9}, {
77 3.7136934974407682`*^9, 3.713693502440379*^9}, {3.71369585042238*^9,
78 3.7136958621256723`*^9}, {3.7136962372383823`*^9, 3.71369625172823*^9}, {
79 3.713699909675071*^9, 3.713699911321764*^9}, {3.713705145967136*^9,
80 3.713705146765848*^9}, {3.713705272059846*^9, 3.713705389130068*^9}, {
81 3.7137054560770063`*^9, 3.7137054590863457`*^9}, {3.713707127878566*^9,
82 3.713707195361743*^9}, 3.713707239591943*^9, 3.713707647051476*^9, {
83 3.713707678231587*^9, 3.713707769699853*^9}, {3.7137079350929337`*^9,
84 3.713707990644376*^9}, {3.713708077620658*^9, 3.7137082291922073`*^9}, {
85 3.713708262381814*^9, 3.713708397505809*^9}, {3.7137086769693527`*^9,
86 3.713708678181336*^9}, {3.713708731659623*^9, 3.7137087637420473`*^9}, {
87 3.7137088432272243`*^9, 3.713708846434181*^9}, {3.713708880202011*^9,
88 3.713708903074548*^9}, {3.713716330266464*^9, 3.713716351175259*^9}, {
89 3.713716434540206*^9, 3.713716525606535*^9}, {3.713716611280066*^9,
90 3.713716646483315*^9}, {3.71371670304727*^9, 3.7137167087341948`*^9}, {
91 3.713716769507305*^9, 3.7137169073458242`*^9}, {3.7137221741594687`*^9,
92 3.713722287147996*^9}, {3.713722317414308*^9, 3.71372256217129*^9}, {
93 3.7137226062439938`*^9, 3.713722633168475*^9}, {3.713722962911133*^9,
94 3.713722973892164*^9}, {3.713723004049493*^9, 3.713723007056541*^9}, {
95 3.7137230748638363`*^9, 3.713723204412738*^9}, {3.713723352532102*^9,
96 3.713723352853589*^9}, {3.7137234296575212`*^9, 3.713723433512855*^9}, {
97 3.713723670492765*^9, 3.713723686101622*^9}, {3.7137237981056843`*^9,
98 3.713723799313353*^9}, {3.713724456612514*^9, 3.713724459841893*^9}, {
99 3.713785021813044*^9, 3.7137850320829763`*^9}, {3.713785506067985*^9,
100 3.713785533348872*^9}, {3.7137856004003353`*^9, 3.7137856533672657`*^9}, {
101 3.713785685385323*^9, 3.7137857186268377`*^9}, {3.71378582659606*^9,
102 3.713786076938623*^9}, {3.7137861087900457`*^9, 3.713786217248049*^9}, {
103 3.713786289465934*^9, 3.7137863597071877`*^9}, {3.713786442650681*^9,
104 3.713786485443578*^9}, 3.713786723240041*^9, {3.713791026836248*^9,
105 3.713791051116667*^9}, {3.7137916305626907`*^9, 3.713791661098472*^9}, {
106 3.713875993628359*^9, 3.713875996727685*^9}, {3.71387605653726*^9,
107 3.71387606768277*^9}, {3.7138762462463713`*^9, 3.713876297245927*^9}, {
108 3.7138763645880938`*^9,
109 3.713876373787952*^9}},ExpressionUUID->"69204edd-9257-4061-bd2c-\
110720bceada484"],
111
112Cell[BoxData["6.349206349206349`*^9"], "Output",
113 CellChangeTimes->{
114 3.713716380335881*^9, {3.713716442147944*^9, 3.713716525851348*^9},
115 3.71371664825263*^9, 3.7137167092665653`*^9, {3.713716786166074*^9,
116 3.713716907857258*^9}, {3.713722191757473*^9, 3.713722234293192*^9}, {
117 3.71372235014289*^9, 3.713722550900547*^9}, {3.7137226070667143`*^9,
118 3.713722633543158*^9}, 3.713722974680765*^9, 3.7137230073468323`*^9, {
119 3.713723078737987*^9, 3.713723205032489*^9}, {3.713723353632893*^9,
120 3.7137233663510733`*^9}, 3.713723466117373*^9, 3.713723544149767*^9,
121 3.713723802685974*^9, 3.713724460245397*^9, 3.7137829555583267`*^9,
122 3.713785033311545*^9, {3.713785574613027*^9, 3.713785600769039*^9}, {
123 3.713785635971496*^9, 3.71378565548559*^9}, 3.71378573699855*^9, {
124 3.713785829361939*^9, 3.713785937477783*^9}, 3.713785978797179*^9, {
125 3.713786037826242*^9, 3.713786059755438*^9}, {3.713786111342845*^9,
126 3.713786217646771*^9}, {3.7137862906383944`*^9, 3.713786325320175*^9},
127 3.713786364492869*^9, 3.713786496119326*^9, 3.713786548948299*^9,
128 3.7137866507004128`*^9, 3.713786747164752*^9, 3.713791662913753*^9,
129 3.713875939481181*^9, {3.713875994409239*^9, 3.713875998971075*^9},
130 3.713876070842866*^9, 3.713876312937078*^9,
131 3.7138763802523317`*^9},ExpressionUUID->"5f72263f-18ae-470d-b99d-\
13270d7cdbdaacd"]
133}, Open ]],
134
135Cell[BoxData[{
136 RowBox[{"Clear", "[",
137 RowBox[{"\[Tau]loc", ",", "Xloc"}], "]"}], "\[IndentingNewLine]",
138 RowBox[{
139 RowBox[{"photonFlux", "=",
140 RowBox[{"2", "*",
141 SuperscriptBox["10", "12"]}]}], ";"}], "\[IndentingNewLine]",
142 RowBox[{
143 RowBox[{"\[Tau]loc", "=",
144 RowBox[{"Table", "[",
145 RowBox[{"0", ",",
146 RowBox[{"{",
147 RowBox[{"i", ",",
148 RowBox[{"cells", "+", "1"}]}], "}"}]}], "]"}]}],
149 ";"}], "\[IndentingNewLine]",
150 RowBox[{
151 RowBox[{"Xloc", "=",
152 RowBox[{"Table", "[",
153 RowBox[{"0", ",",
154 RowBox[{"{",
155 RowBox[{"i", ",",
156 RowBox[{"cells", "+", "1"}]}], "}"}]}], "]"}]}],
157 ";"}], "\[IndentingNewLine]",
158 RowBox[{
159 RowBox[{
160 RowBox[{"RecombinationRate", "[", "c_", "]"}], ":=",
161 RowBox[{
162 SuperscriptBox[
163 RowBox[{"(",
164 RowBox[{"curX", " ", "nH"}], ")"}], "2"], "\[Alpha]Rec"}]}],
165 ";"}], "\[IndentingNewLine]",
166 RowBox[{
167 RowBox[{
168 RowBox[{"IonizationRate", "[", "c_", "]"}], ":=",
169 RowBox[{"photonFlux", " ",
170 SuperscriptBox["\[ExponentialE]",
171 RowBox[{"-",
172 RowBox[{"\[Tau]loc", "[",
173 RowBox[{"[", "c", "]"}], "]"}]}]], "\[Sigma]H", " ",
174 RowBox[{"(",
175 RowBox[{"1", "-", "curX"}], ")"}], "nH"}]}],
176 ";"}], "\[IndentingNewLine]",
177 RowBox[{"Do", "[",
178 RowBox[{
179 RowBox[{
180 RowBox[{"soln", "=",
181 RowBox[{"Solve", "[",
182 RowBox[{
183 RowBox[{
184 RowBox[{
185 RowBox[{"RecombinationRate", "[", "cell", "]"}], "\[Equal]",
186 RowBox[{"IonizationRate", "[", "cell", "]"}]}], "&&",
187 RowBox[{"0", "\[LessEqual]", "curX", "\[LessEqual]", "1"}]}], ",",
188 "curX"}], "]"}]}], ";", "\[IndentingNewLine]",
189 RowBox[{
190 RowBox[{"Xloc", "[",
191 RowBox[{"[", "cell", "]"}], "]"}], "=",
192 RowBox[{"curX", "/.",
193 RowBox[{"soln", "[",
194 RowBox[{"[", "1", "]"}], "]"}]}]}], ";", "\[IndentingNewLine]",
195 RowBox[{
196 RowBox[{"\[Tau]loc", "[",
197 RowBox[{"[",
198 RowBox[{"cell", "+", "1"}], "]"}], "]"}], "=",
199 RowBox[{
200 RowBox[{"\[Tau]loc", "[",
201 RowBox[{"[", "cell", "]"}], "]"}], "+",
202 RowBox[{"\[Sigma]H", " ", "nH", " ", "dx", " ",
203 RowBox[{"(",
204 RowBox[{"1", "-",
205 RowBox[{"Xloc", "[",
206 RowBox[{"[", "cell", "]"}], "]"}]}], ")"}]}]}]}]}], ",",
207 RowBox[{"{",
208 RowBox[{"cell", ",", "1", ",", "cells"}], "}"}]}], "]"}]}], "Input",
209 CellChangeTimes->{{3.713787181483942*^9, 3.7137871831408443`*^9}, {
210 3.7137872424465313`*^9, 3.713787266310635*^9}, {3.713787331346404*^9,
211 3.713787340786622*^9}, {3.713787450698647*^9, 3.713787459236697*^9}, {
212 3.7137888251065187`*^9, 3.713789092606905*^9}, {3.713789123743115*^9,
213 3.713789170454896*^9}, {3.713789440543017*^9, 3.713789464049374*^9}, {
214 3.713789496027947*^9, 3.713789507638431*^9}, {3.713789572442625*^9,
215 3.713789795087261*^9}, {3.7137898669191923`*^9, 3.71378991261971*^9}, {
216 3.713789946695181*^9, 3.713790021766535*^9}, {3.7137902746899014`*^9,
217 3.713790408204462*^9}, 3.71379044152044*^9, {3.713790532146649*^9,
218 3.713790662372389*^9}, {3.713790726308651*^9, 3.7137908149317417`*^9}, {
219 3.7138764604662533`*^9,
220 3.713876463710588*^9}},ExpressionUUID->"926378f1-32fb-44a4-bc18-\
22146ebeef2ba87"],
222
223Cell[CellGroupData[{
224
225Cell[BoxData[{
226 RowBox[{
227 RowBox[{
228 RowBox[{"X", "[",
229 RowBox[{"loc_", ",", "midpt_", ",", "k_"}], "]"}], ":=",
230 RowBox[{"1", "-",
231 FractionBox["1",
232 RowBox[{
233 SuperscriptBox["\[ExponentialE]",
234 RowBox[{
235 RowBox[{"-", "k"}],
236 RowBox[{"(",
237 RowBox[{"loc", "-", " ", "midpt"}], ")"}]}]], "+", "1"}]]}]}],
238 ";"}], "\[IndentingNewLine]",
239 RowBox[{
240 RowBox[{"plot1", "=",
241 RowBox[{"ListPlot", "[",
242 RowBox[{"Xloc", ",",
243 RowBox[{"PlotRange", "\[Rule]",
244 RowBox[{"{",
245 RowBox[{"0", ",", "1"}], "}"}]}]}], "]"}]}],
246 ";"}], "\[IndentingNewLine]",
247 RowBox[{
248 RowBox[{"plot2", "=",
249 RowBox[{"Plot", "[",
250 RowBox[{
251 RowBox[{"X", "[",
252 RowBox[{"cell", ",", "48", ",",
253 RowBox[{"1.55", "*",
254 SuperscriptBox["10",
255 RowBox[{"-", "2"}]]}]}], "]"}], ",",
256 RowBox[{"{",
257 RowBox[{"cell", ",", "1", ",", "256"}], "}"}]}], "]"}]}],
258 ";"}], "\[IndentingNewLine]",
259 RowBox[{"Show", "[",
260 RowBox[{"plot1", ",", "plot2"}], "]"}]}], "Input",
261 CellChangeTimes->{{3.7137908742968473`*^9, 3.713790902172386*^9}, {
262 3.713790934163289*^9, 3.713790947085926*^9}, {3.713791058259803*^9,
263 3.713791284890162*^9}, {3.713791328952856*^9, 3.713791498411188*^9}, {
264 3.713791531479056*^9, 3.7137915320176563`*^9}, {3.7137915705756083`*^9,
265 3.713791607462647*^9}, {3.713791746837451*^9, 3.713791757442541*^9}, {
266 3.713792603606801*^9,
267 3.7137926048563023`*^9}},ExpressionUUID->"3fcfb279-8e67-40ab-b069-\
268a6483b65a7ec"],
269
270Cell[BoxData[
271 GraphicsBox[{{{}, {{},
272 {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
273 AbsoluteThickness[1.6], PointBox[CompressedData["
2741:eJw11mtcjHkfx/FL4d5ti1g6ONxNB0mUVELE16lS0kmJTlPT1FQzE3btbgoX
275VtvBoSKrYk2WlK1WOVQUV2srkijruMmgJbE5F0lu98vvfz2Z13u+85vPg3ky
276xuFxPlINjuMGDOC4/79+fp7P8X95/G6G6sEcegPX59W5Z2xl/gI50d5HFT8w
27762Jvi4fO+AjmEdg9vCOi2pPZAPn9O07bzGQeA7f6wQYbzZlFsP5asblkGLMJ
278bo70/rKy7z7ZDIt/qSpXPWI2h9xw/n5pC7MFmjYPvfqfamZLFDm/lP10mHkS
279JrxXxrdnMFsjYGiDkVEisw0Gu+qtnRbJPAX+o+YVWnkx26LzhlnzAEdmO6wq
280SDM4bspsDy1358IFOsxToZtjd/RY9z2yA7Ze9AgeoGaehl9kJx9YXWCejkdv
281auJnlDHPgIGL5nKzXGZHfPu9/5muzcwz4brxyeM9cuZZ4JReRqZ+zE44Y5Bw
282aJsT82x0J7bU3RzHPAetqdOaBw5hBmwTdPSGd6s/mweG7J3bybWRublwTvfJ
283b65l+1x89d4//8dits/DdxmabqOy2D4POWlHX2cmsn0+nitDdJ5J2D4fTyee
2846rd2Z/sC7E0WTfO1ZfsChLSoLIIM2b4QDVf6Xd04ti+EY0eU6aymu7Q7Y9A3
285dz/+nkvmnfFwg/WEFzK2u+CIbnIB58B2F7wo37z9rgbbXTGuUP/V9stttLti
286m0bGm6F7ydwitLimVklkbF+EYdrqNWn2bHeD7u+/+idzbHdD6yLt3MDGO7S7
287ozRPmjR4D5l3x22bD+GpErYvRsMMnU0d1mxfjJ/vHZpp0ttKuweOLdnVNb2W
288zHsg54ILZ53O9iXQaDveqhHI9iXoTnDorhjHdk9MNlTUej//m3ZP3L86+kLT
289KTLnhenZmklTtpDhhR98qpISPNnnvTBi7NbAEkOy4IXMzFk5Fx/cpntvxJh3
2903fmrmAxvtKse1138nsx7Y+DZuN7SuWTBGxncTyYpX7F7H6Ck/l+fa7fo3gf/
291NTUfOmw/mfeBo/S0U62MLPgga9RW+1W2ZM4X4g6PMr2+m3TvC6X7xPjyWjLv
292i4p3W5b47SALvmi5U977IoDMLcWl6hbXHSbsfikcVCGDJz+9QfdLodTQ17x6
293giwsxdrX+brrNpA5P9hMjtCYvIgMP7g6flncOZzd+yG9QvmwtPU63fvB9MX8
294LUn5ZM4fvUV14piVZPjjOr/WK9SRzPujLffuxKiB7N4fjXYTGjY3XaP7ZSg7
295/GpA+R4yluHd6t4CDQmZXwaRVtNOmRVZWAab0ylZnT1/0X0Abv+WsiH1DzIC
296MG+osbHbNjIfgORiL5eJAWQhAN3X11VNNmX3y3G0dYLr8q6rdL8c31kcKSms
297JPPL0ap976jJFrKwHF9E9XxV50XmVmCnpX581hh2vwLnn5kXZ3a00P0KfMhN
298Cjh7nCysgNXMS+ZjN5K5QEQaJ58p9iAjEA+j3wurR7H7QJg/fnBB8aiZ7gPh
299mO8YpzpO5oJg+OwdtDeREYTiUe3nyzzJfBAStRanZo1l90G4l3h3wMknV+g+
300GP/M+ebQyFNkBKPMt/nBiWQyH4w8rRSzPcvIQjDa7Grqq83ZfQhGVxos9np6
301me5DII9YLypNI/MhWCPK/EU1kSyEoH/jWFfjxia6D4X89fbtUxVkhGLX1IqW
302u0PIfChmh2ZEmJReovtQNCRxmz/6kjkxet8eCozvafxskRg9I4u1MnPJEMPu
303cl6OO8hiMS7VuQ0s+Ocifb8YoiHJ1oVpZJUYxesdurxsyYIYBof1B+fdavhs
304tRiK+WfMczeSuTBsdNjQBUuyKAxXeiIf77x6gfphMJPZZWStI4vDcLPP28PF
305gsyHYbfN/dVHrp6nfhiczasWCBvIQhieRye8SplEVochpHLbvi9u11M/HDHV
306Cz2dksmicKzIfaExYRoZ4ShY8EbV+LCO+uGIjOn4aPEzmQ9Hwtvy9tmuZFU4
307dPsaNIb21lI/HEU3LPuyi8jqcJzKLwlTh5I5CXJKJxy7/zVZJMGfLXtD887/
308SX0JDIpjasesJ4slaMvcbutnT+Yl6O+71u725Bz1JYhcNUar/wBZkGB9eYqg
309DCSrJTg3o0Vv3wgyF4G+1fHaKU1/UD8C7s56mbYpZESgSZSdkr2ALI6A7Z87
31066o5Mh+B1OxNQ/ZX11A/Ap5yfUOnBLIQgZWKLek/zyCrI+AU2jSs5K1AfSm6
311G3St+AqySIoS4xtSnXgypOCPRNn5OZLFUlTNtdQoeXuW+lKENxoFm60nq6Ro
31284sTaQ8iC1LsenS7TbLtDPWlaBUujrTUJ3ORcLTvsQo+UE39SLzL00x8b01G
313JOzFVnM1q6uoHwmEPflt9WIyH4mc2J03Xe6cpn4krpR/aObjyEIk2hvGZ47W
314JKs/3Uuva+jtOUX9KLSGOvQrrMmiKMh8Ez2N6yqpH4WMh5uOWYWSxVFIPRB1
315K/1dBfWjcNDaXOKaRVZFYfa1KQ+9bclCFNIcdJ4euVJO/Sg0DFG/9FpJ5mQ4
3160KOlnDuMLJKhRKeyLv7YSerLYCGNC+rxJ4tl6JiRs6X8/Qnqy5B/Xj+9PI+s
317kkFUI9S8diULsk+/V1eA/MVx6sswwfPD1tG5ZC4aLacvb/+4kCyKxvBv3VJG
318vjxG/Wg0jmvOW7GfLI5Gh+PZtxc9yHw0Oo3aCiM/lFE/GplXbjaOLyEL0eg/
319VLRGW0xWRyM/XP/EsK/JXAyOx2tk29eXUj8GSH8yaXUiGTG4Mc1tbb0tWRyD
320O/dfJk/pPEr9GHAjvg8tPkBWxeB2wybN6UFkIQaPbaatbdIjq2NQNGJQxt8H
321f6d+LNa3dT4c91+yKBYtBy0HfJddQv1YlJ0MeFqrRxbHotvOqnzo7mLqx+LR
322VMNvPPTJqk97/XX7+Jwi6sdipfBqxDYjsjoWvSnLbDbn/0Z9OXKenSvwm0wW
323yZH8wetQf+UR6suhpdXu84MzWSxH77u89qq/Cqkvx6ZrQbvqIsgqOYpq+g7s
3246C6gvhzLvm0PGJ5KVsthMt178CIjMqfAut0nNa1PHqa+Am5nRpadXkKGAqeq
325vRKeduRTX4FZ+fEtp34k8wpUlI61NzYhqxQYxG81H1VziPoKWPw42lEVRlYr
326Pv0fk35ZoknmlCgu/WnG7MMHqa8EzAPWui8mQ4mqgfuSrrz8lfpK3J7kZHY+
327h8wrkdSuNLNaQFYpYWn0QNLbdYD6SliU/NtgmUtWKzFm1RGPSlcyF4fIvUZP
328VT151I/D1EnBZa2HyYgDl5NdEbWcLI7DrE53UydtMh+HivHzTQMEFfXjcC7s
3299Jtja8hCHLTnrWx3n0RWx0G+b0Waw/n91F+J1lvlbx0dyborwdHzP6JT88M=
330
331 "]]}, {}}, {}, {}, {{}, {}}, {{}, {}}}, {{{}, {},
332 TagBox[
333 {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
334 Opacity[1.], LineBox[CompressedData["
3351:eJwVkHs01AkUx8cwdUyvSZPsUWNYWlnJo4cwc68isjESUUnGs5CfaFPUZqV2
336lR2TV7UxW34ztNgdUVhN08OzYlUoG1nUJJNW0/Eor9npj3vuued+7/1+7jUO
337IXzCqRQKxVMTX/KW+rHDNMoHiBXvqvTKVUBGIX1XIKpgTwY3mp2tgPYUI05l
338ykdQatW9UQkVEMzZSuNrjcEPzBHznAwFHKvKz5VRJ+EU/2Z35ykF7Oww7v6L
339SsXahSvpXvEKONvYPBVHm48K5muaHk+j9zsnImgGGBKxMH54rgLU1Sr5bM9K
340dIo0mqmKeg1n5Flnc05aoyBs7ebcjlfQHkIzs09chy2+nPlXua/gWs5kEfnC
341Ae8eNyjnXRqAkEK7HsoAYHqyfuXkdD+0ixKns0o2YYpxiW339n4IaP3xfoy/
342K1Ld85TSa33QNrQ6BJ654fSKuIbd9D6ge14eGzvrgRf1G4I9bvfCBlFJiGSx
343J+4ncjzLT7wEd++IufnnePj2D4LtsqUHBGIV79aK7Uj/YMP4eeIFfFympSsJ
3448sGFSzbPfFX3D1RbqW6l1O5ANvuTnm5SF/z3fTWdwvBDi5Ql/s6c58C64dWY
3457LUTb0xmv7ZZ/gzKcxuTh4r90cCE9l2YsgOcYheJdigDsN46P6+9ox0slPa1
346UrfdGFXIfMa5/hTS8+zjpaI9iOZHF9TnP4GBjboXSnsDcdg3oXWCfAwFKbU3
347L20MwpxHSP2MbbBd7uScdmIfmjre0DOJawXHPI+16XuDUccvXmx86xGoWYMu
348D3T5GCC+d4Y1/gDw6+FzP5XwMe1l+GGT0GawPdbfJg8IQdM13Qcb6hph7+C9
349/ibtUJRdHjcJcm4AnSd0Cx4ZipYHDNdVvauDwPkuBM87DM27WJG2v9wHmzqF
350slUVhkKRjoFbwD2ArPi1LGE46hBSCsPhLsxWjCQNOkZgzfGtleMWcnDgdU5Y
3519UTg409K0/RVMljJKT00kh6JNXP2pUyxaiGzaMjYynI/No93u1z1qwEb8yVX
352uI/2Yx8tS/VbeRWkTVkaNiceQLZdB+8b2k3w9ZC2tBlEITt0aYWTsBI6lcbJ
3533k1R+Pbd+pNq/QoQkCcnCuOjke0lmPCVl0P9c9eSBL0Y9D79Xv7eVAqle6d3
3543b4Tg4yIP1mGTWWwaJ8hnuYfRCGvcsHTmhIoO2o2yqLEIjrK3dKeX4PPr4wk
355NmWaWmjNlauKoO5hasb9bQQGa5dln9omAeUys9RMHoGUVX9/OOgqAUZ4U2Kg
356D4GMambmbq4EgtTzQsf9CbyidSTDwVoCk3a5Gy1CNfOz/KNLmRKwLih+c/6Y
357ps9dnra+RwwFRAvyiwgULk8YZseKoV4Wu97qdwL7RqVWPpFiUOoutpwqJdD7
358W0HB6WAx2Iv9luVe1+ifFs9+9BFDR1fv+2YZgXf7imWDG8RA36T61bqdQLaM
359cdmeKgabzGzhTCeBGN23J2OKhICedWcedhEYF7MFBkZJkBxJOhTWS2D5Bf/o
360nEESsFTb/eKQxq/PzE2vlYTITxJO+LBGf4g7L7GRBIGru53tiIYvcmz45R0S
361ev7NYLWMavZzG+ZUVJCgvXoN89KE5r7OqxyjMhIskp7oRkxqeJzFmQIJCd5N
362CWrbmS//sppWi0hIZOqPqdUa/vNkasJFEv4HkBNGKg==
363 "]]},
364 Annotation[#, "Charting`Private`Tag$28794#1"]& ]}, {}, {}}},
365 AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
366 Axes->{True, True},
367 AxesLabel->{None, None},
368 AxesOrigin->{0., 0},
369 DisplayFunction->Identity,
370 Frame->{{False, False}, {False, False}},
371 FrameLabel->{{None, None}, {None, None}},
372 FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
373 GridLines->{None, None},
374 GridLinesStyle->Directive[
375 GrayLevel[0.5, 0.4]],
376 ImagePadding->All,
377 ImageSize->{913., Automatic},
378 Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
379 (Identity[#]& )[
380 Part[#, 1]],
381 (Identity[#]& )[
382 Part[#, 2]]}& ), "CopiedValueFunction" -> ({
383 (Identity[#]& )[
384 Part[#, 1]],
385 (Identity[#]& )[
386 Part[#, 2]]}& )}},
387 PlotRange->{{0., 257.}, {0, 1}},
388 PlotRangeClipping->True,
389 PlotRangePadding->{{
390 Scaled[0.02],
391 Scaled[0.02]}, {0, 0}},
392 Ticks->{Automatic, Automatic}]], "Output",
393 CellChangeTimes->{3.713792605186795*^9, 3.7138763190906982`*^9,
394 3.713876472036104*^9},ExpressionUUID->"803fca89-9bc0-4c05-bf94-\
3950977d9e14dda"]
396}, Open ]],
397
398Cell[BoxData[
399 RowBox[{" ", "nH"}]], "Input",
400 CellChangeTimes->{{3.713876593631577*^9,
401 3.713876594561634*^9}},ExpressionUUID->"eeab12d1-69f9-4c0b-b893-\
4021866424752e8"]
403},
404WindowSize->{960, 1028},
405WindowMargins->{{Automatic, 0}, {0, Automatic}},
406FrontEndVersion->"11.1 for Linux x86 (64-bit) (April 18, 2017)",
407StyleDefinitions->"Default.nb"
408]
409(* End of Notebook Content *)
410
411(* Internal cache information *)
412(*CellTagsOutline
413CellTagsIndex->{}
414*)
415(*CellTagsIndex
416CellTagsIndex->{}
417*)
418(*NotebookFileOutline
419Notebook[{
420Cell[CellGroupData[{
421Cell[580, 22, 4701, 87, 280, "Input", "ExpressionUUID" -> \
422"69204edd-9257-4061-bd2c-720bceada484"],
423Cell[5284, 111, 1358, 20, 34, "Output", "ExpressionUUID" -> \
424"5f72263f-18ae-470d-b99d-70d7cdbdaacd"]
425}, Open ]],
426Cell[6657, 134, 3244, 86, 221, "Input", "ExpressionUUID" -> \
427"926378f1-32fb-44a4-bc18-46ebeef2ba87"],
428Cell[CellGroupData[{
429Cell[9926, 224, 1531, 43, 131, "Input", "ExpressionUUID" -> \
430"3fcfb279-8e67-40ab-b069-a6483b65a7ec"],
431Cell[11460, 269, 6623, 125, 602, "Output", "ExpressionUUID" -> \
432"803fca89-9bc0-4c05-bf94-0977d9e14dda"]
433}, Open ]],
434Cell[18098, 397, 171, 4, 32, "Input", "ExpressionUUID" -> \
435"eeab12d1-69f9-4c0b-b893-1866424752e8"]
436}
437]
438*)
439