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