(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 1457559, 24096] NotebookOptionsPosition[ 1455251, 24017] NotebookOutlinePosition[ 1455860, 24039] CellTagsIndexPosition[ 1455817, 24036] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"xRange", "=", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1000"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tRange", "=", RowBox[{"{", RowBox[{"t", ",", "0", ",", "200"}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.405075396121993*^9, 3.405075396284467*^9}, { 3.40513870534182*^9, 3.40513871140376*^9}, {3.406292480635222*^9, 3.406292695602671*^9}, {3.40629315163052*^9, 3.406293154001243*^9}, { 3.406293301763631*^9, 3.406293350521651*^9}, {3.406293451589776*^9, 3.406293453636269*^9}, {3.40629351459724*^9, 3.406293516723909*^9}, { 3.406294045198205*^9, 3.406294049541683*^9}, 3.406361624186737*^9}], Cell[CellGroupData[{ Cell["Reaction Term", "Section", CellChangeTimes->{{3.406361638579706*^9, 3.406361640976664*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Q2", "=", RowBox[{"d2", "*", RowBox[{ RowBox[{"(", RowBox[{"IP3", "+", "d1"}], ")"}], "/", RowBox[{"(", RowBox[{"IP3", "+", "d3"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"wInf", "[", "c_", "]"}], ":=", " ", RowBox[{"Q2", "/", RowBox[{"(", RowBox[{"Q2", "+", "c"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"mInf", "[", "c_", "]"}], ":=", " ", RowBox[{"(", RowBox[{ RowBox[{"IP3", "/", RowBox[{"(", RowBox[{"IP3", "+", "d1"}], ")"}]}], "*", RowBox[{"c", "/", RowBox[{"(", RowBox[{"c", "+", "d5"}], ")"}]}]}], ")"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"VIP3R", "[", "c_", "]"}], ":=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"mInf", "[", "c", "]"}], "*", RowBox[{"wInf", "[", "c", "]"}]}], ")"}], "^", "3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f", "[", "c_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"(", " ", RowBox[{"fiVLeak", "+", RowBox[{"fiVIP3R", "*", RowBox[{"VIP3R", "[", "c", "]"}]}]}], ")"}], "*", RowBox[{"(", RowBox[{"cER", "-", "c"}], ")"}]}], "-", RowBox[{"fiVSERCA", "*", RowBox[{ RowBox[{"c", "^", "2"}], "/", RowBox[{"(", RowBox[{ RowBox[{"c", "^", "2"}], "+", RowBox[{"KSerca", "^", "2"}]}], ")"}]}]}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.405075396121993*^9, 3.405075396284467*^9}, { 3.40513870534182*^9, 3.40513871140376*^9}, {3.406292480635222*^9, 3.406292695602671*^9}, {3.40629315163052*^9, 3.406293154001243*^9}, { 3.406293301763631*^9, 3.406293350521651*^9}, {3.406293451589776*^9, 3.406293453636269*^9}, {3.40629351459724*^9, 3.406293516723909*^9}, { 3.406294045198205*^9, 3.406294049541683*^9}, 3.406361624186737*^9}] }, Open ]], Cell[CellGroupData[{ Cell["Parameter Values", "Section", CellChangeTimes->{{3.406361638579706*^9, 3.406361671695335*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"ParamRules", "=", RowBox[{"{", RowBox[{ RowBox[{"fiVLeak", "\[Rule]", " ", "0.004"}], ",", " ", RowBox[{"fiVIP3R", "\[Rule]", " ", "20"}], ",", RowBox[{"fiVSERCA", "\[Rule]", " ", "1.2"}], ",", " ", RowBox[{"cER", "\[Rule]", " ", "1"}], ",", " ", RowBox[{"KSerca", "\[Rule]", " ", "0.15"}], ",", " ", RowBox[{"IP3", "\[Rule]", " ", "0.7"}], ",", " ", RowBox[{"a2", "\[Rule]", " ", "0.2"}], ",", " ", RowBox[{"d1", "\[Rule]", " ", "0.1"}], ",", " ", RowBox[{"d2", "\[Rule]", " ", "2"}], ",", " ", RowBox[{"d3", "\[Rule]", " ", "0.2"}], ",", " ", RowBox[{"d5", "\[Rule]", " ", "0.2"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.406293082091298*^9, 3.406293085577598*^9}, { 3.406293158914969*^9, 3.40629328813798*^9}, {3.406293984699964*^9, 3.40629398480679*^9}, {3.406294494794796*^9, 3.406294494904505*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Plotting Reaction Term ", "Section", CellChangeTimes->{{3.406361638579706*^9, 3.406361689473212*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"f", "[", "c", "]"}], "/.", "ParamRules"}], ",", " ", RowBox[{"{", RowBox[{"c", ",", "0", ",", "0.6"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"{", RowBox[{"Thick", ",", " ", "Red"}], "}"}]}], ",", " ", RowBox[{"Filling", "\[Rule]", " ", "Axis"}], ",", RowBox[{"FillingStyle", "\[Rule]", RowBox[{"{", RowBox[{"LightGray", ",", "LightOrange"}], "}"}]}], ",", " ", RowBox[{"Frame", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"FrameLabel", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"Style", "[", RowBox[{ "\"\<[\!\(\*SuperscriptBox[\(Ca\), \ \(\(2\)\(+\)\)]\)\!\(\*SubscriptBox[\(]\), \(i\)]\)\>\"", ",", "Large"}], "]"}], ",", " ", RowBox[{"Style", "[", RowBox[{ "\"\\"", ",", "Large"}], "]"}]}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", " ", RowBox[{"{", RowBox[{"600", ",", "400"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.406293365930523*^9, 3.406293404067813*^9}, { 3.406293993817085*^9, 3.406293994646769*^9}, {3.406303237901023*^9, 3.406303239814268*^9}, {3.406303341346028*^9, 3.406303380979755*^9}, { 3.406303441244811*^9, 3.406303594704716*^9}}], Cell[BoxData[ GraphicsBox[GraphicsComplexBox[CompressedData[" 1:eJw9mXk8lO/3/6VSoRJCylqWRNmXoiO8i8q+p5ClrNlTqexLkWSZmbKkkn3f qVzKLjQz99iF6C2JiBYh/a734/H7fOef+3E/rpl7rnOu13me15kRcfAycWZm YmKq2MDE9N91aYznSOWZ08eD3ye903fwBgVWVlKWYAw0ins8yxAwRWm075vN RB9CWoe67lRsFGI4V7o1FzyFG7t6kj9vSEIyh40/1uTkQshIyDeRPhL6KcUZ HC5WAvbV+Xs/LaSgE49ZTRQflQPFTGW5XSUBaWoL/kicqAKm8xOr2el+aPrs fpGMvXWgue3XKKU7EhYv6nEqC7wE2fYVV0/NdFB8sIG73g4B9b1LkaJkPsw5 e3o6K72GmD0jzcXMFcBdbbJbp+INWNWfq98lWQdWXKEtFONmsLowt032LwL9 nzCgwGgB9Y8mDspHmsB26tSSTmAbxJ032zoR2wo6d9pNCjk7QDOlgX6ouAOG ywa5zyV2wvJIcOG5S10Q/0lh7JZ8F0Ts1pBN5H4Hi08F2aeru+HRV+/H88+p 4OtucqLd5h3EO1v2zBTRwNEv+ELwTSowWyz9Hgmjw79L6W81NGjA8spsdpMQ AT9WxjVUdtDharVu0/WHBIzqeHwwpdFBi1Ov48cYAXLTCjmpCQRQb2btiVsj gJNT2EfQngGsZYqfdOYJ6Lpt5tEs0gs3dzk4UuoJ4PYtcn++1AsjX37Tz1oS 0Nqgda/sdR80/qgqZummg7GB8Ph4aD9sbefU8ztIBzcnUZUzJgMga1esZW9G A5tXHmiAexCSPX+mTHtSocZryfUBMQi6KGbKe7wHnuRm8EU+HYLpTOW1N8Fd 4HggcTnVcRioVbLcbLs7oUI8f2ROcQSWn777y6LdBiPurIsOyyMgW9ijbenf DF6sOQc31r6HhCMtL2zvvcZXF+fR6FFobFNJlA95CXfJPvQRvTFoXHljNsNZ DaxN70TYdoyDsCqTY355EWz6FOGuETAOmiXDZTmcGeAmPaof2jsOTAtKSXyW dsBbpVc+ofEBNDlEeF9opyNZO+sS88cfQLho23UivAitK0VrjK5/APsJ0uMp mWrEZdSmH3FpAoRXE7W+nWtAKbEFP4+1TYBTqPGJLtYmZMv6VGSX7CRENPxK //dkG5I4mSmtoDAJThWZKlcvt6ED+i1HswWvQbzJhdP9vx3B5stNSYKSCEdz +aV7W66hpz3Fnl/VH8MBm8821znvI6bOMImGvOdwrX12C90vBSUWlpWlTxVA jJUtEtxFRrZDHVuDeMqAyTSl/X1eEtL0SX7meLwSxjup5NFLUYi7XTiWL7gG QjieVztlOYKnZfUet7R6aDxcrnu3OAnko2vFre70AIvSuxOvjtAg3owbye2m Qubeb9vUBOmQVRGfPDpChfO//Qv39dHBMxBF1D+nwXj1hc1WDgTMPM6+3+JF h8NTH2D1NQExPEbNv2QJ2GHJnKiyRIDon4PsF5axzljcTUmrBHjp1l/51siA TY78spWDBHgvM81W3u2FtD2a6cbxBHjyXbtZcK4PGGou53l4CYiLKlPuEOmH LtbgX2eD6WBr/YSNZ6YfEgYcO92baKAlKSLxVG4SBkTVnshYtqEtxLGIPEEX sNW90tat4gn+C6F+YrnxsLZ19/eNOz2REJWm9uRIGpg05tn6/b6LZsh79V5L ZcFVgahiJ5tkZO+54a3UqXxwYg9llAiS0WzgQacMs1KQvJ3VLLiUjDRj3q2c 8aiAkIIE170oFuVHCjX2llYD0xtyfu2AOZKiewhsGqnD+d3MrmwfD5wp7nuE hXuA2mpq+c8sFU7dfTW2u+kd8HkM+UV8okEE15bfxYVUEMwJH2stokNe1WHJ mz40cIAT84paOD/dTZfWj+O6k+JLYC8hoEJmlC62Adf1d1PRvhkC8heLrW82 E3CDyMgL/kOApszcfHQCA6LI27pL/yUg44VUobpFL/y8s/CO+hyvt3D6S0r0 wVma+oi3KgF5oisXX33rAyfpmjNyOXQYIMYDVqv6QThB6kHyLxoYPbyYNYj1 7KIlfW+bYRtysKoeXOUJg9qyl4GyN07DzVFVctL7FAh8G9j0/E4oSpxqoQqI PYFy+Yavjx0eIJ61z5tfSOaA0eYtJyR1Saj5y0UKR1QRULO8HDMLSYiq7lZk OVAGIbzDB97uTkSa56zPBnFUQYic4vUitxC0PKkXP6BLheajWTRmEzo0m/LI 7thOg4/y+oEjq3SQUzP6sDxIgwG49ehKMAFOvndq4zLpUPhAcY8unYDzf19X ddhjHtplZAj8JuDRlucNGtIMGAkt+9v7E+fHMXru3DcGLGb4vzF7S4CWUv8h GdQLtTVuo6L+BCjWjOe5yWN93ds5ZGzXhrpHJXlLBa2g0fGFTBfVC8z+tGxy dowFLeP2H85CjmjUaPLxSu4jYKHMN/x8G4PG1j8Nark/g5WATbfuVyehVv08 S5/APEhb8VHm30xGl++V77f/VALCAj7r14NSUCOXp7jmxgqguqalR0vHo4CP Gecq9bC+HqVWGl5wRXd26grRgnC/e/2+0WDbXWDe18IZ8aUb1mOjQwI6qDCy Es3hdf8d/By5/Z7rLQ38tbUPkZOpcLRH+kwWhQ5tZlO/z1vSgINVN3Yd16m8 kwn9kTgdiuviu8lZBOiqiQ2pf6bD0LjB1mqsnyuePkxx+QQkH84dvIj1lfc4 Stb0GgOGZl9ot2P9KaVDHJN6LzzSnWaPK8V63Ds+f2l7H0i/+iB1RJeAWu+M wJrBPuCLHdWzracD9wOmQ19T+4H7sBHq5qBDP22DYwDWV/OWnwm/Treh7+t8 vwzEb4P6iewOZT8bmP04ZKdwLBlO331xMG/hFlrM7gnvvJgJJrWXWbtuJiA2 x2qSums2DMx3ReROpaBDXdvELKsL4dosx1DzZxIis/HbuAaWwYKb5HuRz4mo 8d+fgV45lRDCFJTKfCscae8rEd98hApPbfgfBanSoTLcxOHnd5zH091uhjgP hgNbKgYRDabzUuk6PgRIbOM7cDCGDrpG7xsCOglYZN32ahzHuZxg1HEd68k2 7LKzPicDKGphOQlYbyx9wofFBxlg9fHKVSrW4yDPLNu3rF5gt7ArvB9KwIv5 EBcWrC9pbvY6sGlDZqk+t//w+AL3UMLhMQVX+CFrWAszWcATFVizOJ+MfsdZ H//IUwA3S7VNX4iSkVje9eq4qlLInJWx/WiZjIqUlVYvWtJhRbJYb6CWAEFR 8tsdewjgJ96EWOD+PxU4W+4xSmC9lR+fwv5goN605GUuAwbsJhX/Yv9wWm5c hcofCeP2SZ2OG/Yj0RSBALXpHBjnz7LeQSGhC2z1bJ3mxcC0JHCzwJqEngyv nF2opUOibW2OL+Z569KN4vJrBNhXuMfvx7zvi5vgPXwS80glSuwy7gfFtzW/ 9OB4datPbHJybENWP5k1XgnqQplmBjn1hDc0b3Tp3vT6DoxcTlHfkXgeaeya Ncn59yG4eU2OcRVHo3PyBk/D2Z5BctkK5bBtElrR+JOevDMPvsu+iz21QEKH zzvQfKJKoLSDvPt+bQoKGTymXt9TDgmnbiXbmd5HQ0mjOft5qqFxt9YBtyPe SLVYKbPmDOa1luOLjdPRcGPMZ6cyrRtM2lPVOmqoIKEwl5sagP2ZV6Al6SUN PGwu6sTGUKHz5Hal3nt0sDjwxkH/NA0aZ2NF+SVxnZTs2erKT4ejL1l/CDwm gAu1X9B9T4eZvndLnycwh/LpRFI6AazRTcw7cD2NWhFUVQ8G+GeOrOyZI4Bp x+mbE4d74bZE5YHcKgLGWltlopj6oIsyvZ9uSAAie4dVdvfBdL2KQFMTHcq9 EiI/xfeD01vEqiBAB/6vxlrauJ4OtBmSvuq2IXbxecoN0SDI3CL5ueSXPRyY PqSz0S0JRjb7PEiyCUKiW0zfCzNlwuziiaqiTQkoirTP7ta2bFjzFYp0LElB ti8vuN8wKAT/bVPXpJjIaGXxm9ziyTJoryycYZxJQiEFLfc++1YCU7GJ7tZf Ecil9oPvXlEq/Fw4X+FzCHM41TtoZpoK+1b3o84xOozmTO7tLKfB6DbUM+WK /YLZmPWvIDroWGjf0W8l4FVCoeasOgF8ZhGhs98JaETtquabGfCy+vzB/BUC jt737GHrZgCLja6teR8BWdE/V+nkXvj+cZaDJYaAdTbaKQb2A80eawXK1m1o IOK+pon4FVDIj+Bst3MHb6eqkJC7WZBd98H5RU4ykr5VFLGlJB+sVGLnXPaT UUyC8b3dkaUw7TJkO4LXd4StycqeoYNRNvXOlkoCpnOOGvJvJ8BAWUlLCp+X bpV+lD+NgJd/V9/p4fNk9880fJrGAO7ez2pf8Hn36bzdOkUJh4QHk7N/zylD Wnu4xo2YHLgZ/K+uiS8JJaxuN5mbK4J28rfd45EkxCzxJqyhkA6ZdeJmi70E PLguKN/gQcDKdc0zijh+1U7tL0pqDHhqKMSu9IOAb48SOtNwPY1/Zv567mIb UhGZbmPidQCp0wKsMlVXIOmC+S6DsTyII2VsT+YmIyO7DcckDpRC+wKHYZ5a CrpdEWoh9wOfi8J1rnPTeP+yHDIRNQQMmfpVRuN4Eo4Y9ay2Z8MBQ9Vc8h4S Wh6uGlUQLAK+T1euz3aTkESks1qvGQEHAu9eG/tFwJoOS+I/Qgxg8n80pYd9 n7p/6giNPwA2AqUuvP4StKZtXm68XQBp7PcePNtHRhFHY8ecvuC4PqpvN8c8 yr1YU7h8Khc+5rQ+N68jIam9l1qyXxWDUUzcDLMkCa159kvlRBCgk2V/7Qvm i79hqI64KQM8ntGcPb4R4OawwLYB++3mmrf7PJ3bkJdYTUy/oCpItLCS/7H2 hjWWujLGfAwkB5exqFEtkYePd+Gf8Icg3rJ+K0QgGvGeDWXJGn4KidxnY6Nl kxC3ssSJ7OFc+E56WjM6SUKJ2kdzV8xLILeJ8s6kNwXZq+8/FVhZDhy3kzR7 Ju6je+63DWzWqkAzm5E5PO2DEm4Vjucr/OcHh7Z/G4wC5RPt23c2d4OV8ieu iWIqnH0n9KPM5R1MD941PllFA5fLC0aRYVTo++h7zBb3GUF2WXKLNg1E3Tc8 SdmP4zb/I2rETYcRnUYJ4zTcd0IfCRn00yHOdveE/wcCQvTCX/GScR/uH6TR cT6vVf3OPHTpPz/YcEbxKwGXThr6m0r2wr6F/a/FcH+4sa++a/dKL2w1rYhK NsX+p2BwS0VbHxjRP+90bacDs92AvFNMP7SHdXVb7qeD+tpvOan/+LI8uGfm VBuiqA25JmReh4yPt/woIQ7Q6Gqm7L4hCVjltHjLOW4gvytS3JSMxxAn02X9 Iv4+0i41bEoffA7tbFJ1+xNT0LXK6B1x/IVgdnq4wGQbGRVOXUivPoz5QjnG uBKUhOx3FFGyrfA8g0IvHS2NREEsq5TsvVQo9a1VGj5Ah9ycv9EfJ6hwTtBg YnmIDtn3jY20CmnA17ocJnaJALWVN/6jAZi/30tyDmKf/NT0NvOiMp5PHPP+ 8mC++NheFx1cx7rRmZoRw3pifEnKWG9lQIbI5dhtA5jTjatL2g96QbZDSlEu FvdN/usqLzBfvFWNxeWt2tBCgfPELVF3OHCNM67MwgNkbk9/zDyXBVJ3rcS5 wpLRyYB+iT1B+fickldqRcjoaMSr5UGPUpA1fubG3ZOMkMT4Gd5/cL8QrDe7 X0ZAjxv7LuEtBFzVOnNk8ywBv03XvZaxz5hpUA5ww/WYTW0VSiExwL5XiDD5 SEBaj+a3hf4w0H2sn8v/WQusJ6qfPLLKAcpxkQedF0gof0nlx3od9tt1Pm1f H5JQRdfz9OfZdJB/beEvysDfo76w1ITz9DNNaLII+5XPLt8ZdfIMiBf7V4wN +xnOPNX6cMwX9k7kaGHfhjISbhJW4rZgbzdblWfoBYXmwZV2BXlAteVJDttO Rsymntlf2EpB98Wdv5I2KcjxqimfyjwdJFU1d2z/hPne08izFcfJ4nwo8gaO x0AvXH7fw2xolvltX8BEQpyMsdB/ZgpBl3uaK3qEhNh1VVaHDAjIfMGaJIf5 UjRZHdTBxwDZbRfnFjFf2qOjij5R/KD7tNL7ICYXGE8+kko/XQDsiryfzwuT kdLgsY5pnKeeuVCrZ7genINmOPbw5QJTgM1fzTwS0uFKcKm9Uwylreu5kSdI 6INs3mzhbbyvo1GecVgPYl/lDuefZcDZXhk7r0U8j3pTVj/hfHCbaSS4OLWh Nenec1t4TcAnzTCWttcbxv8hVsK082DZz+yJzSoJCbGLzuzEvjpXQaBEA8er 06zQ7hJWCKWK83PPf5DQThE+QxobA9gnR3ftx3wfLz5vyJXiDakVVbyfTN2A MzL8w4+JfFjTdZFDuD+RedLTFvpxn5PqeymAnzdWbDRoXZsDZnLBIQ24fxRW n/Wu9cX6jmEeLMXPE1rZJd7L7wwNf7vRCeYrcMP50BV1rnzQ9SQzGfCTkc0D O22mBjxnrhq73fnPj0zGHhM1wXzvG5j//go/L2vYrKS8APyVqEvxvGREpA79 GMec1XEfs1vF+Zk5FaE7EJgLVtlf7pxpJyGegs2uT+5iH1O/vUwT51vXTHcf P+bx8uBsns+lNpQeLiTMLiQJjqn8Etl23kCJ+fBMPToGdJySHnzzN0cRT+QH JU0fwpCGfBDBiEI/Bx6pVjY9BeqW46ct9yQhYfmXnCWvc6F5OqVF6j0+P5Vs LqRRAlstX/6rNZmCqHOkLvOccojh9G9/w52ALln2iATO43kyUM62r8kXqf/7 KSRLog7s9W8E3zONAm7fpbigV93Apyeq7ZZPBTsBtb0aDu/AQuzi34dlNLg0 p2IZFoz9XvrDyYxIzEOmy2WhmjTQnGPtNhAhYJ+kV6n2LjokkwLbpx4R4O4Y zmTMoAOHZD6pehz7a4pN+WASASdFGFOlOB8tk+zbxBxxf4w3ob7CPH44YbzA JtYLN/LbfOzrCDjVppD2/mcvpPl+mvcxJ2CTUvzp8uY+8I/ZOFPbSYd8i5t3 RCP7wcq2tPISntv0N/bkCWAeq7r0DV/GPLZtMjRJHbsGPNN71FTBEeLzXY4I dSfCbK0fh6b7dXQ3Y1azyRPz2CNVsMfoPhLVcP1R9Po5cLR93SoQnoIiS/eN JzPh+vOoVHLC9Rw5zkVQRMvAW9un4m1iEqIiYcHIM3h+Yq66MMoWhSKKtkh7 8FHhZZtXSqwI5q/MqfgPY1QoTwjgjBugw93wGLf1XBrUBnuc1nTCPOlWZKb7 Yt8cmNXM1UTAbqXLk98VsF5ZPaNy8XzQWWTMm4p1lbv7xl07fBXlz2FZbmLA jey1HeOYxxrn173W7vVCzJ9jjW/iCJAT0S8pxjz26CifLLFsQzeRUUVSpiuk pXd8YOn0AHY2Zt4yrSzIqkwnjVxJRqY+Hr9FLuZDrjE90g7zgeu23+GX9qXA J2thgsaTEWXFeJBVC/vW2+46hnie1Tz745voJszdGOnpK9inmDzrc6pvI6De oeLWVVwvS+uD0fFJDDDSOLRpI56Xg5OCnK8GhsFylXEWf81JUPWiEvmaOZBJ 0SzaZ0pCNfNdrVsyiyCkZkbF7xkJxT3mPfTwGR0aNE6RlAnMX4FBRqsj5tqv u2uXMI+1ars6go4wIL/hz4oP5vH+Unb5QMwfYcNeQYZdG6quMQmIELUBq8G/ UgFJXlBnzB3l8iAPkq2PqVzH/bT64cna1L8loDv3xCrQPQX1War/Vpulw1nK tbcOUwTwCD5k6ywioGuv4GMvHI90X5qTdGg2MAw4KruWUpAg1d/UgFoIAwsP 7htj/xN8UHjk/Rn8uWZipg/vh0XhvnfsbgYsnldV1Mb7jd+1lrLY7wv1N457 xKW5wGxCs8ewbAGUnv8VUIP7X7wkd0M+9imSih6qrbgedLZu0j7IlAvJFz7Y vc0kIaOdYkdfexVDCP2iXLM+CZlxJnaXBOH8rz1RtcZ6uK3But9DD/v/0a+5 s5jH6YEefsM4H0zSMqen8DwpmWvx21ncAK7WBZF8Zb1h7oRLV+zBPDB71nR0 DfP2ms3KKC2TgCtuFnoSOF7jxvaHPi6FkFme/ZoL89rFe0E7aSsD1vYe2uqC +XnWJ8tLPdcLnnz2veA64gZCoXwV6+35MG7x+ux/PF5wDnhQjvsoo+OXghJ+ XsMxNza3zBxQHWNMHLpJQoIUJeMXXgRw8+h3++HnFWz/83OG4giVRzVrq22v QNTFRvjndx50zZVsSMN8rRxXUX1Tj+vgusP6fzz+Mia8/ZBqEciuz5ruayWh ukxv8VpKATBFFpfqYH5LfymfeYo5Mn7sydODOJ8djPQb03a5oMmxEvS9kYTu xH3Vy4rGHAobNuDF69Y2+5K3Yx4nf9NznMf+OEanJISbVwve7nv/RN3QGyRX 8tvNMccM4tiHv+H3u9cuqocrFsK0OQv/zCYyUruuZcGRkg8u591EX4qSUeYm 8SDlHuz/bSITbPB+mxQJCk+KPVirfe2M3uwFmdYzExYteUB59FpFeRcZOQVP i5zEc5Rt68G4MPx+gvKoo8OpAMwyXmUWCJLRyB1HhjieQ86eN3l8C3//x81/ a4f5zcGpVUq+faM3tL7U+6Gfg/vLUW9FY/x5ljtfd86t5YPiafd7/81vx9ua cw4PYz8z/d1tI173dCWlKcjkw1n15F95eN5g1oo2OP4az+lLacuReP1x7snw 7K4CoLAfMA7D85Gqvn+WOM6PrHflwhLuV+buyjUVORHA9P9f0hsucMWtVP/f /earqVnbR8b/7/5//4f97/5/v///7/7/AW6HxAM= "], {{ {GrayLevel[0.85], EdgeForm[None], GraphicsGroupBox[{ PolygonBox[{{51, 435, 433, 45, 46, 47, 48, 49, 50, 345, 256, 186, 133, 95, 73, 361, 272, 202, 149, 111, 376, 287, 217, 164, 393, 304, 234, 412, 323, 430}}, VertexColors->None], PolygonBox[{{432, 431, 158, 2, 325, 236, 166, 113, 75, 53, 347, 258, 188, 135, 97, 3, 326, 237, 167, 114, 76, 54, 348, 259, 189, 136, 98, 4, 327, 238, 168, 115, 77, 363, 274, 204, 151, 55, 349, 260, 190, 137, 382, 293, 223, 99, 371, 282, 212, 402, 313, 159, 389, 300, 230, 410, 321, 5, 328, 239, 169, 395, 306, 116, 378, 289, 419, 219, 405, 316, 427, 78, 364, 275, 416, 205, 400, 311, 425, 152, 387, 298, 421, 228, 408, 319, 429, 56, 350, 261, 415, 191, 397, 308, 138, 383, 294, 224, 407, 318, 100, 372, 283, 213, 160, 390, 301, 231, 6, 329, 240, 170, 117, 379, 290, 220, 79, 365, 276, 206, 153, 57, 351, 262, 192, 139, 101, 7, 330, 241, 171, 118, 80, 58, 352, 263, 193, 140, 102, 8, 331, 242, 172, 119, 81}}, VertexColors->None]}]}, {RGBColor[1, 0.9, 0.8], EdgeForm[None], GraphicsGroupBox[{ PolygonBox[{{1, 434, 431, 211, 281, 370, 96, 134, 187, 257, 346, 52, 227, 297, 386, 150, 310, 399, 203, 273, 362, 74, 315, 404, 218, 418, 288, 377, 112, 423, 305, 394, 165, 413, 235, 324}}, VertexColors->None], PolygonBox[{{432, 433, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 72, 94, 132, 185, 255, 344, 31, 71, 93, 131, 184, 254, 343, 30, 70, 92, 130, 183, 253, 342, 29, 110, 148, 201, 271, 360, 69, 91, 129, 182, 252, 341, 28, 109, 147, 200, 270, 359, 68, 157, 210, 280, 369, 90, 128, 181, 251, 340, 27, 233, 303, 392, 163, 216, 286, 375, 108, 226, 296, 385, 146, 309, 398, 199, 269, 358, 67, 320, 409, 229, 422, 299, 388, 156, 426, 312, 401, 209, 417, 279, 368, 89, 428, 317, 406, 222, 420, 292, 381, 127, 424, 307, 396, 180, 414, 250, 339, 26, 322, 411, 232, 302, 391, 162, 314, 403, 215, 285, 374, 107, 225, 295, 384, 145, 198, 268, 357, 66, 155, 208, 278, 367, 88, 221, 291, 380, 126, 179, 249, 338, 25, 161, 214, 284, 373, 106, 144, 197, 267, 356, 65, 154, 207, 277, 366, 87, 125, 178, 248, 337, 24, 105, 143, 196, 266, 355, 64, 86, 124, 177, 247, 336, 23, 104, 142, 195, 265, 354, 63, 85, 123, 176, 246, 335, 22, 103, 141, 194, 264, 353, 62, 84, 122, 175, 245, 334, 21, 61, 83, 121, 174, 244, 333, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 60, 82, 120, 173, 243, 332, 9, 59}}, VertexColors->None]}]}, {}, {}}, {{}, {}, {RGBColor[1, 0, 0], Thickness[Large], LineBox[CompressedData[" 1:eJwl1WWUlWUYBdARkO5GcuguCUFaWkXFpBWxCAFFxURBJAwUAVEBQZAGxe7u 7lawu7vZZ/ljzznPexcz937vM0zh+KnDT96joKBghy/J/r58JVfKG+RC2YEN TDOfKRtwPW+bl8h9GapPlJV4Rl8ku7BMny1b8JHeiZH5WdyjX8TRnGKuxSv6 jXKVLMIA/Wu5nunsT2dGOS/Cvfo8jqEoA519Izfke3IAXRjtvCj36fMZTzEG OftWbuRUDmSSs8o8q19FV8aYi3G/voDz9Cbs1I/Nv9X34lV9sezGdfpc2ZpP 9WtkX/bUB8vv5CYu1jsyQz9LFrKWd8xLZQ+2MMw8WVZhNc+ZL5dd2cRy8xzZ MnfHx+blsjfb2Mc8Nj8/740HzJfKTizUZ8mm7Mrnlr2YoM+QtXlNX50709vw mV6cIfr3cjOncba5Ie/qBzFFr8rz+tV0Y5y5OA/m83McJRjq7Ae5hdM5mO7Z EecleEi/hOMpmX1w9qPcyhkcwk3ZSVkq96//JLcxk+H0oDRlKEs5ylOBilSi MlWomp3wPX6W27P/HEpPqjHMa79kZ3N3HEav7KPzkjycZ8wJVM/zcPZr3mOe E4fTO/vovBSP6JdxIjXy+Z39JndwDkfQJzvnvDSP6os4iZr5/M5+lzdzLkeS X/BqvJAdpC8TzGV4LDvExNyduQ6v62uoledl/kPekp3PezA34j39KKbq1XlR v5Z+uUtzWR7Xr+B8vRkf6JNyr3pd3tCvlN3z/4k+T7bl83wv2S+/U3ne8k+5 St6a+5d7s5FZ+ZyyMet437xM9mQrI8zTZA3W8FLeU34X2MwK84WyFev5JDsq +7Bf9kyW44nci+zMYv0C2ZwP9cnZPb0eb+prma+34wu9dvZC/0velufBSKY7 q8nL+kr6Z0/M5Xkyz4Yp1MmeOPtb3p6fzSgGZE+cV+ApfUnumrrZE2f/yDuY zWgGUi874bV/5Z3MYQyDqJ879dp/8q48F8YymAYU0pBGNKYJTWlGc1rQkla0 Lvj/b0ob2ZZ2tKcDHRnhtfzBuVvMZRxDssfOKvK0vjT7xUxzfd7S17FAb8+X +orcE9uzE+bdqlHAIQ== "]]}}}], AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, Frame->True, FrameLabel->{ FormBox[ StyleBox[ "\"[\\!\\(\\*SuperscriptBox[\\(Ca\\), \ \\(\\(2\\)\\(+\\)\\)]\\)\\!\\(\\*SubscriptBox[\\(]\\), \\(i\\)]\\)\"", Large, StripOnInput -> False], TraditionalForm], FormBox[ StyleBox[ "\"f([\\!\\(\\*SuperscriptBox[\\(Ca\\), \ \\(\\(2\\)\\(+\\)\\)]\\)\\!\\(\\*SubscriptBox[\\(]\\), \\(i\\)]\\))\"", Large, StripOnInput -> False], TraditionalForm]}, ImageSize->{600, 400}, Method->{"AxesInFront" -> True}, PlotRange->{{0, 0.6}, {-0.18284490314672575`, 0.3123497947907735}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.406293395814513*^9, 3.406293404601663*^9}, 3.406293521781127*^9, {3.406293988455219*^9, 3.406293995309438*^9}, 3.406294498866057*^9, 3.406303393339962*^9, {3.406303445337124*^9, 3.406303494208503*^9}, {3.406303555470802*^9, 3.406303595425263*^9}}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztvQWYFMcTNr57t6fI4U7w4G7BCS5BggYLEIIGAgGCBPfgwQJBAkGCWwgQ IMElSIK7u1tw2e6v367pmdm9WeX3/z95vi/zPBx3sztVb1e1VFdV19Ro2b1d 604tu7dv1TJ95a4tu7Rr36pb+kqdu4pboXabLeQTm83+XXobfuc2m/phXI3x w3zn33jjS/m7gw0bNoyHh4fzr776it27d++m6SvZ8SOEJ0iQwObgiRMn5unS pePp06dn27dvLyy/F86WL1+OW/KjhAkT4ss8Y8aMDvlxGO/WrRu32+3sxIkT bODAgezt27e2G+4sbKH4DlCws2fP2opryHAdOHBAfpYkSRKNJFuxYoUkOWvW LHbo0CH2559/eiQpHgMVW0nXe9pfxEY0moeEhPDPPvsM340lgABIiYvnz5+f v//++/6QipT3UrK///6b9+7dm5cuXZqnSpWKJ02aVJKOHz8+z507N2/RogWb N28e++effyQLcTGn08kmTJgAWfPQ0FCeI0cO3rp1a3bw4EHeo0cPqdLIyEjQ ZOPGjbOl9gFbqF9K9caNG75gh2uqv3PnDv/444/lZ6lTp2ZTpkxhN2/eDJEf R7KHDx+ymTNn8hQpUsiviJ7BHj16JNFPnjw5RI4aB1uwYIFsjGAMsHQ7lGfJ koW9fPlS0vKBna1atUpiX7p0qS/sDg07rvv370ux4fOiRYuiI4Rr2C9duoQm SaqiS0pyApH8asmSJUNMVNBDtTmAN2/eHLfS+MALxuhxXbp0CQAvb9KkifxM DDN2/fp13FN4zVLYu3evJIdega+3atXKjPfVq1c63nbt2vmDVwxanilTJl6o UCH/8DrYkSNH5H0AGjp0aKQJp1A2z549O1u3bp0kgx6Pr1arVk0+rXqABs5H DwjDHMDr1KnDM2fOzCMiItjjx499zgcaTN6rVy8d5rVr16IJJh8wYAAGnOyq 06ZNk18RCmNz584F/JBAYIbRPdn058+fs9evX/sNr1ixYvI+xCWu+JoaMcY7 d+4sR9YHH3zAFi9ezN68eRMMIjGG2bNnz/xGhFUA94XAcT+Bpm7RKLZ+/Xpe t25dnjdv3iCQyElPdAG/kajujUlPXIlMXXz//v3+dHEz91KlSvHatWv7zT0s LEze//zzz3E/nvZV0XF4nDhxeIcOHWzmQYcR5AeiB35y16YmzL4GdwckL2/H ixePPXjwIEQ+E8YuX76sM//mm2/k3XfhLrQk72fNmtXgHoft2LFD7x6iH8jV CjO4mH/kPSEyrHTvyp5NmjRJH7NXrlyhQRHJbt++zQYNGoSZnCdLlgxCkKQF WrZt27ZQTdlYZpTJ5tC6KaiJtRZjyGXeu3Dhgly1xBX9DoCFKqRpBMCDBw9O pX21YMGCPCYmhmfLlk38RT+F4YCOzY4ePfq/R0JdUc1ngjXEZ6wixAka0x6L r91evXq1XI5fvHiB2yHGsMGKCzvs3bDSSHI4HJjYpC0AzObhgyWjY8eO8nvC voKVYcYshCqFa8KMMd2mTRv+5ZdfGpgdipgkKK5QgBVNEfas7AD0F/1ER9Eu m3ocj27evFm156EHsav+RdNgDNu4cSM3XZKMMJZkpxBdFfYZBqlqupxNMYcJ i1TQEBS1SUrBiWsSDZbzWrVq8T59+qjnpbZDpdkeKm1Izfqy0b0w+X2FRSxp kkOIaSxYqaoQccbqymbMmMHr168ve2vy5Mn14ZY2bVperlw5OckIq9y8QKBb NWrUiDdo0ABGoFl7pv7yxsSwtqmRYmDyxo0by1547NgxvZHeKdTRviJWXZ4z Z07ZrYUxF8CTuIR9KoUikAf4pOil8kmxsjncbjds2ND9yZry90g2atQonitX Ll1BWodQk5S86Yj1XCg6ixwEwvKy0cwfKmd+PFekSBGbhwcjsTWT81DlypXV EOd9+/aVxNB4T8/hOnnyJAQjV54CBQpgt8LmzJmj91bzc/L6F2x8A7tRWf4X zU6dOiWXFvyrVKkSNgtYcrAI4htSGuI7IfiC0Ib2FflwRZPAvJDRJhmMZCn7 4cOHsx9//LGKfwBgB6HfhPCyZcvSeqJ2WxqA+L4AyHZA49jPC1sfRhG2N1UC wq/3HxjMP/30Ez0dzyfjJUuWyFEKS6xEiRKyH1cOhLO0LiRusUOTdrHYGWhL n0YpPi9fvjwMCAmyffv20ushNqqw/ZQarb5TyQDijYRaGo8fPy433/irKunP J+czZ87gviJx9+5dLM/u+vNFhn/44YeYKuVaIPYFWC6rBAKe5ijMCaIB+FXT n0/G+FwsWvLx6dOny3FQORDOYdAd++WXXyQJeA6EDpX+SAHh7MmTJ9JAEps1 +TVsOkRXD/X+sY+nHZrMtVbLv6r4w5CEBchYwnVheXwqXGN0+PBhNZdrMvL4 iDKDYLGIYYHRjelc+g/E7z6eNm/1sE+OioryxVCJEk4KTEMwRtTjYpF30wec I/iW2ODIb6FlcePGDfX+sY+nvenDwyNhmj4WLVrEu3btCosBn/pgpPTx/fff 848++shFHx4eMZulY8eOlY4Y02jx8bTyWYjBDh+T1Cr2jKLn+2KsRHr69Gn5 8W+//aY/Lsaam17Ut8xEhPbcibh97ONpb3rx8IjSC0Y93JLin1kvXhk5IFxe pkwZWFCwZXzwMqvm6dOnbMyYMVaq8SqaMDmasY+EIVuhQgXTOuRTqhcvXpQf oxNWqVIFczEed1ON2BfLtdE87nLkyBHq/WMfT3tTjYdHTKqRTkzNmvXBSA2Z NWvWqH2ELh6v2Bxs4cKFcvfkpg+v8giV+tBcmtpa6lOEMJfhAzl37pz8WGw+ MJ25mAJhUjtqQWrXrh3v3r27TVHAPYvPK3n/WCkh9vLvndvs2bOxWZVY165d K53pJk144KY67K5duxB9kf6Wbdu2wZbXVe8VqEPuqxB4cFnqvWEldcAhYFaH Vy5h8k9sT3FNnDgRXk1XsyxMdtvChQtjxZcOmwcPHph1YfV5Je8fa7qQcY3o 6GizLrxyw7IHkwauV9jQYiy76cLqcaULbDWbNm0qwwrCGGdTp04168IzUBqG 8EO46sIL1jD266+/yqgEbALxT3Kr7A+3MGzj5Y4aJkXGjBnZypUrlU4WaJrG pXYQ+MFJ6f+TXVV5Tv9rjpJk+GE3fWQLNYdUbREW9+L/d8/vezH/Iiz/9nuy r9n/65D/dch/yb2gO2Sdd7hX9x2e/b8TMNcvDTpjGzYwtngxE5YXpXL4T4Js 7zC5xgorjHOHAw8hQvxOpITxKJdbkBKGgjUpF0qKgAOBES42ok5mtzudvXpF 6F8NAIagsmKFYGm/9iZ5csaePYsKrDFmm0/ssjjPnZsJal/c6tcPd96FGjzZ dvuF12nShJ8+ejRoamHwxnOeNStwfXZz2DCb2BTaooMkNXMmQJ1+lT694/Tx 48GQoi9FMG2jw3mGDEDW6MbYsbZTp069M8nJk4HwyMv33w85deLEu5F0oEdw niYNENa5PmlSKGgFS2r0aCDb/yJXLvupkycDJEU5CVrXELtJLrqroFb16vTp aGNQ1MiGZ0OGANfO5wUKhAVJ5cEDzhMlAp6yV+bO9ULFIKI9y3m/fuC+6Vnx 4miHfDZS/6ovBG/NUrl9m/P48YGi2OVFi3Rq/rfnrSEVzr/+Grh+efLhhxFB Url+nYudnFNs/PNfWrkyOCqcd+oEHEv/qVwZLQqQymuzfC5eFEtBxFseGprj 4tq1QVPT5NO6NXDNeFS/flSQVE6fFrNw2J23iRJlvrBxY3BUOG/UCDimPWzY EC0KkMpLs3wOHkQO2nMWGZn2/NatQVPT5FOzJnCNfdCyZZwgqezaBQqPnXHj Jju3e7cXKgaRl0oq5crh2cH3OnRAO+IENqqem6WyYQMo3XubIEHMmf37dWr+ t+e5ag+Wa7GxF9R63unWLV6QVGi5vv4mWbLoM4cOBUjlqblltFyDWsfb/fqh Ze9Cbe5cUDr/Om1aLNdBU3MgQ4LzzJlBrcXNESNigqQybRooHH+VOXPo6RMn AqTyj7lltASCWr3rEyagZe9CbcwYUPrrRY4cWAWDpiZa+Pgx50mTglqVqzNm JAySysCBoLD9eaFCwBIglUfmliF6HhMDaqUuz5/vg5pB7JFCglmjRw88v+5p 6dKJAxuzOpVr19Sqk+vimjWJg2mPmD86dFCzO9oRIJUHZqkgGTk8/DVzOLJd WLcuaGqafJo0UbiSBknFdZYPjgrniL3R7I4WBUjlnlk+27ebZ/mgqWnyKVsW 1Abe69gxeZBUkORst2NVjn/mwAFrKvnlf8EmiDg8YbhjlgtWiwIFgKXrnV69 IBefWALNdAjzC8jChQBx+XXKlJGnDx/2DCQH/gtDLorYxnbp0gV/BhLA94HG Qfu499+Hgdvq5tChKb2II6BItWQcpvO9pdjRCqL2epJduP5VN3YBBGA9tvOW WerPn6t9V8Pr48ZB6t7bG0CY0T8Asdcp7wCY/yE2HwAcCGNgaWPa0pbaF2N/ 40m+GHM+YABaveVZkSJosTXjElqTcKl4i0fKN8wyvXXLfYdmzcHjLHVDCQhz HYaY3b76SblyaYOkgh1RZCR2RHkvrloVHBXOW7VCi6Y/atAALQqQyjWzfI4d g2PtFQsLy3D+99+DpqbJp1Yt4Jr4oGnTdEFScV2bgqPCeYkSwDHgXqdOaJEX KgaRK2apkMV/802SJPHO/PWXTsN/2+iKag8m0Jw5mWbxZwisPTqVGTOA5+zr dOnCTh87FiCVS+aWGc6lT27QJPcu1IYPN/uXgqbmIEuW/DkVrs6enSkoKsqD suFpiRLAEiCVC+aWXb7MeVSU8qQETU0bFy1bAtf8xzVqZAmSypEjGKcvWXh4 xvObNgVHhfMaNYDjuweffooWBUjlnFk+27aB0iNnvHjKhvRCzSB2TrUHtlbx 4qDR/26nTlkDG106leXL1biIc+bgwayBteeMuT0Yp9mzo/+1uTV4MNrzLtS+ /x64Tr3KkAHjNWhqDvK7pkwJXHWvT5yYPSgqnNOOc/fzfPkwTgOkcsrcsps3 1Xpa8vKCBWhZUNS0cUF7vhVPKlbMGSSVkyfh0cNeL/uFtWuDo8J53bpmj16A VE6Y5bNzJxdzxhNndHTKczt2BE1Nk0/JkkzbW+UOkgqNkFtvEifGSuaFikHk mLk9xrjocKt/f7Qnd2Bj1YVa7HERYLuOqXa5jou8QVFxHxcBUjliblnscREU tVjjIn+QVFzHRXBUMC7QonEPWrRAiwKkcsgsn61b0aJ/nHHipDi3c2fQ1LSV o0gRpq0cBf2mklf+jHbK684dp3PvXiR1Op0LFzqdc+fiSKzTuW0bujwpQVxT FVNytV59kyIFXMD+M6UjbxHkXXj92ulcupT2edBwqVKMxcQgAi0eEfvc1Kk5 r1iR87Ztnc6hQ53ONWto/InrO7MsYcOlS8c0by1k6T+grKpFOMW0bBnnjRsz JgtHiG9FMyY21JqeYNU9fBsv3lNnVNTBl9myPXNGRjIZObfb6eviZ7JkgMvY lCmM0UkMPtYMlQJ9B17kzKlCo16gSqRZFEDsBJo352JPDxIhJKRHAtChF1mz Nr0xalTBS8uXJzyzd2/yszt3hsuIa0TIyePHE5zZtw8fNbkxevTOZ/nzAzfm ZJxQxDFB2JmcN2sG5ZNcRiuWGMQxMSqoWNjviS6DIoDr8WPGxo3DrEmYHaC3 53nevDC4YSwlklgdSc7u2YNbS/6pUgW2ZoiSaYMGIODSD0eoAdm2LRqx6kn5 8pBlYb/Vns5A6HQePsx5wYKcy0x98EQmQoPr48dHnD5yBHSTEUL8mujs3r1d bvfufeJVpkw6QjHRsXr1nM7793WEQ8xaJ2MVm8qsF9av94xUS3FSR55prGj3 cF27ds3UaV1blNYs8zt3OO/Vi7GICMIYiuAudCgbktpoDqb5Sldnzdr6rHBh vTkpUmB+QpKJ3pxBBnXOoR+7fcKDZs2KemuGPJ5r0QyZ5+6pGXTkn/yaclHr 21ctJPffxMQcEgPvg0uLFmWUTQhvc3PQoL9fZM/OtQwWeIUZ04pZ6JjXrwfe B2/jx0989s8/rTHHkf+F4WysrBCzfv165C3DoSXQiN8E9tWrV8s8/RQpUtiy WWBPaVYBvHBI1UieXFPB8ieVKqGPS+Ry82SLVCpAYsPF15jv7HaMTcaWLCEH nLj6mbsSEn8oYtbtTs+eeN66PVEaFpzp+vbbb8UYTwTvYHxA8YxdMNm9m5y+ JNA1T8qWzXdx5cpsBt6oM4cPw5OAGAfhzZULPmu9r/Qx4yUP5rnX772nBpM3 vKHyjLN2Hlr0Iw9gkysO589jTcauFB3k8MusWZGNkcsAi4AH5hoCKjbUbNgw l47dW+kMXrQUKVSIrZg3ocojqyiu8+DBA722h4f+/vTpU3nm2VN/T6r6O9Yf rA5x4qApMJhrXZsyRQUy81B/RxwRkVv4zKhJpUvTNG00hYaM3b7jecGC6Fve mhKKvizPg9NfwY/bxOZxK8YNz5IFIO6/jYlBCFwZ/wWoHYUuLVuGtCpqQ+bM yNTQVfK1uf9cukT9LCSk6OUlS0DCu2pQEqRixYqohoDHlWpQ5MVqRkX2vIsa mF+nRi2+oqTnIpWEJlbsr784r1rVJruqLXTfi9y5E5ylcL6tEIlFhubEFIVR R6JJnJiSlnSxaJpp1gyf//y4evUSXsUBFziGE3qsuDIbzfTncKXFVyybGWOQ dTr37MFaz7UxiXXyo2vTpiF+8QG1Ep0SpiplMKIGw4wZZMaKq5tqJSbxkBAY AunOb97stZUy8oACAhuQXalaGezBT6sGxjP3boRTO3ViYjUHhhdikR1yr317 6K0kNbDGtalTYTGSCtu3Z4KV3ru/UqO0dGkVJ4TarRvoI5QU6sfHciaVP8wn JkhjVJphzpw58gBKnjx5ZBkH8aQsyoJjIpUrV5bncHbv3o2ZLpnG0MehNq8f u0qWzobFYezGjQps0ZSGbOrwMc42TS7zVMkf8vhxmRQh/uGIC+Zuxmzl1ASB 9CWKxh3l2TKlY3vX2Nj1/WKDdONAJ/ywOSLZhV0XedpU1M9mzarsLmGtDo48 hTNQbH+12heRUmIOOXYKFizI/vjjD7BOonEO7jye5/aL/j5/vlpzRQttms3A GOwZbCbu37dVUM1WFh5l8hzkOd/PyPasMje7D+vZnkSXPz9MkVjNJnNHlmyR RxBRjAH1iT4S/CZOnIjCPBI9qk/BhMCJXqx2J0+elMf4cIIU0S7Ip2/fvqqM iKmikrfTXV4Oi1kIR8zcv/9OO0FsBWDXdu+O7ZXsDPpRKTF18ZIlIZErovO8 z3Ys1ztDeymVcPydiv297h8eN5o6xObNoOFBOg5ZzwUlwjp27AhpxNVu4+Qj jnp26dIFDTeXkfJyfMrzqSzXRtc3qxkWcJkyAHudJ0+anW1damqQAw2azRvU IFXXrVvVvSFaMR15urFmzZpCjr/bTOVuZI2p/v374zBzZgsoDRUUYcTwChVk DIgnS1zNOe87m/PaPgmlLUEpxNbNfctDQ2ignjvnEcqyZcvkbALdJdDuoR+g fhrKsGXyhgPXs2eMVaoEks95ZEQO55bFEkdrQ8d2dm3fdl40P4mlZ088FQuP 1a6wEfFxOmF5Yytvs93iSRPFZ6e36n2plcGnAZs2gnhgx/ToUeB8cJ07h90n 6PzKy5eMYae2mPi4TWDTp1f3h34T85jIL+WwlRcrANp6O1pazxSx6Eu9JeVp 0qSx6iE6J1zI/s2SBdQW8lqVwtml3ZJTc++jz5qjW1mkRPhpBaCZAnDmDFQG 5pN4ywboAnpTm1mPlhqeGItRqupm6l0UUXbU1cRfsotKGM3MzYfhQ9b7N6xn B3fu7gOkuk7EjXtb+DS01pkHSHM4fxR3NynQJw4UdIRVdY8njKnHpn8r+TeJ PTC4aWB4lAImZiscLVq08IHDJX3jKY+OKsp+na3Lo5H3AeQJDyx6KzzSuPCE pyXJhXO0xW6/ylMmS8MOrDXhiDXAarnztxpgLVU7cXqifHm1DEezczv0dja0 HmCx6JMFEi3bgmUWp7nFgihPyqMyC0wwlAycPn06tiG5vMHBhZMAtPkaz1s3 gsolnPreR6EHWKHS4YFCvfgozGQVoVhs586dYUXnsIDUSvVIjEuxixHP9GK9 O5qg6OORemPz5h+7Q4ijQUDJChQFxkdRJggFChTA2WvAt/IStDJLBTl1ISEY GCXZyhm6kuoSlA/YmtlObbsBg8gDlLjmQoJymwmlYKuJ6QLWrpCGrDaGcmfi c1UhrV69emoDYLl7bq36ae/eqp8mZUc2SogfG4oLYVf37eP5cpDEhNklrlhQ tTwjWUIQNbpQGWAzKVmrHuxgW7dulZXXYKT//fffVtJrraQHjxPtVdbzD4u5 TKy1DGifs9F9qZunSQOzxQs0qp2CirSo3BytOgrDrHhOOqp69uzpFRKuq1fV fO8y2xIkB5bsGzxZEuri8+fX9QBFlsNETRZIJ64G77vvvpMlpTxpq62CDMdi aOgb7ggtzZZPN/Mfzds1ZbqDxun0yD9TpkyyDKWZP4okp0yZ0gd/srLF7gxt vMZTCGEc3qjLobqhmmxs29JXHLtV+HXItewRT4YMGXQ+Cg/4YArQ8Ug4JhTU ScqWRYs38DIfoJ+aUMiR/hsvW4z6bYcOdXUifnLX7ltKo73SBpzjdCSrL+vR zsy/Nps1mrSBKpV379b31Hqx7VVlgHT+KEOIqrne+ZuOajgc6BFl2bJpujYq G9rAenOep0tNPXPyZDzlEQ9qdGzZssUFDyo4oGC4JzxfkDzUsS7sG5KzQxvc cfRjX7cjmcA1RLEVjziwBqHSlpjeXAYsSkLAKmrVqpXVgP1CiQbdo1w5sNvE SxUJZVf26qKpSCrCsgxTgcTyxx8NraGEyTo7mH1z5swpC208fPhQn9bEvkJ6 MTCtLVy40CskXDgeJiwPhE5Z97bukObzOlVJQth2MuYREryBFcT2SPvHRo8e rRc8h4RQtRWba1StRo25p0+fWkHrpAS7aZPqROXYkqkSUnlDcaXYihm0Vglj k126ZAktqUYLnRll4TRPjNVK3cmQCJJP0eL9PG/2xOzY77pEPiSJQHF/81xZ SUmDBjVyZ2tlK5noU0cQm3/xPFZ+d/rt2YheJPGMGRH+iEU/xtT3fv75Z7mI TJ06VVVxQpnMqlWrSmGLDQsmchgPVsLurChduKDiS43Z5CESTRlD2JDCXZ4o AbV4xQow9ogKpTOR54vuWbZsWQ2VrBJm2svI3jlz5kzPqDRZ0UHf+zxB/NTs r3W6rEoZc+oU3rweyatChSbWqOKKaaCfZF+sWDHlxhKL51V5D3YJBhSqweBv MXisYOk+y5o1mckOMMPJxzYuQI9FlgUTJpEHOKGyEF2HDh3kWx2o5mgEygHL 8s+qmC0Mvho1akCPlvGfr5TyUG+XkrTRdSSeEq7bnR28SD5adkTXFpcHXJEy EAK/nXbJ4YCiiJrppA4vzJ8/X76vgrushV+Z9bZ0KRg+4XGiM7NdK3W9FTOA tWTj+pPeUqVS+55GOjE3YCi1w4xKvtJng8lPW6cUMBTh1OoqWi4L3ZQOKXSw kxfOi+FsguZIxI7/cZsnSaSMpmbWsgqXdaBNtY5lCak2bdqofPdmzZrBFe8L CpyaXJihLnsDDcpU3qwuSUgsG+LyCAWOS9GlZGwWJb/g5RJGegjbt2+fYLFz 505Yk3Lx8jATdDPrDqhDQl6K3U0e9vtCXXeFDd0VYWvn0D4+PBz+Pi/oQtmQ IUP0aGaEthfQCgN7lM7XSlEdO0IRh3iOLLorpbAxD2OOpo7du3dza/4O1OaV tdRRaVnwDTdBcDgcviAgbgEHO9gM551amCF0YMN7koKEiSJWtFgQkisajF5e IToFfKlWGnAJ+CH+TItfcbZ6lq6B/IYGYM2IrU1885wci31KTQPLli2T5y88 NbWnkjbZ0Wd4hrRx2NntJp5ypv2J16tOzf3oo5buvNRWmRY4KtCs4tnYGKL5 PviTJyFfPhdRa/yxtzA2p+fPx+IfahL1qlWrpAOfuZa7RhgaL64QU1chjyg0 DcCfEhkJjjA/dA3kNjQQxi7vQcRF2QF4yiMqzATocGJjhyGoUGGtRvlnmPbi YytUemx+6FA1EnTdmNB8zfp+SbrBK3SoGrUnNNiUy1DhEwrJmes/Qk14pwtE iGq0vlAhap5f+tLG8LZNTaikMWs4WNavb2WNJhTGiprb9WpyiMZXqVJFBTN9 Q9HUtn8/1l4MHExRutpyEqQl/KMKZmPWA6RI9ZYOxUOSBqRPPvkEpwfVcqOh kn8VxE+JrI8aTaKriY+P8GyZ9bkru6GxKmz+RAITLx6MEV1joRaAsNapNxmh ROcmYR6j0h4qTcs3mHgA5CarPkpt8O3TnDaAd2tjgiZdgad5xvdIbaNGtfYg I8jB3aJDUXtZntW0JAvL1A9IpgMEtPLkYpsX6erLStCwPzHv1jxAC5dWAeKY CFshiIuAnqwxzvW5QHtrEkKdVp2qn9IhJd3u4QVy64ZCFkOHmdjulQjAULrS zp2gHAtVHBNbvKhhz5498HF6YUtnyt+TOujKBnQxsZUTMdwLKq7e1p2dmnRx iW2z/KuoR16a4FH/iMw1NEkXfCbi2ZRNHGy4Du7c8cZTX1OteA5QYiX763de srBuQ2cwxOpqg/34I+i+K0+KFidPDprw07nznMEb1zaMLbLB35Gny8Hdhzx+ PPicddmmI9mWYct/MBa2U6faB8tzkJItisva7St5lbI6r7RGOyPYxd0neeb0 JNvBg0E3Fk9V/hpOOBfTLaP84H954lpPUITnKoHcYyL/wR32QN61NalHzITC TvAKG1sozJTin9he4WeoTHJxuWdqyv/iwPYgs8p/+AHtwDSKrqyrITWpHBHs FzwinGaNHTu+8EflWTT6gZ7klqekXfojkoyoLstk3qK+Di6l654RsVKaY6jA bludmB8gAz317SbLIapPiB4mhgXk9CFbOtUdZis2pq8ROaMYSECyfMfT4i7p zYMHQ1qneKZ0WEF1qSYnlcMtjbgLjbqff+7kBabMLTHDfMcz5S4wYa9R6uxg 3uVzd5iGt69aNXw7Fsz/0SnzYUrBOBkVGqoWfAklqaHgymzBJIITNy6cVL4g vcv582GmlmmrPtwD8OLrUkpMUkKE8xxPl0aZSJ3dIakQDzxMir45WqvdV692 k/e8nWA3YaN0BkqShrvHHdu3/IvmJLJ8+WAjecLmMbKt3beMJI9QaqMDZGpJ k/wTGmqD18CItezbB7oB49AyALzjgK7g7rHbl/Hq5Uw4BOMrfxregR49vvLE v0yZMnq0Q/HH/hXGoXf+mj5wwoV8p7XYj2N0fcQQDsTlDG/uP/94xKHtDuGy FTjwM1Tlq3vEMVLJYcIEtBMpX/HY6W2SfzxDHznZ5sXGerNpE1B7xWHzlXCi M0btqmIypDWBt/rExFi6Z3bxQnlIAW3bdnNnaOUmH2kW7KFDcC9hL4fItC7Y OEQfdjEJNm1aDApr+pK8dpRHBZaxhcbmXVKLMsSUhe1c8YxHRdKwXrtWF5Ov kz76USGcwcqc2cWVR/QdrkknLVt+7Y80RpulgXhISAjwuVjoEUT/CzasB0kj ZUoshbHokycqSr7jg5l2tmLakesNUgtgyDVr1syyq41VMmwupxcX6zLMkGEG tme1selfvhy8PGKBy3LlypXypUfFixeXy7P2wkOJRawg3rCQRU+eb91kDDPk jSgbybtRo57uGChVX77jUb0TVnrcEyZMqL06jl4DDFfq7NmzpXF77NgxvIvN MyRNVXPmgC02MC7xo1CC1pqN6kOqQqXXO3diQaNUyiiZkhoRESHfdIcLxjje kVm/fn3EfWSioLJO8Q5JOC3u3r3rRWQmA7CCdIT8yBsa8Se7oUb4bjC7kxoX LsRTHnBGyIgPZqpqZDPIt0bAR8H0FyhKfPoJCyt845VKkSlDJRZ0i89ubHzX 8nIlSKV16vT2ZwiNN7ebPGhmE40R6RZs/ABSCY4a3LgRi7Ra0FesWCG38Mw1 bxMvUvSYt+kCgcwvZjK/3hpST8EO/WYYiTh9yphHKJAmOufOnTtdslmx9xk2 bBgbP3681cw9QX3VsLxysC1L3hoyXs6rfcg1G7CPO3fly4AFg9R4YfOhX6qk YIyQ6tWry1e1CiCeAeh2FkSB8JRKqXlteBwIBF6ucuuWByBRssnnz5+Xb/7D eMV7EmfMmCFPD+JF2gKFFmGLkh0Tb+5ENGvkyJF6wA9zorYFsFztTJCV+QXY bdnI3oD80nXcIFZpHjexoKtse7ij8RIRjGk7DrPiO/SGXNzC0DfdVq8zXoMD wIzJw1QS5iSl0SVLwBZBgZTs4PqXxkRouItq1+6jP+caIZUBcqQk4U3E0Cxs ZbwrC0LFa3PxbjhsizcJ+0EIyyqEoSOhuqOQ0SJes6KGpBMb0p1QYD67ebOf P+N3kln0KL8SLx7IVmdzx0PeTw3RI9fdyJ2ggIhfLKYo1OPH41mkQcZlZ7Y/ NexIeP8IeYsWAwIiiThI4cJM845rJJF0wvWo0cOHAZDUBIFDVOTlLsB+m4fW PzIEgVswnMjao1i3Xyy0c+0omAPEMFfQ+keGX1iFO7joCYMCIEkxNGEAi2c7 s0HdNJJjeZsmRC5HDhhPAZB0MYrggUV+F1p/zxAEko8wx1GQEgf7GfOLxQ9K EJ9/zkzbQCLtgIV9gb+XipAPHDgkAJI0fVCpJQhCIzmLf1KTyAmLSGzaAiCp CQIvFiL/LSwwtP6WIQjsS4h8dDQTe3M84BeLGUoQjRurtCXMJ0TaIROqtOOi OOM6zJKkpDhDNZ98sLJ0IBvT94arywkpFYSzVi0sZEP0531i1MTw889qX6qy Ra4ZLLqzfp2JvLAZtbOW1pjdWMxSYqgtPcSreaXSBmmZ046C1jTe1qwZEQBJ 8npSiatGbMrQKwZadDMVh0GNOaANgLRL+swdnjihSoy8ZLCYy+tVI/IFCqiY YeAsYE1WrowW/MxrVzZYOL5mfTpx3Xtz4sRIf0jPVoLBqWaa7JG4eMFAXdG5 cDKRxSkxYZqLKwDSGuqJE0EDGaMJ2YnNoHuOUCNUd5JlTgcqONU9OmDS2BOX KgXkM3mjWibSo3j7ZjRY3nvP6bx0yS/SPymBoHgLWcZIfjvj6iWjY6Lh8ug6 MARAWkP97bdAfImnSalOyJwi1LBKceCJkA8fPjZg0lgFixYF+Wm8WR2QPuE6 8I20WLEtEfKzZiE5zDOLQ0xn4qmP2E/jTQQRCef6/oqG+Wj9eS+Y55kxo1iH 3Y7sabWHO0biwDEOY7hv2DDeH3G4kEb/oDRhWEcgfcRAn4Qd3QQlUAvo5W9+ sVhgHjjkj8LJEo001q91vFxxIluihAqhBEDaJbSBdEB0DEA+RIKRgWTdvrty ZULApDGTVK1qnmL/du3llFsHpyKlBfjFYqESDJyuNNUi9YxIO9Ky/b8aocZB gyb5T9Kl4DrWnPfYPnn2dp+BGrPJbl4wt9pQoZF+sVhMqFW1pY28dFGMFBPp 7/jnnxDZ9OmRzw0sgZCmMqi0GHdkQ78m0g54UmgvjGoB+/ZN8Z+kJhAUEHJN /dttoIadZjiIKPM7Fgty6ceR1Px8fbfN6v3dIRZQlyrBfvqp2nvCVUkQHUit OcazZiTBNm061Qs0f99Mrb9V3vxqai/QKCWQjETYLLuN2AzXI6BPn3qA5jNC qL11Wf5lRAhdNEjnhOAyVmb1dkODRtIbKv5Q7t8knY4bFq+BQHJYmN48bBUf Xq701a6dMsYRfyFADmT5vJSbHXSnRYt+8CQUD6E+maEegrfreuFNR39z5VIe s62GMHBAyTg4R9F0jxh8vBlWe5cs/eURjalEL52zhKsEUDYTKsw4EJLa8KuZ 3iMq5v0tqRbvWLVZyWqlkhWVpcUYz8q2LyNUjnm8blXDW3/r1gwvaLy9J9Ti LaNe0GhepqZNgQg5+8rvt9HQIDz1Rozx8GE8Egud1ZxnYkH5OvRqKxVTJxby 1A2ldzjkydRZ/pBerWSJjT4KrNntddmMkb8ZqPOzDfMNb8eIEcAQAGkXIxge o/Tsz19Adx2hxu4J0STqRJ99Ntsf0hS0dsA3i+f+5PlzoTP+aqDGwDWcKcIG En0zENI0M9KLumBMa6SxMOKwNakxSxYcwkbzAiBtyrsND4e+lIfpFxKIkX4C f+TRoz/5Q3qtQj1ypNqFQbCrDIHgtIXRO9avBwZr0pLyWoUVllI1uXdDmoby na4wLFSjb0yc+JP+vE+spjPs5ONEOBmklxmYUWZD7Qekr0hc8/wRx3olDpQb EsaFeP4T9v0wjTRmUiQeE1lU/KDjWAGQdqknCgkUZOuld2xR7C5onEZ++HCB Pyw2qN7dRlq5GB2AvMiYdWGjuXfBAEhrB6mpesc2/kF+0ATk+a6bJWMqhZ19 5szCgFjgQqmstFLQ3/NP67mxwLgyevr+/XjALxablHpRtIiqoiHURaQdxjQb IUuZLPafpEtUDYt9CbZqJqDOMVA3YZO0gSQMHzZvHh7wi8UfSq1du4I8Ojac +ETaUYPNGWeYF1OmLPWfpHbgDwVb7Pa9PH9OtfbMNFCbRyqi00AdGAtcd++q 9zSp42QaC7jE1d5Gvh5AS2v0i8UWpU5kwseNCxpd2MCu012DHcYiAZfdmzfL LUlLylvMmHE8iKa98mzx96A2zZgGDDczMsNR6kTD7GsS26aU2b8/nkf+BwzY 713tNYN85cqoFLDSH3FsU+ixe6cMS7QdkRbQnWyFHiUBaRoIjAUuWJw0FSB/ HHQnemoFsVntD4sdigXc0HlkbgYmQ1Wy4TvvrQiMBS5sdqhrwvmACW2c91as 8YfFLsWCyp3geUyKahkc470VgbHAJexyLTN1Iv+soZwkvbdirT8s9qgBhrpO 8CdR1sgo7+gDIG16bQXVmFP5YsO9o1/vD4u9Cv2UKfSsw/ExmzlaIw1z8gTP ksEgS8UQAiDt8jIJCKiac64suDTYNdSlChxQkYHbtzdYspAc9hNmKlxnt78W mGs7Z47WCOI8hvJMYN5WZxLX6897wbxfYYZlhh5vs93nCWKQ6Qug/Q3Mn7Gx /V5rCpHZNGLV3+SPWHQWYAJnPblUUSYjN/tjYV/JIlJZCYbBqqpk0NzvF6u/ FCusYBT7gaccZ7BA/xtqDYKD8HAzTf+IAYHFHwGxwAU2TWWpCKyGWBF6GQJD MBazh7FJW7Nmiz8sDppb0VaeHMFaXo3N+w50exALeDEP8DzZjZScP/8EpsBY kFYQFwUd+CkR3+putAI9AVFDYpMnDyoobvOHxWEzCypND0NH1fvoSizgA1fH AZGRRxWnGfOLxVHFAnVqHbI0NbL5VOrrl8YuCz5eOvcJNh07YnXeERALXMge jYoyl8jqaHReRB/g/jV61fDhqvP6xeq4YoWzLwkTgs4Q3rkVPHvSfUWtQdIL MoGMYy2//QYWuwJigQtpiqmkXlH9CEmE7Qy112c/jMCkYmw3Ll3a4w+Lk4oF ltACBfD8Yl6joor3f24IDKmM8AUbox2TCuVfWbOSnE6a24DwfJ06oICpCXYe MQhHc4ztEop2UVLLHp2OlzacVixgM33zjaqGqpLOWhhiasNGfmMYwmXKoELL fn/EdFaxgNOZoh5zeP2PlFnW1BATDmIYTgJ1Fv7VK7QmMFa4jh9XDkH4uREi Ilbh8PdgSWL6QKRE0r/8YXFesYD19+mnXIYgkeHTsgHiQp8QC/y6gH9cRX1M npSzZw/6w+KiYjFvHowavLT+IYsX90Pnku+xuNYzBIboEQr/GfMupncaiAGw 0ur+XrzIWIECqqTnJNayAeZ18LHVNnjCbkABQENJOIx1/jwoHPaH52VqntOJ jG5kcIudueYlwYivqbESLY3jPLsdkyjcqEYYBRW8nc6j/rC6qiSJCaBwYaUN WKIJ2MnN1YyMOkykY1mbJlgWNAHQkqyFKY75z06rvgvXyNixKkAOTxs8SZhJ q1AfqcrmT8B0ZPSR996DIXvSnRX5XeNLqj5CLbK6Qvr06VFHQB48PXDggJXv 9bpZ78jZqlZNr/fM8+XAtlGa6hUk1Agkd2/hxQsqqIiOI3EB04btlBe8vuIv snAcqlcgJXLx4sUoGGSF96YZL0IfOXMqvJgwYCLLflqWOg9cr0h/NkSLdJI5 c1QGxWlryB7jMsr9j4gECu/Kvyxg3lYwoXxEYVOlUjCxLVG+LltJggmxruEV SpFXhI4Z05z38CG4n/UE00PIRsFEjiCKuegwZdjntuqdKB6CmkvUM5ELiWUr p3PzYrvz6t4ixqBA0K0mmz32D17CeOkCSmb3748Xv+mSDPgtnw7kc+NYB8qg IznZSuf3FV74wvCy4/haziZhxgooNxgGXiwpKGN4madOoeNF5ZIhQ9Te7KIn kfoI/CD/HoUOPWn+oQL74AF8/VT3XgwUARjpwiN4xxZY3fIYBy9g4+KwO+Y/ HWyNGtRNyaFy2RNY5j0etGvXLlmRC2cpUfIkJibGMgLzmJThdGIXhfNU9P4B Dc0PvOnH2C+qXJP3JfZoFTKvxH6eNJ03qQ2XubEWYKuPg/sbN173At5b+AgF rlHlWwgfQWMcTrcE/0TrSfJCTQEUak+cWHvfDCxgWBnYT5dyLp+uHA/pSQHv OfetmcRa1N/P8mYXNnk8XQHK5bRxI/K+brq3wWrWf2Zql0CPECfnH3+s7GcN EA5b/cIqlkIeO9zzytZNRjOs2hMiHofUB9gTh3n2LKixjQ5H6IoXx/4C72kB r9v+oHulxARrE1a9mIRF56QXAUViuYNPfBRr13QVq1S6G+v3ZXHnqpkq4Sqa Jir0V3TfGs7Z48Qq2Rj5O9A90iZgE2qvFWIMZRLEqsR69HA6z5697w9AWU7f FuF0Yk+HOfOzz2BNQ15MDnebPDQsfsa7wlImF+wir/PkSYAaRhxGF84wwQxC NBxTA5JN3Z4Np5eFoCb3yJGUgUUx14f+YJQ/KR0givwke/fiVU8yc5jh5Yxi udVeMC+ZRpF4qT/hbR4wPhBWxQCH6/vAAVj0z/xmTkcKolGFGS++Iq//li1O J1449O23eJuB07lihdO5cqXTefWq0x/CERb36rzDvbrv8GzE/zOA4/93z+97 Mf8iLP/2e9YdUn7nv274Xzf8/+ueTFeL5XXi8jv4YX65C9evxv8f3rDZ/w+m u16Q\ \>"]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Step Initial Condition", "Section", CellChangeTimes->{{3.406361638579706*^9, 3.406361742207674*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"StepInputSize", "=", "0.5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"InitialConditionProfile", "=", RowBox[{"StepInputSize", "*", RowBox[{"(", RowBox[{ RowBox[{"1", "/", "2"}], "+", RowBox[{ RowBox[{"1", "/", "2"}], "*", RowBox[{"Tanh", "[", RowBox[{"50", "*", RowBox[{"Sqrt", "[", RowBox[{"1", "/", "2"}], "]"}], "*", RowBox[{"(", RowBox[{"x", "-", "3"}], ")"}]}], "]"}]}], "-", RowBox[{"(", RowBox[{ RowBox[{"1", "/", "2"}], "+", RowBox[{ RowBox[{"1", "/", "2"}], "*", RowBox[{"Tanh", "[", RowBox[{"50", "*", RowBox[{"Sqrt", "[", RowBox[{"1", "/", "2"}], "]"}], "*", RowBox[{"(", RowBox[{"x", "-", "100"}], ")"}]}], "]"}]}]}], ")"}]}], ")"}]}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.405075078890408*^9, 3.405075099571019*^9}, { 3.405075280201322*^9, 3.405075280380905*^9}, {3.405138720365396*^9, 3.405138769503144*^9}, {3.405138806641839*^9, 3.405138818255853*^9}, 3.406292546382388*^9, {3.406294019810198*^9, 3.406294023785833*^9}, { 3.406294113461681*^9, 3.406294113644303*^9}, 3.406294744834838*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{"InitialConditionProfile", ",", " ", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1000"}], "}"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1000"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", " ", "\[IndentingNewLine]", RowBox[{"Filling", "\[Rule]", " ", "Bottom"}], ",", " ", RowBox[{"AxesLabel", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", "Large", ",", " ", "Brown"}], "]"}], ",", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", " ", "Large", ",", "Brown"}], "]"}]}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.405138822275371*^9, 3.405138826903682*^9}, { 3.405139053394552*^9, 3.405139057448738*^9}, {3.405139100679759*^9, 3.405139118736527*^9}, {3.405139195063642*^9, 3.405139208101695*^9}, { 3.405139252027806*^9, 3.405139279406005*^9}, {3.406294555418806*^9, 3.406294562209322*^9}, 3.406303625962486*^9}], Cell[BoxData[ GraphicsBox[GraphicsComplexBox[CompressedData[" 1:eJxN0wtMU2cUB/BSawOySREZCNGB1QhaGIM1Y9U+GAgBWxaY8rIWKdqiYts7 GJanTJ5FRMEa4iPzwkZJmKJRthWEbsFuOtxsxQywpQgy0YgMKYhAYQyjnu/e pGl++c73P+eetN5ieex+MolEOrz0ef09ap3nh0teckhvn+ge/YtxPIj3RoNc XnWR4/bMMLDzZNzV6tRoMOeGMiafthOMZ+gu9yck8t7lDURaum3hIvD+HTlj YkMKOKTtFcUg2wd21A6RPIck4C0TEbJt5QfArFaPX9nkdPCVyOSNv8tl4I4V tP+cVQqwh1ZjvnoTA68vsmFt7hlgUaBDTZk0E+zsK/hef/FrsJz9hTzdmgUO O1vGKPZXgvWrCuqpJdngIwq/ev7fOeA6Vxc12SkPnff9xKGn5qP7j06Zr+MF KP8enWzfdxSslbUphtWFYJV7wYOcwG/AAffwbyMsyC9sF32mi4+B06PafhB6 F6H9NAYdn+lC7vK5nco/UAye8PzM7hilBDy21sH2/BJyQ+jUmFNUKThQyWGn TSHnGr8rUZ4uA1vOPa45+FE52DWzdL66D9kB+zwZy1KBu/P2rCv1rgD/OLez Pu8P5C8vSD/cIjkOvsV/MLOSVgme9pLcaWpG5koNYp+oE2BRe22c2YoczMwN DT9VBX7a3lE5wjoJNvXzXy5ykddt9Y2YwOlvPciVzj9e1hHGAad6dK+m18WA 5bXnFy+1xIF/GRDwnht2Q94Qf3lIBhvl6yjWyUncBerjfCt8zoh3gQN8qUf2 qZOg/p/ioRHGNnTfkpGoTsP8oL58+eUePTUBPLDa93wrB9Uv1BYJZvFl4LRo cQ0+jf6/I7ihVbUV1edFbxjV7omHvMEAN9kHhLx89vqQQ9haOFe3MzuNV1D/ /hnWbSFhn7k6DYlSN8l9Z0oKDdc1ov5PmsPpwYT+g5yqXj1hXytTsD+jdqN8 C/Xax46E/Hbj9XMK7H0493Cf+zmOksjjGe9mY90KrvnZTVEzoV7pOCZZUzcM 88QefKIRFKJ5njLJCU7EeX6LfO8CYZ6BE4cljwj7ELZ4Tn2isueNG2ldvdqH 3Phnbh0uJsI+sKOnSwj9Wa96BLnYAvSvmm7w6mE4o/fb3n/trD+ax5Rqv7mX cP+rhuygsLr7cH/239jNZ9JRfdpssnWO8Pveq2rUlBHmF8XcId0lzM/QGLzc XKlLLhxu2hvPybqv4Vr1aH4h0+KHEfqbTJxGrmoG+s/7H2Iq17zpX7hxFztp k0TfSaiv7JTeil1YwUuuWFx6HnKD7VQs5ngCb65lh+21/1J9mmNm2EG9TigI uBGK3icJbzo5Ssj7H0YOwbM= "], {{{}, {Hue[0.67, 0.6, 0.6], Opacity[0.2], EdgeForm[None], GraphicsGroupBox[ PolygonBox[{{65, 92, 103, 80, 106, 88, 100, 71, 77, 63, 68, 55, 59, 5, 54, 4, 3, 53, 2, 62, 52, 70, 58, 91, 79, 105, 87, 99}}, VertexColors->None]]}, {}, {}}, {{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[{1, 94, 82, 73, 102, 90, 107, 65, 99, 87, 105, 79, 91, 58, 70, 52, 62, 2, 53, 3, 4, 54, 5, 59, 55, 68, 63, 77, 71, 100, 88, 106, 80, 103, 92, 108, 6, 95, 83, 74, 66, 60, 56, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 96, 84, 75, 67, 61, 57, 97, 85, 76, 69, 64, 98, 86, 78, 72, 101, 89, 81, 104, 93, 109, 51}]}}}], AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{ FormBox[ TagBox[ StyleBox["\"x\"", Large, RGBColor[0.6, 0.4, 0.2]], StyleForm[#, Large, RGBColor[0.6, 0.4, 0.2]]& ], TraditionalForm], FormBox[ TagBox[ StyleBox["\"c(x,0)\"", Large, RGBColor[0.6, 0.4, 0.2]], StyleForm[#, Large, RGBColor[0.6, 0.4, 0.2]]& ], TraditionalForm]}, AxesOrigin->{0, 0}, Method->{"AxesInFront" -> True}, PlotRange->{{0, 1000}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.405075222522332*^9, 3.405075286659236*^9, 3.405075372959513*^9, 3.405138828475061*^9, 3.405139058067793*^9, 3.405139119714719*^9, 3.405139208950967*^9, {3.405139258301105*^9, 3.405139279949971*^9}, { 3.40629402538591*^9, 3.406294054315035*^9}, 3.40629414862294*^9, { 3.406294558153632*^9, 3.406294562665886*^9}, 3.406294746197568*^9, 3.406303627932884*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Setting Diffusion Coefficient and Solve PDE", "Section", CellChangeTimes->{{3.406361638579706*^9, 3.406361764359718*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"CoeffD", "=", "16.0"}], ";"}]], "Input", CellChangeTimes->{3.406294705639785*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"NDSol", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "t"], RowBox[{"u", "[", RowBox[{"t", ",", "x"}], "]"}]}], "\[Equal]", RowBox[{ RowBox[{"CoeffD", "*", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"x", ",", "x"}]], RowBox[{"u", "[", RowBox[{"t", ",", "x"}], "]"}]}]}], "+", RowBox[{"(", RowBox[{ RowBox[{"f", "[", RowBox[{"u", "[", RowBox[{"t", ",", "x"}], "]"}], "]"}], "/.", "ParamRules"}], ")"}]}]}], ",", RowBox[{ RowBox[{"u", "[", RowBox[{"0", ",", "x"}], "]"}], "\[Equal]", "InitialConditionProfile"}], " ", RowBox[{"(*", " ", RowBox[{"initial", " ", "condition"}], " ", "*)"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"u", "[", RowBox[{"t", ",", RowBox[{"xRange", "//", "Last"}]}], "]"}], "\[Equal]", RowBox[{"u", "[", RowBox[{"t", ",", "0"}], "]"}]}]}], " ", RowBox[{"(*", " ", RowBox[{"periodic", " ", "boundary", " ", "condition"}], " ", "*)"}], "}"}], ",", "\[IndentingNewLine]", "u", ",", "\[IndentingNewLine]", "tRange", ",", "xRange", ",", RowBox[{"Method", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"\"\\"", "\[Rule]", RowBox[{"{", "\"\\"", "}"}]}]}], "}"}]}]}], " ", "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.40507523134622*^9, 3.405075255751947*^9}, 3.405075704849061*^9, {3.40513878566262*^9, 3.40513879283292*^9}, 3.405139288064379*^9, {3.405139328615452*^9, 3.405139330545434*^9}, { 3.405139418286599*^9, 3.405139445589445*^9}, {3.406294082840383*^9, 3.406294095295097*^9}, {3.406294748459499*^9, 3.406294750053563*^9}, 3.406303659172774*^9}], Cell[BoxData[ RowBox[{ RowBox[{"NDSolve", "::", "\<\"mxsst\"\>"}], RowBox[{ ":", " "}], "\<\"Using maximum number of grid points \\!\\(10000\\) allowed \ by the MaxPoints or MinStepSize options for independent variable \\!\\(x\\). \ \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/NDSolve/mxsst\\\", \ ButtonNote -> \\\"NDSolve::mxsst\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.406294750958071*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot3D", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"t", ",", "x"}], "]"}], "/.", "NDSol"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"tRange", "//", "Last"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", RowBox[{"xRange", "//", "Last"}]}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", RowBox[{"PlotPoints", "\[Rule]", " ", "100"}], ",", " ", RowBox[{"AxesLabel", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"StyleForm", "[", RowBox[{"\"\