1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 10.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[       158,          7]
13NotebookDataLength[     51108,       1041]
14NotebookOptionsPosition[     48688,        955]
15NotebookOutlinePosition[     49045,        971]
16CellTagsIndexPosition[     49002,        968]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell[BoxData[
24 RowBox[{"SetDirectory", "[",
25  RowBox[{"NotebookDirectory", "[", "]"}], "]"}]], "Input",
26 CellChangeTimes->{{3.660861388735448*^9, 3.660861437338705*^9},
27   3.6788887777601137`*^9}],
28
29Cell[BoxData["\<\"/Users/sandrews/SSA/software/Smoldyn/trunk/examples/S7_\
30surfaces/surfacediffuse\"\>"], "Output",
31 CellChangeTimes->{
32  3.660861439140813*^9, 3.660864234844933*^9, 3.660864284584375*^9,
33   3.6608650350472393`*^9, 3.660865096205949*^9, 3.6608651398801737`*^9,
34   3.6608654954230947`*^9, 3.6608655277567387`*^9, 3.6611075827902327`*^9,
35   3.66110792705791*^9, 3.661108437278463*^9, 3.661108526271861*^9, {
36   3.661108635909975*^9, 3.6611086646319447`*^9}, {3.6611087260259323`*^9,
37   3.661108740145318*^9}, {3.661109730652629*^9, 3.6611097462013273`*^9}, {
38   3.66110980654559*^9, 3.6611098476378813`*^9}, 3.661110721470619*^9, {
39   3.661153760115095*^9, 3.6611537824119587`*^9}, {3.6611538422247953`*^9,
40   3.661153869715589*^9}, 3.661154420555513*^9, 3.6611544797574797`*^9, {
41   3.6611545316576242`*^9, 3.661154604459819*^9}, 3.661154808938043*^9, {
42   3.661173580925185*^9, 3.661173714180011*^9}, 3.661173804943112*^9,
43   3.661173845147468*^9, 3.661173881470805*^9, 3.6611739205264378`*^9,
44   3.661223969204103*^9, 3.6612240024286947`*^9, 3.661224090186607*^9, {
45   3.661224161116886*^9, 3.661224185907539*^9}, 3.6612242190617943`*^9,
46   3.661224425801344*^9, 3.661224780263822*^9, 3.6612250824917088`*^9,
47   3.678888778366929*^9}]
48}, Open  ]],
49
50Cell[CellGroupData[{
51
52Cell[BoxData[{
53 RowBox[{"timestep", "=", "10"}], "\[IndentingNewLine]",
54 RowBox[{"diffconst", "=", "10"}], "\[IndentingNewLine]",
55 RowBox[{"binlow", "=", "0"}], "\[IndentingNewLine]",
56 RowBox[{"binhigh", "=", "100"}], "\[IndentingNewLine]",
57 RowBox[{"binwidth", "=", "0.5"}]}], "Input",
58 CellChangeTimes->{{3.6611075990872917`*^9, 3.6611076381188173`*^9}, {
59  3.661107727171955*^9, 3.661107736299667*^9}}],
60
61Cell[BoxData["10"], "Output",
62 CellChangeTimes->{
63  3.6611076397678757`*^9, 3.661107739145013*^9, 3.661107927124093*^9,
64   3.661108437313519*^9, 3.661108526356292*^9, {3.661108635947103*^9,
65   3.661108664711038*^9}, {3.6611087261036158`*^9, 3.6611087402005453`*^9}, {
66   3.66110973073179*^9, 3.661109746284515*^9}, {3.661109806626606*^9,
67   3.661109847716003*^9}, 3.661110721535725*^9, {3.6611537601629553`*^9,
68   3.661153782467613*^9}, {3.661153842304234*^9, 3.661153869791154*^9},
69   3.66115442063766*^9, 3.661154479822852*^9, {3.661154531723963*^9,
70   3.661154604527279*^9}, 3.661154809021887*^9, {3.6611735810048313`*^9,
71   3.661173714259088*^9}, 3.661173805022574*^9, 3.66117384522873*^9,
72   3.661173881551258*^9, 3.6611739206080236`*^9, 3.661223969300562*^9,
73   3.661224002494129*^9, 3.661224090270707*^9, {3.661224161184846*^9,
74   3.6612241859748077`*^9}, 3.661224219144142*^9, 3.661224425885468*^9,
75   3.661224780331427*^9, 3.661225082559111*^9}],
76
77Cell[BoxData["10"], "Output",
78 CellChangeTimes->{
79  3.6611076397678757`*^9, 3.661107739145013*^9, 3.661107927124093*^9,
80   3.661108437313519*^9, 3.661108526356292*^9, {3.661108635947103*^9,
81   3.661108664711038*^9}, {3.6611087261036158`*^9, 3.6611087402005453`*^9}, {
82   3.66110973073179*^9, 3.661109746284515*^9}, {3.661109806626606*^9,
83   3.661109847716003*^9}, 3.661110721535725*^9, {3.6611537601629553`*^9,
84   3.661153782467613*^9}, {3.661153842304234*^9, 3.661153869791154*^9},
85   3.66115442063766*^9, 3.661154479822852*^9, {3.661154531723963*^9,
86   3.661154604527279*^9}, 3.661154809021887*^9, {3.6611735810048313`*^9,
87   3.661173714259088*^9}, 3.661173805022574*^9, 3.66117384522873*^9,
88   3.661173881551258*^9, 3.6611739206080236`*^9, 3.661223969300562*^9,
89   3.661224002494129*^9, 3.661224090270707*^9, {3.661224161184846*^9,
90   3.6612241859748077`*^9}, 3.661224219144142*^9, 3.661224425885468*^9,
91   3.661224780331427*^9, 3.661225082562871*^9}],
92
93Cell[BoxData["0"], "Output",
94 CellChangeTimes->{
95  3.6611076397678757`*^9, 3.661107739145013*^9, 3.661107927124093*^9,
96   3.661108437313519*^9, 3.661108526356292*^9, {3.661108635947103*^9,
97   3.661108664711038*^9}, {3.6611087261036158`*^9, 3.6611087402005453`*^9}, {
98   3.66110973073179*^9, 3.661109746284515*^9}, {3.661109806626606*^9,
99   3.661109847716003*^9}, 3.661110721535725*^9, {3.6611537601629553`*^9,
100   3.661153782467613*^9}, {3.661153842304234*^9, 3.661153869791154*^9},
101   3.66115442063766*^9, 3.661154479822852*^9, {3.661154531723963*^9,
102   3.661154604527279*^9}, 3.661154809021887*^9, {3.6611735810048313`*^9,
103   3.661173714259088*^9}, 3.661173805022574*^9, 3.66117384522873*^9,
104   3.661173881551258*^9, 3.6611739206080236`*^9, 3.661223969300562*^9,
105   3.661224002494129*^9, 3.661224090270707*^9, {3.661224161184846*^9,
106   3.6612241859748077`*^9}, 3.661224219144142*^9, 3.661224425885468*^9,
107   3.661224780331427*^9, 3.661225082565363*^9}],
108
109Cell[BoxData["100"], "Output",
110 CellChangeTimes->{
111  3.6611076397678757`*^9, 3.661107739145013*^9, 3.661107927124093*^9,
112   3.661108437313519*^9, 3.661108526356292*^9, {3.661108635947103*^9,
113   3.661108664711038*^9}, {3.6611087261036158`*^9, 3.6611087402005453`*^9}, {
114   3.66110973073179*^9, 3.661109746284515*^9}, {3.661109806626606*^9,
115   3.661109847716003*^9}, 3.661110721535725*^9, {3.6611537601629553`*^9,
116   3.661153782467613*^9}, {3.661153842304234*^9, 3.661153869791154*^9},
117   3.66115442063766*^9, 3.661154479822852*^9, {3.661154531723963*^9,
118   3.661154604527279*^9}, 3.661154809021887*^9, {3.6611735810048313`*^9,
119   3.661173714259088*^9}, 3.661173805022574*^9, 3.66117384522873*^9,
120   3.661173881551258*^9, 3.6611739206080236`*^9, 3.661223969300562*^9,
121   3.661224002494129*^9, 3.661224090270707*^9, {3.661224161184846*^9,
122   3.6612241859748077`*^9}, 3.661224219144142*^9, 3.661224425885468*^9,
123   3.661224780331427*^9, 3.661225082567481*^9}],
124
125Cell[BoxData["0.5`"], "Output",
126 CellChangeTimes->{
127  3.6611076397678757`*^9, 3.661107739145013*^9, 3.661107927124093*^9,
128   3.661108437313519*^9, 3.661108526356292*^9, {3.661108635947103*^9,
129   3.661108664711038*^9}, {3.6611087261036158`*^9, 3.6611087402005453`*^9}, {
130   3.66110973073179*^9, 3.661109746284515*^9}, {3.661109806626606*^9,
131   3.661109847716003*^9}, 3.661110721535725*^9, {3.6611537601629553`*^9,
132   3.661153782467613*^9}, {3.661153842304234*^9, 3.661153869791154*^9},
133   3.66115442063766*^9, 3.661154479822852*^9, {3.661154531723963*^9,
134   3.661154604527279*^9}, 3.661154809021887*^9, {3.6611735810048313`*^9,
135   3.661173714259088*^9}, 3.661173805022574*^9, 3.66117384522873*^9,
136   3.661173881551258*^9, 3.6611739206080236`*^9, 3.661223969300562*^9,
137   3.661224002494129*^9, 3.661224090270707*^9, {3.661224161184846*^9,
138   3.6612241859748077`*^9}, 3.661224219144142*^9, 3.661224425885468*^9,
139   3.661224780331427*^9, 3.661225082569577*^9}]
140}, Open  ]],
141
142Cell[CellGroupData[{
143
144Cell[BoxData[
145 RowBox[{"bendpos", "=", "30"}]], "Input",
146 CellChangeTimes->{{3.661110490137539*^9, 3.661110493057027*^9}}],
147
148Cell[BoxData["30"], "Output",
149 CellChangeTimes->{
150  3.661110498569437*^9, 3.661110721582816*^9, {3.6611537602081957`*^9,
151   3.661153782512244*^9}, {3.661153842354282*^9, 3.661153869841504*^9},
152   3.6611544206894007`*^9, 3.6611544798770933`*^9, {3.6611545317781553`*^9,
153   3.661154604581534*^9}, 3.66115480907236*^9, {3.661173581054943*^9,
154   3.6611737143081636`*^9}, 3.661173805073143*^9, 3.6611738452794857`*^9,
155   3.6611738816019783`*^9, 3.661173920656692*^9, 3.661223969360214*^9,
156   3.661224002551424*^9, 3.661224090325699*^9, {3.661224161240304*^9,
157   3.661224186031185*^9}, 3.661224219199957*^9, 3.661224425940599*^9,
158   3.661224780387721*^9, 3.661225082615237*^9}]
159}, Open  ]],
160
161Cell[BoxData[
162 RowBox[{"(*", " ",
163  RowBox[{
164  "Everything", " ", "below", " ", "here", " ", "should", " ", "take", " ",
165   "care", " ", "of", " ", "itself"}], " ", "*)"}]], "Input",
166 CellChangeTimes->{{3.661107643372525*^9, 3.661107653918783*^9}}],
167
168Cell[BoxData[
169 RowBox[{"(*", " ",
170  RowBox[{"Define", " ", "functions"}], " ", "*)"}]], "Input",
171 CellChangeTimes->{{3.661107671773728*^9, 3.661107680216037*^9}}],
172
173Cell[BoxData[
174 RowBox[{
175  RowBox[{"gauss", "[",
176   RowBox[{"x_", ",", "mu_", ",", "sigma_"}], "]"}], ":=",
177  RowBox[{
178   RowBox[{"1", "/",
179    RowBox[{"(",
180     RowBox[{"sigma", "*",
181      RowBox[{"Sqrt", "[",
182       RowBox[{"2", "*", "Pi"}], "]"}]}], ")"}]}], "*",
183   RowBox[{"Exp", "[",
184    RowBox[{
185     RowBox[{"-",
186      RowBox[{
187       RowBox[{"(",
188        RowBox[{"x", "-", "mu"}], ")"}], "^", "2"}]}], "/",
189     RowBox[{"(",
190      RowBox[{"2", "*",
191       RowBox[{"sigma", "^", "2"}]}], ")"}]}], "]"}]}]}]], "Input",
192 CellChangeTimes->{{3.660864315305353*^9, 3.660864402104426*^9}}],
193
194Cell[BoxData[
195 RowBox[{"(*", " ",
196  RowBox[{"Load", " ", "and", " ", "look", " ", "at", " ", "data"}], " ",
197  "*)"}]], "Input",
198 CellChangeTimes->{{3.661107683797237*^9, 3.661107689696024*^9}}],
199
200Cell[BoxData[
201 RowBox[{
202  RowBox[{"simdata", "=",
203   RowBox[{"Import", "[", "\"\<simple2cout.txt\>\"", "]"}]}], ";"}]], "Input",\
204
205 CellChangeTimes->{{3.660861462009038*^9, 3.6608615084983997`*^9}, {
206   3.660863565188951*^9, 3.6608635660204144`*^9}, {3.660863601090904*^9,
207   3.66086360284127*^9}, {3.660863643565474*^9, 3.660863646377057*^9}, {
208   3.6608638330074987`*^9, 3.660863833863155*^9}, 3.6608639590084248`*^9,
209   3.6611105171738577`*^9, {3.661154407781514*^9, 3.661154407858136*^9}}],
210
211Cell[CellGroupData[{
212
213Cell[BoxData["simdata"], "Input",
214 CellChangeTimes->{{3.661110523760148*^9, 3.661110524496921*^9}}],
215
216Cell[BoxData["\<\"10 5368 5452 5838 6094 6401 6830 6808 7253 7636 7925 8087 \
2178557 8930 9036 9568 9547 9932 10411 10763 10678 11261 11070 11565 12029 11790 \
21812309 12795 12748 12733 13014 13442 13531 13606 13598 14043 14169 13794 14175 \
21914003 14052 14179 13886 14106 13683 13997 13942 13553 13724 13484 13407 13205 \
22013164 12614 12735 12458 11928 12012 11638 11082 11181 10872 10864 10170 9894 \
2219324 9552 9031 8882 8638 8205 8064 7598 7248 6775 6832 6383 6029 5884 5740 \
2225292 5050 4678 4696 4373 3991 3983 3724 3540 3413 2939 2934 2714 2429 2420 \
2232234 1913 1954 1746 1638 1824 1277 1415 1256 1183 1028 975 720 1003 741 787 \
224474 753 479 524 514 244 457 249 234 226 257 264 241 232 265 0 263 0 218 0 0 \
225222 0 0 0 255 245 0 0 0\\n10 0 0 0 231 260 0 0 0 259 0 0 234 0 226 0 209 246 \
226236 211 240 262 236 245 527 235 478 481 490 745 465 731 711 950 722 1019 1001 \
2271267 1217 1490 1254 1617 1737 1714 1858 1959 2217 2492 2385 2694 2867 2966 \
2283458 3474 3583 4041 3977 4425 4663 4641 5064\"\>"], "Output",
229 CellChangeTimes->{
230  3.6611105248552647`*^9, 3.6611107216888113`*^9, {3.6611537603115673`*^9,
231   3.661153782632226*^9}, {3.6611538424725637`*^9, 3.6611538699620037`*^9},
232   3.661154420883669*^9, 3.661154480035396*^9, {3.661154531920308*^9,
233   3.661154604721693*^9}, 3.66115480919103*^9, {3.661173581170137*^9,
234   3.661173714382021*^9}, 3.661173805191244*^9, 3.661173845397778*^9,
235   3.66117388171874*^9, 3.661173920730735*^9, 3.66122396956768*^9,
236   3.66122400269109*^9, 3.661224090465967*^9, {3.661224161381133*^9,
237   3.66122418617059*^9}, 3.661224219341011*^9, 3.6612244260812387`*^9,
238   3.661224780527691*^9, 3.6612250827541933`*^9}]
239}, Open  ]],
240
241Cell[CellGroupData[{
242
243Cell[BoxData[
244 RowBox[{"simdata2", "=",
245  RowBox[{"StringReplace", "[",
246   RowBox[{"simdata", ",",
247    RowBox[{"\"\< \>\"", "\[Rule]", "\"\<,\>\""}]}], "]"}]}]], "Input",
248 CellChangeTimes->{{3.660863844544029*^9, 3.660863872547986*^9},
249   3.660863961751296*^9, 3.661110551566107*^9}],
250
251Cell[BoxData["\<\"10,5368,5452,5838,6094,6401,6830,6808,7253,7636,7925,8087,\
2528557,8930,9036,9568,9547,9932,10411,10763,10678,11261,11070,11565,12029,11790,\
25312309,12795,12748,12733,13014,13442,13531,13606,13598,14043,14169,13794,14175,\
25414003,14052,14179,13886,14106,13683,13997,13942,13553,13724,13484,13407,13205,\
25513164,12614,12735,12458,11928,12012,11638,11082,11181,10872,10864,10170,9894,\
2569324,9552,9031,8882,8638,8205,8064,7598,7248,6775,6832,6383,6029,5884,5740,\
2575292,5050,4678,4696,4373,3991,3983,3724,3540,3413,2939,2934,2714,2429,2420,\
2582234,1913,1954,1746,1638,1824,1277,1415,1256,1183,1028,975,720,1003,741,787,\
259474,753,479,524,514,244,457,249,234,226,257,264,241,232,265,0,263,0,218,0,0,\
260222,0,0,0,255,245,0,0,0\\n10,0,0,0,231,260,0,0,0,259,0,0,234,0,226,0,209,246,\
261236,211,240,262,236,245,527,235,478,481,490,745,465,731,711,950,722,1019,1001,\
2621267,1217,1490,1254,1617,1737,1714,1858,1959,2217,2492,2385,2694,2867,2966,\
2633458,3474,3583,4041,3977,4425,4663,4641,5064\"\>"], "Output",
264 CellChangeTimes->{
265  3.6611105522063313`*^9, 3.661110721736672*^9, {3.661153760343286*^9,
266   3.661153782678008*^9}, {3.6611538425207253`*^9, 3.661153870008266*^9},
267   3.661154420920793*^9, 3.661154480078126*^9, {3.661154531962288*^9,
268   3.661154604766026*^9}, 3.661154809243697*^9, {3.661173581205278*^9,
269   3.661173714426399*^9}, 3.6611738052389393`*^9, 3.661173845445801*^9,
270   3.661173881741942*^9, 3.661173920774186*^9, 3.661223969611335*^9,
271   3.661224002734375*^9, 3.661224090509625*^9, {3.6612241614240828`*^9,
272   3.661224186214877*^9}, 3.661224219383051*^9, 3.6612244261249647`*^9,
273   3.661224780571287*^9, 3.661225082798772*^9}]
274}, Open  ]],
275
276Cell[CellGroupData[{
277
278Cell[BoxData[
279 RowBox[{"simdata3", "=",
280  RowBox[{"ImportString", "[",
281   RowBox[{"simdata2", ",", "\"\<CSV\>\""}], "]"}]}]], "Input",
282 CellChangeTimes->{{3.660863881721616*^9, 3.6608639034920473`*^9},
283   3.660863963551323*^9, 3.661110561500711*^9}],
284
285Cell[BoxData[
286 RowBox[{"{",
287  RowBox[{
288   RowBox[{"{",
289    RowBox[{
290    "10", ",", "5368", ",", "5452", ",", "5838", ",", "6094", ",", "6401",
291     ",", "6830", ",", "6808", ",", "7253", ",", "7636", ",", "7925", ",",
292     "8087", ",", "8557", ",", "8930", ",", "9036", ",", "9568", ",", "9547",
293     ",", "9932", ",", "10411", ",", "10763", ",", "10678", ",", "11261", ",",
294      "11070", ",", "11565", ",", "12029", ",", "11790", ",", "12309", ",",
295     "12795", ",", "12748", ",", "12733", ",", "13014", ",", "13442", ",",
296     "13531", ",", "13606", ",", "13598", ",", "14043", ",", "14169", ",",
297     "13794", ",", "14175", ",", "14003", ",", "14052", ",", "14179", ",",
298     "13886", ",", "14106", ",", "13683", ",", "13997", ",", "13942", ",",
299     "13553", ",", "13724", ",", "13484", ",", "13407", ",", "13205", ",",
300     "13164", ",", "12614", ",", "12735", ",", "12458", ",", "11928", ",",
301     "12012", ",", "11638", ",", "11082", ",", "11181", ",", "10872", ",",
302     "10864", ",", "10170", ",", "9894", ",", "9324", ",", "9552", ",",
303     "9031", ",", "8882", ",", "8638", ",", "8205", ",", "8064", ",", "7598",
304     ",", "7248", ",", "6775", ",", "6832", ",", "6383", ",", "6029", ",",
305     "5884", ",", "5740", ",", "5292", ",", "5050", ",", "4678", ",", "4696",
306     ",", "4373", ",", "3991", ",", "3983", ",", "3724", ",", "3540", ",",
307     "3413", ",", "2939", ",", "2934", ",", "2714", ",", "2429", ",", "2420",
308     ",", "2234", ",", "1913", ",", "1954", ",", "1746", ",", "1638", ",",
309     "1824", ",", "1277", ",", "1415", ",", "1256", ",", "1183", ",", "1028",
310     ",", "975", ",", "720", ",", "1003", ",", "741", ",", "787", ",", "474",
311     ",", "753", ",", "479", ",", "524", ",", "514", ",", "244", ",", "457",
312     ",", "249", ",", "234", ",", "226", ",", "257", ",", "264", ",", "241",
313     ",", "232", ",", "265", ",", "0", ",", "263", ",", "0", ",", "218", ",",
314     "0", ",", "0", ",", "222", ",", "0", ",", "0", ",", "0", ",", "255", ",",
315      "245", ",", "0", ",", "0", ",", "0"}], "}"}], ",",
316   RowBox[{"{",
317    RowBox[{
318    "10", ",", "0", ",", "0", ",", "0", ",", "231", ",", "260", ",", "0", ",",
319      "0", ",", "0", ",", "259", ",", "0", ",", "0", ",", "234", ",", "0",
320     ",", "226", ",", "0", ",", "209", ",", "246", ",", "236", ",", "211",
321     ",", "240", ",", "262", ",", "236", ",", "245", ",", "527", ",", "235",
322     ",", "478", ",", "481", ",", "490", ",", "745", ",", "465", ",", "731",
323     ",", "711", ",", "950", ",", "722", ",", "1019", ",", "1001", ",",
324     "1267", ",", "1217", ",", "1490", ",", "1254", ",", "1617", ",", "1737",
325     ",", "1714", ",", "1858", ",", "1959", ",", "2217", ",", "2492", ",",
326     "2385", ",", "2694", ",", "2867", ",", "2966", ",", "3458", ",", "3474",
327     ",", "3583", ",", "4041", ",", "3977", ",", "4425", ",", "4663", ",",
328     "4641", ",", "5064"}], "}"}]}], "}"}]], "Output",
329 CellChangeTimes->{
330  3.661110562313175*^9, 3.6611107217855988`*^9, {3.661153760391265*^9,
331   3.661153782728606*^9}, {3.661153842571579*^9, 3.661153870059094*^9},
332   3.661154420970058*^9, 3.661154480129277*^9, {3.661154532013538*^9,
333   3.661154604817532*^9}, 3.661154809289579*^9, {3.661173581256136*^9,
334   3.661173714475307*^9}, 3.6611738052898817`*^9, 3.6611738454969378`*^9,
335   3.661173881784297*^9, 3.6611739208234663`*^9, 3.661223969659273*^9,
336   3.6612240027841473`*^9, 3.661224090561145*^9, {3.6612241614739113`*^9,
337   3.6612241862649937`*^9}, 3.661224219433659*^9, 3.6612244261752777`*^9,
338   3.661224780622611*^9, 3.6612250828486347`*^9}]
339}, Open  ]],
340
341Cell[CellGroupData[{
342
343Cell[BoxData[
344 RowBox[{"simdata4", "=",
345  RowBox[{"Join", "[",
346   RowBox[{
347    RowBox[{"Drop", "[",
348     RowBox[{
349      RowBox[{"simdata3", "[",
350       RowBox[{"[", "2", "]"}], "]"}], ",", "1"}], "]"}], ",",
351    RowBox[{"Drop", "[",
352     RowBox[{
353      RowBox[{"simdata3", "[",
354       RowBox[{"[", "1", "]"}], "]"}], ",", "1"}], "]"}]}], "]"}]}]], "Input",\
355
356 CellChangeTimes->{{3.6608639152067137`*^9, 3.660863929371982*^9}, {
357  3.660864081742077*^9, 3.660864081826098*^9}, {3.661110656925498*^9,
358  3.661110676244894*^9}}],
359
360Cell[BoxData[
361 RowBox[{"{",
362  RowBox[{
363  "0", ",", "0", ",", "0", ",", "231", ",", "260", ",", "0", ",", "0", ",",
364   "0", ",", "259", ",", "0", ",", "0", ",", "234", ",", "0", ",", "226", ",",
365    "0", ",", "209", ",", "246", ",", "236", ",", "211", ",", "240", ",",
366   "262", ",", "236", ",", "245", ",", "527", ",", "235", ",", "478", ",",
367   "481", ",", "490", ",", "745", ",", "465", ",", "731", ",", "711", ",",
368   "950", ",", "722", ",", "1019", ",", "1001", ",", "1267", ",", "1217", ",",
369    "1490", ",", "1254", ",", "1617", ",", "1737", ",", "1714", ",", "1858",
370   ",", "1959", ",", "2217", ",", "2492", ",", "2385", ",", "2694", ",",
371   "2867", ",", "2966", ",", "3458", ",", "3474", ",", "3583", ",", "4041",
372   ",", "3977", ",", "4425", ",", "4663", ",", "4641", ",", "5064", ",",
373   "5368", ",", "5452", ",", "5838", ",", "6094", ",", "6401", ",", "6830",
374   ",", "6808", ",", "7253", ",", "7636", ",", "7925", ",", "8087", ",",
375   "8557", ",", "8930", ",", "9036", ",", "9568", ",", "9547", ",", "9932",
376   ",", "10411", ",", "10763", ",", "10678", ",", "11261", ",", "11070", ",",
377   "11565", ",", "12029", ",", "11790", ",", "12309", ",", "12795", ",",
378   "12748", ",", "12733", ",", "13014", ",", "13442", ",", "13531", ",",
379   "13606", ",", "13598", ",", "14043", ",", "14169", ",", "13794", ",",
380   "14175", ",", "14003", ",", "14052", ",", "14179", ",", "13886", ",",
381   "14106", ",", "13683", ",", "13997", ",", "13942", ",", "13553", ",",
382   "13724", ",", "13484", ",", "13407", ",", "13205", ",", "13164", ",",
383   "12614", ",", "12735", ",", "12458", ",", "11928", ",", "12012", ",",
384   "11638", ",", "11082", ",", "11181", ",", "10872", ",", "10864", ",",
385   "10170", ",", "9894", ",", "9324", ",", "9552", ",", "9031", ",", "8882",
386   ",", "8638", ",", "8205", ",", "8064", ",", "7598", ",", "7248", ",",
387   "6775", ",", "6832", ",", "6383", ",", "6029", ",", "5884", ",", "5740",
388   ",", "5292", ",", "5050", ",", "4678", ",", "4696", ",", "4373", ",",
389   "3991", ",", "3983", ",", "3724", ",", "3540", ",", "3413", ",", "2939",
390   ",", "2934", ",", "2714", ",", "2429", ",", "2420", ",", "2234", ",",
391   "1913", ",", "1954", ",", "1746", ",", "1638", ",", "1824", ",", "1277",
392   ",", "1415", ",", "1256", ",", "1183", ",", "1028", ",", "975", ",", "720",
393    ",", "1003", ",", "741", ",", "787", ",", "474", ",", "753", ",", "479",
394   ",", "524", ",", "514", ",", "244", ",", "457", ",", "249", ",", "234",
395   ",", "226", ",", "257", ",", "264", ",", "241", ",", "232", ",", "265",
396   ",", "0", ",", "263", ",", "0", ",", "218", ",", "0", ",", "0", ",", "222",
397    ",", "0", ",", "0", ",", "0", ",", "255", ",", "245", ",", "0", ",", "0",
398   ",", "0"}], "}"}]], "Output",
399 CellChangeTimes->{
400  3.660863930323296*^9, {3.660864063968863*^9, 3.660864086715946*^9},
401   3.660864234921925*^9, 3.660864284709508*^9, 3.660865035179634*^9,
402   3.660865096334959*^9, 3.660865140008404*^9, 3.660865495568335*^9,
403   3.660865527889092*^9, 3.6611075830063257`*^9, 3.661107927347301*^9,
404   3.661108437430786*^9, 3.661108526510558*^9, {3.661108636061977*^9,
405   3.661108664862555*^9}, {3.6611087262645483`*^9, 3.661108740363573*^9}, {
406   3.661109730888185*^9, 3.661109746440199*^9}, {3.661109806761413*^9,
407   3.6611098478868628`*^9}, 3.6611106774851847`*^9, 3.661110721837788*^9, {
408   3.6611537604282312`*^9, 3.661153782778905*^9}, {3.661153842621489*^9,
409   3.66115387010841*^9}, 3.661154421023695*^9, 3.6611544801787043`*^9, {
410   3.661154532063583*^9, 3.661154604864925*^9}, 3.661154809339509*^9, {
411   3.661173581310166*^9, 3.661173714530191*^9}, 3.661173805340065*^9,
412   3.6611738455466557`*^9, 3.661173881834208*^9, 3.661173920878592*^9,
413   3.6612239697106113`*^9, 3.6612240028355207`*^9, 3.6612240906099157`*^9, {
414   3.661224161523966*^9, 3.6612241863149643`*^9}, 3.6612242194846973`*^9,
415   3.6612244262264957`*^9, 3.661224780671605*^9, 3.661225082899951*^9}]
416}, Open  ]],
417
418Cell[CellGroupData[{
419
420Cell[BoxData[
421 RowBox[{"Length", "[", "simdata4", "]"}]], "Input",
422 CellChangeTimes->{{3.661110694144629*^9, 3.66111069883423*^9}}],
423
424Cell[BoxData["200"], "Output",
425 CellChangeTimes->{{3.661110699243342*^9, 3.661110721886673*^9}, {
426   3.661153760476221*^9, 3.661153782830151*^9}, {3.661153842671913*^9,
427   3.6611538701596518`*^9}, 3.661154421073723*^9, 3.661154480229265*^9, {
428   3.661154532114209*^9, 3.661154604916913*^9}, 3.661154809388966*^9, {
429   3.661173581355053*^9, 3.66117371457561*^9}, 3.661173805389761*^9,
430   3.661173845595358*^9, 3.661173881884508*^9, 3.661173920924575*^9,
431   3.661223969762288*^9, 3.6612240028856153`*^9, 3.661224090660366*^9, {
432   3.6612241615741367`*^9, 3.6612241863669367`*^9}, 3.661224219533532*^9,
433   3.6612244262752743`*^9, 3.661224780721755*^9, 3.6612250829486*^9}]
434}, Open  ]],
435
436Cell[CellGroupData[{
437
438Cell[BoxData[
439 RowBox[{"nmolec", "=",
440  RowBox[{"Total", "[", "simdata4", "]"}]}]], "Input",
441 CellChangeTimes->{{3.660864954122253*^9, 3.66086499765872*^9}}],
442
443Cell[BoxData["1000000"], "Output",
444 CellChangeTimes->{
445  3.660864998833362*^9, 3.660865035214733*^9, 3.6608650963755903`*^9,
446   3.660865140047003*^9, 3.660865495587974*^9, 3.6608655279246063`*^9,
447   3.661107583052704*^9, 3.661107927398151*^9, 3.661108437449894*^9,
448   3.661108526561778*^9, {3.661108636080645*^9, 3.6611086649122458`*^9}, {
449   3.661108726309012*^9, 3.661108740411521*^9}, {3.6611097309385567`*^9,
450   3.6611097464900093`*^9}, {3.6611098068160257`*^9, 3.661109847934825*^9}, {
451   3.661110702625884*^9, 3.661110721937879*^9}, {3.661153760525971*^9,
452   3.661153782879583*^9}, {3.661153842703539*^9, 3.661153870191695*^9},
453   3.661154421123344*^9, 3.6611544802804127`*^9, {3.661154532164153*^9,
454   3.661154604967559*^9}, 3.661154809425294*^9, {3.661173581406*^9,
455   3.6611737146247597`*^9}, 3.6611738054398413`*^9, 3.661173845628528*^9,
456   3.661173881916006*^9, 3.661173920974317*^9, 3.661223969813113*^9,
457   3.661224002935402*^9, 3.661224090710704*^9, {3.66122416162426*^9,
458   3.661224186415546*^9}, 3.6612242195848207`*^9, 3.661224426325164*^9,
459   3.661224780772645*^9, 3.661225082999001*^9}]
460}, Open  ]],
461
462Cell[BoxData[
463 RowBox[{
464  RowBox[{"xvector", "=",
465   RowBox[{"Table", "[",
466    RowBox[{"i", ",",
467     RowBox[{"{",
468      RowBox[{"i", ",",
469       RowBox[{"binlow", "+",
470        RowBox[{"binwidth", "/", "2"}]}], ",",
471       RowBox[{"binhigh", "-",
472        RowBox[{"binwidth", "/", "2"}]}], ",", "binwidth"}], "}"}]}], "]"}]}],
473   ";"}]], "Input",
474 CellChangeTimes->{{3.660863989428227*^9, 3.6608640224666843`*^9}, {
475  3.66086410371815*^9, 3.6608641478083897`*^9}, {3.6611077570821867`*^9,
476  3.661107762473649*^9}, {3.6611078836257277`*^9, 3.6611079035136337`*^9}}],
477
478Cell[BoxData[
479 RowBox[{
480  RowBox[{"simdata5", "=",
481   RowBox[{"Transpose", "[",
482    RowBox[{"{",
483     RowBox[{"xvector", ",", "simdata4"}], "}"}], "]"}]}], ";"}]], "Input",
484 CellChangeTimes->{{3.6608641596188*^9, 3.6608641833119926`*^9}}],
485
486Cell[CellGroupData[{
487
488Cell[BoxData[
489 RowBox[{"sigma", "=",
490  RowBox[{"Sqrt", "[",
491   RowBox[{"2", "*", "diffconst", "*", "timestep"}], "]"}]}]], "Input",
492 CellChangeTimes->{{3.6608644506570587`*^9, 3.660864457870208*^9},
493   3.6608649175533648`*^9, {3.660865092376251*^9, 3.660865092474855*^9}, {
494   3.660865135842045*^9, 3.6608651361812983`*^9}, {3.66110791247908*^9,
495   3.661107921033339*^9}}],
496
497Cell[BoxData[
498 RowBox[{"10", " ",
499  SqrtBox["2"]}]], "Output",
500 CellChangeTimes->{
501  3.660864458729279*^9, 3.66086491906544*^9, 3.660865035282187*^9,
502   3.6608650964589252`*^9, 3.6608651401142817`*^9, 3.660865495652907*^9,
503   3.660865528006706*^9, 3.661107583168273*^9, {3.6611079223547487`*^9,
504   3.661107927491674*^9}, 3.661108437532053*^9, 3.661108526640664*^9, {
505   3.661108636161117*^9, 3.6611086649828053`*^9}, {3.661108726392659*^9,
506   3.66110874049479*^9}, {3.6611097310209503`*^9, 3.661109746569604*^9}, {
507   3.661109806898305*^9, 3.661109848021763*^9}, 3.661110722014998*^9, {
508   3.661153760592411*^9, 3.6611537829460278`*^9}, {3.661153842760065*^9,
509   3.6611538702466793`*^9}, 3.661154421216859*^9, 3.6611544803991623`*^9, {
510   3.66115453228211*^9, 3.66115460508428*^9}, 3.661154809505584*^9, {
511   3.661173581489431*^9, 3.6611737146921797`*^9}, 3.661173805523896*^9,
512   3.661173845697287*^9, 3.661173881954224*^9, 3.661173921041184*^9,
513   3.661223969925725*^9, 3.6612240030543118`*^9, 3.661224090828182*^9, {
514   3.6612241617417507`*^9, 3.661224186533325*^9}, 3.66122421970347*^9,
515   3.6612244264449987`*^9, 3.661224780891039*^9, 3.6612250831172943`*^9}]
516}, Open  ]],
517
518Cell[CellGroupData[{
519
520Cell[BoxData[
521 RowBox[{"Show", "[",
522  RowBox[{
523   RowBox[{"{",
524    RowBox[{
525     RowBox[{"ListPlot", "[", "simdata5", "]"}], ",",
526     RowBox[{"Plot", "[",
527      RowBox[{
528       RowBox[{"nmolec", "*", "binwidth", "*",
529        RowBox[{"gauss", "[",
530         RowBox[{"x", ",", "50", ",", "sigma"}], "]"}]}], ",",
531       RowBox[{"{",
532        RowBox[{"x", ",", "binlow", ",", "binhigh"}], "}"}], ",",
533       RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "}"}], ",",
534   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
535 CellChangeTimes->{{3.660863937053011*^9, 3.6608639513266993`*^9}, {
536  3.6608641879635687`*^9, 3.6608641880479*^9}, {3.660864311465605*^9,
537  3.660864313050109*^9}, {3.660864517827593*^9, 3.660864582236044*^9}, {
538  3.660864619503529*^9, 3.6608646655904837`*^9}, {3.6608647032550917`*^9,
539  3.660864721064691*^9}, {3.660864924891444*^9, 3.660864925083404*^9}, {
540  3.6608650038562393`*^9, 3.6608650059230328`*^9}, {3.661107941133526*^9,
541  3.6611079647069187`*^9}}],
542
543Cell[BoxData[
544 GraphicsBox[{{{}, {{},
545     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.01388888888888889],
546      AbsoluteThickness[1.6], PointBox[CompressedData["
5471:eJxd1l2IFXUcxvFRzEIiJkNSenG21LR8GdN0S8vR1F1zq9ParuvdIciL1JxC
548bC9M/le1N5uDgm5pNe6yKmY6am1pL560dLOQuRAJwjwXEV4ETmARYVLyfH/I
549zgF5+Pj8f/OfPS8z0/DS+taXh3ueN/j/vxupV77AG/K6XPLVkq/J9bWRPFwZ
550xXhkNHT9qJLvKPlOObT5u0r9mJLvZr91/P+4Un8PvZ3ffaV+vBy8onSBsvYq
551biA5vntAWbX1D3I8Wz+BnvN3E0vzk5SprX9I2enwZOZt/RRl/U38sNLfjB/h
552eOapyua38TTW2/x05U9v4RnKRhyFyqxbWcOJ9TOZ30KPm3H0qPLHbfR4BI5m
553sX4HPe6xfrZy8bv0eMJO+seUg+/R49vep5+j3PUBPW7dTT9X2dlHj3t76RuV
554t/fT49V76B9XTt1Lj0d8RP8E54tr+Lr185RbD9Hj0Tiar1x2mB7PPUL/pLIB
5551/DgUfqnlH9/Qo9bPqVfoDw3QG/+TBlEymHHlFV88LgyxT24jtu+ZH6h8vxX
556zOM/v2Ye7zjBvOwGaswvUr72DfN4D05xdpJ52fXj4Gnl6lPMy67tW+bxLd8x
557z/qf6YPF6q/TV2XvQ5zSbz/NPN58hvklyltxVXajBpnH/+I663/DwVL1v+Cq
5587K35nnm88CzzrF+HgyblNlzFW3DaxPofmMcHcNCs/MPmZXeIPsWrbJ71r9v8
559MvWHcVX2xts87rZ51h+z839Gbrd52Ztn8/Sd9vfTn7P55cp2O3/ZzbC/H5/l
560/ayzfgAHLco+e/9ld8k+P/o2+/xwi33+zyo32uePv+D7kuJ++/7J7gIOnlPO
561t+8fbrTvH77I97uOV+DgeeWck8zjfTiVXR+/pzr9KuxXlJf4PYayG4kr9I7f
562b4wP8ntPcIXrQYbf4HqR4wwX+Mrn7P+CspvrUYivcb2q4LU4xge43iX4ONfD
563DC/lepnjDlzgMVxv/VbldK7HIb4fV/DYj9kf7+P6nuD9+9kf/8X9I8e/4wLP
5645n7jr1D+08v++DKu4A3cz2L8K/e/BHdx/8twM/fPHPdwfy2w28X+Lyqvcj8O
5658b3b2R+n9DHetJX9sZ+wP+58h/2x43miwB08P/htysk8v4S4p4v9ccrzTIy7
566WJ/ggj7DVZ6vcuzjAjuex/x2jrOJ/XGwnv2xPW/GZp4vE+zz/Jdhh3M7HvsV
5672ON4/kqOw/oQ26uCY/q41Cc4WsP+pT4vucA11vsdQ/uw5ErJMa7z/iTYnm+z
5680vq85OKm/wN8k8zc
569       "]]}, {}}, {}}, {{}, {},
570    {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
571     1.], LineBox[CompressedData["
5721:eJw1mXc41f/7x636pFQICRWyMsuKyHlaxzkHWcc5x0yyV/ZKMpKSVRGVjEKi
573zIg2ZZ4io4GiRVpCKCLf8/vj99f7elz367rf7+u+n/d4XW8J18PW7hxsbGw9
5747Gxs//eUzFi+mr+tWS/2prB94l4dmARd4/b3KiC49+m/lw7UgdV7xQw5ryrC
575Y4rMqZ2lOqBd+Gm12ushQUhbtI62SReMOUu1Oc8ewopr32e7GV3YXWUv/u45
576SnhnTQ45+1gP9ot2jhOePwlS48o8jQ76CG479uDc7ilCcIpzu+5RfTRm9O1X
577q5oi7B9zfLaUrw8DycjDweXTBEkRr127PuqDYfK4eqrwFyEqqfBGQIABjp+x
578U/uR9psQRPv3Iz/TEEPSJ7TGPVcIljtculJWE7EsdXxUaGqFUFUcePKeGhG8
579/PqPd+iygRLM/tngIBGqfAQf6hM2NKesbbS4T0TkRp3G+gF2+LNTNaKiTMDF
580o0aNnONEkkZAYgoHGWKrdqQua66B7yyHdruDGSq5rowKHlkDC9VB3RtZZiBw
581SagpP1yDuZK8+T/PzODKuX3I2YQbD781Sz/VN0cZu6jcI/paRLWeeuWrvB9q
582//ifJETy4PfFKuNNGyxRcoT9H6cTL14ruUyryVojYbY9tTuGF6soy0X+JGu4
583+KeLXsjjhZLQeasRb2uIHRDVVn7Di/COa9/cb1rjrIFGKMOBD70lgr1Ce2xw
584jNv76w07fpwc6W0h2lJhl/P8pS1NAG7c/sNKD2nQ5M1xFw9ncSyv28RnGgRO
585Oc1+zRYAMbE08xMvHc+jv/LFvRCA7v2fhlcP0UF05jIvpwoief9fRSMeBtSk
586tR7/sxZCiUKIrIePHdbXFVRdsxBGUcnF8eeOjggXWsjtCxCG7cSlntg0R4xE
587WScspwmjSzf6dMQDR1Rjla3NU2FcSFNdbSrhBOtun8V/pC2o5TMeP/7VCblf
588NIxpBiIQG/x2Ju/0AUiKPx3i0hBDNaOBfCHpIE4nSj9RoYqh4s9/1WxXD2J2
589PPamfYgYsmturs16dBCtN3fHVdeI4bDH1BHC0kF46eRIOypvRdK+p8FVYa64
590SXMNrJPdBonZ68elgw4hj7jEsXu9OO5cef+KM9UdPAdiLwpvE0eG0rFu+1p3
591xESwqbKpiGNoc47K5Gt3OJZxHeyxFMeX/PRApowHxNauf+ifJQ6/mLYmr1YP
5925D3beqRcVAI8ay63hW3wwmUq4dcOeUlcufjqlFSfDzb4N6es05EEPWSP0nU2
593X8QmGUr+MpVE1NB0R6iKL5wbTKxa/CXx4GaaxYc0X2wTtqxyqZFEzHepMPH9
594fsgfdvG9rLUD0dFr2KPe+CPfNeGjkIkUPMS0BLMvBCL3d/LNtzQp5I6/kXhU
595F4hzqWkRxR5SeK/sr76xOxAnG3LXqZ6QQlX8bR4RziAEr61SM2+VwkwHYXE8
596IAjE2uHERENpHKwPuxFvHoxJdg2pKYIM9r4yzByVDcXE+b2TDRYy0Krhrx8m
597huKjAhqPHpDB8OFNO9g9QvHa1tRs3TEZ7PKtSekoDkXL9YOhMg9kYFA54Lxf
598KgznrdKfOOrKQrNwWcFJIRz7ij67dWrJYRN9381W+0jM5rVYfSDLgW3FlKWc
599SJTnXtb7ay+HcQ/BGKHySAhl2ggrHZVDXEfDZuPFSEzGPmRmNsths5WDmeHl
600KBQ65qgyKDtxOdDcbOFHNDiETdg+O8jj28a1q9juHUXjJskfK37yyN/mH/3t
60161EEbFweFI6VR2VcyYKgSCyGV9fVUQrlwbM0JbwzKha35rZ5Vn6Sx8W/8Qv7
6029x6De//8szB/BYyOlFt4u8WhI730EtcxRTB12GnGPAl4stJV/yxDEdSyJR0f
6032QQ8CvzZc75QETlqd07UGCTgtrUW184WRUxvVHh6LSoBpUJdfuarlMAhaFGi
604+SUBifk/9p0/rQQCm4z8Ynci9CrV38nmKsP4gH7ySn0S9m63W5wqU0b/ulUV
6051S+SoJl5VOBOkzKkHw7anpxLglJwG8lsWBlRosfbKzROQEyTURO4TQUprnpP
606226fwML9IwlNxSoI1xBf8XmcjCGPg2dLGLsQXtuaZL9wCped6hujXHZh7wHl
607iGLRFLhQuUfNvXbBxr+naJNeCsb1axTmI3ahqewIv2tiCmbE2FuJObsgKWEU
6082LLxNNb1F/7+/GIXRH5uuHhdJRW6hHeO8ja7MZ/b/+RFajpWNNQT/9nvhmL5
609gMhyXTpaFE9e73PdjfwEPXnz4XSQRHbNRwfvRvRqjHnKZ8B27lgG88xu9Bho
610DKd0ZSCgYnuL3/PdGEjwzLM4lYmCzQdkq8xV4ZlAMHd9fwaPmrceUKSpgnPD
611ngtHF87gne+b89edVTHYU/Khlu8sJB/ZrS4+rIoP8VN1vgZnUepFHcs9o4pI
612tkPawiVnUdlEvprwQhWtbXP/+R0+h3sO6ttpTmqInEjTiRDOxlAB9+ZlP3Uk
613SEpptD/MBdasnfgSqo5oSXm5F29yURq4tulljDrU6sUOLC/kIlh/nUN1qjpE
614SpK78tQvYN1HnoJDN9RBuxs4LVxxAboyvLJd39SxJZ6syHf5IvJvCO3J8dHA
6154J5Gd2p2Hg417aDv9tKE+5KpdpdfIfzedOjNH9Zk9RFhZb/YQoSu+MvcjdDE
616mgRBPrHMQpwwbpwzStaEvW5GbUldIcp7zbIYpZqI0d/z2uVvIWa+hPUe+6SJ
617Z7U1WtFpRYgT6TDtdtmDUaX+QZt7V3D5iB987bTQMrrTkmFcAi7rGO3Bg1q4
618rP51g4tjCXzkUlVNfLSgXLWyNSqkBFovKqR2HNFC5gdVzdGiEvQrf1sznKeF
619PaOBKUdWSsD9wbuX8k4LdZPwmb9fijCSp6u8pzYKp1V6/5mUwVzQ9fiX0L3A
620sRfcAZkVMEW78H4/XdybeFLtql6DhwaD//1O08PphTixeK96zMzzPO+zBjiU
621w8fHfzVBUdDd6LKZPm75ZJEkQh5AomDT+CkHAxjGbB8MaWzG/emFKm+GIbz1
622rEQdOJ5AoFrzmq+PEcJF9uVetWnDATaPOwUBRjhUeX6TrmMbyi2zn/UHG+Gz
623JnG2x60NhOlfv3SOGEFShNr2JawNPqo14Ek1wsvnNn2fctrw6JbC0I2bRjj/
624QvDO9FAb/JrE109OGeEpTaxjl1M7HresDQmMMMb8fYPcy5Yd2DUzZO14xBjg
625T3tIp3cgX6JClXTMGGKmfUY8zh2IjDOd2Z5sjKySi9Zuvh1Q1ksN6jlvjI/k
626DW97kjqQ27Q+ULnBGAompp0ZTR0IqOb1/z5rDPuLrpWpYp14M/rO7PUfY1xu
627nhhNluwEeWON4pMlY1yfc2+OkeuEVIDV90tcRKS0FB+zUe/EoMJZX1MBIjgp
628/Xa3TDthdG2TTzlrbyr5HL+FENUJkXwhT69gIlatX/Sj9XRCUzavpyyciMqx
6293+3fBzphXS2h9SWaiL2y71aihzpxukWR2yeRiF9H3X7GjnViadygwjeLCLrj
630BSWpv50YUTn8M6CBCJlgccNRqS4sNs7Sq+6w7P5RnVM7uyBkEP3o5wMiiiSF
631BxaUu2BBTToT2E7Ew4uCB2e0uvAo8pJa8GsipsY8OmLMunC1uT0ibJEIF/Mk
632+5kglp1i/q7+HxGzutEPosK78Ka/jzTPYYLQJV3an+guCIy/3RKxzgSDz7Q3
633DCd24cTa2buRYibQ8HXeqpndBS8bcbYYPRMU/J5+qNfQBcWxyJMJCSbYNqC5
634K2a6C77mBfejT5hguxLXzrLZLpTXt84Ep5jAquUm+9PfXZA7we906KwJ4q5M
635EP8ud2GH7I3dRldMMCY3+vw3NxMi3qNDq1pM8NNSL9pDnAm73lW8/1pNsGj/
636UWLjDiZytRWN5ztNcDF4Y2yVNBNCa6OqxntNsPROmfJKngm+Cv7j7e9NcKUy
637a7JFnYk1P4yVTrGTMESN8V5jwoSJrZ9r/CoSXtdYGCqQmThx/2xOFDcJZj9+
6383DQ2ZYIrfZTdh48EYmiF0kELJlaUo15SJEjYrNJapEpnQi+nYJ2hNAn2hjZP
639+eyYiFlphc5OEnZ6NEVO2DOx2MNfobCbhDMUecsYZybmA28c49EngXdMgrnP
640nQmNwb5bXMYkFIfee/Hag4lQ/YUvSySWv3hNW18vJmb4iNQfliSI7vl0I9yX
641icnaUbluFxJm0otzZIKYUBJd7dzmRsJkUbzE0WAm/BIVzz3wIiFhzkiQGcLE
642V5uo5cpAEgyTLmqYhzMxNsvflxFHQrr3Wu2v0UyE3znNH3qchPcH8krnjzCx
643Oo7LhnGShLTTf7oWY5iQWTffL55JwvYxD/WvsUw0PA8QWJVFgnTS7YgXx1jx
644O/+Z+iWHhBB9i4DGOCa8JAZf1BSQ8OPcjmD3BCb+jFsJnb9KAvXO0FGVRCZO
6453uiiRV8joWV7t/Y0i6/vuffKoIqEl3k1Dx2SmNi7rC4sW0dCYTLjJNcJJrpa
646bjLW3Sah2u76v2IWfzUvGOx/QEKRPffM02QmojcJizS2kFD3YF0A9SQTawcz
6477fPaSFgdnp/dz+KL+dyX4rpI8GDb70U5xYSCW8KwWzfr/f6VY40svrfzryi5
648j5V/Z4f121OYMPsZ4qj0kmVf9WckhsVvbn3P4xsi4ayMin0fi/2i3d/OvSXh
6490ea62O2nmVgijGwdek9CVcs+MzcWp66iOz8YY+W3OLOlkMVizJ78K19Y8X8a
650+n6AxTczSaMnfpBwbfROGVsqSz+05u2+0ySMC+kIS7G4W3Svi8UcCSqEFyoE
651Fju/ry1UWyDB807gTysWT5YqvN+8TML+Y0s0RxbH+hVLLLGRcXHLYQ9nFm9Q
6523er6josMjvK2rQwW5//OvvJkDRnG42uOklmscn/DxzIeMiS3GB5XZfGjhOQd
653abxk6EqmqG1isSWJzS1IgIz/whdTvrG+9936qGJbYTLStjam3GVxYP/0J20x
654MigNk6qJLGa74CO9TZyMl0G9iQYsznT+6M4hRUa54KWjC6x4iUs5lo7LknF3
655MU28jMU1XwbGuxTIEFk/67efxQZV5rJVKmQkMrcd+s7KR19om+c5NTJGiJ6c
656CSyeWbk94aBDxvuMYcMsVj4TWnftBIEMhu3oJ14W85++7i1lyPLXVr07iaUH
657NaFLX79TyBB/8vA5laWXJ8ObFHr3k9HWzCdVzdIXtSjVt96ajLdke2lOFocq
658xn4/ak+G3weKeiZLj1wz84quzmTUNRsadrH0m3X7sD/RlQze9gGOpXgm6g0O
659Tm7wIcPlYd15Ikvv8wyjqYIoMkaXm2QjWfWiERnc9PooGTvE2DnsWPUUmlOY
660wJdAxvSjaLFdrHr79WJJIDGFjAfjnJubI5iYsq7X9rhEht31EpMAVv2qBH/k
661KCggY2h/IImdVd8BZ/ievrrKyi9fRFbKYSa+9/g7k2+w7L3xAnF+THwxk0lU
662vE9G2HTnh3lWP5Hzo5q5N5Ox/T9fX103Vj2eThDMbyXj8OtPZpGuTIx3jlzb
6632M3Kj8sviX5WP/pAzHk6M0JG/Gi5qSCrn73BGqEmNgquzGpHPjNmQtRFc3Sa
664iwK9aw4nogyZsD/mVibPTYGsmAqvmD4Tg/cf7c3jo2DdP7YiI10mXu6NOnBM
665koLhQAE2KVUmetS/lhkbUUCtixbZJsbSJ1U4OJZEwX1ik4zaFib2hxB1bptR
666YHAnYlhPiImnNVeeydlSYC62TNnLy0SHssPMOg8KfOfZ85o5mWiRe6rTl0zB
667zAD7rMCXLkjJXuptOU2BWMgXqb4x1nyS9vGsy6CAbf8B1+QPXSBLrjmXlUPB
6680u31G4eGu/Bc1Ogr7RoFYRElSeLdXRjdcC93uI2Cph/c1eo1XViarZj7tMoU
66938bkrpxhzUfnX9GnX3CbQlbSv4Xmz5qX02SJtvWmaDDgcBLw7kLi5Geza4Ks
670LUz5LG+USxfWTuwo9pY2hXz/85GJ/V0QHr5kPWlkCrmiS7nu8l3QaD5d9TvR
671FBkkYufUcCeIW3i7dU+aQkq9Lc35ZSdoQVnf41NNsSJ4SfTJ805ESlzeyZNt
672Cm9KS3l4ayfuxFdeFS81xW0luX3RlZ2AQW8OuYN1nosrhu1YJ8weC8Vd4jHD
673nnYrRQ2RThxqLbIkZJvhhd7G+rq9HWiLjowOvGCGw+W39bzVOiC3y6K46LIZ
674oqLFjggrduD7xaXfnKVmUPttJue2tQOhgfSCjgYzbAq9LFKz0o7jousnrV6b
675wUO6tkD2cTtKgyJPu4ma44moPNcd43Z83WrRdqrIHG38DdZCOm0QJqi7tV/a
676jwOXzY9yb30C9qeHVWQrLDBYKC9UW9aMjswgvn1US/ztzLOyLH2Ae2Phc57P
677LXHlulzXdYU7SA+UndpBt2Ld07cZkbPrUTf6K3ZwwAqfqk6ol9yoQcNxXylx
678J2vcvS+rPjdYgfUbVpxc3rAWXul1ldGvSvAmJd3NZsQa2Y3M66+aS1D+3zZf
6794jtrlLYvHTW8UQIi+75IxU/WOJy8/bpJfAliZ6PP/vlmDQ/+4UF2xRJMDs23
680Zf61RtDXhEWd48Xovjap0ixiA5mU0czbhKvI1B9hl7SzwcafvhMKHwqh6dth
681RnCwQUKjh39JdyGGs2pzHJ1sMHJ1/WmVu4WQmTihlHPQBmU7Vx0JyyrEvXQV
682Ox4fG3yLyMtuNynExHBc9Xy0DQLZDIpWXAuA8B0HmHk2aP1cyq0seRlT5V73
683Qt/Z4Eu3ZdC3f7mQNbr6w/+DDcRN1zm9eZ8L57dvt3l+skEa1wW7z09ywdxo
684E283YYO5Dz02pqdZ963QfUS9KRtcy5qZIm/JhZMe//P/2KjYRNDuFdqbg67e
685ux8vbqdi6dNiSHZ6Nth95wWyJKgYWJn5WBWaDS2u3cS0HVSUsn3v+GCfjRLN
6860rJjslS0FDrdOyObjbhLmf5uKlT8a9fu5W/Jwh53999KBCpksreMZP89h+I/
687G9Y1O7PO92rcIsSdRa1ecdwnFypya8k8bp5n8ei49vx/h6hYXTa/+er+s3jL
6886/Z+vycVknw1HEFbz0JoZ9PtN4ep+HSgeBfl3hmctDvk9ieOCutnNns+/MuE
689353b91WuUhHGzjFMVc5ANJu5mk0Jy1/04atcnBk4SfxQFn6NCrpPmfnYq3SU
6909PFk3a+gwio8vF4rIR0jXw/6mN6i4rnzQP+OwTRYivJs9mylomDynuW3jFRo
691xLgE5n+m4j0SDhHEU/B9Kv/b7S9UDEly84f8PYUr7m89er9RESdIjux9eQq8
692FvZOXFNUCFG3f1iffgrfJagUnz9UDHMqeMb/O4nidhMpTW5b8OywlNb4nAyB
693TSqvn8nb4r1zL9u7l0lgJvvbfFa0xbuRjwS520lIWKroZlOxRYWxysz5nCTW
694Xi/XqqZmi49b63Z+ZiThaaNk7QUdW1z0lt1j8fY4jjsLpbqb2aJpw4lj4j8S
6958ev6MmHZ3xa9Z6b9R2UToJueWjRz2BZP5srjkzck4HiwKOdEkC0OLPKHmM/F
696Q0BHu60vzBZkp4Ek1cfx0HwaalZ21BavhcLTnF3iETn5jWGTZgv142N02cI4
697LKkOBV2/YYtPFPf/MrxjYbTZuz+/0hZOoqyddV8sUv/+Vs+qtsVQXFhnCn8s
698tj4R/BN7yxZjYfNcDveOQs/WKpZ6zxbbvp6qCuY/iriIjpQVpi36z2Xk2HQc
699Ade9hqu232wR3Xmt9YRdFKTMC8lhP2wxElfXw7c7CkYjp35m/bSF/Or7rgNr
700opDI5qQz8MsW6eyj6iuNkeAw5uy3XrJFUGzfmjHRSLB1W3BYrqeB0CP34f63
701cCyNTrhQVGj4Wqkdv7E2FGJBfat9dtMgsu5BsUhmKHQ57t04pUaDV8dOcVpA
702KI5IZfzp2EPDdd0IAyeFUCx6aWSagAae9st7Sq+F4M90/EMjKxoG5ndG/SoL
703xi9O0W2EEBrk4wW41zwLRDGXbfneMBrOMc0iztcGgrYqQ1MzgoZZ9gQzcm4g
704GldzWigdoSFTLO+kqFsgYri/HRNNpEGUK6ZyePEwuDbeff/7LA0Br/yczf0D
705ICDiUFJdS4OF4KLwzeO+aBXJ3n3jFg3kAgcOCsMX4aI996810JCawbkgreiL
70612KGL/Pv0FC2Z37swoAP8rYr/JfWQkPRM1qtt5wPpKT/enn30eC2Re2I97AX
707NHZdUpScoWHf2ZADC34eSD8Q+2DTLA3+n4xGSkkemEg/aLlqnoaQQ3FxJ6U8
708kPdDNnRigQaF+kD1z2/cwVFx624lOx1yOt0p6yzd8Vy6m6LDR0dpq+fOOwQ3
709+Imye1N30aG63MDDY+iKNsqnRWNVOkJu8yZ1SLlCPLo9dY86HVF719fWrHZF
710/+u0ahEtOoa4hCnczIPYe17kzzsCHVYUkQPttgexhk892d+CjsKVN0oWIS64
711utqz5EQAHf2Pjfh3X3TCwIfiI1cC6ejmMkpvO+SE1Q8/WD0IpmNfTUZ1vZIT
712vCKc/82F07GeXWsyv8URihO2dPdjrPMTMjkZUw5o6DTiNsqkY63j6Yqjjvbo
713TJX0Z6+lQ0zebZIURGfpwcVw6y06/q5z15IzpkPROH+LdgMd+dQ4B+stdKQv
714b2kLvEPHpuAXYUmsuFsH8G1/10JHsto5AQdhGt5YsPU+6Kdj4OVJD49eKn7y
715j6jHzNExsfazoXy4NWqd/aZDftNRP6ct94dqjdDyxZu+C3TYRiRvWVKzxoL+
716ZlmHZTqaQod+Fs9YgTPIcsveVQxwJyoLBwdbYXNPy/JvAQb2xJPfNh6xBOH0
717tbYQdQbKjhbzLrmbg+OVRqKvJgO6X7yJ/1RYc17yCeGQFgNc0vnl9otmIN95
71812ity8Bu1ajv2zPNYPNlS8VuIwYWSJuWnj00hadJasZPGwbOa244bqJCQQZn
719EMM3hIFbXvdHyKYmMOL4UX45jIE0nRtmaSIm+MPmvdwTwYBnqdSQ0lciXP65
720XFGPYeAe6fDbxBQiVBcsvy8fZ0A8yC46uMcYLyd3xWWeZ+D9J67kAW8jpPy4
7210deSy4BzZsNQkZ4RCN/lpOcuMhBAdHrdsskIZV/Eu+wKGNgfmHBm+YEhoj/x
722CuwoY3FSsRdxiyG2DU1dq29iQECN807FG330vfZb/HyXgYJu7tdBDfpIfjVh
723JvKAgckx9uWETH1MDbybjm1hYNWqI/uPG+ujpee5DonJAEe7zNvZWsCjtbpn
7246A3L//eOLfs0CRB5oii5fpSBHZq8GrUzeuhpKQslvGdggJZfdLNKD1qPCrcU
725jzHQ7nyhVk1BD2vvZh7yn2Rx891RV7l9eNDE01AwxcA24eHSqQldhDSeXNM3
726w0DnGIff8nVdvKmPu6n5m4HN74wPZirpIvPW0orXAgMyu3SeCUzpwLgu0vrS
727XwY2Ln+pFKjTwWLNbPGzZQb6lUpkM8J1UFkd+HtlhYH//7/5P+zhHyc=
728      "]]}}},
729  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
730  Axes->{True, True},
731  AxesLabel->{None, None},
732  AxesOrigin->{0, 0},
733  DisplayFunction->Identity,
734  Frame->{{False, False}, {False, False}},
735  FrameLabel->{{None, None}, {None, None}},
736  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
737  GridLines->{None, None},
738  GridLinesStyle->Directive[
739    GrayLevel[0.5, 0.4]],
740  Method->{},
741  PlotRange->All,
742  PlotRangeClipping->True,
743  PlotRangePadding->{{
744     Scaled[0.02],
745     Scaled[0.02]}, {
746     Scaled[0.02],
747     Scaled[0.05]}},
748  Ticks->{Automatic, Automatic}]], "Output",
749 CellChangeTimes->{
750  3.660863951837199*^9, 3.66086418988871*^9, 3.660864234984974*^9,
751   3.6608642847860107`*^9, {3.660864563230816*^9, 3.660864583696065*^9}, {
752   3.6608646208487167`*^9, 3.660864666081502*^9}, {3.66086470383584*^9,
753   3.660864721673966*^9}, 3.6608649263607283`*^9, {3.660865011654098*^9,
754   3.660865035413919*^9}, 3.6608650965792933`*^9, 3.660865140220324*^9,
755   3.660865495727442*^9, 3.6608655280829067`*^9, 3.6611075833364477`*^9,
756   3.6611079275712767`*^9, 3.66110796841026*^9, 3.661108437603694*^9,
757   3.661108526714575*^9, {3.661108636236356*^9, 3.661108665059026*^9}, {
758   3.661108726467654*^9, 3.6611087405687923`*^9}, {3.661109731094899*^9,
759   3.661109746644104*^9}, {3.661109806974688*^9, 3.6611098480946217`*^9},
760   3.661110722089058*^9, {3.661153760672658*^9, 3.661153783021706*^9}, {
761   3.661153842835071*^9, 3.661153870322942*^9}, 3.661154421387599*^9,
762   3.6611544804747353`*^9, {3.661154532362813*^9, 3.661154605164646*^9},
763   3.661154809583425*^9, {3.661173581577107*^9, 3.661173714766621*^9},
764   3.661173805598506*^9, 3.661173845774172*^9, 3.661173882031178*^9,
765   3.6611739211174717`*^9, 3.661223970094832*^9, 3.661224003132491*^9,
766   3.6612240909079523`*^9, {3.661224161822092*^9, 3.661224186608893*^9},
767   3.661224219783635*^9, 3.661224426524351*^9, 3.661224780968639*^9,
768   3.6612250831963177`*^9}]
769}, Open  ]],
770
771Cell[BoxData[
772 RowBox[{
773  RowBox[{"residuals", "=",
774   RowBox[{"Table", "[",
775    RowBox[{
776     RowBox[{
777      RowBox[{"simdata4", "[",
778       RowBox[{"[", "i", "]"}], "]"}], "-",
779      RowBox[{"nmolec", "*", "binwidth", "*",
780       RowBox[{"gauss", "[",
781        RowBox[{
782         RowBox[{"xvector", "[",
783          RowBox[{"[", "i", "]"}], "]"}], ",", "50", ",", "sigma"}],
784        "]"}]}]}], ",",
785     RowBox[{"{",
786      RowBox[{"i", ",", "1", ",",
787       RowBox[{"Length", "[", "simdata4", "]"}]}], "}"}]}], "]"}]}],
788  ";"}]], "Input",
789 CellChangeTimes->{{3.661107996422367*^9, 3.661108004591732*^9}, {
790  3.661108051967945*^9, 3.661108120407748*^9}, {3.661108153731009*^9,
791  3.6611082099449053`*^9}}],
792
793Cell[BoxData[
794 RowBox[{
795  RowBox[{"residuals2", "=",
796   RowBox[{"Transpose", "[",
797    RowBox[{"{",
798     RowBox[{"xvector", ",", "residuals"}], "}"}], "]"}]}], ";"}]], "Input",
799 CellChangeTimes->{{3.661108235523424*^9, 3.6611082543762217`*^9}}],
800
801Cell[CellGroupData[{
802
803Cell[BoxData[
804 RowBox[{"ListPlot", "[",
805  RowBox[{"residuals2", ",",
806   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
807 CellChangeTimes->{{3.661108263102892*^9, 3.661108273467305*^9}, {
808  3.661224819324401*^9, 3.661224826739822*^9}}],
809
810Cell[BoxData[
811 GraphicsBox[{{}, {{},
812    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.01388888888888889],
813     AbsoluteThickness[1.6], PointBox[CompressedData["
8141:eJxVV3lQlOcZ/9QBrbG6RkS0oosHIlFcDpWbB1iQa0/2EoysyykUXRWmxWjy
815YRqqlNBVEyTKJKt4cajgHacNWwcrTgysxIFaTbIqsTQqbox4ANq+7ffbP2Bm
816h/nN7/29z/sc7/M+n49pozp3LMdxHez3v//Cnz3m7Oc/fNs/FGkTcH/MxFNX
817G2fnEfCzmI7TftdDdsUBj8QcnvOfmr8fzicBj6WSD09nj99fCOxOy8WLP920
818OwXrJ9Ks+C83ld6VAU8m94XDJ8M7lMBTKVmar/d8ngf9NCqsbz631UsLfjpd
819qPXsqjigB55Bm//g7XNsshHrZ1JA7+lh75J3wf+GJAu+aDjSoAfvTU3b028G
820JZnAzyWp/telB5asEnheTD3irJEerRLYh9KeTZMd/zoWeB7ZvtScb/PwFfT8
821fHKYio+WXPcEXkBc1rY3lR7zgReSM2hk7L47CmBfknyT9lZHBvzhF5H9Ly17
822tXk52N+PpPN/WHPnYA74xcSFX15lVcSD9ycu8lb7zjMDbQJ+h+w+96885OKx
823fgkd77qk+qTftd9Sap43YpN87dovgERXAlJ2DUrBLyOjV4FH/TSpwJOEOo7N
824lh+oQDxtErpY/Tqi7F/QUyDZxmj69hxJAR9IfGbCVFHYKvBBxN+sqfeqdOkZ
825fpTzXa18IfhgmpCvfVFUjPPZgkniYa781Yo88CEkWmd8VbVFDT6EjLf/aMp4
826BEzLiS+7WVs0WS2sty0nR9SW+n1T5dCvIKusVD74sQH8CuLGTB8c/sof+pVk
827nzEypaPLdb6VZK6uG1vQZoQ+lLif7l+sMM0Az7ChLDtM9yZG4MPIeG+J/7WX
828GuwfRtY38g+fD6DeKZy4TsXT/sFo6MPJ+n1tRKCXCvtHkKX4lLqkJRd8BBkN
829hTen7l4DPpI4n9vvJe5fDT6SbGfH3D9VDExR5NBvX9HTjnzYooizLJ5bECUC
830H018SF1spzIJ+miyyt0C6upd/sUQt/vUpyc6g6BnWJZnd7oh/2IiLvpthbo9
831QMBGIhqWtV9rzBL2sxIZj+53r0/NF3gHEd9aHT15MBZ6dk92hNyL74J9I8MX
832087GmlG/1ljid7R/8tB7NfQMux+yNHfjvonjiFdd+rY8ch30ccQl+ZXPfJYm
833rLcyvro/7Ok5xM/BcOv1F60S3Gcxs6N/672hW6k4fzyRh81e8gHybWV40vHT
834uc0F0MeTLeVGY7hHLvRS4h/39XiWuexLWX7GjJvVvwX22b35prvprw3Iv0NK
835tqClgy9zsb84gawT802HL5fCfgJxDeNrq7eg31kTiD/4pE7shXpxMP7E3jlt
836oYmwn0hc7I0HvbUboU8k25P1yyY9hX9WxiuS2lefy4SeYZlv34ZyxFvM+pjH
837nM/1L32hZ7jxwaWFueugZ/hJ+PxK2XroVxHv/+9tfwvZAH0S8aUfvRucux56
838tu+RKd6vFK78M94to+76Vpf/jM+qW1krTcb5k4m3B/ad8U5H/JLJNnapPO+7
839AthPJi5lbvjTuWmwz/hLl09+dgjxFbO+ou1+tWO9CnqGvy9N9MhE/VtTBH/c
840XPWXQnxU5KRBawbOn0qcPSduwJODPlWI1yD6h5XhWENc1wXUo4PhtoCZfqvw
8413ojTiB4evrpMngf/04R81ZqgZ/idkTOhFTro2frpne1JCYU4v4z43Iyffv7T
842GtiXEVdQs/axDO+BVSbUkx/8dTD+/N7MzD4Z7MuZf6UN/mvDoJeT2F1zK+jo
843ZtiXE9/7wVrP2EzoGW76KOr1b/EeihVkMzwYNkz4Hc6vwH3Jhl4h3KcSvEcO
844hutLjo0rgz8iJdE4VZJuF/yVsHv5/HGzameggJUMB0wrFkUiXmaGJ3TvnBK1
845SMAWpdAPStFPWpRCvygqErBdSdbXvvdypuO+ONn6P7fv9WjFfCFSEe8bXPAy
846LkPAElYHG5WF6sAI2FeRsXXoRvYa+GNm/L6r27LXRgjrLQyHFjlT7yLeLSqh
847Xy5A/7Mz/T/5n3t7XPbZ+rjzL6pCkuG/WujHU9GfJGoS3157pN+BfCoZ75Pf
8483rkI9WRWE//L3sitUQb4rxbeAzn6U4taeC8uIl92hn/8+EdDGOrbyfRDoScP
849/gP5FqWTY3aYm/19HeynE2U5Y8QxCthPJ56vHNqehffFnE72o+HHhtzRHyzp
850JL5ctaZgAPNKC7uHvxwP7J2NfNiZfo+jtEa+UsDOdJJ6GaLvV22CfY3wHuuy
851YV9Dos7Ggd5ZqA+lhmy7ZLvPGCJhX0OcOrPa5If+Y2H4QuU88eZgxJ/h8crd
852qTXzBN6uEeaJRei3To0wb/hifhRpiV9nVVXJF8C+lnaaZhxp34r4KrXCPNOE
85385m1mHcQb4tWmIf4JNjXkujAozuBbgmwryXO++58+zHkz6kV5q1y7CfSkbVs
8543aKmNMyLEh3mNfQ7pY5+P/P98tQNmFfNOrLN+KJt3pVExF9HFvm15iuf+SL+
855DN8ofNs+hHnLrhPmyVRgp46aY796eaLJVX96OtRUcUlWg34n0Y+eV5V68trT
856dbO7BPVr1o+edy166psQESf3Rr226EfPy3b96HnaqRfmbZ0J9g2j53GJYfS8
857rjSMnufNBmHeH8J7YjHQ/78HcoBbDDTqe8FuoFHfE04Dub43/guA86fu
858      "]]}, {}}, {}},
859
860  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
861  Axes->{True, True},
862  AxesLabel->{None, None},
863  AxesOrigin->{0, 0},
864  DisplayFunction->Identity,
865  Frame->{{False, False}, {False, False}},
866  FrameLabel->{{None, None}, {None, None}},
867  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
868  GridLines->{None, None},
869  GridLinesStyle->Directive[
870    GrayLevel[0.5, 0.4]],
871  Method->{},
872  PlotRange->{{0, 99.75}, {-368.52172460536167`, 298.3874142915902}},
873  PlotRangeClipping->True,
874  PlotRangePadding->{{
875     Scaled[0.02],
876     Scaled[0.02]}, {
877     Scaled[0.05],
878     Scaled[0.05]}},
879  Ticks->{Automatic, Automatic}]], "Output",
880 CellChangeTimes->{
881  3.661108274612281*^9, 3.6611084376751842`*^9, 3.6611085267765303`*^9, {
882   3.6611086363252697`*^9, 3.661108665111405*^9}, {3.6611087265203733`*^9,
883   3.66110874063727*^9}, {3.661109731163712*^9, 3.6611097467059517`*^9}, {
884   3.661109807055561*^9, 3.661109848161718*^9}, 3.6611107221317453`*^9, {
885   3.66115376075438*^9, 3.661153783077216*^9}, {3.661153842884677*^9,
886   3.6611538703745747`*^9}, 3.661154421464151*^9, 3.661154480592321*^9, {
887   3.661154532446109*^9, 3.661154605279621*^9}, 3.661154809669935*^9, {
888   3.661173581648666*^9, 3.661173714828924*^9}, 3.661173805660392*^9,
889   3.66117384586515*^9, 3.661173882115933*^9, 3.6611739212054996`*^9,
890   3.661223970175799*^9, 3.661224003248719*^9, 3.661224091024156*^9, {
891   3.661224161940206*^9, 3.6612241867287703`*^9}, 3.661224219900621*^9,
892   3.661224426639511*^9, 3.66122478108354*^9, 3.661224828087525*^9,
893   3.6612250833130817`*^9}]
894}, Open  ]],
895
896Cell[CellGroupData[{
897
898Cell[BoxData[
899 RowBox[{"Min", "[", "residuals", "]"}]], "Input",
900 CellChangeTimes->{{3.66115468011102*^9, 3.66115469553443*^9}}],
901
902Cell[BoxData[
903 RowBox[{"-", "368.52172460536167`"}]], "Output",
904 CellChangeTimes->{
905  3.661154696741271*^9, 3.6611548096906548`*^9, {3.661173581690281*^9,
906   3.6611737148576736`*^9}, 3.6611738057071657`*^9, 3.661173845894189*^9,
907   3.661173882134103*^9, 3.6611739212240963`*^9, 3.6612239702293158`*^9,
908   3.661224003281253*^9, 3.6612240910564203`*^9, {3.6612241619712152`*^9,
909   3.661224186761483*^9}, 3.661224219928742*^9, 3.6612244266719*^9,
910   3.661224781119219*^9, 3.661225083346334*^9}]
911}, Open  ]],
912
913Cell[CellGroupData[{
914
915Cell[BoxData[
916 RowBox[{"Position", "[",
917  RowBox[{"residuals", ",",
918   RowBox[{"Min", "[", "residuals", "]"}]}], "]"}]], "Input",
919 CellChangeTimes->{{3.6611547020972357`*^9, 3.661154719103921*^9}}],
920
921Cell[BoxData[
922 RowBox[{"{",
923  RowBox[{"{", "125", "}"}], "}"}]], "Output",
924 CellChangeTimes->{
925  3.661154719677985*^9, 3.661154809721857*^9, {3.6611735817404346`*^9,
926   3.661173714878026*^9}, 3.6611738057618313`*^9, 3.661173845916881*^9,
927   3.661173882170998*^9, 3.661173921257353*^9, 3.661223970370145*^9,
928   3.6612240033198643`*^9, 3.66122409109451*^9, {3.661224162008938*^9,
929   3.661224186799014*^9}, 3.661224219969226*^9, 3.661224426711216*^9,
930   3.6612247811736803`*^9, 3.6612250833840313`*^9}]
931}, Open  ]],
932
933Cell[CellGroupData[{
934
935Cell[BoxData[
936 RowBox[{"xvector", "[",
937  RowBox[{"[",
938   RowBox[{"First", "[",
939    RowBox[{"First", "[",
940     RowBox[{"Position", "[",
941      RowBox[{"residuals", ",",
942       RowBox[{"Min", "[", "residuals", "]"}]}], "]"}], "]"}], "]"}], "]"}],
943  "]"}]], "Input",
944 CellChangeTimes->{{3.661154729523085*^9, 3.6611547634745407`*^9}}],
945
946Cell[BoxData["62.25`"], "Output",
947 CellChangeTimes->{
948  3.661154764779113*^9, 3.661154809741576*^9, {3.6611735817934093`*^9,
949   3.661173714908766*^9}, 3.661173805806725*^9, 3.661173845962756*^9,
950   3.661173882217852*^9, 3.66117392130685*^9, 3.661223970412418*^9,
951   3.661224003370822*^9, 3.661224091146294*^9, {3.661224162059444*^9,
952   3.6612241868500433`*^9}, 3.661224220019877*^9, 3.661224426761208*^9,
953   3.661224781223856*^9, 3.6612250834341297`*^9}]
954}, Open  ]]
955},
956WindowSize->{808, 752},
957WindowMargins->{{Automatic, 14}, {Automatic, 48}},
958FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (December 4, \
9592014)",
960StyleDefinitions->"Default.nb"
961]
962(* End of Notebook Content *)
963
964(* Internal cache information *)
965(*CellTagsOutline
966CellTagsIndex->{}
967*)
968(*CellTagsIndex
969CellTagsIndex->{}
970*)
971(*NotebookFileOutline
972Notebook[{
973Cell[CellGroupData[{
974Cell[580, 22, 198, 4, 28, "Input"],
975Cell[781, 28, 1260, 18, 28, "Output"]
976}, Open  ]],
977Cell[CellGroupData[{
978Cell[2078, 51, 409, 7, 97, "Input"],
979Cell[2490, 60, 965, 14, 28, "Output"],
980Cell[3458, 76, 965, 14, 28, "Output"],
981Cell[4426, 92, 964, 14, 28, "Output"],
982Cell[5393, 108, 966, 14, 28, "Output"],
983Cell[6362, 124, 967, 14, 28, "Output"]
984}, Open  ]],
985Cell[CellGroupData[{
986Cell[7366, 143, 122, 2, 28, "Input"],
987Cell[7491, 147, 679, 10, 28, "Output"]
988}, Open  ]],
989Cell[8185, 160, 250, 5, 28, "Input"],
990Cell[8438, 167, 163, 3, 28, "Input"],
991Cell[8604, 172, 603, 19, 28, "Input"],
992Cell[9210, 193, 195, 4, 28, "Input"],
993Cell[9408, 199, 496, 9, 28, "Input"],
994Cell[CellGroupData[{
995Cell[9929, 212, 99, 1, 28, "Input"],
996Cell[10031, 215, 1650, 22, 216, "Output"]
997}, Open  ]],
998Cell[CellGroupData[{
999Cell[11718, 242, 287, 6, 28, "Input"],
1000Cell[12008, 250, 1651, 22, 216, "Output"]
1001}, Open  ]],
1002Cell[CellGroupData[{
1003Cell[13696, 277, 252, 5, 28, "Input"],
1004Cell[13951, 284, 3566, 53, 250, "Output"]
1005}, Open  ]],
1006Cell[CellGroupData[{
1007Cell[17554, 342, 531, 15, 28, "Input"],
1008Cell[18088, 359, 3918, 55, 250, "Output"]
1009}, Open  ]],
1010Cell[CellGroupData[{
1011Cell[22043, 419, 131, 2, 28, "Input"],
1012Cell[22177, 423, 673, 9, 28, "Output"]
1013}, Open  ]],
1014Cell[CellGroupData[{
1015Cell[22887, 437, 158, 3, 28, "Input"],
1016Cell[23048, 442, 1117, 16, 28, "Output"]
1017}, Open  ]],
1018Cell[24180, 461, 572, 14, 28, "Input"],
1019Cell[24755, 477, 241, 6, 28, "Input"],
1020Cell[CellGroupData[{
1021Cell[25021, 487, 376, 7, 28, "Input"],
1022Cell[25400, 496, 1173, 18, 33, "Output"]
1023}, Open  ]],
1024Cell[CellGroupData[{
1025Cell[26610, 519, 1004, 21, 46, "Input"],
1026Cell[27617, 542, 13105, 225, 233, "Output"]
1027}, Open  ]],
1028Cell[40737, 770, 714, 20, 46, "Input"],
1029Cell[41454, 792, 246, 6, 28, "Input"],
1030Cell[CellGroupData[{
1031Cell[41725, 802, 248, 5, 28, "Input"],
1032Cell[41976, 809, 4452, 83, 230, "Output"]
1033}, Open  ]],
1034Cell[CellGroupData[{
1035Cell[46465, 897, 128, 2, 28, "Input"],
1036Cell[46596, 901, 496, 8, 28, "Output"]
1037}, Open  ]],
1038Cell[CellGroupData[{
1039Cell[47129, 914, 199, 4, 28, "Input"],
1040Cell[47331, 920, 505, 9, 28, "Output"]
1041}, Open  ]],
1042Cell[CellGroupData[{
1043Cell[47873, 934, 337, 9, 28, "Input"],
1044Cell[48213, 945, 459, 7, 28, "Output"]
1045}, Open  ]]
1046}
1047]
1048*)
1049
1050(* End of internal cache information *)
1051