(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    122693,       2380]
NotebookOptionsPosition[    120707,       2339]
NotebookOutlinePosition[    121080,       2355]
CellTagsIndexPosition[    121037,       2352]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Problem1.", "Section",
 CellChangeTimes->{{3.8126996046459904`*^9, 
  3.8126996080004005`*^9}},ExpressionUUID->"fd76de14-de24-4a76-939d-\
182b44b1f370"],

Cell[TextData[{
 "\:0414\:0430 \:0441\:0435 \:043d\:0430\:043f\:0438\:0448\:0435 \:0444\:0443\
\:043d\:0446\:043a\:0438\:044f, \:043a\:043e\:044f\:0442\:043e \:0440\:0435\
\:0430\:043b\:0438\:0437\:0438\:0440\:0430 3-\:0435\:0442\:0430\:043f\:0435\
\:043d \:043c\:0435\:0442\:043e\:0434 \:043d\:0430 \:0420\:0443\:043d\:0433\
\:0435-\:041a\:0443\:0442\:0430 \:0437\:0430 \:0440\:0435\:0448\:0430\:0432\
\:0430\:043d\:0435 \:043d\:0430 \:043f\:0440\:043e\:0438\:0437\:0432\:043e\
\:043b\:043d\:0430 \:0437\:0430\:0434\:0430\:0447\:0430 \:043d\:0430 \:041a\
\:043e\:0448\:0438 \:0437\:0430 \:041e\:0414\:0423 \:043e\:0442 \:043f\:044a\
\:0440\:0432\:0438 \:0440\:0435\:0434, \:043a\:0430\:0442\:043e \:0441\:0435 \
\:0438\:0437\:043f\:043e\:043b\:0437\:0432\:0430 \:0441\:043b\:0435\:0434\
\:043d\:0430\:0442\:0430 \:0442\:0430\:0431\:043b\:0438\:0446\:0430 \:043d\
\:0430 Butcher:\n",
 Cell[BoxData[
  GraphicsBox[
   TagBox[RasterBox[CompressedData["
1:eJztnQdYFMf7x++XWBK7YkMNoolGSTAa28+IFStGBYyJvVcURYwR7CWx994l
Yolib9g1aqLYO0YRjDXWqH/BaAz8v/+bv/NM9m6POW6RveX9PE/y4N3s3s67
73znfWdnZot16uvf7R2TydT/PfzPv+PAWsHBHUOb5cI/mvfp37N7n65dGvb5
tmv3rsH/7fQuPsz9H5NpD/77v7+TCIIgCIIgCIIgCIIgCIIgCIIgCIIgCIIg
CIIgiLQm8Q1pfSEEQRBOwP/8z/8cPXqUNJMgCEIGaOZ///tf0kyCIAgZSDMJ
giDkIc0kCIKQhzSTIAhCHtJMgiAIeUgzCYIg5CHNJAiCkIc0kyAIQh7STIIg
CHlIMwmCIOQhzSQIgpCHNJMgCF2h842DdK6ZOreeniHTOQJZ7+3Dbf7PG/R5
F3SomYkCOreeDiHTOQJZL01gFn716lVCQsKRI0fmzJkzbty4lStXXr9+/eXL
l/wupPVl/j9600zmrn/99detW7d27NgxZsyYqVOnrl+/Hsb8+++/mfXS+hr1
C4wDK7148eLkyZNhYWEjR45csmTJuXPnYE+9OZ7eYMZBC71///6ePXvGjh3L
HQ9tmayXerAm/+DBA19f37x585r+jY+Pz+LFi/ktSOuL/T90pZmsr1m4cKGn
p6fCdLly5erevTsMqx/T6Q2Y5eHDh3C8/PnzK6xXt25dOB7vstP6SnUH62t2
7txZuXLljBkzWjre7du3X79+TdZLDWDVmJiYmjVrwtr58uVbvnz51atX0dLR
c9WvX5/dhdGjR+tHNvWjmayXnzRpErOSn5/fr7/+CtNdunQJgXrmzJnxIS71
/PnzOjGdroBBrl27VqNGDWa9H3/8EY4HCT18+LC/vz/78PvvvyfZtIT11IsW
LcqaNSusVKVKFbTWe/fuoSGj/TLToS0fOnQIskmm0xbY89SpU5999hmMXLp0
aSREooXRkXXq1Indgrlz5+rE/vrRTNhn6NChzD6hoaFwY/GqDhw4wL7y8vJC
7kkNXwSmOH36NHM8YOl4ffv2ZV/NmTMH/yTTicCXRowYwezTq1cvdCuifWBM
HnPevXuXHE9b0My5Kq5fv97StkeOHGHffvDBB5cvX9aD/XWimbiACxcuvPfe
ezBOnTp1FH7LCnDbQlGp4YvAGoGBgcw4q1evtrQMPPPzzz9njoe4XQ+OpxNg
h4sXLzLHy549+61btywtM3PmTGbbNm3aWHomkWJgSeRBzLYwvlXbwlerV6/O
ygQEBOgh1NSJZqLVDxs2jFlm0KBBVq8H3RDv8eHnaX7NOgF2uH37tqurK3M8
RZDJy/Tv35+ZrkePHtTjcNAkeXbj7e1ttTd59OgRE1UQHh5OPY5WQAA7duzI
h9zj4+OtGnb06NGsTMWKFa2691tGJ5oJU3h4eDDL/Pzzz1avBy39/fffZ2Wm
TZtGrsuAHaZMmcLMUq1aNTWn2rlzJytTsmTJFy9ekOkYiG2446FtWjULzFWv
Xj2evFOPoxUwftGiRZPtyletWsWDpT/++CPNja8HzcSvR0ZGiqNGatfz6aef
sjJNmjRRDHimW2CHxo0bc81Us8muXbu4hZEQUY/DOHr0KDcL2qaaTXr27Kmr
UMcYREdHc+OrZZdJQncPNmzYkOauqwfNhBFgMW4WG2LYoEEDViZbtmwULDHQ
hHPmzMnMUr9+fTWb8LF0MHz4cD2MC6U5sMDcuXO5WX755Rc1mwwePJgX+/PP
P8l0miAGkDY0U+zu9fAsQw+aCSMgQOJmsWGT9u3b82IITcl1wY4dO7hN2rVr
p2aT2NhYXszLy4ui9CTzeBpybW6WmJgYGc3cvn07mc5xYEMxUpLUzM6dO6e5
6+pBM1++fGkSsKGZAwYM4MUiIiLIdYHoeN9++62aTeLj43mxcuXK0fPfJIvO
+tmzZzKauXTp0jRPDw0AbNi7d29u1bFjx8poZuPGjdPcdfWgmffv35fUTIXr
kt/CAmKkhMxFzSawKi9WsmRJGpRLMg8Fly1blpvFhk1Exxs/fjyNbDgObCim
jfPmzZPRzMqVK5NmAuREKdDMmTNnUneP6iMft1czCxUqRKPBSeYEp1ixYtws
Nhqj6Hj9+/cnzXQc2NDPz49bdc2aNTKa6e7unubdvR408+TJkzKaqRgAQXdP
mgkLiI5nY1AIMRUvljNnTtLMJLNmiptCSGomOinSTMeBDevWrcutunnzZhnN
zJIlC2lmkrAuMlnNDAkJ4cWGDx9OmgkL1KlTR0YzxUHjjBkzkmYmmW3C9jGw
SzPr169Pmuk4sGHDhg25VcPDw0kz5eHrp5LVzICAANJMEViAz78y2czNSTMt
gU3YvhykmW8f2FAcVpo5c6aMZhYtWpTGM8HVq1clNVMcNKbcPMmsmW3btrVX
M3PlykWamWS2ibu7OzeL5DMgGJw003Hgunw9r8m86ZbMnPaKFSuSZoInT55I
amazZs1kOqb0AywQHBxsr2YWKlQozRMcPfDq1atKlSpxs9joR0TNhMFJMx0H
mjlu3DhuVcl1QI0aNaL5mUnm5pwjRw4ZzRTz0CVLlpDfgrFjx3KboOOW0cwS
JUqkeWetB9D6vvzyS24WyfmZY8aMoQRHExYvXiyjmeKqjQ4dOtA6oCSz6374
4YfcLGr9CBzVy8uLF7MxOSFdgb6D26R9+/ZqNhH3gSxbtmyad9Z6QNxXB1y7
dk3NJuKEjUWLFpHpNGHLli3cqjZSJL7zs8m8aiPNOyw9aCZcF9fAzaIWAsFW
pUuX5sVoCRtj69at3CbIXNRscv36dV7siy++oOwyyexR4sqyo0ePqiU43bp1
48XQ0sl0mvD06VNu1b59+6qJ4YwZM3gxPayY1oNm4teR73CzqO2BgGZeoEAB
ViZbtmyUXTJgh+zZszOzVKtWTc3xzp49yy38ww8/pHlnrRPEDbU2bNig1ll/
9dVXvBhaOplOE5Ds8CcUjRs3tjqejE/4zrp58+bVw7NLPWgmuHLlCvfJY8eO
Wb0eWDhDhgysTJ8+fdJ8WEMnoCvhr66wsfnGnj17uIVhbTIdA20wX758zCwT
J05U66xr167Nyvj5+dGwhlYkCvtKubu7WzWs2GF17NhRD8bXiWYiWOKWWbhw
odUo6NSpU7zV79q1K82vWT/wCWx58uRR26d9+PDhrAx6dgrROeJbutq2bWvV
8dBOCxYsyMrMnj2bQnQNuXz5Mm/UZ86csbQtjP/RRx+xAmFhYXqwvE40E7aC
VDLLwIctY0j8c86cOaxAhQoVqNWLwBqwCTPOypUrLR0P9uSTatS6pHTL0qVL
mWVKlSpl9UVU4tre6OhoMp2GwBX5qyvatGkDR+XOmWhm/Pjx7NvOnTvrpNXr
RDOTzA92v/76axgna9asUVFRzHTcevwtYCZ97NWsK2AK2ERMzxWON3/+fPYt
LEwzMxXAXMzxwLhx4xSOh1bMVw1MnjyZnp1pzosXL/ieCRs3buSvk4apL126
xD739vZOSEjQieX1o5lwzps3b5YoUQImqlWr1qFDh3ing2bOJ8iNHDlSD2Ma
egMmGjVqFDPR2rVr8U/ueLCqm5sbPkeOg7/JdApgkFu3bjHHM5lfz8cdD562
adMm9nlQUJBO4hzjce/evXLlyjE7165dOyQkZPHixXywrmzZsnFxcfqxvH40
M8k82L59+3b2CkXg7+8P0/Xr149HmIgH9PDgTJ9AJ7/55htmqPLly8NusB5/
LlmwYMEdO3ZQfG6VRPMbqQoXLiw63qBBgxDesE/QfvUT5xgSdE9Lly4Vdzoy
mTfHnj59ut4yI11pZpLZe9GnIKrkyslo1aoV5JQiTNsgQEKQ2ahRI9F0UEs0
/9jYWBJM29y5cwe5eZkyZUTr+fj4REREkOO9BVhGibtw/Pjx/fv3X7ly5fnz
5zocDNGbZjJgKJjr0qVL+/btgwHv379PTisJyyjheMeOHYP1kGniFuvQ8fQJ
upX4+Hi0VrRZON7t27f58FpaXxqhF/SpmQRBEPqENJMgCEIe0kyCIAh5SDMJ
giDkIc0kCIKQhzSTIAhCHtJMgiAIeUgzCYIg5CHNJAiCkIc0kyAIQh7STIIg
CHlIMwmCIOQhzSQIgpCHNJMgCEIe0kyCIAh5SDMJgiDkIc0kCIKQhzSTIAhC
HtJMgiAIeUgzCYIg5CHNJAiCkIc0kyAIQh7STIIgCHlIMwmC0D+JAml7JU6n
mfoxHUEQb41Xr14NGjTo5cuXad7wnU4z//nnn6ioqPDwcPzhRJdNEESKQUuf
OHGiyWSKjIxM81bvdJqJjqZ48eKwnh56HIIgrKJtPoggs1q1au7u7n/99ZeN
EyZaw/FfV/AWNFPDi8dJEGFCMMePH//69es0t55uSc91dxAynePAaEgD/zbj
uAFxhh07dqDVt2/fXq3VszuFb6GuCQkJCKhem2HZqLY3MbU1E2d+ZUaTVBpG
8PPzg/WOHz9uw3QAN4v9LjomZsDUsJ7eUFT/pZl0UnfHEdsdMx33W7KeDNxQ
aHQzZ87UKpXGLejWrRvOtnTpUrVWjzu1ffv2jh07li1bFiXx/6pVq4aGhuJD
XIy2/p9KmsmtB8339fVFXK1Jj3P27FkYpEKFCjCR5dmYw8NEgwYNaty4cbFi
xXLlygXTtW7des6cOdHR0VpJtz5hnfv9+/c3bNgwePDgihUrvv/++5988kmH
Dh0OHz5s7Lo7DutofvvttyVLlrRr187FxcXV1RVNY/r06S9evOD9Tlpfpu4Q
A3KYCG62du3ar7/+2mQG8aHjRnv48KGrGfxhtdWjyUMe2S9mzJixXLly+fLl
M73Bx8cnLi4ON1erO6ihZorWg+lwnaiIt7c3LlsTzcThI0aMwNmGDx9uWX32
o+jXvLy8mK1y5syJ382cOTO3Hvo++L8hnZ81+X379pUsWRI1Ra1LlCgB2SxY
sCCre3BwsFHr7jis3U2ZMoXZKnfu3OiXYcn33nuPBS0//vgjwk6yniXwOjS6
sLCwkSNHIsyrWbOmScBxzcThCxYswKkQalpt9bhxLVu2RIG6devu3bv32bNn
LME8cOBA4cKF+ZXs3r3b9miePBpqJuyDLgb9cr9+/ZBB89aqlWbiDOXLl8fZ
zpw5Y3kq2HPTpk3s5wIDA5G8s9zqyZMnCOnFTseQ0gF/QHjJutdSpUqh+jAX
i7oRc7K6t2jRwpB1dxDW7oKCgpiVOnfuzAPL06dPly5dmn2O/tpqdpPOQSsz
/RsetGiimXBj5Iw41bZt2yxPhW/79+/PfvT58+cK3+aCALJly3bz5k1NnF9D
zURGLJoub9682mrm9u3bcSoY0PJU+CdUokiRIigAQyl8G/6PbItfTK9evQzm
/KjLunXr3n33XWbqhIQEsXb4e8uWLazuOpnhpivgTs2aNWP2WbZsmaJZ4dvq
1auzbw8dOkQ9jgK4E/SqadOmiAPHjh175MiRkJAQDTXz/PnzOA+CJattFjrp
5ubGf44PXbJv0fALFSrEv0Usp0moqa1mIhNv3bo1lP+nn366ffu2hpqJY7t3
745TIVC3PA9MMXXqVP5zoaGhCt9GbiXq+ZUrV4zk+XAnxJCsakjPrXbH/v7+
rAC6Hmr4HNhh//79zDJoX1YbJrI8VgAt9/Hjx2Q6BSyj4U4lxk4OaiaOhQ7j
PEOGDLF6HqaoYsMXVRF/s8yUgeRXk6cqGmomm13Arxnup6FmxsfHFy1aFLnn
3bt3rWpCmzZtROvFxcWJxW7duiV+Gx4ebhjPR0WQQrIgE61eLYzkfQrUlUJN
DjynZ8+ezDK+vr5WzQJP5iNjs2bNoh7HNhpqJixfo0YNnAfhq9XzREdHi+16
7ty5igSzXLly/FurKWoKSL25RhpqJg5cvXo1zgNhtBpd40P0RKL1YEyxWGxs
rPjt/PnzDeP2qDv3UthZbdhh586drEz+/PkVyXt65sGDB/wZ66RJk6yaBV0M
ck9Wpl27dgYb2NEcrTQTB0ZFRdlWD9yaihUrst/y8PA4e/asYlxFfAzUo0eP
9KOZ6NnZgOSKFSvUTrJ7927+c23btlWEUuvWrRM1E/fCMG6PmlapUoXbWa1e
u3bt4tXfsmULBUtJ5lYZGRnJzbJx40Y1m/Tr14+VKVWq1IsXL8h0NtBKM+Gi
Q4cOxUl++OEHGyc5derU7NmzFy9efO/ePYXIoGlkzJiRX8yMGTP09txcgYaa
iWAA0dEHH3zw7NkztZMg3Nq3b9+ECRPwf8UCK/w9evRofjEI122vwHIuxBAa
2YdavU6ePMmL9e3bV5MO19lBq0R75GZRWygB+NwD8Pvvv5PpbKCVZkLxKlWq
hJMookdLEs0Tsy2fC69cuTI1Wr3+NRNHLVmyBCfp2rWr7ejoHzNsgqj4OezJ
RkV4d6PVTC09sHDhQl61zp07q9VLHNFFXEoJZpLZMZo3b87NcuPGDRnNtBGO
EkkaaSaOYpljvXr17G2tTAHg4d988w2/EqSoOpzTrkArzYTF2FSQDRs22HsG
Zr2DBw/yK/Hz83v+/LlhfB4VCQ4O5rWDu6pVDfbnxcqUKWOkSDvFwEVZJMOw
YRNRM+fNm0cjGzYQNTPFayfR6nv37s2sLXmGxDfgWMSo4mVAMDV87vnWNDPF
gc2lS5cyZsxYtmxZyTMkCsB6yEn5tGQIpsG0AnVB+M3tHBoaKqOZxYoVM5gd
UgbaEfcNk82dskTNHDVqlJHyFM3RRDOfPHnywQcfFClS5NGjR5JngD4gGZ86
deqQIUMqVKjALqBXr17Hjh3T9n7pXDNxCBtxGjFihOThiAEuXLiAdB4HIjjP
nTs3DndxcRk9erTxhAKVFRMQSc3Mnz8/PchIMmumOOdZUjMDAgJIM23geG6O
QxYvXsxMLR/S85khHMRImzdvRqvXdrsAnWsmWjqyp0yZMkEG5TUTP6ewHoyP
eFWxTMAAoDo+Pj68mjZyc/F2ZMmShTQzyayZ2bNnt1czmzVrRk/QbOC4ZsK8
bHLX3r175Q+HMEI2EdniRyMiItgGaKBBgwb4p4Z71OhZMxPf7JvXvHlzu7x0
uZnw8PD169fPnDnz/fffx0ny5csH5YyKijLSxmioC1/ZZ1szFeuCSTOTLGwi
qZlabdJlVERbpUwzz549+5///KdmzZqOPKnEPeKyaTIvRtBqlxU9ayYqyNZL
2piWKcPWrVtdXV3ZlbD9KwyTW8FEfK20vGZSnMmATfLnz2+vZvr7+xvGf1ID
BzUz8c32ZVOnTnXQyLhNbF81xg8//KDJdBE9a+ajR48KFSrk4eHh+DgkdJJF
m0w2T5w4YYwkHVXo0aOHvZqJkNt4Q7spADb55JNPuFkkn5ujHyfNtIGDmgmh
+Pzzz11cXDSZB8t2PeKsXLnS8YavW81E4UWLFjEdcDwVgpNDe/n1VKpUyRhr
rlEFtlaCYeMZkKiZRYsWJc1MMrtorVq1uFlsxN6iDgwZMsQYHW4q4YhmovC+
fftwYKtWrWRafeK/sSwQEREhaiYaiONiolvNhMp99dVXOPbgwYP2TjGyWl58
VvLuu+/eu3fPGG4/Y8YMXq/AwEC1SkEkeTF0H8boMhwEzUecdfDgwQM1m4jP
NaZPn06aaQNHNBOG/fbbb03mTflkBBMq8cKM2ljlH3/8IWpmvXr1HI8WdKuZ
iMyzZctWpUoVmdadaN6lPMGM2ib2ffv2Fa23e/duY3g+272EAQVQqxFuNC+m
9nKQ9AZaXK9evbhZLl68qGaTAQMG8GKrVq0i09nAEc1EY/f09HR1dZXZcw+F
p0yZAm0pUKBAcHDw3bt3LVu04jFfjhw5HF/Sok/NRMnZs2ebzNMyk1U29ni9
a9euefPm9fLymjx5stXwW7yVYMyYMcYYlWLvSGIg01SrUUxMDC+GAN4YdXcQ
WADews1idetRVqxTp068mI1l6USSA5qJkmxf8Q4dOsj4J/o4/mzXZH5SbCky
d+7cEVu9JqvO9amZfEN7Gf9kL03jv1WwYEGrAYNCMxcsWGCMOBOGrVy5Mrez
WqVgSV739evXG6DimgBX4WYJCwuzahaYlO8FV7x4cRoKtk2KNRM6GRgYiKOQ
OskcpWjRiJcsc9ITJ06IZRCOOp5h6VMzL1++nDFjRkRNMkchHzf9G8vHx/in
4tUb+/fvN4bno0XzYAluo/Ygg+8Fh0SGJhpxIIB8K2+1pAZl+CqJPn360ORM
26RYM2HnEiVKFCtWTHJ/V8UGj3B+RXSKv8PDw8Uymmw2rkPN5BnTuHHjZA6B
qcXJq2Dz5s2KA9EWFBu537592zCeL+7wb/XZVqKwLL1nz540mMlBK+Njlb6+
vlZTwqdPn2bKlImVkQyB0hWKx6+iZkZGRooviFd7PstOsnbtWpN5KpdkAojI
iv+Qh4fHqVOnLFu9+CIwCJ0mOUIqvasXiKOv0Ez+YCtZ66FF16tXD0edO3dO
5sL44CdD8U4QVuC3337jbg+GDRtmpAE9eALaO6va1KlTFS6Hv69fv+7i4sIK
UGKuAAZhlnF1dbVcXZsobMNYsmRJw0y30ArW0tnDa1gPIaKomWPGjLl161Z8
fDy+Yv9Xsx6i91atWpnsCU3RhNlv5cmTB3rLF5Vzli9fLoZJR48e1eTeaaiZ
7M2wDBgwLCyMXy3C5j179sBoMCm+hZHVJAsfsiGIBg0ayCdBqEWdOnVM5pcg
8zdycnCe9u3b84vx8fEx3kyba9eueXp6msxvdYedRQuIL78YPXo0BZkK4B6j
Ro1i9mnatKmic4fHspdf4P+GGc/REHiX+HJe2wwcOFDNgBcuXEBUU6VKFbuG
PuDM/KU2nTt3PnDgwJUrV5AXQB4h14qf1uH7zcXXBCSL2pAa39B+2rRp8leF
klw08H8EBidPnrx79+6DBw82btwojmR6e3vHxMQYz/PZNvVsrdMnn3xy6NAh
9np32JmtRDO9eXO38eruOLBVly5dmJVCQkLQ6cB0aLzI/thI5rvvvrthwwYa
ybREoZnwwMKFC3/22WeVK1cuVapU/vz5xbdLwLZqkRJr9ePHj7fXwujU0N7F
XRcU9O7dG7Guhnll6mlmjhw53N3dK1SoULFixWLFiiH+kdFMeC8OgeWhgXZd
FQrHxcUFBwfzNZIKPvzww+nTpxtpw2EFaNGrVq364osvWH3d3NzgumwrPGSd
6DgowrQBHO+7776D08JcuXLlguk++uij9957z2ROTFavXm2kwRxtYbk562U4
lp+gmFqng5KffvopjJ+yeIaND8yePbtt27ZQj6xZs6Kx+/v7I33YvHmz5m6v
7bt62bCGaKtXZkTTsczd6jBv4puNjPz8/FLmojj5sWPHgoKCkNoXLVo0W7Zs
VatWDQgImDdvHnJ2w4sGcx4k4OXLl0cXDydEXI3OHa5o+Lo7Dpzz9OnTCDg9
PDzQ85YuXbp9+/Zr1qxhGzCS9VIJGHbFihUm86sPHbQze10415xUmkyYes/N
UwDf0H7hwoWOXBIbxBOtp5MKvh3gKqg+Km6jeyLUgLnYmHyqtjuCAwuzVdLO
8nRSV5oZHx+PLN7FxeXu3bs6uSSCIFKVuLg45IOenp7O8mRWP5qJa1i1apXJ
vHKKOneCSA+gmc+cOdNkXofiLK1eP5rJZwSl4P2SBEE4I69fv27YsKFJu8mT
bwH9aOYff/yRL18+hOj0tIIg0glsLU/9+vWdqNXrRDNxAXPnzjXZ3GmcIAgj
gZY+a9YstHq0fWdJzJN0o5noZdDXwHqWi0YJgjAkaPXNmzcvWrTow4cPnajV
60QzDx06xKZlprd5QQSRbnn58mWePHnatWvnXK1eJ5r5ww8/QDMjIiLS/EoI
gng7IB9ftmyZjVeK6BOdaOZff/21ZMkSJxoHJggifaITzSQIgnAKSDMJgiDk
Ic0kCIKQhzSTIAhCHtJMgiAIeUgzCYIg5CHNJAiCkIc0kyAIQh7STIIgCHlI
MwmCIOQhzSQIgpCHNJMgCEIe0kyCIAh5SDMJgiDkIc0kCIKQhzSTIAhCHtJM
giD0T6JA2l4JaWa6Qj+ORxDywF137twZGhr68uXLNHdd0sx0BVwuJCTk8uXL
dMcJJ+Kff/75/vvvTSbT8ePH09x1STPTD7jLCxcuhOPNmjWL7jihOYkqOH5m
9PXVqlUrV66c7Tgz9S5ARJ+a+Xbqrk9Sr+6vXr1q1KhR1qxZbceZ6dn4SSrV
T+uL0jvMSq9fv4am/fXXXy/M4I+///4bIaKDNsSxkZGR6Ov79OnDzqZ2AfBw
/Cg0Db+ekJCAi8ElOX4BCnSombgYVJNVnxk/lequN1jtFHXHH/inJo73+++/
Z8mSpU2bNmqvk2Y/gW9h8Pj4eGZ5/Hp6MH7Sv9udtq3e2MAy8BNkzf369fP3
9/fw8IC+ffTRRz4+PgEBAevXr8e3alonA9yvW7duOOfGjRvV/BY3a+3atZ06
dfr8889Nb6hRo8bQoUNPnDjBHdixiv4/utJM5pnPnz9ftGhR//79vb29c+TI
UaBAgQoVKkyYMOH8+fPcgdP6SrWHOR5u/aBBg1q0aOHp6Sk6HrwF7deRxotj
58+fj3OuW7fOhuNFRER06dIFLpEtWzbR8dAiDGz8pDf2X7ly5eDBgxs0aJAz
Z044XqVKleB4586d491WWl+m7mC9THBwMFcqLy+v7NmzmwSaN2+OqC/F3vvg
wQNXV1e0CKvv8sYnuEEtW7bED+XLl6927doIR5HI81/PkyfP5s2bNbyDutJM
VOrmzZvwWFbZr776Cveiffv2n376Kftk6tSpDvZZ+oRFd8uXL4dvsJripufO
nVt0PPShZ86cSbFwwW5NmjQpU6aM1REhfHLy5EnmeEjeK1eu3LlzZ8QMCscz
pPGTzNVHo0aVWWUbN27ct29f9B1ogOwT+CHvs9L6YvUFHBK9PExUunRphDq3
bt1iSfGBAwcqVqwoOjAyF7UExwYov2DBAhyOX7Hqt+joq1atyu4R1Iw1EHD4
8OFSpUrxXx85cmQKft0q+tFMXMPFixfRWlkdt23bxuuIPoK339DQUKvdjVOD
3GH79u2onbu7O3O8V2aQViCV5vfdxcVlz549KZhugfKnT59mjoffUhyOf0IQ
vvjiCxRo1arVpUuXuONduXKlXr16ouMZz/ioTmxsLJIaVsdly5bxOsIO3P41
a9Y0apeRYmAKJCBMMC2faN+7d+/LL7/kzoN0CfpmrwFxC9CF4fCoqCjLA3Gn
AgMD8S2898aNG4qTw+czZ87MLwDJmia3TyeayTKj8uXLs9qdPXtWcUkwHZdN
yIuRXBcViYmJQb0yZcqEqinqBbP079+f33eE3Pfv37e3+tDJESNG4PCjR49a
HgjbMsfz9fU9ePCgQlQvX74splobNmwwkvGTzO2OD4JFRkYqqgZr9OvXj327
dOlSyx4nPYOQkgnapEmTrMbh+/fvF0PNGTNm2BvsQQpwIDpuqwc+e/asUKFC
/PzwVfEacLOKFSvGv+3du7cmoaZONBM1/emnn1jV2rVrZ9Uz586dywp4eXnh
BqX5NWsF6oKs2WQei7DqeLhHYpbRsWNHe4M9nBb9ERzP6oHPnz/nYwKWqog/
qlev7siv6xlUZPXq1axqCJastqlTp06xAghajh07Zpi6O86dO3cyZszIjNOl
Sxc4ksJ78U/m2wy4MUJNeQOi5KhRo3DglClTrB6F2yFqMpJQ8QLwt+i6uBJN
psTrRDNRl+bNm7OqhYeHW72e27dv8+pv3rzZMNEOvCh//vysXn5+fteuXbOs
GosSOYhL7XK8NWvWMMezarQTJ06IJ+/Tp49COoYOHcq/LV68uB7WYmgFaoqu
ilWta9euVuuFMuxZMBgzZoxhHM9xVq1aJXpO9+7dLaMdxJZimUuXLslbD70z
1Ald1ZUrV6wedevWLfHkK1euVMSZn3zyCf8WIbEmaYJONDM6OppX7eLFi2rX
w+Pwvn37ajWim+YgoVA4nmXVDh06JJaxTOFtgLN98803cLzffvvN6lHwAXd3
d35yyIJ4AfgDl8S/zZ07t5HizOvXr/OqzZ49W00zW7Rowcr4+voaqfoOsnTp
UtEtO3XqZNmhKELBsLAwSeuh2N69e3FI06ZN1bQO3XeVKlXYmZEmKKQDNy5f
vnz8pwMDAw2jmfh1nncDG2EMj7Q9PT0Nk54rNBOOZ3lnb9y4IZZB3y1f93Pn
zmXIkAGOp9bL8IF0ULJkScWYJ1pBs2bN+E8jxzeMaKAWGzZs4FXbuXOnWr0G
Dx7Myri6utqVXRqbe/fulS1b1nZX/vjxY9F11bJsS9AKgoKCcMi8efNsHAJv
3LhxI9T7wYMHls8CxJ9GFGqY5+Ywjvh02EaT5GkUuHr1qjFcF3dWrL7VJ1yw
iXj3x48fL99ZDxkyBIegV7LteJGRkfBniLNCsfFVjRo1+E8PGzbMMMkpKjJy
5EhetTNnziSrmeD8+fPGqL7jwFXQ0XTr1q1z585r1qyx2ikfPHhQdN3Tp09L
Wg/tokSJEshrkICnwOA4ZNKkSfx3K1eurNWYkh40E63Szc2N185G0h0QEMCL
bdq0yTCue+rUKbRKf39/caKLSEJCguh4u3fvlqw7zoacBY538+bNlDneunXr
+O96eHg8evTIMGYXBzPB3bt3ZTQzIiLCML2GJkA51WaM45PZs2dz05UpU0Yy
SkcZNljasmVLu6yd+AYopDhNbsWKFVrdMj1o5rNnz0RBsKGZoutqaAQ9wFYs
qrmH+JgGiTZyIknHW79+PQ5p0aKFXSM53PHQFho2bMh/evny5UaSC9SuUqVK
vHY2RntEx5s5cybNOJIEhvLz8+OmCwwMlHwMwceQEb7Kp1SJ5qXHuI+XL18W
fzcsLEzDxx960EyEQCnQzPnz5xup/doAdWQDOwz5LAP26dKli8nOgXfmeLgL
cXFxouMtXbrUMCOZDJgRQTivoA2rio43bNgw0kwZ+EMcBkx97do1SbtduXIl
S5YsOETe5SCVvXr1QrImLpz09vbetWuXtvdLD5p58eJFGc3Eh2yVFmPq1Knp
QTMTzasaRfts2LBBstZPnz4tUqRIgQIFHj58KO94wcHBSIj4uhjueIaZqMCB
SIozom1opuh4PXr0IM2UAXLn4+PD7YaIUdJuKDN+/Hgc8t1338nbeefOnaZ/
U7BgwYEDB8bGxmq7VYIeNPPXX3+V1MzQ0FBebPTo0elBM/l6Xh7kSGoXyqxc
udJkXiMgbyg+f1vheNHR0cbbpwIiKS5xktRMRDLG6z40B/aJiIjgRoMLyUeM
KOnl5YWjjhw5Im9nNqyakJCAfh+pK08N4MD4W8N1r3rQTEUHYUMzxRQVbmyw
JmwJasfWoTN8fX3ln/3Bhdq2bYujfvrpJ7usBMdj26DFx8cvW7aM/TS0BWm+
wfapUMxFkczNkfqRZtoGxoHz8HEPCKC83/KMvm7duo7YWZGWItzSynv1oJls
+wgZzezatSsv9v333xup/VqSaN6zhW+tU65cObumpN67dy9fvnxubm7Pnz93
xEpcNgFEOAXbLOgWNGS+AkteMxFnUm5uA+a3LFA0mYcy7JrkA9uy2cJ2TUK2
imLlplavE9KDZoqLIm1rJtusjGHs8Uw2WYJvS2iv46HkkiVLTOrrAe26EnES
Y+/evQ0TbcKk4to6yefmVlcIEgw214ILJmTK3lgRelu8eHEXF5fff//dcSOz
ycmcFStWOO66etBMuO67774ro5lirzF//nyj+m2ieTtNnlngD3sfWONwtnhH
bYdhuxBDTQA1NoZowKq1atXi9bIxdVDUTDRDY3QZmsMEk21Dihxn4cKFot8m
Ctg4AxsFRXQk42OJ/8aywKZNm0TXHTBggONzP/SgmahFwYIFeb3UKoUPeaJq
Mtz8TA4TTPb4xmTeZSsFjnfx4sWMGTN6enrK9PLJOp5iyTAahTG26YBxvv76
a14vG3PaxZGxadOmGaDumpNofh8Qsj+TeVpRVFSU5YIy9sYQG9bDIZ06dTKZ
dzGV8Vuc888//zxx4oTakBHuqei6mmyJpgfNhOuWKVMmWc2ETfgGmyb1l4M4
NYnmuZH8uU94eLjCGnCqBDM26s5fbIroSNLx4uPjIYxqjqdYMoy8yRhrrmHM
Xr168XrZWBQpPnxctWqVAequLUwwmd9CMCFiliaaM2cOiwFsWO/JkyeFChVy
dXWFyyVrZDhht27d8ubNi9M2atQIcYKl9z5//tz0bxx3XT1oJmoqdvdqMQxu
irj9jiGX/cIU27ZtYxVEWmHpA+ydPgMHDrRRd7685eTJk8maCNYODQ0tWrQo
64UVu2kx4GYKx9u1a5cBjI8qsLiIsW/fPrUEh81AYMhYNb2RrGAmvpn0cvjw
YTXr4fPly5ejTIcOHWRGP3755RfRJ+fOnWspHXFxcWIZb29v2/GGDHrQTLB1
61ZeL7XFAghHc+TIwcpUq1bNYGtSkoR3kmbIkOHOnTtW3YY53vTp0204HjuJ
j4+PzIgQOmLxVX09e/a0TF4UO4Qg63/69KkxjH/lyhVer8WLF6t11nxuNjTB
GOMSGgJrHD16lBkHHQp/V5oIdIbZ0MbWOvxNIjLrJVFgzJgxoltandGkmPs9
fPhwx+eJ6UQz0U5LlCjB6mX5fgGGODQxa9Ysg43Doy5wJ1Qtf/78XDAtYY6n
FhElmRt4jx49TMntoMVR7MyJTMdSMxXzGBs0aGAY3UBF+K5Nil2+OTAIHzv6
7rvvjPH8SytgCrb2+eOPP4ZyIiVhb3n+6w0IgQ4cOLBjxw5mwGfPnqlZDw0c
7odcUiYURIEtW7aIbmm5PkuRR5g02thHJ5oJ8Q8JCeF9gaXrJgrvIMiWLZsm
8xD0Q6L5LWbM8Y4fP87eaC86HlIM9Jjc8Sy3CuQ8fPiQjQjdv39fxkQ4lehU
nTp1UnTEiRaLN62+iM1JEbeDq127ttUgBM7G627Xbs+GB6ZQ+I9t3NzcbMyK
QS/PpE/SwlBjfmZfX19FX8+GWMUNE5o1a2aY+ZkMaAXbVLlOnTqWMwDxTz45
0/LlC04NE0xxg51kHU/t1uPDRYsWoUy3bt0k43Dk3XzbWDie5Qg5zjN27Fj+
6+zuGMb4SYLjoS8+d+6cwm74m1e/fv36xhsRcoR79+55enrKa2a9evXUelu+
OF1+qBx3is8Bg9Lytb0M/JA4Qww4PpLJ0I9mQgbZkzWTeS05X1bP4Clk9erV
jfHQlgMBbN26tbzjoeWqdRl8WibSFkkToRh7rmR6s1BCMSzAhuUZ2bNnt7Ex
r5OCxsVfE4A8XWGB58+fZ8qUyWROAdTeSpNuEadgydC3b1+1vv7IkSMm8wuR
7QqHUJgLI3QMmfjevXuh5PB/sU2VKlXq1KlTRlpvzkEA07t3b1bN9u3b79u3
D10Dgn/u0ogETpw4YbCRTJ5xSxIUFKTWWcfGxmbNmvWzzz6zKxzC2djcJFC1
atXJkydHRkYiIY2Kivr222/57+KriIgIw2TlIjAXd7y2bdui3bGxkXXr1rFF
05kzZ969e7ch6+4I0Mzq9oD+16oB0aIHDBhgMu9Naq+F2bstxIk3IkjKZs+e
rVWEydCVZiaZ2y86Dr5DV8mSJd977z32t5+fH9qsfi5VK6CZdjleSEiIWmfN
3vEHT7bXSuivN23axF/9qaBw4cLjxo2D4xmstxKB4w0fPpw/iITj8UkayBnl
9zFLV7y2H6s2RPcEyxcpUsTGQL1t4MCHDx+G5AYGBjZo0KBDhw6TJk1Cmv/w
4UPNb5zeNDPJfCPi4uKgDOjxvb29fX19g4ODEXMadShJK8eDfRo1amSycwct
ETjeL7/8MmfOHBi8adOmkNAxY8Zs3779/v37RhpAVgM9wpMnT5DfderUqU6d
Ok2aNOnVq9fatWsNM0lAn8C2bHFunz59HNc33ETWRlKvf9ehZnJemdF2v1Cj
kmheL/mf//xH7eGvXbwFx9MzzPHSbfXfMrAzG4T/+eefncLaetZMQh60bjbF
d/LkyXQ3CSciOjo6U6ZMNWvWdJZ4njTTGMDfqlSpkj179tjYWLqbhLPAp3Kh
r3eWqJ400wDg9u3Zs8ckvYMWQeiEv//+m/X1MTExzuK3pJkG4PWbra1pyx3C
uWDrwdHXO9ETXtJMA/D48WM3N7cSJUo4y4gQQSSZ8yM2MVht3qY+Ic10dvhC
3aCgIGcZESKIJPP8hLp16xYrVgwq5ER+S5rp7CAxr1evHjRz7969dB8JJ4Lt
Mdi9e3fnyo9IM50daCYcr3Xr1gbbOoMwPEiLbty44VyCmUSaaQhiYmJIMAni
7UCaSRAEIQ9pJkEQhDykmQRBEPKQZhIEQchDmkkQBCEPaSZBEImENM+ePatc
ubLam2EJgtAVqaSZ/oQ0TZo0yZMnj5+fX1pfCEEQtli0aFHqaWYMIc3Zs2fL
li179erVtL4QgiCSIfU0k5CHxjMJgiDkIc0kCIKQhzSTIAhCHtJMgiAIeUgz
CYIg5CHNJAiCkIc0kyAIQh5n1MzUXhNBEFYhxyOSnFAzcamvzdD7woi3CXO8
v//+mxwvneN0mgmn7dOnT0BAgBO9EJkwAJDK0NDQhg0bOt37awhtcTrNvHnz
Zs6cOf38/Ox13bew0QFhYOBv3t7eHh4e9r56iRzPYLwFzdTQZ3CGKVOmmEwm
e/crYIkVQlOEqYZ33be8JY6uSKUdgXD4wYMH4XjffvutXbk5SsLlXplJJ7fA
8KSqZjJfjYyM/OOPPyBZjv8K3K9atWp58uS5c+eOXX4Lj/3qq6+qVq2KizGw
6zKDo7IvXryIiYnZv39/bGws4iI+CpdO6n7ixIkdZqKjoxEfOj76jTMEBQVl
ypTp1KlTdp0Hlm/evLmXlxfyegMbP12hoWZa9uxw12XLlqF3htA5HuDh8D17
9uBsLVq0sOtsaC8zZ840mTGwZjLFGDRokJ+fX9GiRU1vcHV1RYM9duyYgR+c
oVIJCQnh4eGVK1fOli2bSaB06dI9evR49OiRI73G06dP3d3d0e3aFS6iJHc8
3BdDWj4doqFmwifZE2341fHjxydOnFi3bl3mMJpoJs7fu3dvnA06bJffRkVF
ubi4GFszUSnEkz179uQtdMmSJRs2bJg2bZqHhwc+KVy48MaNGw05NIEaIbZs
27Ytqvnpp582a9YMdggICIB+cuX86KOPLl68mLJkB4f8+OOPOMmMGTPscjz0
U8iJSDMNhoaaCa8IDAysVasW9xOOJpqJS0VfX6RIkT///FPyVCz0+vLLL/mV
GFIzWUjfqlUrVseffvpJDIfwFW4x+2r37t3GCzVR2Y4dO6J2nTt3fvjwIfc0
1HTlypWiH545cyYF1ccJocM4/MKFC/LH4ijR8UgzDYO2msk9pEKFCgMHDtRQ
M3EsS/O7dOkiHy2ggYwbN05sNYbUTFRz5MiRvG1a2icuLo5lrD4+PgjJjCSb
qAh3vJCQEEXt8MeIESP43a9YsWIKZgpdunQpU6ZMVatWlfdhFBs/frzoeKSZ
hkFDzdyxY8eePXuQAT19+pQNxWuomdCB5s2b41TIN+X99ujRozlz5vT09DS2
ZiIrL1OmDKvguXPnLCsIGWnfvj0rMGvWLE2ex+kE1AUdAb+/LVu2VKji0qVL
Re26cuWKXXXn6icvemw4SOF4pJmGIfWem0M2NdTM33//PXv27B4eHpJxAsvK
69evj18PDw83sGaiOryCsLPaQ4rly5ezMnXr1rV3hqGegWbmypVLVEV01mLt
2HNDzv79++2qO+xZq1atd955B3m95IFw9SZNmigcjzTTMDiFZuLA2bNn4zwD
BgyQzCtRbMiQITike/fu27ZtM7BmwrBt2rThdlar3a5du7gRkKobxgjQTPZk
kNGgQQNFjxAWFiZq5rFjx+x9joOjkONIBucoM23aNBzSo0eP7du3k2YaD6fQ
TD6cfvDgQUm/3bp1q8k8z+Tu3bs7duwwsGYi8HZzc2O1q127tlrtdu/ezY2A
9NxIQ5pr166tVKkSu91Dhw4Vxc1yXNGumb18oHj16tWSjnfhwoWCBQviSvBD
O3fuJM00Hk6hmdHR0RkyZLCReIqgwP379z///HP8LhJS/NPYmnn06FFeuxYt
WqjV7uLFi7wYMkcjTTpic40QSOP/imgQosceqTMaNmxoV8XZlAP4kuSIEELc
li1bcscjzTQk+tdMHDVx4kScZOzYsZJZOUvW+BN2Y2umGEf16dNHrXaPHz/m
xcqXL58eNppABY8cOSIGmYhI7UrMf/75Zxw1atQombAcBebNm8ccj3k7aaYh
0b9m4jzVq1fPnDkzok0Zv920aRN+sWTJkjdv3mTlDayZqEvfvn1lGibsz4uV
KlXKSI+BRBIFUGVfX1/ROHZNGEDhfv36vfPOO+fPn5dxvN9++83NzQ2Od+PG
DVaeNNOQvDXNTNkeBTiEDcQ1b948WdVNNK+FYfO3w8LCeGFja2anTp3s1cwP
PvhA8XDZGCQKG2L8+uuvfn5+omXs7bWfPHni7u7u4+Mj47o4eevWrRWOR5pp
SN6aZqYssEFO1KNHD5N52kayh6Pw4MGDUbhjx45iRGFgzUSV2bTVZBumeDvy
5MljSM2ENdq2bfvxxx9nz56dV7ZZs2Znz561VzBRePHixTh8+vTpMkHm/Pnz
meOJAkuaaUh0Hmc+ePDA1dX1ww8/jI+PTzbIPHDgQIYMGVBYMZfG2JpZr149
mYYpLjHInDmzUTWzmpny5ctnzZqVVbZRo0ahoaHPnz+3a48ONlXj/fffj4mJ
SdbxWFYOx4uNjRULk2YaEj1rJsrPmTMHhwcFBSU7EoXzM/Ww3FrT2JpZu3Zt
XjsbG469fPmSF3vnnXcMqZlJ5kFI8NIM+tDSpUuzKnt5eR07dkzyzRR8Wqav
r2+yASp+jmXlcDzFyUkzDYmeNRPeWKdOHRz+888/J9vXT5gwASXbtWtn6eTG
1swWLVrYq5lIXY2qmQpOnz7NZbNUqVJRUVEySTqsGhISgkPmzZuXrOOxrByO
Z+nhpJmGRM+ayWZ61KxZ0/ax+ApNw8XFxd3d/erVq5YlDayZqEtgYKBMwxQ1
s2DBgkZ9bm4Jok1eccgmAshkQ82EhASUzJ07961bt2w7HsvK4XhWl7GTZhoS
3WomCrO+fuLEibYTc7YVNkouWLCAj1mJQCf5lWzfvl1RRou6piWjRo2yN84s
Xrx4+tFMOE/OnDl53ZFH2647vlq1apXJvN2HbaflWTlCTb6fs4jYWePWGMzx
0i261UwoYdmyZTNkyHDp0iXbB7IxT3j4tm3bII87LNiyZQu/kk2bNrEyEE+E
AQZwXTaPmtGlSxe1GkEoeDFPT0/DvJ7Gsou0LFCjRg1ed+Qjtp8nQv2++eYb
k8R6yblz5zLH27p1K3e8SDPs782bN/PfHThwIPuc/d8Yxk+f6FYz2bTMRo0a
JXuguG+nXUA2DeC669ev5zXy9fVVq9HDhw95scqVKxtj7WSieUJmdHT0pEmT
0GMip7aad7PdWjgoaaPuV69ezZIlS+nSpZNdKpVix7N9AYTO0admomRQUBAO
nDlzZrJHIeupnhz8Sry8vHAx+AT/N4Zm3r17V7SzWo0gBbxYQECAMfboQEy4
cOHCvHnz8npZfbbVr18/UbImTJigVnf+MPG7775L1kTQTNteh9th1fFIM50a
fWomCqOjz5o1q2LCm1X+sQmalRgPwF3Z0BPDAK7Ld2tkdlar1K+//sqNcODA
AQNUPMk84MD3wQNubm6PHz9WVA3/ZKOOMtWHMSFuKANzOeh4QHQ89OwGc7x0
iw41E8XYCKS/v78mb15ji4MYxogtRRLfrFhhwYzVjDLxzVvAQJkyZQzzAAhR
ZY4cOUQ9tJwXoZjCmjNnTvi8WvX37dtnMm/LrMl4r+h49K5ew5B6mik+qIVm
ym+kAydn709csGCB4xeGMxhbM8GdO3f4XuW//PKLZRiDCIcvvrZ3qwo9A81k
M3g5165dU1QN6icupfzyyy/V9DDxzcQtu94vaQPR8WiukWFIvfebi4v1EP8g
trGcjGH1PDiwWLFiefPmtWt7WLUrUaRIiGCNN+UDGhgcHMwq2KFDB4Uk4m+e
mLu6utq1UbnOSRS2dYKPrVmzRtFf4G/x1ZMFChTgmw5ZAhdljmd7Wmayl8RR
5ObGc7z0iYaauXPnzq1bt27btm27GTiJGADgn99//z0K4CuUUQs78SE8H+Xb
tm2b4nAIzrn9DfBb9oJvRoMGDUaOHLlixQpciTHmGjHi4+P5vmf79+8Xmyek
gO+7u2rVKsMEmQw4kr+/P6vdsGHDEEOKdT98+LD4tqAjR46oVR8fRkREmFSW
kskTGRnJHU98+Ojm5hYSErJs2TI4Hs01cmo01Exx6niyqK3dg0uzBg7lTPFV
sb0akr0Gg+Xpd+/eZW+fzJw584QJE06fPo1Pdu3a1bRpU1bfiRMnGmZapgg6
Bd5f1K5dG/n1vHnzVq9ejdj7448/Zp+7u7tfvnzZxkAu39Fd8jUWashMQKLn
5k6Ntu/qDTUDtxlsE5RR08yHDx8WLFiwRIkSjrRu+L/llQwyI16DwdZRgmvX
rqFqhQoVUjRSHx+ftWvXGizCFEG0uX79+kaNGlnVqICAgISEBNse9ejRI+Z4
Dj4gU7iZpePhb4N11ukNbccz5V8KabUkPmFvCQwKCnJwPkayhxt1ygeEMS4u
buXKlZMnTx45cuTcuXM3btwISTFkZRVAFdFroCucNWsW6r506dL9+/dfv35d
ZrdqNq9AZgct2yQ7VmlUx0s/pN5z8xQA32a79Ozbt08nl0Q4KXY9Z4FOMsfb
u3cvOR5hG/1oJq4BQUK2bNkqVapkjJV9hLMQGxvLHC89vFqOcBD9aCYSFqST
6OtHjBihh+sh0glwNuZ4w4cPp6yZSBb9aOarV69q1qwJ1z1x4oQerodIJyCp
YY53/PhxcjwiWXSimbiAgwcPsrkiBn68S+iQQ4cOwfG8vb0NORGL0BydaCZ0
ki3omDp1appfDJF+gLMNGzYMjjdlyhRyPEIGnWjm48eP3dzcHFy2RhD2wkaE
8uXLd/PmTXI8QgadaOaSJUvQ13ft2pUG4Ym3SUxMDHM8mqpBSKITzRw7dmyh
QoWioqLS/EqIdMWJEyegmeHh4eR4hCQ60UyEl3zjo7S9EiJdwfYwoceOhDw6
0UyCIAingDSTIAhCHtJMgiAIeUgzCYIg5CHNJAiCkIc0kyAIQh7STIIgCHlI
MwmCIOSBZl6+fJk0kyAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAI
giAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAI
giAIgiCIFPC/g3x97Q==
     "], {{0, 273}, {444, 0}}, {0, 255},
     ColorFunction->RGBColor],
    BoxForm`ImageTag["Byte", ColorSpace -> "RGB", Interleaving -> True],
    Selectable->False],
   DefaultBaseStyle->"ImageGraphics",
   ImageSize->{301., Automatic},
   ImageSizeRaw->{444, 273},
   PlotRange->{{0, 444}, {0, 273}}]],ExpressionUUID->
  "5929afdc-9e06-4617-9ea4-fdfd8854b226"]
}], "Text",
 CellChangeTimes->{{3.8126994908663597`*^9, 3.8126995746972456`*^9}, {
  3.812699612871337*^9, 3.8126996503556414`*^9}, {3.812699750815791*^9, 
  3.812699759035096*^9}, {3.812706340955819*^9, 3.8127063464210787`*^9}, {
  3.812706377860961*^9, 3.8127063993182735`*^9}, {3.813035102390889*^9, 
  3.8130351037352543`*^9}, {3.813035214902547*^9, 3.8130352690049067`*^9}, {
  3.8130353906817045`*^9, 3.8130353946027346`*^9}},
 FontSize->18,
 Background->GrayLevel[
  0.85],ExpressionUUID->"87dcf015-e835-4d58-a9f0-e1cc62932ea9"],

Cell[BoxData[
 RowBox[{
  RowBox[{"RK3", "[", 
   RowBox[{"f_", ",", "n_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", 
  RowBox[{"(", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"h", "=", 
     FractionBox[
      RowBox[{"T", "-", "t0"}], "n"]}], ";", "\[IndentingNewLine]", 
    RowBox[{"solution", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"n", "+", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"solution", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", 
     RowBox[{"{", 
      RowBox[{"t0", ",", "u0"}], "}"}]}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", 
       RowBox[{"Length", "[", "solution", "]"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"k1", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "]"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"k2", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", 
            RowBox[{
             FractionBox["1", "3"], " ", "h"}]}], ",", " ", 
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
            RowBox[{
             FractionBox["1", "3"], "k1"}]}]}], "]"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"k3", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", 
            RowBox[{
             FractionBox["2", "3"], " ", "h"}]}], ",", " ", 
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
            RowBox[{
             FractionBox["2", "3"], "k2"}]}]}], "]"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"solution", "[", 
         RowBox[{"[", 
          RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", "h"}], ",", 
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
           RowBox[{
            FractionBox["1", "4"], "k1"}], "+", 
           RowBox[{
            FractionBox["3", "4"], "k3"}]}]}], "}"}]}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "solution"}], 
   "\[IndentingNewLine]", ")"}]}]], "Input",
 CellChangeTimes->{{3.813053956565979*^9, 
  3.8130539600634584`*^9}},ExpressionUUID->"1b56c63c-48a1-429d-a6c3-\
7377e0ae88f0"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Problem2.", "Section",
 CellChangeTimes->{{3.812706327392376*^9, 
  3.812706331687019*^9}},ExpressionUUID->"cd97b77b-6eda-4198-be3b-\
ddb0d0e4e616"],

Cell[TextData[{
 "\:0414\:0430 \:0441\:0435 \:043d\:0430\:043f\:0438\:0448\:0435 \:0444\:0443\
\:043d\:0446\:043a\:0438\:044f, \:043a\:043e\:044f\:0442\:043e \:0440\:0435\
\:0430\:043b\:0438\:0437\:0438\:0440\:0430 4-\:0435\:0442\:0430\:043f\:0435\
\:043d \:043c\:0435\:0442\:043e\:0434 \:043d\:0430 \:0420\:0443\:043d\:0433\
\:0435-\:041a\:0443\:0442\:0430 \:0437\:0430 \:0440\:0435\:0448\:0430\:0432\
\:0430\:043d\:0435 \:043d\:0430 \:043f\:0440\:043e\:0438\:0437\:0432\:043e\
\:043b\:043d\:0430 \:0437\:0430\:0434\:0430\:0447\:0430 \:043d\:0430 \:041a\
\:043e\:0448\:0438 \:0437\:0430 \:041e\:0414\:0423 \:043e\:0442 \:043f\:044a\
\:0440\:0432\:0438 \:0440\:0435\:0434, \:043a\:0430\:0442\:043e \:0441\:0435 \
\:0438\:0437\:043f\:043e\:043b\:0437\:0432\:0430 \:0441\:043b\:0435\:0434\
\:043d\:0430\:0442\:0430 \:0442\:0430\:0431\:043b\:0438\:0446\:0430 \:043d\
\:0430 Butcher:\n",
 Cell[BoxData[
  GraphicsBox[
   TagBox[RasterBox[CompressedData["
1:eJzt3f9vVPe95/Hs7i/74/4JgIRjnAJJ7pUaKSyrZQApNkS1m0oYchUb0gJp
sbO9t2RviO3uxnalC4mCHfaGVrq2Idc2zl7sthSiBjsq2FcKzPBDcdRkHCWM
qxDGKmFQCWPF9sd37OPzZcbz7cycMzPnPc+HTpuUzpk5b+bzfs2Zcz7nzNr9
zXU/+s+PPPLI3//X2H/VNb7yP3/608Z//P5/i/2PHzT9/eGDTT988Zmmf/jh
wR/+9Kn9/yX2h7/7T4880hv7z9K/LwIAAAAAAAAAAAAAAAAAAO9TFsXeFgBw
0pq167SFfAMgDPkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkG
QCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3
AFKRbwCkIt8ASEW+AZCKfAMgFfkGwC0qjUK8PvkGwHlaht0Njp5tPbh7x8aV
nNm0bdfB184M//HLhYWFQqQc+QbAYbE0mYv439mrx9rqZdOet8bvfKulnIsb
Qr4BcJaaC4+27UiRbOayfv/ArW/dTTjyDYCTlAr21Jk55jt0+kooEl0OmNmZ
qctv733SjLhtHf6HbgYc+QbAOUpN9dWa4dbhf7AqWr75uPsHG/Tk2X7ypou7
cOQbAMeo6HhbpZ4qjQPTKXLl6/ebH9Mz8PDwPdcCjnwD4BQVudSkR0pF20Q0
daxMv/t9PXxqez53K+DINwBOiVxo1iOlsv1a2lD57EyNHj4vnLvjUsCRbwAc
ovydxpfT5ouR9KFy44TxFfUnF++RbwBKWmigXo+s3X3BDJkyM/Ki/uANx6+7
swNHvgFwhgq0G2dOj1zKsPu2uBj4hTFRpPH8HfINQAkLjzQa+dYZyBgp0wP7
9Pxp6P+SfANQwr46Z+RbxRuZ823mvf1G/mQ4GZEr8g2AMyz51jAUJt8AyEG+
AZCKfAMglSXfXjiX+YQo+QbAMyz59vTbNzNGSnio0bjY4Xjm0xG5IN8AOOP+
+83G/JAOf8YZu6H+veaV+LfJNwCl7EaHkW+vjj7MkCnK32nM7814MVeOyDcA
DpkZbNQja8/ArQw7cOHhA/qDq074uf4UQGkLvGncuLLlcvodOBU4XmXjYq7c
kG8AnDL7YZseKY++/lG6HTg11VdtZzJJbsg3AI6xnGJ4+q3Udx6P+42Gut4v
uH8vgNI3O/Fz4ytqQ38o2S8ALv0u6qUm48vpSyMufTldJN8AOOuzM3XGWdRt
7VfvWX7KWfvng8mu54zbYO44HkhzF/N8kW8AHDUbeHOnESxrHq8/+ftP7zxc
TrloOPhBd/0TxrSQdb5ON9ONfAPguPnpkR8/bkZciqXiwGBojt+vB+A18/cD
p/dtShlum+tPToRdDrdF8g2Ae76eGnu37dCzO/Wg2+zbvbfpjUG/O1djrUa+
AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCry
DYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKR
bwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCK
fAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU
5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwBy
oOLkvOKqNVUGtjbSwXzzRL0AHKBUqL9xOToaB27b6kJ189TW5RUrW65E49bU
+jka8g/9suVgTfVTlXo6bdq26+Br/zxwPfRwYWHBVtM7lm8eqRdAvmINNz3Y
sNKMNvtdBXt3ayu2jlrbPfavc+Hxk3s36omUdPnOvq6rd7Smz+rVnMk379QL
IB9Ka/Yqowft9bu6rQfFkUsRcz2lHkx2PVeZptPNZUND/+fZtnz++eategHk
avnr1Cd9lma32+8qcqFZW7FhKGysplR4+ICl2Z+oazkbCEX03Z1vZqYu/6q5
ZoOl5Vsu/zWrjs8v37xXL4BcLPd6cOgVX+IOhp1+V5GLR7S1anqnzHaPjrVW
GE/o6xiPJHvC+elLP9tpvO761vGHWTR87vnmzXoB2KMdBr890fX85mRfoOz0
e3SsRVtrd1/QbPdQb63l2aZTP9v8x93bjUcePP+XzA2fS755uV4AdsQ6PdB7
rC7+MHilz7c1h35X1zq0/ZYtpybNdYwjVGvXVbRNRNM+2f3fHDY248hv77mQ
b96uF4AdNzoSdmCqao7+OhR6b7/9fjdnSrRfsxxq//dO48lbxtK3++LiZ2dq
9Ac/8y+fupBv3q4XgB1x/b7x+e7x5e6eyaHf1WTXFm2VuJkS96+faT20q+bZ
XdufXGfNgeQsOz9rOgNu55vn6gVgx0q/VzzV2HXFPAOYQ7+rqb5qbZVjGfda
Uj6Hmugwwsd6RjKVnPPNo/UCsONGd93Lp6+EonNxvWW/31V4SJv/v67pQtLz
hVlQ0fE2Y1pFZeadn5zyzdP1Asif7X43Z0rYvb7JeAalrF/WqlrHs9gryuX8
aTJeqRdA/uz2u0o6U8IGpeZCA5Y5sbU9WT1NsfKtWPUCyJ/dfo+Ota5ZPVMi
a0pF/Z07zHOavhP+7HZmipVvxaoXQP5s9rvyd2q7Ilu77Ld7bIXQUKM54X/t
juyfpEj5VrR6AeTPXr8bMyWqOvxztlp16UqCUH9csx+9EMr+KYqTb8WrF0D+
bPV7jjMllq8AvXmqztLsmw/+m71mL0q+FbFeAPmz0+8q2FNje6aEWro92mib
5RhUTnsyxci3YtYLIH82+l2FBxq0Rx4avpvlzH+lHgR7rXcQqqppH8tldmsR
8q2o9QLIX/b9riKXmrRHNgxm07FL39K+GDnqW2dp9sbeT3I8fVj4fCtuvQDy
l3W/K2OmRHUWM7hiD4gEumut95P0tY7msSNT8Hwrcr0A8pdtv5tXGGWeKbF0
6vBC3P0kKw70BR/k1eyFzrdi1wsgf9n2+1ygvSqrmRKrm93XNha2N7kiiULn
W7HrBZC/LPtdTXZv0fZM0t/bR6notQ5Ls1fWnpp05AhUgfOt6PUCyF92/W7O
lEh7I8eEH6iqbOgJOtXshc234tcLIH9Z9bv5QwPNF1NPBEv4MSlfp9/4VeNs
pN/OguZbCdQLIH/Z9Lt5e590MyVU9Irlx6TW1v/TxfGrf5gYv5LF8ic3fn/B
2/UCyF8W/W7+9Gft2dTz8FVk+CXL7Ahbiyv3J/d2vQDyl7nfVXT0mPYA609/
rmLMhvV6vpVGvQDyl7nf5290PLb8gC3dN9O05VfnGmXkW2nUC6B8OJVvAFBq
yDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCpyDcAUpFvAKQi3wBI
Rb4BkIp8AyAV+QZAKvINgFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZA
KvINgFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCpyDcA
UpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCpyDcAUpFvAKQi3wBIRb4B
kIp8AyAV+QZAKvINgFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8A5ADFSfnFVet
qTKwtZGFzjfHNrxoPPG2Au5SKtTfuBwdjQO3bQ1PdfPU1uUVK1uuROPW1AZ6
NOQf+mXLwZrqpyr1dNq0bdfB1/554Hro4cLCgq1uKFy+aRt/Nzh6tvXg7h0b
rVt+ZviPX+obXuJ97JG3FXBRbCRODzasjFKbjaCCvbu1FVtHrX0Q+9e58PjJ
vXoyJF++s6/r6h2tG7J6tQLl29LGR/zvpNn4TXveGr/zrY0tLwLvvK2AS5TW
BVXG4LTXCOq23kFHLkXM9ZR6MNn1XGWaFjCXDQ39n2fbC4XJt1gLj7btyLjl
6/cP3Pq2RJvYW28r4ILl7xmf9Fm6wG4jqMiFZm3FhqGwsZpS4eEDli54oq7l
bCAU0fcDvpmZuvyr5poNll5oufzXrFqhEPmmVLCnztw236HTV/Rtn41t+dt7
nzT/urZ1+B+WXAt7720FHLbcBMGhV3yJn7x2GkFFLh7R1qrpnTL7IDrWWmHm
Q8d4JNkTzk9f+tlOc1+odTyboHA/35Sa6qu1bLz/waqX+ubj7h8YXbz95M1S
amFvvq2AY7Tjw7cnup7fnOybhZ1GiI61aGvt7guafRDqrbU823TqZ5v/uHu7
8ciD5/+SuRNczzcVHW+rzLzxX7/f/Ji+5YeH75VCC3v5bQUcEmuBQO+xuvjj
w5U+39YcGkFd69A+0LecmjTXMQ7drF1X0TYRTftk939z2NiMI7+9V/R8U5FL
Tdlt/PS739c3prbn86K3sLffVsAhNzoSPtmrao7+OhR6b7/9RjCnELRfsxyD
/vdO48lbxtL3weLiZ2dq9Ac/8y+fFj3fjONOCUUlYdnyF87dKXbAefttBRwS
1wgbn+8eXx72Mzk0gprs2qKtEjeF4P71M62HdtU8u2v7k+syRMRi3F7Bms5A
sfNN+TuNL6fNF5MeXzLdOGF8Rf3JxWLvo3j7bQUcstIIFU81dl0xT43l0Ahq
qq9aW+VYxo/zlM+hJjqMrrSeqkvF5XwLDdTrXWk99JTczMiL+oM3HL9e5B04
b7+tgENudNe9fPpKKDoXN+jsN4IKD2kT49c1Xciwo5P6OawH8zN9H1zmbr6p
QLux1xE37yu5wC/M4+3n7xQ53zz9tgKust0I5hQCuxf+GM+glPVbTFXreBa7
C+7mW3ik0c63qumBffr2NPR/WYpd7JW3FXCV3UZQSacQ2KDUXGjAMlm0tier
p3E33746Z+RbxRuZ883yl5bFIali8MrbCrjKbiNEx1q1x8dNIciaUlF/p+UC
KN8Jf3af8gXLt2yOGsnLt2K9rYCrbDaCcZ5xa5f9PoitEBpqNGfCr92R/ZOQ
b7Z45W0FXGWvEYwpBFUd/jlbY3hpin2oP64Ljl4IZf8U5JstXnlbAVfZaoQc
pxAsNUH05qk6SxdsPvhv9rqgYPn2wrnMJ0SF5VsR31bAVXYaQQV7amxPIVBq
1U2HcvmIL1i+Pf32zYwvYUylWLO28nhJTmL1ytsKuMpGI6jwQIP2yEPDd7Oc
Eq/Ug2Cv9dY6VTXtY7lM+3Q33+6/32xsYYc/44zdUP9eW0fvC88rbyvgquwb
wbwCvWEwm6G89PXli5GjvnWWLmjs/STH82ru5pv1KqdXRx9m+j7n7zSKyngx
V3F45W0FXJV1IyhjCkF1FlObYg+IBLprrTda9LWO5vEJ73K+zQw26tu5Z+BW
hh248PABo7VP+Euysb3ytgKuyrYRzEtvMk8hWDqndiHuRosVB/qCq28XaYfL
+bYYeNO4cWXL5fQ7cCpw3GjwLC7mKgqvvK2Aq7JthLlAu9bUmaYQrO4CX9tY
2N6sgyTczrfZD9v0l3j09Y/S7cCZZxtL+BJyr7ytgKuybAQ12b1F+8hOf3mm
UtFrHZYuqKw9NenIoRm38816iuHpt1LfeTzuNxrqer8o0Qb3ytsKuCq7RjCn
EKS9w2HCLzdVNvQEneoC1/NtcXbi58ZX1Ib+ULJfgIr9wd1LTUaBL42U5pfT
Re+8rYCrsmoE8w786U4XJvzKkq/Tb/zcbzbSb6f7+bZ079k6YxdlW/vVe5af
ctb+GfcbeTuOB0q3y73ytgKuyqYRzPvepJtCoKJXLL+ytLb+ny6OX/3DxPiV
LJY/Ff/3F5bNBt40fwFqzeP1J3//6R2tmaPh4Afd9U9Yjj51lnC6eedtBVyV
RSOYv4lZezb1BHUVGX7JMm3A1lL8+5Pr5qdHfvx4xg2uODAYKu1D6155WwFX
ZW4EFR09pj3A+puYq1h+f8rD+RYzfz9wet+mlFu7uf7kROmfN/TK2wq4KnMj
zN/o0H5OZUt3uiszLZdwejvfNF9Pjb3bdujZnXrQbfbt3tv0xqC/JK/GWs0r
byuAxcLnGwAUCvkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkG
QCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3
AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+
AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCry
DYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKR
bwAKRMUpwAsWOt9UGoV4/QIpkzITlEudyI1SDybatywHTmdAWr5pg/5ucPRs
68HdOzauvO6mbbsOvnZm+I9fLiwsiGiLMikziSKMXnhHbEREx//v1pXAEZZv
sWefi/jf2av3++pl0563xu98q7W/ixvirjIpM4nijF54w9IHetTfucNsBFn5
pubCo207UrS8uazfP3DrWw+3fpmUmah4oxcesPShv6ovJOWbUsGeOrM036HT
V0KR6PILzs5MXX5775Nm4ds6/A892vllUmaCoo5elDTtSEx44vhzlYmf8nLy
Tampvlqz6zv8D1a91Dcfd/9gg74l20/e9OK+TZmUaVECoxclShsbc5GbQ69U
VyX7FiMm31R0vM0Y/40D0yle5+v3mx/Taz88fM9znV8mZWpKZvSiNMXGRvCD
7oan4j/4fDt84vJNRS416S9R0TYRTf0y0+9+X9+Y2p7PPdb5ZVKmpnRGL0rT
zHv7Ez7yNj7/S39kokNcvkUuNOsvUdl+Le2LfHamRt+YF87d8Vbnl0mZmtIZ
vShNcSOkqubo0GRkLjYibojLN+XvND7lmy9G0r/IjRPGd7efXLznpQ4pkzJX
lM7oRWnSR8jm2mODN812kJdvoYF6vaLdfcEMrzEz8qL+4A3Hr3tpz6ZMylxR
OqMXpWnm/MuNvxi5eTdhFIjLNxVoNyo6cinDfs3iYuAXxleexvN3vNMiZVKm
rnRGLzxFXL6FRxrtVDQ9sE/fnob+L73TImVSZgbkG9ITl29fnTMav+KNzBVZ
D+xkOEpfUsqkzAzIN6QnOd8ahsJiG79MysyAfEN65Js3G79MysyAfEN65Js3
G79MysyAfEN6kvPthXOZzxR6tfHLpMwMyDekJznfnn77ZsaXCA816ttTedxD
LVImZWZAviE9cfl2//1mo6IOf8aprKH+vcbEsIHb3mmRMikzA/IN6YnLN2tF
r44+zPAayt+pPzjzVU4lpUzKTI98Q3ry8m1msFGvaM/ArQx7NuHhA/qDq074
vdQhZVJmeuQb0pOXb4uBN407OrZcTr9nowLHjZuJZXGVU0kpkzLTIt+QnsB8
m/2wTX+JR1//KN2ejZrqq7Yzy6KklEmZaZFvSE9gvlmPvT/9Vupbcsf9eEFd
7xde648yKTMd8g3pScy3xdmJnxvf3Rr6Q8l+Gi/2B3cvNRnf2l4a8eC3tjIp
Mw3yDemJzLelO9bWGXVta796z/Ibx9o/H0x2mb9RsuN4IM3tvUtYmZSZEvmG
9ITm2+Js4M2dxgutebz+5O8/vfNwuf2j4eAH3fVPGPMl1vk6vdv2ZVJmKuQb
0pOab4uL89MjP37c7P0US8WBwdCclzujTMpMjnxDenLzLWb+fuD0vk0pu35z
/cmJsICuL5MykyDfkJ7ofNN8PTX2btuhZ3fqCbDZt3tv0xuDfjmXKS0rkzLj
kG8oOYXONwAoFPINgFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvIN
gFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCpyDcAUpFv
AKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkG1CGVHkw8m1hYaHY2wIgOcfzzWh8
FhYWliIu5BsLC4vUhXxjYWGRuriRb2WCv0MAUpFvAKQi3wBIRb4BkIp8AyAV
+QZAKvINgFTkGwCpyDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCp
yDcAUpFvAKQi3wBIRb4BkIp8AyAV+QZAKvINgFTkGwCpyDcAUpFvAKRyMN9U
nJxXTLamSi/PDbfHrTJLp8Ikm5PziuVa5upHFL/WcuRYvikV6m9cfqrGgdu2
nkzdPLV1ecXKlivRxDW1sRAN+YdOND1f43tiZWvXf7d6z8vtPb//9N63CwsL
BRwubpSp1Xg3OHq29eDumr+t0mp89G+eqT/yf3ovB++tVFjIjqDMTGt6adCW
MWfyLbby9GDDylPZHCoq2LtbW7F1NGGkxP7XXHj85N6N+kYmWR6vf8e/0hq5
b3+2m+pCmVnU+J3vtZz/vIAdQZmZ1/XOoC1v+eeb0sZJlfH22Rsq6rY+xo5c
iljXi/37g2DvgcqUg8RcNv3o/99ye7C4U2asxsmu57KpcfvrHz0sQENQZlar
e2TQIr98W94N/6TPMk7sDhUVudCsrdgwFI4bKXOhAcs4qahuPv1BMKx9Vs7P
znz24S+bdj1qvmhtT9C9weJWmUqFh6298MTerliND1b+/9lbN/r+0VLjhob+
kKsNQZnZPoMXBi2W5J5vy8MkOPSKL/GDyc5QUZGLR7S1anqn4r7O6AdGlhbf
P4yE5pI85/TAfnO0HB6+58ZYcbFMFR1rrTB64UBf8EGSJ5z5sG27/pj1P7no
So2LlClr0EKXS75ph09vT3Q9vznZjredoRIda9HW2t0XtI6UB2MtVUZHDCYd
J8tmBhuN12363T1Hh4rbZapQb202zzY79toG/WEN/V863Q+UKWnQIp79fIsN
kkDvsbr4w6eVPt/WHIaKutahfbJvOTVpHSmhs3X6s9XFfUSuMvth25q1G/62
uqZ6967TAQeHivtlfnWu0diraZtYdRLOtFzjyiNfG33ocIhTppxBi0T28+1G
R8JnX1XN0V+HQu/ttz9UzJPs7dfiRorxiZ++I9zkfplfT42996sTrx2uqfnu
0QuRNM8VayhjMxIO+OSNMiUNWiTKM982Pt89vjwwZnIYKmqya4u2SvxJ9vBI
g/78LWPFGinul5klZR7Nji3tE87+fVCmpEGLRDnnW8VTjV1XzA/ZHIaKmuqr
1lY5Fjcg1ITxOd44rL+CSs1+0SVRZlZij58LDTQYvdl8Md0eUA4oU9KgRaIc
8q277uXTV0LR+MOn9oeKCg+tnGxqiv/aEurfqz9Vh18fD5Gpsd5jh6p9KweH
13+3+kevujrp3fUy066kmwuPtu0w3iNf3AEfR1CmpEGLRPbzLTnbQ8U8yZ7w
eOXv1D/HGwbvxIZCONU5r6XlO/tOXy/gTHDnykz1eK3fo5HbwfGhjvonzErT
npJzGGWm2GxPDtpyVqx8UylOssedPX9h4Nbng/HzMJMt247+LlSg0eJcmSlY
DuNYls31JyfCher6RcpMwaODtpwVK9+iY63a47ckfhmxDJVdtd9bGSeV1f+7
z/9FZOX7xezM1OVfNddssIyW9ut/LcRgca7MFIzDO9ble00FviibMpPy6KAt
Z0XKt9j+vHYZy9auxIYIW44za8uOlrFkswXmpwdfNEfL0x3+h+6PFefKTLHC
9c5Un/jrnzna/6dCXJu5SJkpttqjg7acFSffjJPsVR3+xO8jiVOVGvpDKZ/u
mxsd241HHjz/F9fHinNlJjf756npu/p5uW9WfeJvWL5isdRinDJLe9CWs6Lk
W9q5BFN9z1qGypbum2lb5v5vDhsPfuHcHbfHinNlZm1++tLPdlq+1HT98dvS
inHKLPFBW86KkW8q2FOjPTjZXALr1XlZHM+5/36zMa5edfiqntWcK9OO+am+
OvPv5MhvXb8mmzJX8fCgLWdFyDdlHKw4NHx39SPjhkoWU8E/7v7v+uNfGLhT
Oh2RoUx7Zq+2PVrAjqDMRF4etOWs8PmmIpeatEc2DCa9yHDidXOoxF3il+KV
LUNreWplPmVk4mCZ9lg/8fcM3KJMJ5TJoC1nBc83ZZxkr+5JPlVqqneXsVWv
f5TxO4rl1FWWk6/y4GCZKlH6V7YewS6dGKfMFaU8aMtZofNNRcfbMswlsN4n
Z8+/3sowVlTguDGX8sCIs3edWM3BMkOXThw/dqh66ZdWWi5n/CamAu1G42/p
vkmZTiiTQVvOCp1vc4F27Z1NM5fA+iUl0w1OzTvhZz+5NA8Oljn5//6H/lSP
xj7x03eEmh6sN/5OEu757wLKjOPxQVvOCpxvarJ7y/LDKjrT3Nfv/siPjS8p
tT2fp24KZb2pYGUWxz3y5WSZ1zuNY+nrW8fTzfNUKthjlOn4vdGSoMy47fL4
oC1nhc038yR7+nNM1vNo6/cP30l6yU7sDyx3hM446cgRTpY5bzmJtjb1jFaV
8EtPDpymzIgyrVvm9UFbzgqab+Y9TjPd4Gt+2nrCfVvHuHFZ4tJa2j/jfnKu
Mt/JV9lxtsyZ85afGtnWfvWepUbjyHx4rMVn/lXUpji+7SzKNHl/0JazQuab
edghm7kEfz7X+Jg5WtY/c7Tno1v3Hi4PmGg4+EF3/E11CnSQ1uEy469GXLOh
9vVz/lt/XWmK6O3g6DvN1VXWMgt07yDKNAgYtOWsgPlm3oC69mzqq/Osbp07
9KQ5HlItqX5yzg3Ol/nNVJ+191MvG18qqRujUWYKpTdoy1nh8k1FR49pD6hJ
/+tCce5/PGj98d/EpfTuGJZDmfP3A6f3bUrTDlU1R4cmI5TptDIZtOWscPk2
f6ND23W3P69p9s83Bt98ee+zO/Xu2Ozbfajl7FjQ/UPQCVwsc3Y68N4bzX+3
a7v+6V/xVE39/+ro/SBYyJbXUOYKEYO2nDmVbwBQasg3AFKRbwCkIt8ASEW+
AZCKfAMgFfkGQCryDYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCry
DYBU5BsAqcg3AFKRbwCkIt8ASEW+AZCKfAMgFfkGQCryDRBDIZ6RbwsLC8Xe
FgB5uQ8AQkUAQKg/AIBQ/wFiAmBf
     "], {{0, 437}, {416, 0}}, {0, 255},
     ColorFunction->RGBColor],
    BoxForm`ImageTag["Byte", ColorSpace -> "RGB", Interleaving -> True],
    Selectable->False],
   DefaultBaseStyle->"ImageGraphics",
   ImageSize->{233., Automatic},
   ImageSizeRaw->{416, 437},
   PlotRange->{{0, 416}, {0, 437}}]],ExpressionUUID->
  "81818d1f-f4bb-418c-bd0a-440a66d5754f"]
}], "Text",
 CellChangeTimes->{{3.8126994908663597`*^9, 3.8126995746972456`*^9}, {
  3.812699612871337*^9, 3.8126996503556414`*^9}, {3.812699750815791*^9, 
  3.812699759035096*^9}, {3.813035073034068*^9, 3.8130350839243917`*^9}, {
  3.813035258654602*^9, 3.813035278845145*^9}, {3.813035373601288*^9, 
  3.8130353868555107`*^9}},
 FontSize->18,
 Background->GrayLevel[
  0.85],ExpressionUUID->"f084202f-33da-4926-824d-86be44fbdf50"],

Cell[BoxData[
 RowBox[{
  RowBox[{"RK4", "[", 
   RowBox[{"f_", ",", "n_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", 
  RowBox[{"(", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"h", "=", 
     FractionBox[
      RowBox[{"T", "-", "t0"}], "n"]}], ";", "\[IndentingNewLine]", 
    RowBox[{"solution", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"n", "+", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"solution", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", 
     RowBox[{"{", 
      RowBox[{"t0", ",", "u0"}], "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", 
       RowBox[{"Length", "[", "solution", "]"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"k1", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "]"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"k2", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", 
            RowBox[{"0.5", " ", "h"}]}], ",", " ", 
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
            RowBox[{"0.5", "k1"}]}]}], "]"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"k3", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", 
            RowBox[{"0.5", " ", "h"}]}], ",", " ", 
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
            RowBox[{"0.5", "k2"}]}]}], "]"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"k4", "=", 
        RowBox[{"h", " ", 
         RowBox[{"f", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", "h"}], ",", 
           RowBox[{
            RowBox[{"solution", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", "k3"}]}], "]"}]}]}],
        ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"solution", "[", 
         RowBox[{"[", 
          RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "+", "h"}], ",", 
          RowBox[{
           RowBox[{"solution", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "+", 
           RowBox[{
            FractionBox["1", "6"], "k1"}], "+", 
           RowBox[{
            FractionBox["1", "3"], "k2"}], "+", 
           RowBox[{
            FractionBox["1", "3"], "k3"}], "+", 
           RowBox[{
            FractionBox["1", "6"], "k4"}]}]}], "}"}]}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "solution"}], 
   "\[IndentingNewLine]", ")"}]}]], "Input",
 CellChangeTimes->{{3.812699653823759*^9, 3.8127000856273637`*^9}, {
  3.8127002519425097`*^9, 3.81270025218266*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"274b5909-0cb0-4260-b141-d58a5d1a171f"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Problem3.", "Section",
 CellChangeTimes->{{3.8127010649899545`*^9, 3.812701068818507*^9}, {
  3.8127064205417128`*^9, 
  3.812706420648537*^9}},ExpressionUUID->"995c2255-176b-4e1a-aab4-\
2ba8523595de"],

Cell[TextData[{
 "\:0414\:0430\:0434\:0435\:043d\:043e \:0435 \:043b\:043e\:0433\:0438\:0441\
\:0442\:0438\:0447\:043d\:043e\:0442\:043e \:0443\:0440\:0430\:0432\:043d\
\:0435\:043d\:0438\:0435\n\n\t",
 Cell[BoxData[
  FormBox[
   FractionBox["du", "dt"], TraditionalForm]],
  FontSize->24,ExpressionUUID->"f1607318-304a-4f2c-851d-b1a48c5880b0"],
 StyleBox["=10u(1-u), 0<t<=1,\n\tu(0)=0.01.",
  FontSize->24],
 "\n\t\n\:0441 \:0442\:043e\:0447\:043d\:043e \:0440\:0435\:0448\:0435\:043d\
\:0438\:0435\n\t\n\t",
 StyleBox["u(t) = ",
  FontSize->24],
 Cell[BoxData[
  FormBox[
   FractionBox["0.01", 
    RowBox[{"0.01", " ", "+", " ", 
     RowBox[{"0.99", 
      SuperscriptBox["e", 
       RowBox[{
        RowBox[{"-", "10"}], "t"}]]}]}]], TraditionalForm]],
  FontSize->24,ExpressionUUID->"ec0f1487-0f55-4588-a77d-00b5dfa38b42"],
 "\n\n\t\:0430) \:0414\:0430 \:0441\:0435 \:0440\:0435\:0448\:0438 \:0443\
\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0435\:0442\:043e \:0441 3-\:0435\
\:0442\:0430\:043f\:043d\:0438\:044f \:043c\:0435\:0442\:043e\:0434 \:043d\
\:0430 \:0420\:0443\:043d\:0433\:0435-\:041a\:0443\:0442\:0430 \:043e\:0442 \
\:0417\:0430\:0434.1 \:043f\:0440\:0438 n=10,100,1000;\n\tb) \:0414\:0430 \
\:0441\:0435 \:0440\:0435\:0448\:0438 \:0443\:0440\:0430\:0432\:043d\:0435\
\:043d\:0438\:0435\:0442\:043e \:0441 4-\:0435\:0442\:0430\:043f\:043d\:0438\
\:044f \:043c\:0435\:0442\:043e\:0434 \:043d\:0430 \:0420\:0443\:043d\:0433\
\:0435-\:041a\:0443\:0442\:0430 \:043e\:0442 \:0417\:0430\:0434.2 \
\:043f\:0440\:0438 n=10,100,1000;\n\t\n\:0414\:0430 \:0441\:0435 \:0438\:0437\
\:043e\:0431\:0440\:0430\:0437\:0430\:044f\:0442 \:0432 \:0435\:0434\:043d\
\:0430 \:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442\:043d\:0430 \
\:0441\:0438\:0441\:0442\:0435\:043c\:0430 \:0433\:0440\:0430\:0444\:0438\
\:043a\:0438\:0442\:0435 \:043d\:0430 \:0442\:043e\:0447\:043d\:043e\:0442\
\:043e \:0438 \:043f\:0440\:0438\:0431\:043b\:0438\:0436\:0435\:043d\:043e\
\:0442\:043e \:0440\:0435\:0448\:0435\:043d\:0438\:0435 \:0432\:044a\:0432 \
\:0432\:0441\:0438\:0447\:043a\:0438 \:0441\:043b\:0443\:0447\:0430\:0438.\n\
\:0414\:0430 \:0441\:0435 \:043d\:0430\:0447\:0435\:0442\:0430\:044f\:0442 \
\:0433\:0440\:0430\:0444\:0438\:043a\:0438 \:043d\:0430 \:0430\:0431\:0441\
\:043e\:043b\:044e\:0442\:043d\:0438\:0442\:0435 \:0433\:0440\:0435\:0448\
\:043a\:0438 \:0432\:044a\:0432 \:0432\:0441\:0438\:0447\:043a\:0438 \:0441\
\:043b\:0443\:0447\:0430\:0438."
}], "Text",
 CellChangeTimes->{
  3.812700893965129*^9, {3.812700946173775*^9, 3.8127009624543037`*^9}, {
   3.81270114351902*^9, 3.8127012259974203`*^9}, {3.812706427722856*^9, 
   3.8127065532706647`*^9}},
 FontSize->18,
 Background->GrayLevel[
  0.85],ExpressionUUID->"eb5c5bca-f5fa-4910-a0f0-d971673b8378"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"F", "[", 
    RowBox[{"t_", ",", "u_"}], "]"}], ":=", 
   RowBox[{"10", "u", 
    RowBox[{"(", 
     RowBox[{"1", "-", "u"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.8127001366949544`*^9, 3.8127001524561615`*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"f3597dcb-7738-47d2-82ab-7769297cd647"],

Cell[BoxData[
 RowBox[{
  RowBox[{"exactSol", "[", "t_", "]"}], ":=", 
  FormBox[
   FractionBox["0.01", 
    RowBox[{"0.01", " ", "+", " ", 
     RowBox[{"0.99", 
      SuperscriptBox["E", 
       RowBox[{
        RowBox[{"-", "10"}], "t"}]]}]}]],
   TraditionalForm]}]], "Input",
 CellChangeTimes->{{3.8127009690982356`*^9, 3.812701003186534*^9}},
 CellLabel->"In[55]:=",ExpressionUUID->"8464c760-38b2-4e23-8e97-ed2c0348e6ae"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"solRK310", "=", 
   RowBox[{"RK3", "[", 
    RowBox[{"f", ",", "10", ",", "t0", ",", "T", ",", "u0"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"solRK3100", "=", 
   RowBox[{"RK3", "[", 
    RowBox[{"f", ",", "100", ",", "t0", ",", "T", ",", "u0"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"solRK31000", "=", 
   RowBox[{"RK3", "[", 
    RowBox[{"f", ",", "1000", ",", "t0", ",", "T", ",", "u0"}], "]"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{"exactSolution", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListPlot", "[", "solRK310", "]"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"p2", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{"exactSolution", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListPlot", "[", "solRK3100", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p3", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{"exactSolution", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListPlot", "[", "solRK31000", "]"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"GraphicsRow", "[", 
  RowBox[{"{", 
   RowBox[{"p1", ",", "p2", ",", "p3"}], "}"}], "]"}]}], "Input",
 CellChangeTimes->{{3.8130528010859394`*^9, 3.8130528551848*^9}, {
  3.81305292752722*^9, 3.8130529520137634`*^9}, {3.8130541535302796`*^9, 
  3.8130541658816657`*^9}},
 CellLabel->"In[64]:=",ExpressionUUID->"3a6a925f-56b8-4a77-a813-c11c80d05961"],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[{{{{}, {}, 
         TagBox[
          {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], 
           LineBox[CompressedData["
1:eJwV1Hk8VPsbB3DbjVtRQjWWK+m2iEkq7X3qqihKli4VspYkrZarhBqSLUlI
1iyFjLIkS2SdKLImZEkZM2NmvslFxeV3fn+c13m9X+ec53yf7/M8XzX786ZO
YiIiIi3U9f/7QSdOy2vu0V2ZLFrYZ28GJPQUL/TTNqLGYt9oazADG5Ji4/Np
e/GjyTqxOY4Bw8bOvASaORzlL/z4VsKAhePkq0iaI0Q+GT6vm2Kgd90tvTDa
FfjJJ7g3+QYgWcE4MJjGwKcY0dzIu4Hwl8pfdZsWhcMOlUeX9QVBQ/9uUN6F
KMhe7p+VHQtCa8B5bg8rCkZ9TQoKUrexQlwzm+5xH5Y7sl4Z69xG3X9p9PbW
aJRvcQuJvXUb88eiN6mGPkBKNMlnbApGbK+3XuFMIhreSKyiJ4XA8lhK7rfR
NChoN7xYnhMO0mM+95FiOpzrCsW1a8MRaCPlZKaXjiK9C4PGveHIdziv+CIq
HbX3voaypO9Axm1nwNXNGUiNG1GOPX8HtTe6js/xeYzn3nv2OClFQOep7Bxl
qSy03fn1gM+NwLxZfyv9P3KxOcX497H0SMSesbDR3J+LDaIpXjsKIrGiXdNW
1i0XXtMF6RFVkdiV2enQ8yoX+na6lsf6I1GcGV1YRH8G1W9bf7rS7kFgxi8z
bHoG9SZ52bQ793D0Sczby9J5eL3Tr+s2Iwp/mgg51SEF2F0peUkyMhqKvvrb
TqcU4GVEh2pRWjQW5iSHzC0qwPZl4LsXReOXlBndZLAA7kNBC9R6o9H8uuhy
35ZCEHftN/Q1MfDW9p/5MVQIDQNlG52aGDQtkJen7y5Ct0l1Nn3OA3g2bUXs
v8WIO9GaV1r2EDL0wqF6qRJM+aqslW59iLQw7dAp5RJ8v3g2wm34IVqMVnXZ
7CsBg9592G5RPNa+lb/8Z3QJ9ihIfZA8E49+FsnI21yK3xI1pu1pCdCvTJdp
9C6DVMJFxqPARCwpkO0TFa1AjtauNd8bkmEVOmX+dXEFdl2/v6n7YzIeOQ69
rdOsgInfFbtGdjI0FYqLQywrMBhdv6VDLAV/edjeV3hWgXU8buPZ7SnQjbfI
iFr3GmZmDl9KmSmws8mwftL3GoExxUfUEx/hxaDeu/c7qmAbxMmdTknDBwP/
0TD9KvwR8OxmSl4aJpjli41Mq5BjY7rBrDoNm69utXtzugodYtsP9n5Nw0u5
dROvI6uwN9/UPGx1Oor3Ki7L41ThS798p2NeOsoyvl2OiqpGS/wnP53mDFS5
JCgd49cA6RdCT+pmIp979qzNRA0qOC2GGgcykea8rdRBpBa7VsqbilllIuBU
53E3+VqcCDtj3u2fCQN72bibO2pxzt+WtrYpE++PBSxlhtbiabnlj7/PZOGT
wTkFca06OPUstpPLysb4yp0Lcs6xsLnnfoeEBRO7Bv9iBXuxELt0c5HBeSYC
Ewx8nW+y4D2woC7xFhOL5c2J+gMWjvwnGuf/kgld0bPv42pYeJcoba+uRPVx
T+ydIKU3OPmksfjzl1xM3/l3gWP9GzTFXTNqqnwGsZ/ZC5X+bMD8FPcNlr55
0J21fWHGfIf5wmSPCZ0XeMfY2FSn8R6el728DD1LoBox8jNT2IzjqdP7HrMr
MDLZ39eh3Yr1lbrXY4erkMYrbZAJb0NSsv/6w3F1kFF5TL5+bUdOad5Ln+oG
hJ4pPJW84gM0taxdB5a9x8Tu7K/Ovp0I8j1yrpiKc0WnpmBx9UewlgcZ6va1
o59ZE9Ki1g3XKNazyKlOTLm0P1zj3YMjtJmVSzf1YO0eyZlTpZ9w9JiKvpFx
H/QfqZ1+vqQPlVp/H43rHoCnysnln636UW+RoH3VdRCa/KPzOg8N4Fe65akZ
3hfILfgmVfNrAGROqvpL9yFIVDzMVY//jIPMo943Gtk4u7G063fDQXB9nBz5
yznIPmCfXD4xiHcs7j6pI1yMjyguZEd/wdClvCt1B3goNNMXj9r7FaNDkoes
N4/A1yAgqYD9FQxtydTgyRGo9R78mBowBEtZGuNAHB/bNjOiaDpsvFjbo8/W
FsDqj/c/RdvYEIavVil/LIAM7VJhaQcb6i5iO2IyBaiUU7jo/pGNhuTd4uez
BVj5+wkOp5eNJUFLwhVzBSBj7A/vOWxUPLbPtnshwI36mfyEGTbirJwLn9UI
8OQK3W2rxjDuWjkxLAYEGG8IG7zgP4y5p3g7x+SE+DSvw3n65jA0uv+9a68g
RJWRMrkVOAzpBp2h5sVCRDRlTSWGDOOzXG5IFk0IjVaWfOP9YRxyP801URXC
pktUf00WFa9p6bmLGkLUDbtnf24dRvAWy4JJCBErYXPFdAUHe+k2fnuchZhx
WLKpZSUHH8ju9efPCOFY3TxuvIYDmc/BMvEuQmjf0PM8ROdgz6ztkTFXId7M
rvY22MJBSob9bMxFISZ/jfnuMuIgfu/payXeQvw9ejtE4woHHm93XjkdKoRc
X+EjsRoODq2zZV7LESJ8vo9PUR0Hsd9WvgpnCjF3+15L13oOnD3CJ5NyhRCL
bZX+0MRB2+Go0YrnQnw3IV6ZXRxM6BTPnSgUorl2tfERwsHCDBPT7eVChOY8
nEpQ5OKMNVOL3iTEbz7+5lsvcpH/rZnVNiJEwSj/Z8VlLhjeMmoZfCq/U5ZJ
+z24+NKb2u8pEKLGeB3P7CoX36a9w2lECMbyXl+3AC72/W1cfvS7EBJvtmSn
PuAiTcuj8vkPan1y30Rlqrj4Pvlz2ZQEgUimTe7gIh72tqkMcpQI6vidbkkK
PCgeH/0rQ5kgVNuEbrWUh/QS0XkOKgRLXurldKjwcHXycWv3HwR01ursN6t5
WOTyyKhajcBq6HsGcxcPAw40Y69VBMVqtxK9XXjgWa5IcNchuBSXG7aoiger
J4eH2/cTFE3H3aqq4eGoivMTV32CaZtA/0ssHuSuzRBxA4Jb6tbure94aFsv
L7/+AEH807nWkZ08LLYeLw4wJKgtd9JcJOChofbIu6VHCJZ+UXoru3QE6rXt
1ziWBK/WBknJuo1AfHNh6+4zBLrtSxojL1BzJG1yr4Dys2uP78pfHoHExvzp
VS4EjxrrFJd6jeDkPFaD9FmCoPO/aaneGIF2/7wTHa4Epvk3TLSiRyDd+s8N
8wsEw9t84g6UU3OqMx0j70kge/CS5g0ZPpQ91uY/uUmQXLlg3jZZPuw0F4Ut
ZBBob83hjsrx4SQaIetF+fBqToYdjY/7o7ZK+wMIgufYqO1ZwceYYdnWvkAC
ieoDi0W28aE72dkzc5tgYruayHUnPhR23/QVjyD4pNXc9k8ZH67vRCZXPSRI
GpH1Ea/gY+H15x0elO0zzVaGVfLRzfTwrKXMXdHplVzHR2L/uIR9PBVPqU+F
1czHiQfVRrEJ1Prn8k/LD/ERs0mu6GcSwX625DRTWgC96Odr76ZR+5WIFV9t
BMjKm9OllEMw5j0maWcnwFCzT8gJyroWj0d6HQQIeOov/ZBy6YIF+Z3OAsj2
VfjRmFR/+fXteXtJgDIrsdXyuQQ99j4n8wIFcP6cEyvynGDOqtI4P6YA6OCP
vyggsM7VlVWZEeAQaeEllxK8OP5DplmEmmP3p+LdlGXmlMy/KS6EyPVdgkVl
BK+tdkpxpYTw5i/fwKC8fO7e2QLqXDKIUQp3fEXAdjQRGK0RQrflU4pyBYGb
omv9NXMhcrxV9c5WEfgFpPh/yhLieF6+evAbAmlPyYN21Nw3cze4FVKOcz63
iE3NuabYipgByvmGW1NJgRA3XUoYuvVUfNmWGrEKIa6KB7cMUDZMFJFa00ad
i/W1+uveEiwuOhnuPiWE4no3s4JGghyOykMZI4LUzOxFjq1UvyUE2PYdJkgJ
n9cTQtnRVPAn04Qg4hrXL5+yclnZs0MWBDVl953E2qh+uXOiLtSeYNyZfT+R
srPug9G5/xB0jL3Lbm6n8mfIG0hmUHW9yb6k3Ekws+WadOcT6jvr28t3U+4S
fGnNyKbqWE7Pd6AcYZFvvZ+qw+G7I55ZlGc1TC8HUPsYw5ON1P1I9VvLnURx
6j+qfnzVA10ELwMnHdo6CJ5vkeC7UI7afnJNKvXe+MnyuFDKhun0gr/6CNSn
rye/p1zi1VjvzyU4/Z4Ym3cTRGttumPCJ/jjR9cGd8oXB+PN1QjBl7P0qfuU
Vxu59r/+l8r3ntfGTsoSou1pEZMEDL2kjEnKA4XbXWx/Uc83RU4v6SEoc0ld
p/0fde4UnNHdQjlGdd747CzB/Imd5paU/wcrpKuc
            "]]},
          Annotation[#, "Charting`Private`Tag$26694#1"]& ]}, {}}, {{}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.012833333333333334`], AbsoluteThickness[1.6], 
          PointBox[{{0., 0.01}, {0.1, 0.0262729185613093}, {0.2, 
           0.06739169520001892}, {0.3, 0.16292082784235823`}, {0.4, 
           0.3458319391572666}, {0.5, 0.5898375307923677}, {0.6, 
           0.7956474125088959}, {0.7, 0.9150849161063017}, {0.8, 
           0.9687114610576664}, {0.9, 0.9891582429569918}, {1., 
           0.9963362178932972}}]}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.012833333333333334`], AbsoluteThickness[1.6]}, {}}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.012833333333333334`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0, 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{11.052083333333343`, 1.552083333333428}, {
       6.264878372495332, 0.5000000000000142}},
      Method->{
       "DefaultBoundaryStyle" -> Automatic, 
        "DefaultGraphicsInteraction" -> {
         "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
          "Effects" -> {
           "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
            "Droplines" -> {
             "freeformCursorMode" -> True, 
              "placement" -> {"x" -> "All", "y" -> "None"}}}}, 
        "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> 
        None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0, 1}, {0., 0.9955255170204409}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {193.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{{{}, {}, 
         TagBox[
          {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], 
           LineBox[CompressedData["
1:eJwV1Hk8VPsbB3DbjVtRQjWWK+m2iEkq7X3qqihKli4VspYkrZarhBqSLUlI
1iyFjLIkS2SdKLImZEkZM2NmvslFxeV3fn+c13m9X+ec53yf7/M8XzX786ZO
YiIiIi3U9f/7QSdOy2vu0V2ZLFrYZ28GJPQUL/TTNqLGYt9oazADG5Ji4/Np
e/GjyTqxOY4Bw8bOvASaORzlL/z4VsKAhePkq0iaI0Q+GT6vm2Kgd90tvTDa
FfjJJ7g3+QYgWcE4MJjGwKcY0dzIu4Hwl8pfdZsWhcMOlUeX9QVBQ/9uUN6F
KMhe7p+VHQtCa8B5bg8rCkZ9TQoKUrexQlwzm+5xH5Y7sl4Z69xG3X9p9PbW
aJRvcQuJvXUb88eiN6mGPkBKNMlnbApGbK+3XuFMIhreSKyiJ4XA8lhK7rfR
NChoN7xYnhMO0mM+95FiOpzrCsW1a8MRaCPlZKaXjiK9C4PGveHIdziv+CIq
HbX3voaypO9Axm1nwNXNGUiNG1GOPX8HtTe6js/xeYzn3nv2OClFQOep7Bxl
qSy03fn1gM+NwLxZfyv9P3KxOcX497H0SMSesbDR3J+LDaIpXjsKIrGiXdNW
1i0XXtMF6RFVkdiV2enQ8yoX+na6lsf6I1GcGV1YRH8G1W9bf7rS7kFgxi8z
bHoG9SZ52bQ793D0Sczby9J5eL3Tr+s2Iwp/mgg51SEF2F0peUkyMhqKvvrb
TqcU4GVEh2pRWjQW5iSHzC0qwPZl4LsXReOXlBndZLAA7kNBC9R6o9H8uuhy
35ZCEHftN/Q1MfDW9p/5MVQIDQNlG52aGDQtkJen7y5Ct0l1Nn3OA3g2bUXs
v8WIO9GaV1r2EDL0wqF6qRJM+aqslW59iLQw7dAp5RJ8v3g2wm34IVqMVnXZ
7CsBg9592G5RPNa+lb/8Z3QJ9ihIfZA8E49+FsnI21yK3xI1pu1pCdCvTJdp
9C6DVMJFxqPARCwpkO0TFa1AjtauNd8bkmEVOmX+dXEFdl2/v6n7YzIeOQ69
rdOsgInfFbtGdjI0FYqLQywrMBhdv6VDLAV/edjeV3hWgXU8buPZ7SnQjbfI
iFr3GmZmDl9KmSmws8mwftL3GoExxUfUEx/hxaDeu/c7qmAbxMmdTknDBwP/
0TD9KvwR8OxmSl4aJpjli41Mq5BjY7rBrDoNm69utXtzugodYtsP9n5Nw0u5
dROvI6uwN9/UPGx1Oor3Ki7L41ThS798p2NeOsoyvl2OiqpGS/wnP53mDFS5
JCgd49cA6RdCT+pmIp979qzNRA0qOC2GGgcykea8rdRBpBa7VsqbilllIuBU
53E3+VqcCDtj3u2fCQN72bibO2pxzt+WtrYpE++PBSxlhtbiabnlj7/PZOGT
wTkFca06OPUstpPLysb4yp0Lcs6xsLnnfoeEBRO7Bv9iBXuxELt0c5HBeSYC
Ewx8nW+y4D2woC7xFhOL5c2J+gMWjvwnGuf/kgld0bPv42pYeJcoba+uRPVx
T+ydIKU3OPmksfjzl1xM3/l3gWP9GzTFXTNqqnwGsZ/ZC5X+bMD8FPcNlr55
0J21fWHGfIf5wmSPCZ0XeMfY2FSn8R6el728DD1LoBox8jNT2IzjqdP7HrMr
MDLZ39eh3Yr1lbrXY4erkMYrbZAJb0NSsv/6w3F1kFF5TL5+bUdOad5Ln+oG
hJ4pPJW84gM0taxdB5a9x8Tu7K/Ovp0I8j1yrpiKc0WnpmBx9UewlgcZ6va1
o59ZE9Ki1g3XKNazyKlOTLm0P1zj3YMjtJmVSzf1YO0eyZlTpZ9w9JiKvpFx
H/QfqZ1+vqQPlVp/H43rHoCnysnln636UW+RoH3VdRCa/KPzOg8N4Fe65akZ
3hfILfgmVfNrAGROqvpL9yFIVDzMVY//jIPMo943Gtk4u7G063fDQXB9nBz5
yznIPmCfXD4xiHcs7j6pI1yMjyguZEd/wdClvCt1B3goNNMXj9r7FaNDkoes
N4/A1yAgqYD9FQxtydTgyRGo9R78mBowBEtZGuNAHB/bNjOiaDpsvFjbo8/W
FsDqj/c/RdvYEIavVil/LIAM7VJhaQcb6i5iO2IyBaiUU7jo/pGNhuTd4uez
BVj5+wkOp5eNJUFLwhVzBSBj7A/vOWxUPLbPtnshwI36mfyEGTbirJwLn9UI
8OQK3W2rxjDuWjkxLAYEGG8IG7zgP4y5p3g7x+SE+DSvw3n65jA0uv+9a68g
RJWRMrkVOAzpBp2h5sVCRDRlTSWGDOOzXG5IFk0IjVaWfOP9YRxyP801URXC
pktUf00WFa9p6bmLGkLUDbtnf24dRvAWy4JJCBErYXPFdAUHe+k2fnuchZhx
WLKpZSUHH8ju9efPCOFY3TxuvIYDmc/BMvEuQmjf0PM8ROdgz6ztkTFXId7M
rvY22MJBSob9bMxFISZ/jfnuMuIgfu/payXeQvw9ejtE4woHHm93XjkdKoRc
X+EjsRoODq2zZV7LESJ8vo9PUR0Hsd9WvgpnCjF3+15L13oOnD3CJ5NyhRCL
bZX+0MRB2+Go0YrnQnw3IV6ZXRxM6BTPnSgUorl2tfERwsHCDBPT7eVChOY8
nEpQ5OKMNVOL3iTEbz7+5lsvcpH/rZnVNiJEwSj/Z8VlLhjeMmoZfCq/U5ZJ
+z24+NKb2u8pEKLGeB3P7CoX36a9w2lECMbyXl+3AC72/W1cfvS7EBJvtmSn
PuAiTcuj8vkPan1y30Rlqrj4Pvlz2ZQEgUimTe7gIh72tqkMcpQI6vidbkkK
PCgeH/0rQ5kgVNuEbrWUh/QS0XkOKgRLXurldKjwcHXycWv3HwR01ursN6t5
WOTyyKhajcBq6HsGcxcPAw40Y69VBMVqtxK9XXjgWa5IcNchuBSXG7aoiger
J4eH2/cTFE3H3aqq4eGoivMTV32CaZtA/0ssHuSuzRBxA4Jb6tbure94aFsv
L7/+AEH807nWkZ08LLYeLw4wJKgtd9JcJOChofbIu6VHCJZ+UXoru3QE6rXt
1ziWBK/WBknJuo1AfHNh6+4zBLrtSxojL1BzJG1yr4Dys2uP78pfHoHExvzp
VS4EjxrrFJd6jeDkPFaD9FmCoPO/aaneGIF2/7wTHa4Epvk3TLSiRyDd+s8N
8wsEw9t84g6UU3OqMx0j70kge/CS5g0ZPpQ91uY/uUmQXLlg3jZZPuw0F4Ut
ZBBob83hjsrx4SQaIetF+fBqToYdjY/7o7ZK+wMIgufYqO1ZwceYYdnWvkAC
ieoDi0W28aE72dkzc5tgYruayHUnPhR23/QVjyD4pNXc9k8ZH67vRCZXPSRI
GpH1Ea/gY+H15x0elO0zzVaGVfLRzfTwrKXMXdHplVzHR2L/uIR9PBVPqU+F
1czHiQfVRrEJ1Prn8k/LD/ERs0mu6GcSwX625DRTWgC96Odr76ZR+5WIFV9t
BMjKm9OllEMw5j0maWcnwFCzT8gJyroWj0d6HQQIeOov/ZBy6YIF+Z3OAsj2
VfjRmFR/+fXteXtJgDIrsdXyuQQ99j4n8wIFcP6cEyvynGDOqtI4P6YA6OCP
vyggsM7VlVWZEeAQaeEllxK8OP5DplmEmmP3p+LdlGXmlMy/KS6EyPVdgkVl
BK+tdkpxpYTw5i/fwKC8fO7e2QLqXDKIUQp3fEXAdjQRGK0RQrflU4pyBYGb
omv9NXMhcrxV9c5WEfgFpPh/yhLieF6+evAbAmlPyYN21Nw3cze4FVKOcz63
iE3NuabYipgByvmGW1NJgRA3XUoYuvVUfNmWGrEKIa6KB7cMUDZMFJFa00ad
i/W1+uveEiwuOhnuPiWE4no3s4JGghyOykMZI4LUzOxFjq1UvyUE2PYdJkgJ
n9cTQtnRVPAn04Qg4hrXL5+yclnZs0MWBDVl953E2qh+uXOiLtSeYNyZfT+R
srPug9G5/xB0jL3Lbm6n8mfIG0hmUHW9yb6k3Ekws+WadOcT6jvr28t3U+4S
fGnNyKbqWE7Pd6AcYZFvvZ+qw+G7I55ZlGc1TC8HUPsYw5ON1P1I9VvLnURx
6j+qfnzVA10ELwMnHdo6CJ5vkeC7UI7afnJNKvXe+MnyuFDKhun0gr/6CNSn
rye/p1zi1VjvzyU4/Z4Ym3cTRGttumPCJ/jjR9cGd8oXB+PN1QjBl7P0qfuU
Vxu59r/+l8r3ntfGTsoSou1pEZMEDL2kjEnKA4XbXWx/Uc83RU4v6SEoc0ld
p/0fde4UnNHdQjlGdd747CzB/Imd5paU/wcrpKuc
            "]]},
          Annotation[#, "Charting`Private`Tag$26778#1"]& ]}, {}}, {{}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.009166666666666668], AbsoluteThickness[1.6], 
          PointBox[CompressedData["
1:eJw1kg1Q03UYxwdKAYeVppRgeiqmxjoSDixwfDPD7HgpTeRV4abYiaJgIrkl
pIyhEegRqKQslJdMgUBA5U5gbCii4AkDhrzMjcHY+9/UqXiZsfn43P3vf/+X
3/N8v9/Ps5C7Z0OCPYvFSpi6rPfXxZ9d/ZWKJ8DrOyvlUqXnGSE9F8B3pdqR
ZzmCq/OzDOuvilBZ7pa2/dcc+l6GfWsz7Gd9n4ffi6z1F6Krolveu3Gc/q/G
3Gf8xOxD+TCut76phaOIE9m7o4DO14NTsNAvPPMEnOWlyRznK4iYFjWjtecU
9WuE/YFzn/luPo3lsccklmXXMBIm8A5YJqL+zTjV/cZLzZkSlE2dlpeKsXA5
P9cr5BzNa0WeS92GpCVl8LGVFKsn//ELZ1fQ/DYE3hH2D207D6ll2dSE69ii
4v3oK7tAetqRJeLqevlV2GAbeBMmfHGwfkUN6evAx1v3eO1qr4WKZ214C4Ic
c2fxoTrS24lE7+YIPbcBKVZ5yV04Hj1tkf3uK6T/DhyWNJuKyxvxisZdJEy/
kLbf5Rr5uYuLo8Nh/uXNWGcb0I1BydqHnzeKyV83/J0bNPEBElinxR7rwRBn
TckHg1LyK0OXoJqfeuI6bHbkMgTv5UQEpLST/178VBsic9nZgQVWOfP7YJzY
sXGd4Dbl0QdRm6/btKYubLdVP3ZdLZ7IcL5L+chh92iJvSyoG9Zus6vlYF3X
HX0xpeNVXgOoc8jbH6aV4Ym1nWUA3wZGnb0f2kf53UOcd0pbYUs/bPFwBuH7
zlaXRasGKM9BrEx3C3eQ3INNXtYQPL0+cS38ZojyHYb7gTS9nXoYXZ3WGoa3
/42sOJ6C8h4B+5aHvt9JCZvd2QpEsoNTU/qVlL8CHnnZfz8/q4KtXex9XAtW
hP2wc5R4KBF/8q1tmhVq8qPEvJiokaTHauKjxBzXP5P9Lo2RPyW6/pDfPpw4
TrxU4M5/aqmapyG/Kpx/kD3AvqkhfipMV1oe5e2aIP+jeBhXH5LupCWeo+g7
OVjRLNJSHqNY7KX3ZLN1xFeNE09bihprdZSPGn6TpW67vfTEW43NfO44t1xP
eY0hN7SpVfyugfiPYVblwNEmnoHyG0NwgF1T0aCB9mEcw0yNj8jHSHmO48Nk
LstNaKT9GEfNuQ5BfI+R8tXgS0l+S/v7JtoXDXIj9bKKKBPlrYHTL7c7Nhaa
aH8m4F4aJ1jaaaL8JyC0z9Hv/c9E+zSBNJ4D5wjbTDy0WN/2vKBhk5l4aLGt
01UYdNBMPLSI/+ilIr3ETDy0kCzOyCwTm4mHDoltjrGWETPx0GHdp0WbKp+Z
iYcOjeE/p7JmMsRDD6PLC63HUoZ46LFGMWPF1wEM8dCjodY/qCiUIR4GPG3g
iTlbGOJhQAy/8u3YJIZ4GCB1jz/rymeIhxGrfR7Nyc1miIcRUWt89wznM8TD
iMLfJgMXFDPEw4Tv8hf1JlUwxMOE+n97EhTVDPEwYaMgZ1XWZYZ4mKGaOfPe
viaGeJhxWek4t1nKEA8zYthBUckdDPFgcPF0e2HRHYZ4MEjPFGavlTHEg4Hw
yZviw3KGeDzAY6lOEjnE4H8rniDe
           "]]}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.009166666666666668], AbsoluteThickness[1.6]}, {}}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.009166666666666668], AbsoluteThickness[1.6]}, {}}}, {{}, {}}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0, 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{11.052083333333343`, 1.552083333333428}, {
       6.264878372495332, 0.5000000000000142}},
      Method->{
       "DefaultBoundaryStyle" -> Automatic, 
        "DefaultGraphicsInteraction" -> {
         "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
          "Effects" -> {
           "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
            "Droplines" -> {
             "freeformCursorMode" -> True, 
              "placement" -> {"x" -> "All", "y" -> "None"}}}}, 
        "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> 
        None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0, 1}, {0., 0.9955255170204409}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {580.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{{{}, {}, 
         TagBox[
          {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], 
           LineBox[CompressedData["
1:eJwV1Hk8VPsbB3DbjVtRQjWWK+m2iEkq7X3qqihKli4VspYkrZarhBqSLUlI
1iyFjLIkS2SdKLImZEkZM2NmvslFxeV3fn+c13m9X+ec53yf7/M8XzX786ZO
YiIiIi3U9f/7QSdOy2vu0V2ZLFrYZ28GJPQUL/TTNqLGYt9oazADG5Ji4/Np
e/GjyTqxOY4Bw8bOvASaORzlL/z4VsKAhePkq0iaI0Q+GT6vm2Kgd90tvTDa
FfjJJ7g3+QYgWcE4MJjGwKcY0dzIu4Hwl8pfdZsWhcMOlUeX9QVBQ/9uUN6F
KMhe7p+VHQtCa8B5bg8rCkZ9TQoKUrexQlwzm+5xH5Y7sl4Z69xG3X9p9PbW
aJRvcQuJvXUb88eiN6mGPkBKNMlnbApGbK+3XuFMIhreSKyiJ4XA8lhK7rfR
NChoN7xYnhMO0mM+95FiOpzrCsW1a8MRaCPlZKaXjiK9C4PGveHIdziv+CIq
HbX3voaypO9Axm1nwNXNGUiNG1GOPX8HtTe6js/xeYzn3nv2OClFQOep7Bxl
qSy03fn1gM+NwLxZfyv9P3KxOcX497H0SMSesbDR3J+LDaIpXjsKIrGiXdNW
1i0XXtMF6RFVkdiV2enQ8yoX+na6lsf6I1GcGV1YRH8G1W9bf7rS7kFgxi8z
bHoG9SZ52bQ793D0Sczby9J5eL3Tr+s2Iwp/mgg51SEF2F0peUkyMhqKvvrb
TqcU4GVEh2pRWjQW5iSHzC0qwPZl4LsXReOXlBndZLAA7kNBC9R6o9H8uuhy
35ZCEHftN/Q1MfDW9p/5MVQIDQNlG52aGDQtkJen7y5Ct0l1Nn3OA3g2bUXs
v8WIO9GaV1r2EDL0wqF6qRJM+aqslW59iLQw7dAp5RJ8v3g2wm34IVqMVnXZ
7CsBg9592G5RPNa+lb/8Z3QJ9ihIfZA8E49+FsnI21yK3xI1pu1pCdCvTJdp
9C6DVMJFxqPARCwpkO0TFa1AjtauNd8bkmEVOmX+dXEFdl2/v6n7YzIeOQ69
rdOsgInfFbtGdjI0FYqLQywrMBhdv6VDLAV/edjeV3hWgXU8buPZ7SnQjbfI
iFr3GmZmDl9KmSmws8mwftL3GoExxUfUEx/hxaDeu/c7qmAbxMmdTknDBwP/
0TD9KvwR8OxmSl4aJpjli41Mq5BjY7rBrDoNm69utXtzugodYtsP9n5Nw0u5
dROvI6uwN9/UPGx1Oor3Ki7L41ThS798p2NeOsoyvl2OiqpGS/wnP53mDFS5
JCgd49cA6RdCT+pmIp979qzNRA0qOC2GGgcykea8rdRBpBa7VsqbilllIuBU
53E3+VqcCDtj3u2fCQN72bibO2pxzt+WtrYpE++PBSxlhtbiabnlj7/PZOGT
wTkFca06OPUstpPLysb4yp0Lcs6xsLnnfoeEBRO7Bv9iBXuxELt0c5HBeSYC
Ewx8nW+y4D2woC7xFhOL5c2J+gMWjvwnGuf/kgld0bPv42pYeJcoba+uRPVx
T+ydIKU3OPmksfjzl1xM3/l3gWP9GzTFXTNqqnwGsZ/ZC5X+bMD8FPcNlr55
0J21fWHGfIf5wmSPCZ0XeMfY2FSn8R6el728DD1LoBox8jNT2IzjqdP7HrMr
MDLZ39eh3Yr1lbrXY4erkMYrbZAJb0NSsv/6w3F1kFF5TL5+bUdOad5Ln+oG
hJ4pPJW84gM0taxdB5a9x8Tu7K/Ovp0I8j1yrpiKc0WnpmBx9UewlgcZ6va1
o59ZE9Ki1g3XKNazyKlOTLm0P1zj3YMjtJmVSzf1YO0eyZlTpZ9w9JiKvpFx
H/QfqZ1+vqQPlVp/H43rHoCnysnln636UW+RoH3VdRCa/KPzOg8N4Fe65akZ
3hfILfgmVfNrAGROqvpL9yFIVDzMVY//jIPMo943Gtk4u7G063fDQXB9nBz5
yznIPmCfXD4xiHcs7j6pI1yMjyguZEd/wdClvCt1B3goNNMXj9r7FaNDkoes
N4/A1yAgqYD9FQxtydTgyRGo9R78mBowBEtZGuNAHB/bNjOiaDpsvFjbo8/W
FsDqj/c/RdvYEIavVil/LIAM7VJhaQcb6i5iO2IyBaiUU7jo/pGNhuTd4uez
BVj5+wkOp5eNJUFLwhVzBSBj7A/vOWxUPLbPtnshwI36mfyEGTbirJwLn9UI
8OQK3W2rxjDuWjkxLAYEGG8IG7zgP4y5p3g7x+SE+DSvw3n65jA0uv+9a68g
RJWRMrkVOAzpBp2h5sVCRDRlTSWGDOOzXG5IFk0IjVaWfOP9YRxyP801URXC
pktUf00WFa9p6bmLGkLUDbtnf24dRvAWy4JJCBErYXPFdAUHe+k2fnuchZhx
WLKpZSUHH8ju9efPCOFY3TxuvIYDmc/BMvEuQmjf0PM8ROdgz6ztkTFXId7M
rvY22MJBSob9bMxFISZ/jfnuMuIgfu/payXeQvw9ejtE4woHHm93XjkdKoRc
X+EjsRoODq2zZV7LESJ8vo9PUR0Hsd9WvgpnCjF3+15L13oOnD3CJ5NyhRCL
bZX+0MRB2+Go0YrnQnw3IV6ZXRxM6BTPnSgUorl2tfERwsHCDBPT7eVChOY8
nEpQ5OKMNVOL3iTEbz7+5lsvcpH/rZnVNiJEwSj/Z8VlLhjeMmoZfCq/U5ZJ
+z24+NKb2u8pEKLGeB3P7CoX36a9w2lECMbyXl+3AC72/W1cfvS7EBJvtmSn
PuAiTcuj8vkPan1y30Rlqrj4Pvlz2ZQEgUimTe7gIh72tqkMcpQI6vidbkkK
PCgeH/0rQ5kgVNuEbrWUh/QS0XkOKgRLXurldKjwcHXycWv3HwR01ursN6t5
WOTyyKhajcBq6HsGcxcPAw40Y69VBMVqtxK9XXjgWa5IcNchuBSXG7aoiger
J4eH2/cTFE3H3aqq4eGoivMTV32CaZtA/0ssHuSuzRBxA4Jb6tbure94aFsv
L7/+AEH807nWkZ08LLYeLw4wJKgtd9JcJOChofbIu6VHCJZ+UXoru3QE6rXt
1ziWBK/WBknJuo1AfHNh6+4zBLrtSxojL1BzJG1yr4Dys2uP78pfHoHExvzp
VS4EjxrrFJd6jeDkPFaD9FmCoPO/aaneGIF2/7wTHa4Epvk3TLSiRyDd+s8N
8wsEw9t84g6UU3OqMx0j70kge/CS5g0ZPpQ91uY/uUmQXLlg3jZZPuw0F4Ut
ZBBob83hjsrx4SQaIetF+fBqToYdjY/7o7ZK+wMIgufYqO1ZwceYYdnWvkAC
ieoDi0W28aE72dkzc5tgYruayHUnPhR23/QVjyD4pNXc9k8ZH67vRCZXPSRI
GpH1Ea/gY+H15x0elO0zzVaGVfLRzfTwrKXMXdHplVzHR2L/uIR9PBVPqU+F
1czHiQfVRrEJ1Prn8k/LD/ERs0mu6GcSwX625DRTWgC96Odr76ZR+5WIFV9t
BMjKm9OllEMw5j0maWcnwFCzT8gJyroWj0d6HQQIeOov/ZBy6YIF+Z3OAsj2
VfjRmFR/+fXteXtJgDIrsdXyuQQ99j4n8wIFcP6cEyvynGDOqtI4P6YA6OCP
vyggsM7VlVWZEeAQaeEllxK8OP5DplmEmmP3p+LdlGXmlMy/KS6EyPVdgkVl
BK+tdkpxpYTw5i/fwKC8fO7e2QLqXDKIUQp3fEXAdjQRGK0RQrflU4pyBYGb
omv9NXMhcrxV9c5WEfgFpPh/yhLieF6+evAbAmlPyYN21Nw3cze4FVKOcz63
iE3NuabYipgByvmGW1NJgRA3XUoYuvVUfNmWGrEKIa6KB7cMUDZMFJFa00ad
i/W1+uveEiwuOhnuPiWE4no3s4JGghyOykMZI4LUzOxFjq1UvyUE2PYdJkgJ
n9cTQtnRVPAn04Qg4hrXL5+yclnZs0MWBDVl953E2qh+uXOiLtSeYNyZfT+R
srPug9G5/xB0jL3Lbm6n8mfIG0hmUHW9yb6k3Ekws+WadOcT6jvr28t3U+4S
fGnNyKbqWE7Pd6AcYZFvvZ+qw+G7I55ZlGc1TC8HUPsYw5ON1P1I9VvLnURx
6j+qfnzVA10ELwMnHdo6CJ5vkeC7UI7afnJNKvXe+MnyuFDKhun0gr/6CNSn
rye/p1zi1VjvzyU4/Z4Ym3cTRGttumPCJ/jjR9cGd8oXB+PN1QjBl7P0qfuU
Vxu59r/+l8r3ntfGTsoSou1pEZMEDL2kjEnKA4XbXWx/Uc83RU4v6SEoc0ld
p/0fde4UnNHdQjlGdd747CzB/Imd5paU/wcrpKuc
            "]]},
          Annotation[#, "Charting`Private`Tag$26862#1"]& ]}, {}}, {{}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.0055000000000000005`], AbsoluteThickness[1.6], 
          PointBox[CompressedData["
1:eJw113dcjW0cBnAjKyublJG9R7J1mVnZm8xkvbaMZKaJknYJLSKltCiV9t57
d05nzwclm/cpv9M/PueTnvPc6/pe98hDZzabdGjXrp2kY7t2rf+qfsz7v17J
vWqJXy8/Fa87tBGTNduPUw9VfT6A+bVFN6N5lvhx40ScqfMxbKyzdBMNtqLf
n8XUSebhT9Zb4d9zTNE7bPCwOksr+v9XEGGfPTbgvRXMJGNndA2+BuHtntpf
Pqv+/hYct5puLZpojQUNG+zdT1tgyYbXo+YbW0P1XtMS7PnTva0Ru7pHh0k2
1vC62qs5usyanm8LhxtG4Xm9bTAvMvNSwtO7+Hm72uLSGhv6vvvo4zs+N9zK
Bu+GWck2vXNARfdJ3+4k2dD3O+Jzgoap7LcNvFNHhJ0/9hBWK3aHSefb0vs4
YZwlZ8SdK7ao+d72hSi/X7MwOtqW3s8FDeaSG1ZfbDF0Wt3ozrdd4bV3gMnn
WXb0vm74rVV/9o+pHXYfbnsihFcXjXsVZUfv7wG1XqEaP1rs4OXZ9oNta+fE
yOfepfF4Ybj7iTBL87uoyv+6hffsEdbuM7VLSbhL4/NG17nDaoI63MMQte29
poQ/xoveFVi48h6N9wkebI/dZnb/Hnb+eyBG6bsnGBffo/H7wPv71+G/B9+H
4Ry7JQcEPvDqMmHUygP3aT58cfDO4ovLXtyHfFPrE33hdX1MRtPH+zQ/frhu
5J63a7497p0U/Xz5yQ+Ll+ktuWxpT/PlD79bVvXrCuwxyWYgOwR/XFP7bdKg
6UDzF4CRtZ29pxx1QI7viuim3wG4EtV99fxIB5rPZ3A/ZjavS4cHaJ2twWrP
kWyTvcht4wOa3+fw/WI/UfL0AdQr/M8uUg/E1dunK9t/fEDzHYiA5dp9z41w
RFDrcDVe4EPhxa8rVjvS/L/AqJEVDrfOO2JN6+MGvkSX+Obckd6OtB4v4bjZ
/8XEdEe0TbdWEPq41mp6f3Sk9QlCgHbKsMtDH6J19gp1XuHzuudv4gwe0nq9
wii9ziWnzj/EBKMHyV/GB+PDJXfLe08e0voFw1x3ZfbAnIfIYp+mOS0ECU36
QRrfHtJ6hiDTWTz56hgnHHuouAa919gx47PR1i1OtL6v8VJz80qf205gX44d
cijCE9b83BfmROsdilNPTt992OCErezmYZcUQxDQdV5vZ1r/MCzWnnh/O5wR
wM5ehX8Y9vp4TZGdcab9EAa13h+7fPNxRvPv0iuSsW+w6Y7t3BvFzrQ/3sCq
YULDDTUXLG/bAG/wJLnTsJ+zXWi/vEHMtfghzcddwC5e64AwuOOR7yceu9D+
CYcv757wQJEL/r1eOFznzOlZ0cmV9lM4Drsyj/Lnu0K37ScCe9onpqw660r7
KwJ/aue/M3juijuJ2ewOikBj6fXQrFpX2m8RSMrVayzu54YSdvXZOcHjLL97
xmvdaP9F4n3VI8m1O24Y1fbASEz54bRwRJwb7cdIZH33yFj/xQ0X2hYgCnq7
DhZ2m+ZO+zMKjeMMnm4+7o6UL+PZFY7Cy7q1NycGuNN+jUb29DEJTg3u6Pvv
BaEpnh18f6gH7d9odH75yHvITg+0vV1mNOYnjRPMd/Wg/fwWnrKT9Z+KPRD+
mAlMX/sWBhoJPPTxpP39Fmf/ztwydqMnOrQN+C2ejA22f/7Ak/b7O+gsLm9J
KvDE5n8bBmd+TflipuFF+/8dOjzr+Llgkxf8Wodb+g5W727sSnP2ovMQA8lN
X6895V743Pq4nTHo6W04zmXIIzofMXhzb6np2b2PsLR1+9fEwHDCpO+ffB/R
eYnF5RHu0zVFj9C2fffFwnl3Rw/pZG86P7HQP1QbeviCN7hXWxckFivib/+0
i/Wm8/Qe0QtXtzfq8Bgz2jb4e/gstVhTu+Yxna/3GP/7inEPl8e47bqM3TLv
8dDNrEZR95jOWxzid5Ustxj/BEU6bQuMguyaWTkXntD5i8PViE/Lsj88wcjW
9FfEwUHnVs/bPZ7SeYzHzn3tp3/a9RTnWo/H2Xjsy2q/TuvFUzqf8egUo9NN
7etTJLWublM8yp/06rxtlA+d1wTwTh1qZ7vWBxr/DiTeaVplDzX1ofObgC7F
fxb2eeyDtun7zv6eI55xPN2HzvMHdPqotBz90Qdtr3ftA7Lnbuy3XNOXzvcH
8EI47zOX++Kf5olYpK2+P+KML533RJiszKtXe+SLF+lr2ROdiCWK5GEJ6b50
/hMRV6Klxvnsi43saWT3KHrFP7/433A/ygP273dk7Dlo6IdvrfqkJsJzxG2f
ZDM/yodE+OpfmOYY6AcfzdYDkQRbw+ov2WV+lBdJsBz6sPGCmj9WtS1wEj4d
7SO8p+tP+ZGE7Numn/sZ++NjYOsLJuG4+MLPHs7+lCdJGF11vOVyij88Wt8u
OwnWhZuqtzX7U74kIeF35tPAMQFo3VyreySj/YteS07vCKC8ScY5K7XYILsA
tKZJ8bpkrOk4reu+uADKn2SIbkRNtGMC0Lpb2C2D1WcO6UwZ9YzyKBl5N/w/
rdnxDG3DLUxGiIudN+/eM8on9u/zMsd/S3yGtsf1TcEck3eP7VqeUV6lIMZu
s5rn5Odoi6MtKagaNevoOOPnlF8pGGvrXazn9Rz/6kMKXgw6tTWt6DnlWQo2
5o3/Vd4tEG1xUJECnXdN1f8tDaR8S8Gzg04drcwDcbMtcFKxcMTxu6OiAinv
UvF5q7nVKiYQ49oOUCqMXE70+DLhBeVfKnrHXpw70eQFCtrqRSo2Ci+Plvi8
oDxMxR0TkXxm3Qtcbj0OdalQ3CoN7ar5kvIxFTvniu//t+Mlhrdu52Fp0OYf
embk+pLyMg26HzS16kpeIqN1+x1IQ/TOqwOUfYMoP9NgLb/x1nFzEM60BWQa
rvpFT0h3CqI8TcPrCH3XByVBGNR24NNw/WLSaGX/V5Sv6fBZbarG2/4Kba83
Oh0dLgw4dMbzFeVtOtLyJSddal/hSNtPOn6fjd22ZUQw5W86fEcXrw85HIzW
8sQmLObm/7jn/zKY8jgd8WFXl85hgtGajlck6Zjv0uflcb0QyucMeOR9+brg
Wgj2tT0wA6PGFR0LTQmhvM7Awpfxs9O6v0bntoDKQMzSpU9ubH1N+Z0Bk40b
6mofv0bravR/nYGJHRfockWvKc8zsCYgnnN/Zija2h2TgVlvFeO410Mp3zNx
6NsCw4asUPxpPW4zMmHkPWf18u5hlPeZOKu+cQQzIwzPW7fz+Uz0XK5Z8HFn
GOV/Jr4HK76vvhWG9f8GjIkRKSt+BoaRB5l4v19vZ6fCMLS0LseXTNiUXqk4
9i2MfMjCyOh9J8aOZN1vHe6cLEhub3inv+YNeZGFloJVx8MvvIFB2wRmwf5W
9RiLx2/IjyzwS0e5v854A2UbD+zfny88MvfzG/IkCxlbBmCEdjjcWuP9ZxaE
3iHZR1aFky/Z0DGuiOxhGo62eF6UDcWOjS/7+oSTN9kYfmDzVrPccLTF641s
TOrveGjJ93DyJxtXJsy4/d/YCPyLv2y4pubv/bYlgjzKhnrkmefi2xH/4qVD
Dm45W3ZbFBZBPuXgTzlPr6U+Am3Hd1kOmuPTqvv3iiSvclA3JTHQc1Ek2o6H
VQ6mmCYbWp2KJL9ysIrRuFr5mO0dbRsmB5ufLOY450eSZzmwdZ85/d3fSJS1
Lm/XXIgO9x5qMCOKfMsFY7B0s75xFK63vWEupu81tfF1jSLvcnEj0ubk2cwo
jGnbgLkwm+8U8epnFPmXi6Vulb22TotGXm7bDwYJM6efMI4mD3Px4uCXHKV7
NC62BWoeBJxezo250eRjHvqvbJhr0OEttNs2dB46Fy87qzn3LXmZh2Hrxsn3
n36LtLbrVR76jL+2f9Czt+RnHsafm319Se1bnGoL6DwUhK/4Vt/vHXmah/WR
Dnc+rn2Htrjqn4/L2oKiS5bvyNd8TPmt5nY+/h3a4mB7Ph7MvG/Bb3lH3ubj
1bum1XnTY/Dv+sX+vtT2yZT/YsjffHiPOTKz8/MY9Ph34PDN/eTbHdwY8jgf
HectbNDRjkXbdtEsQLRi8p5ju2LJ5wLkqK+sHe8Wi7blMCrAs/wF/UxKYsnr
AtTG7Q3T7vMe/65nBdhR42G8ZcN78rsAzYUnajs7vEfb4xoK4COozpqb9548
L4SEd65a0iOOznshmmpmxGmuiyPfC1GeYa+Vax9H578QJn0OvviTH0feF2LP
d2H3CI14yoNCfG9481myOZ78L4TmGb8uAa7xlA+FKH5W38ytjKc+UIinaZnG
AVoJlBeFWDay5ZfsQAL1g0LsyxtzLOZZAuVHIUQz5xj2lCVQXyiE150RG8TT
P1CeFMJGQzBg7eUP1B8KwbG4t3w222P+5UshHv/shrIOidQnihDxe/bM0rGJ
lDdFcJgzx3Dw2kTqF0Uo9W+s9T2TSPlThNOPzHTPuSRS3yhCstfTJ9YxiZRH
RSj8tu1kY30i9Y8ibHVskViqJVE+FUF7z9TDpycmUR8pQujfr2OfbkyivCrC
5U1u1wddTqJ+UoTAFIfwisdJlF9FaMr3781JTaK+UoQTawrqpsuTKM+K8PDt
QpPsfsnUX4qgbrFMEbwgmfKtCMqZB+MqjZOpzxRjsofu4PX3kynviqHFOb+i
Z1Qy9ZtiuL6U2QysT6b8K0ZGVzOdE11SqO8UQ/v2n+1dZ6RQHhbjWvPBI592
p1D/KcZ4WyvHMVYplI/F8Ly5pP+z0BTqQ8UI7Llv2tnqFMrLYqzXsh9m0ymV
+lExLnQ/rC2dnkr5WYx2R7cd8TRKpb5UjBiXoJmutqmUp+x4rl6PrYtMpf5U
jBFDw2dd4KZSvpagafrH8m290qhPlWDle+8yqwVplLclML9bYtzxeBr1qxJs
ST0Qke+WRvlbgo3Gsc2i1DTqWyWY/Ob79g1NaZTHJfgcV9qlk0469a8S6OW5
63bblE75XAJR/rBORrfSqY+VYENS/rPfoemU1yUIGv5gkbIhnfpZCdRWp3ye
opFB+V2CYuEyadziDOprJajevHPlo3MZlOcl6PY7dFKWXwb1txIkFu15a1Ca
QfleigW+GV37d86kPleKQNuIhTPmZlLel+Lh3imnn5zIpH5XCn2tu9F7HmdS
/pciaNf+KccLM6nvlSI+5po4vWMWeVAK0xOCn6as0//6XykeBWSZXfgvi3wo
RfXMW+eTn2ZRHyzF2783mw+WZpEXpVj1alePzd2yqR+WYueEN6lO+tnkRyn8
j80ZO9o0m/piKb53DJjTOSibPCnFr9+vf8/lZFN/LMViW3XzuIE55EsZHMuN
whzX5VCfLMOAv/O8oyxzyJsyzBZrL5gUl0P9sgyS5UdsPjXlkD9l2OhjaNZ1
ci71zTIYm3E6nz2cSx6V4eaZ9qPHP86l/lkG/1u7kqeX55JPZZixe0eRXe88
6qNl6BQ/Xl93dR55VYaYL5FNk+/kUT8tg+C4e/Wl+DzyqwzzF2hWaXzLo75a
hlXtZ5T8mJlPnpWhyOtC0KzT+dRfy3Fh0BD9hJf55Fs53GM/n3AX5FOfLceJ
wUe7J44sIO/KkbXu1Nc5+wqo35bjxxZbtb9eBeRfOX6dTe40oLKA+m45wg+V
CuLUC8nDcqy4NsrJbUIh9d9yfL798JT7qkLysRz6D7yNPhwtpD5cDk+++dJe
NoXkZTmO/9Fub/m8kPpxObbFjzUbkV5IfpZjpY/Cjsfm+r++XI4HGtf65HUq
Ik/LYZ14T147poj6cwXGalu19DUoIl8r0FJ8vduZI0XUpysQc2TOn4/WReRt
BWY3agd7BBZRv67AQ1cDiUlmEflbATvF6bs7JUXUtyvQcf/PQ2fVi8njCiw5
1rI8ZFIx9e8KjH7Ru31vNif/+VyB8u6Pjriyufivj1dAr5/jssWOxeR1BfZP
tjvcO7yY+nkFgsemuHQoLSa/K3A5XuvJkJZi6usV6LRg4Oqtg0vI80pU6BYZ
hc8vIc8rEb9b46nu3hLyvBItp+aUVtwsIc8r4bto7YcnfiXkeSUKcz6NtE0r
Ic8rcdP6fJyzuIQ8r0TIo+h1Sd1LyfNKcI6udtaYVkqeV2L1352rb28uJc8r
kV4UPbHfpVLyvBK/j1d+TPMsJc8r8X5l5RbP+FLyvBIFZZfU7bml5HklnMyH
Zvt0KiPPK7Egu8ygdEIZeV6Jmen7dMetLyPPq3DTUHOJ5/ky8rwKl1wGj5jo
XkaeVyE5b9DNqvdl5HkVONVW3QM5ZeR5FUaPtzjn1KmcPK+C0dOGPY8mlpPn
VRj0bqxp0oZy8rwKAQHuh9QulpPnVUh9NDnf2KucPK+CNryNGj6Uk+dV2DfM
KPySoJw8r8KIHf7nJ3SvIM+rcHTTKXydXkGeV+GDx7WYhu0V5HkVhEV7bzWy
++yf51Xg5p6Y+devgjyvxsmIvJN6WRXkeTX8tljmWLL75J/n1YiwNpbIB1SS
59V4EKF54NTCSvK8Goan5mR2Mq4kz6sxfMqEmCi7SvK8GmHjVmRdC6skz6vh
tCjbfndFJXnOfl9ngdf6P5XkeTW6fIu4s3NMFXleDfeAtdWXDavI82qsG/R+
WvCFKvK8GtYrNfp/8aoiz6sR+dxo4ObkKvK8GlMyYlKTJVXkeQ3Shq4PNuhT
TZ7XIL5g0666udXkeQ2SovvtsD5QTZ7X4NrWP9OW2FaT5zXY8dP7VE/2vf55
XoPf7X09pBXV5HkNJvQWza38W02e12Ct+bWysnE15HkN7DMCmhs31JDnNXgx
oWHo38s15Dn795PlBRN8asjzGlw8VOJyKLOGPGf/nS5RvPhYQ57XYJV/6d5f
g2vJ8xrsP+B6xGhJLXleA3cfO7vM47XkeS3cFpuuXOxUS57XorHBYkRKbC15
Xos+GltdNvJqyfNazHD+PVLcvY48r0XCi76r7GbVkee1sD4WZKm7t448r0X+
xCPbhFZ15HkttH9fnuP7uo48r4VS52WgcUUdeV6LEN8DOlPa1ZPntVhU8XLU
7/H15HktTmouGFeyqZ48r4Xt4WsvXl+tJ89rcb6jg84D/3ryvBYXtmzQMM2t
J89r0e+WidzoSz15XgeDtUtnrB7WQJ6z46ieuGvuygbyvA6O1vr1E842kOd1
yIy8Pl/bs4E8r8PD4EJRn+QG8rwOZaa9DnSRNZDndfCovO3wpxuHPK/D5sSF
hmbDOOR5HTZOHfa8aSaHPK8Dx7zPrlMrOeR5HYLs1XWEezjkeR2cl3RL2HuW
Q57XYU2NlrTEkkOe12HimaPGqzw55HkdBq8dKYoN4ZDn9ZAZec+clMwhz+sx
yG1lR89yDnlej6mG4QPVZBzyvB4c1wcTT/7lkOf1GBYQ86uwH5c8r8efgXpr
Z47nkuf16L9+VumDhVzyvB7FEZrbJRu55Hk9gjveeqBvwiXP6/Hwri4emHHJ
c3YdD2J0jT2XPK+Hq6/s90g/Lnlej78Pk+4ZR3PJ83psyHewf5rNJc/r0f7I
sZSyei55Xo/FzSGKTk1c8rwBaaE98qZ1aSTPG2B53kJny9BG8rwBZe8OhJ6Z
1kieN6DUfYOe1bJG8rwBgUEbLznvaCTPG2DSkjTL679G8rwBe1IP6HrebCTP
G3DkvKGeo3Mjed4AnXnd/twIbCTP2eeJLDaZvG8kzxvwzKV389KCRvK8AbbL
zwcN4jWS5w24Y1k3j9fSSJ43YM6PL0efqfPI8wZ8EJX+2DuMR55z8NdfPabH
TB6ddw44nX2Pv1nBI9/ZfZSlXbp2F4/OPwcWVZaJNSd55D0H33x2fz1wi0d5
wEGvxc+X1DjzyH/297FNx9YE8igfONinuX3O61ge9QEOQm4cv9w5n0d5wYHW
2G3Nm7k86gcczDPJM3Nu5lF+cCB6OKA0owuf+gL795z2hYwmn/KEg8Wzvs/v
PpVP/YGDKQ0tjUOW8ClfOLi7477b0K186hMc+NU2jtY4yqe84eBmr/i9X834
1C842HT5Sp+i+3zKHw48Ovfu8/gpn/oGB11vDJ66J5xPecTBuMKc+d3T+NQ/
OHA82049pIJP+cRBg+u0M4ulfOojHCw7dVE/7Ref8oqDUDgZLugtoH7CQYSt
7LTfSAHlFwdNf+6b/9AVUF/hYJvDhNXLDASUZxxUL9MKvbZTQP2FgxmTjls9
PyGgfONg4LCEgMRrAuozXNywuiDMdhBQ3nHxrIf96DQfAfUb9hwZxM4LDRdQ
/nFhvnRHe7tUAfUdLoZqntqxuVxAechFxNHwwd3FAuo/XMyIvqcT8V1A+ciF
Haf7ljXdhdSHuOhyuNm2UEtIecnF7GdR/iumCqkfcXFu3Jk7LyGk/ORi2Vum
66+NQupLXITaRA1YdEhIecrF3aXKeycvCKk/cdGxx+X1dy2FlK9cXDerXezm
KqQ+xYXEfPdKx+dCylsuYhbVGlx5K6R+xYV388/R6zOFlL9cTA3skNmnSkh9
i4sEoWhgikRIeczFruubvh/8IaT+xcUaDecdSnUR5TMX096GaB4bKqI+xgXP
ZsPkgkkiymsudg8quTRqoYj6GRdBjL7gsKGI8puL4z+XGzsbiaivcfFIuKTu
9UkR5TkXWy4uXPz2moj6GxerCvlmwfdFlO+NsMtwvuLgLaI+14hb3ddM2hss
orxnPxuF3h4cJ6J+14heb/qcTMwRUf43YsiWwzVba0TU9xqxrVk/qUwqIg8a
8SWmtveyHyLqf41gMDTWu5uYfGhEDxeziMbBYuqDjZCPH/Sx/3gxedEIg6Tw
w7pzxNQPG2G75WVnfQMx+dGIocLmnFnbxNQXG2F+eNKrQYfF5Ekj9p66/kh4
Xkz9sREDNKNcfG+LyZdGXP0w68EqRzH1yUZcevHlTs0TMXnTiJvZ1id3h4ip
XzZimmn3ZWnvxeRPIya0zOuglS2mvtkIe4Pdz/ZVismjRox5uny8vVBM/bMR
jgMP2AQ2i8mnRljP8v0Q0kFCfbQR927b5j/VkJBX7PP7jgy/PkxC/ZSdz55/
j62cLCG/2OfN1P/8a56E+mojCjbfWvdkpYQ8a8SaYtnVSdsk1F/ZXNWeaB5w
SEK+8RA47MiGbmcl1GfZnA9d833XdQl5x0Pnbu3N3e5KqN/yMGvqtIoEdwn5
x4PNjhaNkgAJ9V0e4nd9mlTyRkIe8nC1u+7kDwkS6r88lL9YPtA9R0I+8jBe
bzSzu1JCfZiHRTVL49QFEvKShymGd+48/yShfsyDUmvu2ql/JOQnD9U6Opp+
6lLqyzysLTj9tf0gKXnKw5W70dJ1o6TUn3lI6vnzp9U0KfnKg238M92gBVLq
0zxkdrD2il0pJW95uBtQMO/tFin1ax72jxwy1Ge/lPzloX233asu/ielvs3D
+u2SnNmXpeQxO79Z5mF8Cyn1bx4Uy3Ta33SQks88zNRZldXJS0p9nAd9rzf9
zZ5JyWse+g6+ylSGSamf83DLu++RMXFS8puHL1qLrfZlSKmv8+DDH7LbslhK
nvPR6WuXz251UvKcD23TpN1uYil5zsfoP2v8LJqk5DkfE0d+Kd/zR0qe82ER
t7uLTjcZec7HlNP7l5b1k5HnfGxpN9X94jAZec4Hr3L2gI4TZOQ5Hzf+m5p7
XVdGnvMREtI9V7BIRp7zYWPfecyCVTLynA/OCTfR9c0y8pyPDs5nR782kpHn
fAz9YC3POSIjz/nY59ppU8VZGXnOx/fJW4wLr8rIcz5ELj313lrKyHM+3u3t
VnnPQUae87Fp7MZ96z1k5Dkf6ZOGV//xlZHnfBjWCvY+fiUjz/kwOrX814Qo
GXnOR6+L6pkBCTLynI94a73Snpky8pztJcWMrkmRjDznQ5nt8COoWkaes+uT
8ZBdZhl5zofuGx/1DgoZec7H0/C75v1aZOQ5H7X85X79/srIcz7CO0ucOnSV
k+d8XJ8VeZCjISfPBZgc0GlUyBA5eS6A4PBe2TEdOXkuwG6dPQX9JsnJcwEy
ZHG8EF05eS5A2J+9i+YslJPnAux01fsUulxOnrO94t22/oPXyclzAaIHSd+f
2SYnzwUwq5/9O3qvnDwXwMTNXaEwkZPnAizrfMOv/2k5eS7A80HnZky6JCfP
BZB/6B4244acPBdAXfFTd7y1nDwXYMql9iW9HOTkuQDTNRR+Alc5eS4Ad8XT
mODHcvJcwO7/LlMOP5OT5wLEZen06xEiJ8/ZHrX25+2ASDl5LoBv9fMHU+Lk
5Dk7n0sNNgemyMlzAXSWq3H65MjJcwGOhfTfcKpYTp6zzzsZERVTJSfPBfge
NGriV46cPBcgRPE6eYxYTp4L8GZWgOsKRk6eC9BOZ0fM9hY5eS5AkM7w5bt+
y8lzdnyPthuuU1OQ50JEDt0tntldQZ4LcdreZqp6XwV5LoT7uUV6JYMV5LkQ
D4ZU93UYriDPhUgVxtbNH6sgz9meY78guGqygjwXonF3lMsJXQV5LsSumzdC
FPMU5LkQK/dIehgvVpDnbM8qXJaVbaAgz4W45FL+bfQ6BXnO9rh7/DfntijI
cyE2jgxSf7NLQZ4Lcfu50RDefgV5LkRQ5ITmLkcU5LkQ3x/OfjP8pII8F+KD
59MjE88ryPPWXnd54vgrCvJciGqbjL5DbijIcyF6HbHX+31HQZ4LcWxF3ItS
OwV5LsSCywssnj5QkOfs/A0UVux1VZDn7PxZ3/3Q+5GCPGd/3yRZH+WjIM+F
KKx4+2j9cwV5LsRWC8/3Na8U5Dn7fb5jUo3eKMhzITwKeNlF0QrynB3vKEPB
/DgFeS7E8ZT08R5JCvJcBD2zV0GSdAV5LoKR4ZJb03IV5LkIE7QSEo4XKchz
EfwWGZh7lCvIcxGahsqz39coyHMRLtaEJpdwFOS5COv/3rjMESjIcxGe9pn1
hytVkOcizFrv9l8loyDPRYjZsKs0pVlBnouw+4PWuoDvCvJchLLOJ6VmfxTk
Odtb1RISl3dUkuci6Kaa8tW6KslzEabMqjse20NJnosQzZ9qYtJHSZ6L8DLl
nrTTQCV5zvbq2r89vTWV5LkIYUNv88YNV5LnIizXYG68HKUkz0UoiOzbPGK8
kjwXwXb2k+0Ok5XkuQjqiZNjP01Xkuci5K1YOnOtnpI8F+HF1rXFXvOU5Dnb
cxu9QziLlOS5CHuGm1ZoLVWS52zPNmw5tMFASZ6LEOEw3fjKGiV5LkKH7EEy
9/VK8lyE2Cu2PV9vVpLnYky3NaiL3a4kz8VI1ft9MmG3kjwXo0KhX/Bun5I8
F+PCZKuBQYeU5LkYzvunb3M6oiTPxeiRY/P43Aklec728tMx7VedVpLnYiyR
VjkPOK8kz8WYe5JvUn1RSZ6LEeYaa+dmpiTPxeh5fHSfNdeV5LkYxrMLe365
pSTPxWiw2/zQ3VJJnovhM3RsyAxbJXnO9ujEV+dT7inJczHu7pvYZPhASZ6L
oazNR76TkjwXo5128PGVbkryXIxAZeztd55K8lyMvV2rPUY+VpLnYvz4m/fB
wkdJnosRvXxJl1p/JXkuxkbN1JtTApXkuRhbXXPnXQ5Skudi3PCIXhEToiTP
xZhVkhLRFKYkz8XonTTUZWykkjwX43hQTsumt0ryXIxvz9JlF2OV5LkYRo8E
Fk7xSvJcjFcGopzARCV5LoZi4dnyqBQleS6B9Sf18Lh0JXnOfo4bdzouS0me
S1DL1dKMylWS5xIYe45Ne16gJM8lSBt09MbDYiV5LsFsD+HWi2VK8lyCL0ec
d2yuVJLnErjIl7iOr1GS5+y9ITZY61udkjyX4Kqr6a9EjpI8l+CtzYhld3hK
8lyC61Ltv/pCJXkuwdbLifrNYiV5LsHIHUs1/GVK8py9t8xMvb1WqSTPJfix
ZJ2v4qOSPJfg3YGca3ZNSvJcgi3a/XSGtyjJcwkG5jU8ef1NSZ5L0PtpSbu5
P5XkuQT5oZe3vf+tJM8lCNJeGDS3HUOeS3Bf+0qPsA4MeS7BXI0TdiM7MeS5
BP085ky178KQ5xIcDFV0aerGkOcStCTcGrWlB0Oes/ekqFz7kF4MeS7B0WGu
Wzr0YchzCYYtsLq9qR9DnrP3HM3pox8NYMhzKULDNFY1DGLIcykmNtl919Zk
yHMpZjPC5Tu0GPJcCgh2z743jCHPpYg61L46ZgRDnktx36h+AU+HIc+l6Lym
/miXMQx5LkUJyo+PHceQ51Kc6m65bvEEhjyXInO9w7BtkxjyXIpF3pbiw1MY
8lyKSKeG6DPTGPJciiE/2rlenMGQ51J8a9B1uKTLkOdSHIsMDjmvx5DnUgw+
avfn+ByGPJfisWXRPaN5DHkuhXzLpQNrFjDkOXtPa95so7uIIc+lSDSb3H0w
GPKcvbf1DG/5upghz6Xo1f361pKlDHnOjk+rz7Sg5Qx5LkXqnlrHawYMeS6F
7tk5lmtWMeQ5e08LNunXfw1DnrPvN811VdVahjyXQqddwxSvdQx5LsW7qSvz
tm9gyHMprjbFTe+9iSHPpTDP0zyUspkhz2Xo4Df2vwtbGfJchrqah9uGb2fI
cxmS7n6fmLGDIc9leGshbD6xiyHPZfAwz4tT38OQ5zI8nN77/nMjhjyXIeLm
t//09zHkuQy2dskmxfsZ8lyG14VOFocOMuQ5ey/ZZZ6jPMSQ5zI8yLq6+vJh
hjyXwdvmWqdfJgx5LoOwemfX60cZ8lyGn+plu38cY8hzGTbc9O9keoIhz2Wo
rDfuI/2PIc/Ze+M0f2ujUwx5LsOtlppTOacZ8lyG3dsvFcw+y5DnMkimZUY8
OceQ5zL0XzJgVMcLDHkuw6HlD8YcNmXIcxn2aGxPSbrIkOcy1EdfVtO6zJDn
MhwL+Pjx/BWGPJdh1YdXjulmDHkuQ461hWygOUOey7Bp5qJuh68x5LkMWwqd
PoVcZ8hz9p5rMiG4+QZDnsvwXt1t8dxbDHkuR0jPvdFXbjPkuRwRPT/3j7Zg
yHM5Dgx9Z/zxDkOey1Fd3j5onBVDnsux2qEPs8eaIc/l0DMaMt/ehiHP2XvD
igVO720Z8lyO8p12f4R2DHkux7TsTta97zHkuRxr+/lP07vPkOdyLNFY3mGn
PUOey2GhHad2xYEhz9l7mV7BItcHDHkuh6jnzMhQR4Y8l8NV7+7pjIcMeS7H
uhmnL9Q6MeS5HGM8TbKUzgx5zj5fT3LpjwtDnsvxd6ixRXc3hjyXY07u1M8D
3BnyXI6Thp652h4MeS5HLwPuoFGeDHkuh47OnOoxXgx5Lof1mPghYx8x5Lkc
A+xtakd5M+S5HIEvvccPf8yQ5+w972dzx8FPGPJcjpkbzU17PWXIczkkVT1u
dfBhyHM5Jm0wm9zMfv7nuRwd1O7c5Pky5Lkcnp0abhT6MeS5Akv19Ka/92fI
c7YnG4/w9A9gyHMFep+Y+sHuGUOeKyBdyQk89ZwhzxXYMT9o54ZAhjxn73UV
TeVTXzDkOXuPmDh4fI+XDHmuwI0sk00i9vM/zxVY3DRsZ2IQQ54roB9otsT9
FUOeK3Dyw9O+J4MZ8lyBsG6hhfohDHmuwMSDMbd6v2bIcwU2ZaSOq2c///Nc
gU/XcrKDQhnyXIGfmnGnL4Yx5LkCM1+ZDdZ/w5Dn7D0lsSZXLZwhzxXY5vja
MYv9/M9zBWq6Rhy7H8GQ5+y9SvZk17pIhjxX4BFX61iPKJXnCsg9cj2yolSe
K/C6fOlHy2iV5wpsfdjtgv5blefsPa6Lz5iWtyrPFfAov9oj+J3KcwWUfdUm
H4hRea7A+b2PbPvGqjxnx/fDdWxKrMpzdjxSffXz71Wes/eKX+VzhsepPFdi
pcglIjtO5bkSn/Y9u20ar/JciUsTJz3XSlB5roT7PIxJSVB5roRt7McOxz6o
PGd7z7QVy7onqjxXIvrOXGlwospztgeNjP9pmKTynL2HdIozlyWpPFdCv7z/
adtkledKmLY8Lh+VovKcvZd5jY6JT1F5rsS++pNDtqeqPFfiaPmMn/JUledK
dAscc8giTeW5En/dPm0amK7yXAm7UbPzX6SrPGd7ff971fMyVJ4rMczh5tWs
DJXnSvwe5/VmR6bKc7aH/9p5h5+p8lyJxR9mNJ3JUnmuREjuQ7UfWSrPlZCi
PtYiW+U5ey9UMx2inqPynO3dwsBRjjkqz9le+iq9tn+uynO2pw8dsMojV+W5
Eu21C49p5qk8V2JqwuSlj/JUnjPINtxZoZmv8pzBfpMrkzzzVZ4zWNvp1ZKB
BSrP2XF7d9Z2KlB5zuCstldij0KV56xLfU5Nsy5Uec6gOuDmyT+FKs/ZXI/P
vnqxSOU5m8NdFu+VFak8Z/d1l+yBB4pVnjN4OsEwqKRY5TmDn6V+mgYlKs8Z
KK74n3hbovKcge0TrcfjSlWeM9j1rSzSrVTlOYMV5x0j1cpUnrPndEnnp+fK
VJ6zbjdJL9WVqTxnv/e6TH9Vucpz9jnrXn9/U67ynHVtZ8lLzQqV56xjY5s3
WVSoPGfgLvZvFleoPGdz4NYZl/WVKs8ZXHGq042oVHnOnivX02UDq1SeM5jP
8TY3q1J5zmBqVO34miqV5wxSt3o3LKhWec6Oo9tyv0fVKs/Z99WbfeFHtcrz
jyg9VLB1Zw2D/wFfaalh
           "]]}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.0055000000000000005`], AbsoluteThickness[1.6]}, {}}, {
          {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
           0.0055000000000000005`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0, 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
         Charting`ScaledFrameTicks[{Identity, Identity}]}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{11.052083333333343`, 1.552083333333428}, {
       6.264878372495332, 0.5000000000000142}},
      Method->{
       "DefaultBoundaryStyle" -> Automatic, 
        "DefaultGraphicsInteraction" -> {
         "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
          "Effects" -> {
           "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
            "Droplines" -> {
             "freeformCursorMode" -> True, 
              "placement" -> {"x" -> "All", "y" -> "None"}}}}, 
        "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> 
        None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0, 1}, {0., 0.9955255170204409}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {967.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}]}, {}},
  ContentSelectable->True,
  ImageSize->{
    UpTo[600], 
    UpTo[360]},
  PlotRangePadding->{6, 5}]], "Output",
 CellChangeTimes->{{3.8130541631866355`*^9, 3.813054166972516*^9}},
 CellLabel->"Out[70]=",ExpressionUUID->"65d4e4ef-a5ef-461b-96ad-2d035da53501"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"absError", "[", 
   RowBox[{"approxSolution_", ",", "exactSolution_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"exactSolution", "[", 
    RowBox[{"approxSolution", "[", 
     RowBox[{"[", 
      RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}], "-", 
   RowBox[{"approxSolution", "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", "2"}], "]"}], "]"}]}]}], "\n", 
 RowBox[{"GraphicsRow", "[", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK310", ",", "exactSolution"}], "]"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK3100", ",", "exactSolution"}], "]"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK31000", ",", "exactSolution"}], "]"}], "]"}]}], "}"}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.8130531630599203`*^9, 3.8130531911121416`*^9}, {
   3.8130532559461966`*^9, 3.8130532690711155`*^9}, {3.813053311132979*^9, 
   3.8130533474193425`*^9}, 3.8130541928613296`*^9},
 CellLabel->"In[71]:=",ExpressionUUID->"60afd7b7-e7bc-4158-80c0-baf5b2b32e9e"],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.012833333333333334`], AbsoluteThickness[1.6], 
         PointBox[{{1., 0.}, {2., 0.0004507124280859233}, {3., 
          0.0020614644563615447`}, {4., 0.005743960864461828}, {5., 
          0.009629047979380245}, {6., 0.01002207102066699}, {7., 
          0.007309740261387199}, {8., 0.002113767033160796}, {
          9., -0.0008547566534250484}, {10., -0.0012283462227195008`}, {
          11., -0.0008106999637826906}}]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.012833333333333334`], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.012833333333333334`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{43.80631249999999, 6.509312500000021}, {
       6.165799472760057, 6.632386334297166}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 11.}, {-0.0012283462227195008`, 0.01002207102066699}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {193.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.009166666666666668], AbsoluteThickness[1.6], 
         PointBox[CompressedData["
1:eJxV1GlMVFcYBuArGNRqRVwSf5AKiJU2iDhAtSxyRBARijriMhYbg0qwlOIe
tOrcaFBxHLEk6rDYsqRGiStxQVrpLY0oSx2D4FDbihYUFVFGWxQV9fuY9/zw
JjeTJ+95T+6558z1TErXr3BSFGUF3fzruLrClfcvoSi/XPnk7vKpDg4ku3gb
nm6Ah5GHu1pbt8EjydHN7kUmeDS5smjT4xzYnWxZXdpqgT3IwjKl9kfYi/x/
nn9LCexNrtsXln4E/ph8K2T1T8dhH6GoVQHrak7Dn5LzLl32Pgv7kk26V53l
sB+5tTfL7yLsT/PdnfjmuQZP4vX6/HXod1hH470bl5ir4QChaDPiztTWwIE0
foku2/AHHCQUYRy38Ydr8GeUJ/YYf26AJ1O/X/jS/k3wFJrfxfIswwZ/Tn17
duzom3Aw9SPN+2//DYdQ3unhZG2BQ8mlT5rT78BhNP+9oU0urfBUytWelUfb
4HByjdU54R4saHzD2xyn+w6rZDE8ufDUA+TTeD+KUwwdyNnOFb3PHyGPIM+N
Stv5GDlZC1i2rl8X8umUh2zrLoNVsuZl1i2wI4/k/ehJegCrZM21zjX1KfIo
6i/Vu92EVbJIHWAOeoZ8Bjk+pHMzrJLV3slNJ2UeTfNZO0ZclTlZWWg8aJP5
TLLXoOX1Miersc7Zx2QeQ/mFmSJD5uzMxPKJMp9F47deXnxdPh9Z+Sq/Ypl8
/ljKL9YdviPXx742KjBerj+O7J5bVSrfD1kpaHtlf4L8C/JY+4n18v2yhwzr
bpPvP56c2LA1Uu4Pe0LJ+Fy5f7PJMfHN7e3I2e0RXwbJ/Z9DTtUH78D5EOxZ
Lwtv/Yvx7BvGqAicL42try5LkedvLrlriFM6zqdg666W5P+JPntxyzH7DfTZ
H2zyXduIvp68Kma6J/4fgn1m5ZYXVvTZu6q/GVCPPrt+b0LaFfTnkY9Uho+/
hD673PDdpCr02Zb/9rn+ij7bZ5VbbAX6CeQUz5rB59BnNx+1RZehzy789reH
+N5o7NtFaaPk92g+Wfj3vC5Gn6193/3RIfT7nKdk7EefrRZsrzWjv4AdGJqD
76fo84c2jzXos5U90b7T0O9zznl9cqWjv5CtJY3Z4LBgqwc8BhodVtkitHHE
doc1tuae+U8m+os4D/769U702WpR3a4s9NlKss6yG/0+Xz8cZkLfwPPZOrJg
wRZ+502mynewHo68
          "]]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.009166666666666668], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.009166666666666668], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{43.80631249999999, 6.509312500000021}, {
       6.165799472760057, 6.632386334297166}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 101.}, {-2.0523193378085836`*^-7, 
       0.00001892295400351962}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {580.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.004583333333333334], AbsoluteThickness[1.6], 
         PointBox[CompressedData["
1:eJxV23l8TNf7B/ArlLFP1TJFGaqEWqYEYx9BxD6WEqoyNCoIRouOtWOrqG3s
4YuOpYQq87VUpMHEOtaOPTQYe2xfI6Viaf1e9/N8zh8/f7Sv9+vcc5/z3JO5
M/ee51QZNKrH4AhN0x7k1zT9//Iv3Er7//9s+n9qXfe0EBp0uw2F1tJG3bbE
wGa6NNpDff9Lm3Sbz93ZTVfU7e0+M5026/ZvanaArgqXLHCQrqY7NHXsYbo6
xmPPPUpH6vZ4Jxyna+m2rCxxiq6N4xO2naHr6k4r9WuQtmD8+Qufoz/T7Ss3
/DxdX7cj4fwFuoHu8D+tL9FRyO/Vb5fphroDw+pfoRvpjk3cdZVurDunYIts
2qrb2v/UNboJzjdl0A26qe7Omd+H6GaIf6D7Tbq57qiYGrfoFrrHxGu36Za6
TeWuKbfSneHad4e26S7gWXtX7Ia3Js6+x/bWulNLjL3PdvjOvIQctkfrTs7t
84DtcDC220O2t8H5V3d8xHY49U2Hx2xvi/GP7vKE7XBikc//x/Z2urULjqds
h10jrGG2x+iutOwg7YYnPez6jO3tdTcafY12w28tI3PZHovxW/L9xXZ467hl
tNZBt7VI3edsh825x2ito+6Ehgkv2A6/vRrxN9s76S72cD3thn1D279ke2fd
tQc+od3wrstL89jeBdfrgu0V2+Gk+P/RWldcr+/XvGY7bG7e/Q3bu+m2rC74
lu1w7S37ac2u+9Sw8f+IbbD/XqN/eTyc88lL2g8bLenv2L+77mUFbVpL6Q+v
GrSBdsMmS5F8Yj9c6dvRtNYD8epdpW1w2WFtI9gfjqnuo/3wlMSP8rN/T933
oubSNjhjwVvaDbedPLIA+8P7892itV66J1Xp8x77w1mh07Qbrt02piD7w1e7
+Wntc929Dc0LsT9c5Ju9tBt2LrQa2B82Dt9La711J+VvXpj94bLxftoNe6fF
FGF/OPWbM7TWR3egQVxR9oczTtym3XD15qOLsT+cPlcrzv5xultmLqRtcNL1
aiXYHy7/YC/thw337CXZvy/Ge/UBbYOLD442sj888ZqH9sOBuBCt9dM9/Irl
ffaHswdOo93w8ucXaD9cYUlkKfb/Qvfd6Mm0DZ6c/zzthr/OqvkB+8O1Dk+l
tf66exy/Stvg4Q+jSrM/nL+Wh/bDqT8+prUvdf9VomMZ9ofH7U2l3XBgrqEs
+8PuOUNpbYDupmknaRs8t0y9cuwPV9y4mPbD7qRXtBavu99wh4n94a4/B2g3
/G/F+h+yP1zw4ipac+guFTSUF5thv3EcbYMzVtyhHfDGEb0q8Pxw5IIjtBe+
pVkrMh485cgvdAhOvWH+iPEH6o7osZw2w3Url6zE+HB2h2TaAW+4EFGZ8eH6
B6fQXulf5i3thysktjAzPvxywmRaG6R7Qtv9tBl+fuQdbYN3a9FVGB8u/PcM
2g3bNgVoL7y+TPGqjA9Pie5Bh+CK9VJo7SvdB7Ov02b4987VP2Z8uIN7JO2A
d7n20G64aeN81RgfrryvE+2H1xRdTofgvI9v01qC7sGFLJ8wPlx5z2TaBsdF
naIdCTKeCtUZHzavGk574eKzM2g/3CK2RA3Ghy9fddDaYN0FonfRZvj+FEMk
48OrFn1JO+DoKTtpN/wsukhNxocP3xxI++GjX6TTIbjYjg9qMf7XulfcHUGb
4e9eBmgbvDmn2qeMD9dLn0q74bdJN2gvbIhoWZvx4bYTVtMheOWlf2htiO6n
5eLrMD5co1UmbYPLd6lWl/Fhb3Qy7YbnmJ/QXnjd/R71GB/elrKXDsGNVxot
jJ+o+/S1trQZ/mb5eNoGTzq2jXbArUbeod2wYW75zxgfrlalO+2Hs+sk0yF4
zI4DtDYU13P7S9oMx1Sz1Gd8+O/3h9IO+FvXOtoNP+mTTXvh5b6yDRgfXjyt
Ox2CC5+ZS2vDcH0XB2gzfOtsgSjGh8f82Jp2wA8zptBuuMvQDNoLpy54Tfvh
fJYmDRkf7tPRRWvDda+7voc2w9cev6Rt8Isx1kaMD18aN552w1/lptNeeNr9
t7QfLtK7VWPGh8+2mEZrSbo3rjlCm+HersJWxofXn+hCO+COSxfRbrhQ1mXa
Cycv+qgJ48NTD39Fh+C5I7bQ2gjdPec9o83wsI+bNGV8eGi9qbQDvug7Trvh
ZptKNWN8uGyp/rQfvh3+mQ7BxduEaW0kxm9o1pzx4bqtfqBt8J37Z2kHPDlf
pRaMD3eaM4z2wgtn7KH98IrHBVoyPuw52IPWRune8d5a2gxX2/2UtsGlgi1b
MT489/P5tBve1vY67YXbdC1nY3zd7mLbWtEhtNdvnUhrTt11Ci2kjfD4d3tp
s273jJq3aQvaD88o1prjQ/saYyPajvapp+JpB3xiz2zaieOvn91Ju9Fetvx1
2oP20wsM0cwH7b3rN6B9cES+AbQfLv8mmQ7CjyrvokPwnpE36DCcea9IG+Y/
GtdrViPaCLu6DaLN8OKW82kLHOiZTttg2/x7tB2OfFKqLfOHN4xqRTvhC6WT
aDd85c8U2gOfPHSE9sKZJ3NpH3wqt3I75g+/bdyFDsJfrpxAh+AXplQ6LOPf
cZHWvkH8r/PHMH+4UtRntBle/mE8bYEblp1H2+Cnkb/Tdnib/QHtgAfPLdee
+cPPr7aj3XD7ZmNoD9xs6zraC/+39lnaB0///R3th+f1qRvL/OGf3vWnQ7Bn
xxw6DNcflU5r3+ruEvWANsLbNVMH5g9XOh9DW+BBW8bSNrjxzA20Hf560Hna
AW+2RXRk/nBGpc9oN9znTTztga0X59NeuPgv+2gfvHDiY9oPj2pXoRPzhwcU
6kiH4EJ+Fx2GK4zcRGtjdLc2XqKNcK2fC3Rm/vD4TxvQFvjK2oG0Df5fQQ9t
hxP67qcdcvzyx7QTfuEv34X5wwPOx9Ie+L+nx9FeePm2DbQP/mXMOdoPL6uk
dWX+8MstdegQPNr0BR2G1w1OprWxuqMW7qaN8JEVt2gz/GpyyW7MH3Y2b07b
4AvnE2m79I9eSjvg/83KpJ3itU9ot27378eK2Zk/vHRTVdqL4/8cZaV98M3K
XWk/jtf2fEUH4flR4+kQfHf5fDqM/v1D62ltnO6uxffSRthU6Qxtht8vdZu2
wKtzXtI2uOzaYt2ZP5zdtArtgHvsbkg74T9LdqTd8IUOA2gPvGXQN7QXPtD3
B9oHT667kvbDnUK/0kF42chMOgQHsi7QYbhXpRxa+073hpZvaCNcp0mJHswf
Hli8Cm2B96Q3oG3wOlsMbYe9a+JoB1zjwjDaCe8OTaLd8LFD82kPfGOSl/bC
3xfdQfvgWqMO0X546KYLdBA+sPMuHYL3L/6bDsOPYgr1ZP4u3VmnytFGuPQn
kbQZjupmpS3w8Y6xtA1OKhdH2+Hfdg2hHXCbqt/RTnhr/A+0G16ctJT2wCti
N9BeeNBfO2gfj8+k/fD5nX/QQfh04Bodguv88ogOwz/1e6XyH6/78qWCvZg/
PLVyadoMfx5VhbbARcrUpW1w/8ymtB1+Y21PO+BNE3vSTtg8O552w58OHE57
4O/zf0d74ZxR02gfXH3jPNoPF92YQgfhnqPW0yH4bL5tdBju8WUarU3QvXHq
QZU/nDLylMof/vfjSyp/ePPaGyp/eOz9HJU/bHnxTOUPp514rfKH7w7J/znz
l3initIeeFPeB7QXzn1UgfbBX2z8mPbDB2t8Sgfh95316RBc6fsmdBi+3MtG
axN1V3kWQxvh+z270Ga46OSetAWemNSXtsHmqvG0HX60OoF2wJk3htJOePn9
kSp/uN/ub1X+sNbBpfKHR66bpPKHvZlulT88feMMlT9cqnuyyh9ukTlH5Q+/
/Hu+yn+S7qa5C1X+cMTOJSp/uEWT5Sp/+MGMFSp/+N+l/1H5w+OSVqv84e6G
n1T+8IyhXpU/XMKzVuUPP5q4TuUPV66zXuUPr1qv7IeH3lAOyvlDyiE4d6Ny
GN5aX1mbrDttuopnhE0r1XjM8LmxarwWOFRW5WOD7ZNVvnbYvFldDwc8YKW6
Xk64SHd1Pd1w7dPqenvgEyUWqfzhvDILVP7wtltqPv3wv2PUfAfl+FPq7yEE
r3uo/l7CcIkz6u9Jm6Lb9p36ezPCne+pv0cz3LbCKJU/HFV2mMofrntZ/b3b
4ZgB6vPggJN/UZ8XJ2w8oD5Pbjg3RX3ePHB0k/Yqf9jkVZ9XH7zotPo8++H9
merzHoT9k9T9IAT7/lX3izC8o31Flf/3urP7lFb5w/Z6xVT+cKMz6n5lgbdE
vVH3P/jCoFx1/4MvxT1Q9z85f9mQuv/BRVar+6kbXvRI3W89cErEIXX/hzvf
UvdrH3z9R3U/98Put+p+H4QTGq9Q9384o9l8df+H/YWnq/u/G/msVd83Bnhb
/iT1fQCXsDhoExxRo5f6foDP5ajvu0g4bWwz9X0B555W35dWOf6Z+j61wWtu
qu/bWHjR6kLqekq86q/5fR0HN5/4WP2+gJesvk4nwl3nBdXvDTjQ8SDtgntc
2ql+f8Axlp/pZLhA32Xq9wh81T6LToHrmlzq9wk8aHsinQpfKNtX/V6B73Xt
QKfJ+eOaqN8vcEr9mnQANtwwqd8z8P4vDXSWxN/+Uv2+gw+ev0fnwHknL6rf
e3DTlMN0nlzPJjtpbaruE9vW0gbY+HoBbYQ7VphCm+DXpZLU72PYe7MvHQkH
ZrVXv5fhe+81pK3wj/2rqt/PcPO5JelYeMbSt+p5As51PaDjZHyNL6nnCzgp
eJBOhFe1266eN2DL8v/QLjj9yCz1/AGfO/stnQwXyYhXzyNwwZmd6BS4X53G
6vlE4u2uSqfCzkol1PMK7Ep8xeelNLmeS+6o5zd47ro/6AAcWpyunufgh4k/
01lw56oe9XwHb/FPoHPEbQer5z347bZudB7cUmuqnv+m6S7fpBptgJv2K6Ge
B+FVgTCfV02w5jtMm+E7w1PoSHjrP8PV8zKc9nUr2gonpJZSz8/w20N3+bwd
C4fS02g77Jo/h46DS9sG0A7YdNZCJ8JZ0flpJ5yx7CLfF7jgyDOb1PsHOPHh
eDpZ8nncSb2PgPMufUSnyHg3P+X7DC/ca0gmnSo2LqZ9cNzGBDoNPlyrEe2H
c9YUogOST4Er6n0OvOuLLXQWbNg0Ub3fgcfc70znwG0rVVLve+BA56d8X5Qn
4xnjp7XpumOXLaQNcHDnINoo7Wca0CY4714B2gxHvrvI912RcGnTJvW+DC7W
wEVb4ebdO6j3Z/Dbb8vTsbBr1SO+j7PDOScz6Dh4Rv75tAOe2zaeToRNHgvt
hJPu56NdsLfTefX+UfLZt4FOlvO3GqfeR0r/s+3pFLk+Yz5U7yfhQK1HfN+Z
KuN/nkH7ZPzn5tNpcO1jDtoPx52rTwdg94sCdBBO+ewy39dmwbtmbVbvd2U8
eRPpHPjxzK7qfa9c33pV6Dzxy7/4PlmboTv7+lHaAM+9u4I2wnnFRtAmWOtt
o82wJ/MDOhJOtt9X78vhQP7faStcMVu9T7fB/a+q9+2xsOOdeh9vhx93KkrH
waf86n2+Az48QL3vT4SXVJ9NO2W8H6r1Ahe8q5FaT3DDTrdab0iG7/x9Ta1H
wL6Var0iRfJNUusZXslnhFrvSJXrtaoh7ZPxvi1Kp8HV5t1S6zOwt6tabwnI
9WrtoYOwZcgQOgtOPNBSrefIeDuWpXOkvXQ214PCcNqbOXSexOvcnNZm4voY
nnC9yQAH666hjXBeZjfaBFsO5lPrVXCiZRfXuyJhe/EhtAX29S1PW+G4smfU
+hocajONjoUNtxvRdtj5/BHX6+LgVNda2iHnc/ahE6V/dnHaCdsyDnO90AV7
3p+k1g9hV3YDOlnG/9Fjrkd64JQLG+gUuT4RA2gvbE0tR6fK9TlyluuhPtjc
by6dJucf0p72w8YnEXRA4j3az/XYoPT/aiKdJfPR06rWcyXe/hdc782R863b
SYfhtHej6Tw49rKF1n7QnRMZ5nqzAU55vp02wq6mTtoEe15b1Ho1HIrK5Xp3
JOx8sJO2wHbTONoqzrSq9XU53503XH+PhbOmHqDtcOSK6XQcHFcrlnbI+RoW
pxMln73nWA/ghLW0FNoFO+rFq/oBOb5ydTpZ2mc/YT2CR/IfvptOgU1HJ6t6
B+mfEkOnwqm3S9I+OX7jFdZLpMnxN9bTfjjgGUkH5PqlNaGDEr/3e3QWbLm4
TtVzwP6cZnQObL5/kfUgYTnfH046Dw56i9HaLN3efptZb2KAg29iaKO0T7/L
ehUTbHs2gzbDvjaf0JGwecJR1sNYYM/iRNoKGxcUpW2we/h21tvEwpbIXrQd
Dh96xfqdONjfwks75PxL2tOJsOPYU9YHOWU8F1Nol4w3I1rVD8n5Jz1h/VGy
5F96Be2RfKe3o1Pg0NFcVe8kvu6lU+X8x7vRPsnvh3esn0qDtdI+2i8ePZAO
wPZVH9BBuf5LjrJeK0vGFzdB1XNJ+/26dI60t7rD+rCwXE/HSjpPzt++O60l
o/25gTbA9sGZrEczwp4VE2iTeEGUqleDzTFPWf8WCWsZW2gL7H/6NW2V9tDH
tA32zrnJertY2PbsJ9ou7SXj6TjYcbcS7ZDjb+1nvV8i7P6qNe2U4zsdYf2g
S7y8I+2W/u3Osv4wWcbXoy/tkfGn32L9YorkO3kE7ZV4y16x/jFV2gvMon3S
fqgMnSb5nP2Z9ZR+Ob52Yzog4wkeZ31mUPz7l3SWnO9hLus7Q3K+L2bTOXCo
aBU6LP3z9rJeNE/6V+ql6kdn4/iJT1l/aoC9xefSRthxtiZtgm2HA6xvNcP+
24l0pLTXKUJbxGu2sl7WCmtRdtom8R89Z/1trPjQStou8dNa03Gw+9gD1vc6
5HxPF9GJMv5PW9BO6T8xh/XDLvGfS1R9MWzuEE0ny3gPhlmf7BG39dIpMr5T
dtor5+sbQafK9Tg9lPXQPhlfTpD102ky/tdNVD21tOf7mfXXAen/8n06KM5y
s347S/zTM9Z3h8TtE+gc8cks1oeHxTW60Xni3sdYX679CPdsTRtgd/l9rE83
wraNTWkTrP2l6tvN0v6mKR0pTt/H+niLnK92NG2V9i4BVW8vNtnpWHHyFdbn
22V8KxLoOGmPfcb6foeMZ8H3dKLEG1aSdooverl/wCXHn2qg9heIOwa4/yBZ
3GkA7ZF4J19wP0OKnO/ofNor7daadKq0N/mR+yN84puHuX8iTfxlhNpPIfFS
o7n/IiBOm8n9GUHxzJPcz5El/Q2l6ZC0N43n/o8caS/1K/eHhMWz/+F+kjw5
foVd7S+ZA1s3cv+JQRz/L/erGGF3ob60SVzzN+53Mcvx+8vQkdK+/zvuj7GI
q2ZzP41Vjr8bTduk/Z9fuB8nVtoTytF2cemZ3L8TJ8cXeMH9Pg6xZQidKMdX
3sP9Qk5xsw3cT+QSr1uq9huJO8/hfqRkcY2Z3K/kEdeczv1MKeJ2M7nfySse
MYf7oVLFy5Zyv5RPvHs991OliTN2c7+VX7zqJPdjBcRd73K/VlB8OILOEj/9
mPu9QuLjsdwPliNuM5r7xcLijU25nyxPfHWF2m82Fy6icT+aQVxxBPerGcUR
N7M4/+J1X3C/m1l8L/si51+clsD9chZxgb+4v84q3jXrLOefx1flfrxY8Z1Z
pzn/4q2XTnD+xSmfBTj/4hkpRzj/YlvRQ5x/sW+2n/Mvjq2RwfkXh7vv4fyL
TYt3cP7Z/8IvnH+xtn0d558eupjzT8/vz/kX+99bsA9OEwdbeMV+sXNhqjgg
tg/wiYPivPhd4iyxoWuaOMTrF5khzmH+4QPisNh1LlOcx+u1/pBYmwcnDD8i
Noijah0TG8WlQwGxSVx73gmxWbzqs1PiSPGY06fFFnG24w+xVfx6bFBsEyc2
PiuOFZ96Tts5nm3nxHHirV+dFzvEbctcECeKCx6inWLDiItil7jHB5fEbnF4
N50sft7rstgjHhamU8T9k7PEXuZX8Yo4VZz7K+3jeJtd5fyLKx6l/eLYLn9y
/sXNz9JBcYEe2Zx/8a4/6JDY3vEa5198J5MOiyc0us75F1dPpbX5cE7ZG5x/
8aVptFH89yPaJK5RM8T5F7/uR0eKv5lDW8Sr0mmreFwObRMXLXOT8y/ub6Pt
YucwOk7ceTHtEOem04ni4Tdpp/j3grc4/+I7tWi3+HEXOlmcNYr2iLd66BTx
KB/tFZv/oFPFxx/TPvGwwrc5/+Iin9B+8a+t6ID48750UFz4GzpLfGI2HRKv
9tI54lm/0WHxjyfpPOZzg9YWwOFc2iDu994dzr/4UVnaJN5SgzaL/9OYjhTv
i6Et4vKf01bxtkG0TTx9FB0rXjHxzr7/A82Omxw=
          "]]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.004583333333333334], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.004583333333333334], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{43.80631249999999, 6.509312500000021}, {
       6.165799472760057, 6.632386334297166}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 1001.}, {-1.4903200895588498`*^-10, 
       2.0153955171409166`*^-8}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.05], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {967.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}]}, {}},
  ContentSelectable->True,
  ImageSize->{
    UpTo[600], 
    UpTo[360]},
  PlotRangePadding->{6, 5}]], "Output",
 CellChangeTimes->{3.8127014063061647`*^9, 3.813054193884497*^9},
 CellLabel->"Out[72]=",ExpressionUUID->"8afb6170-cee0-4697-ade9-386237506070"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{"exactSolution", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListPlot", "[", "solRK410", "]"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"p2", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Plot", "[", 
      RowBox[{
       RowBox[{"exactSolution", "[", "t", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListPlot", "[", "solRK4100", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"p3", "=", 
    RowBox[{"Show", "[", 
     RowBox[{
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"exactSolution", "[", "t", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", 
        RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"ListPlot", "[", "solRK41000", "]"}]}], "]"}]}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{"GraphicsRow", "[", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK410", ",", "exactSolution"}], "]"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK4100", ",", "exactSolution"}], "]"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"absError", "[", 
      RowBox[{"solRK41000", ",", "exactSolution"}], "]"}], "]"}]}], "}"}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.81305285826786*^9, 3.813052884308997*^9}, 
   3.81305291494473*^9, {3.813052968816243*^9, 3.8130529733224645`*^9}, {
   3.813053020203862*^9, 3.8130530609129214`*^9}, {3.8130536551368885`*^9, 
   3.8130536679536676`*^9}, {3.8130542651316147`*^9, 3.8130542896238565`*^9}},
 CellLabel->"In[81]:=",ExpressionUUID->"f0a62991-e83a-4f24-8ae5-33f6ddfd941b"],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.012833333333333334`], AbsoluteThickness[1.6], 
         PointBox[{{1., 0.}, {2., 0.00008930518423820588}, {3., 
          0.0004215034563744835}, {4., 0.001253408785784449}, {5., 
          0.0023924796971642737`}, {6., 0.002667709008317698}, {7., 
          0.002031192758975231}, {8., 0.0015552019939377715`}, {9., 
          0.001071970252280674}, {10., 0.0006212055422888385}, {11., 
          0.0003186279870531772}}]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.012833333333333334`], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.012833333333333334`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{42.78545833333334, 6.487458333333393}, {
       13.159202025497564`, 0.4999999999999858}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 11.}, {0, 0.002667709008317698}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.02], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {193.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.009166666666666668], AbsoluteThickness[1.6], 
         PointBox[CompressedData["
1:eJxV1GtMVEcYBuBjtdWgEowVMUHBYg2xTWUBFbzAEcFW0JZCAUlIWUXAGIRV
WRQby6SgtRRkK3KzgJs27T8SYLFeomFNpH+8BmNrRM0KJksb0DXGBFRK5x3e
+eFJTjZP3vn2fGfOzCzdWZKW/45hGPnyxu/U5Ys33r5MeZ8XabPjpjgLfq9q
9hI6QFoM3bgaTr+PvCVzWSQdJO2eDPTG0MHSzvigApMOxfiCdW8S6Q8w3r0g
OZleJm091heZSi9H/a4wVzodDt8f6M6iV8B/VMTn0B/Dd3pcVvoTaXOgvzWP
jkDu5x0soC3Sqd2jtXvoSDw/88VvxXSUdMCtspZ9dDTGX0haUkqvQv5w4kkZ
vVraYbnrK6fXSNuelSceoWPQz+gLm6BjMZ+e4Pzv6LXS90IvfVtFr4M7oi4f
o9ej/7wVIT/QG+A/S3/5kY7D84aMxBN0PP6/eKjCQZvo15c7/NOUBXy8fv/2
euYbpcNTWq+fYg5bA0fMRuYJ6P+of1cTc+XQvKAW5pvwvJKF+08zh23jmy/8
zDxR2vNRyD+tzOHbl96daGOehPm92utrZw47q5vcZ5hvRj/23mgnczhiUWMp
bXwq3dm38IzOYevZSZfOP8PzDx/RFrAx3dKu8y3Sjraldp3DqW2dsTpPxvv6
zRtgPwL2/DX/Q91fCvr9uzhD96/s22TT77cV/dg6yvX7w87+druen23o58q1
Qj1/cOhwQ7qe38+RN01b38wcttZ0hvF7GF9g/MzBOfxeAvY8Lxvn9zRSUf9y
7F9+bxM2M38dPMnxsPP3osdcH27YWhUwrNfPl9LinOVpHeths+36Yq43oWwp
ya1hvXJHrqua9Wmorz23iOvXVA4rOvk962GzMmMx17sbNuZHX9T7IR2e6Nld
yXpYXLGHc/8I2AyZOc795Vbjxx65Klj/FXwju4H70YRF196Gb1ivXD7Zzf3r
VuPrR0YOsj4DPjQ3gfvdVPau6eJ5IGDhNWMPsF7ljf/d0+dHJnxtS52N9cpz
J8/yvBHKB2aUFbFeudeeo8+nLPhmVOFu1iu/XnmK55lQztn7eBfrlXumbdPn
33Z4deDdHaxXfuU9yPNSKB93FH7NeuWEWH99vmbDzZc92axXHot5kBX3PwnC
VXw=
          "]]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.009166666666666668], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.009166666666666668], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{42.78545833333334, 6.487458333333393}, {
       13.159202025497564`, 0.4999999999999858}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 101.}, {0, 4.977917835713441*^-7}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.02], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {580.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}], InsetBox[
     GraphicsBox[{{}, {
        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
         0.004583333333333334], AbsoluteThickness[1.6], 
         PointBox[CompressedData["
1:eJxd232UE+XZx/FREGNFGxAlUsQBKSdVSiNuMQrolK4SLEJUkMibEQFXu0JE
XqJSGBExAmpEpAtsNSLCQqmmsMKKFKeINeJL04oaXxsp1WARU0UaONg+55nf
9/7jefhDz+dcc93Xfc3OSTJz39Nz0vRrp5xoWdaBdpb1v//Xv+rl1v/95/j/
3XzCYDHg2/2oEw7Kp/0Ad5HreuGQb+doH9zdt/fOBdhWfHIE91L+gf64t+rf
81PcR8eHojis+K5L8fm+7fRg3Nd36hwH91P8/p/hiO/4N0PwhYo31OP+ms/f
r8AX+c5PGYrr5EMx/FPVm3cVHuA7eeZwfLHqbbkaRxVPjMSXKH9BHF/qOzH+
GjzQd+2Sa/EgHf+D6/Bg323WKHyZ7+IB48t9V98djR3fgcL1sivbfxhD/Geq
vzVBXG7ZcgPxIb7Dz48lLte2jyP+c9/WK+OJy+m3JxCv9+18NpG47H53I/Er
fK9dmiQu7+1+E/ErfY9/DrtypH4S8aG+J3+IXTkw62biMd+hTpOJy7nnsDXM
d1N8CnE5fBhbV/mOrZpKXK4OuYX4L9TfIezKkeYG4sN9l4ffSlyOWrcRv9p3
fCt25cj0XxIf4bt0fiNxOX4AWyN9extvJy4Hp08jHvedGDBdduSYleJ4+eBU
7Mld/oyta3zvjt5BvhxYi125FJxBvhybj61rfY/6CjuyddOd5BPfiz05MWwm
+depnocduTE6i3x5+WbsyTP7zSZ/lPrZhB05c8Ec8uXCJuzJpZ+kyR/t+41W
7MjLB95Fvuzsxp5cGnE3+dernw+xIxdvvYd82TmGPblt6Vzyx/i2e/6KfDm3
Dbty3/g88uXdX2Ar4dt9YD75ct8zXPLl/M3Yk4c8j60bfH998r3kyzvHYVdu
yWNPzndYQP5Y359MwI48YCt25a3fv498eept2Brnu/5P2JGvP28h+fIT92JP
7vUptsb73vez+8mX9z+NXTl88iLy5U2/xNYE3/P+gh35iegD5Msdn8KeXPxe
hvyJvg/Pwo7csA+78qD4g+TLqZewdaNvK7KYfPn4U9iVJ5+5hHy5bjG2kr7n
WktlW65LY0du+BdOyl0aH2J8ub6Cc3J1ysPUk7v/A5flj6Y8Qv2bfPetYFu2
BmSpL29P46Qc3YFdeZH1KPXldVdgT16xBJflsX/F1iTN9+xl1JcnTsKOvO63
OCnv+ha78k7nMerLy5diT77qfVyW9/dZTv2bfY+fhW15x27syCee+Tj15f5T
sSvH2nBOHn7qCurLl96Iy3K3VmxN9n3olF9TX96axI7c2IaTco9OTdSXd9+G
c/LUV7Ann9hzJfXl5l9ha4rvAR9iW957ySrqyzNW4qTc5Rh25e3jVlNfTu3E
nhzu1Ux9+etF2Jqq+R/Etrzuut9QX171Ik4S7/0E9eVdD+Oc/PVR7MmXTn2S
+nLz29i6xXfHnjnqy6UEduSLHsVJ+Zw92JVXtnuK+vLawdiTh6RxWZ69BVsN
vgd9hW15xQVrqC8vasBJueM67Mrh/Tgnf9DraerLvSfhstxhDbZu1fj7sC03
n7eW+nJ8Ck7KufXYJf8LnJOD/Z6hvhyZgctydSu2bvN97XFsy6OGrKO+XM3g
pNy/iF05GFpPfXlZEnvy1g24LM/9Blu/9L1vcAv15UMZ7MjNe3GSuL2B+vIn
t+OcPO9F7Mk7T9lIfXltAluNvvu2YFueWMOOPGDYb6kv51dhV373IM7JzZdv
or7c5TFclvt8jq3bfX828HfUl696FDvyxM9xUj79smepLzc+jnNy6kvsyfaV
z1FfnvYktqbJR7Etn3x+nvrypjE4Ka+5H7vyV1twTl6yD3vy/E6/p778Rwdb
030nUtiWnRx2ZLeIk3LHEzZTX/7nhTgn97gZe/LTy3FZTv8JWynfv6nhoHzW
BVuYj/zBBByRD2exI9+wG8flTjWclM/o2yqn5FuT2JUDj+OsfPQ1nJMH/Bfn
Za/uefqTl92Gi3JrDpflHu/hqvzuaVvp/w65Hgdley625fwWHJGX/hM78sbz
ttG/3GU8Tsq7luOU3PomduUvOrTRvzzVwTm5z904L/duxR7HH8JFuRJ+gf6Z
3824ip/A1gzVfx8H5eSZ2+lfPusaHJHbP4Qdud9rOC4vPelF+pfDQ3BKPjwP
u/LxF3FWjh7FOXntgB30LzszsSd33oyL8ulVXJYv7fcH+pebGrF1p+/uv8VB
+Y0D2JbXhXfSv7zpFuzIb63DcbnHZzgpL/vhS/Qv956CXfndtTjL+PtxTh7R
2aN/+fr+2JNvvRYX5SUzcFnetgxX5a82Y2um735v46B89zfYlt8844/0L/+o
Djvyg6NwXP56Jk7KEx7HKcZ7Hrty/bs4K79wBOfk/l130b/ccjH25F4JXJRX
pnFZ/v5KXJXnvYCtWb6/fB8H5euPYVve3u1l+pfPGIgd+ZZxOC7n78FJ+evV
OCX33oFdOf4RzsqNx3FOvrv7bvqX5wzCnjx1PC7K9XNxWe7cjKty4UVszVa9
D3FQPn4M2/KMbq/Qv/zuJdiR7RtwXB6Rxkk5+Wucksdvxa7c/x2clavf4Jz8
cOc/0b8cvBB7zH8kLsott+OyvHMJrsqtG7A1x7f7Kg7Kff6BbXntia/Sv3z8
XOxw/GAcl/uOxUnZmoNT8sbHsCv3zuOsPO0NnJMXV3BentG+QP/Us3FR3jQQ
l6k/Blfl0AxspX0fW4qD8tr12JZDu3BEjn+EHXnEERyXOwdfo3951Y9wSi4P
wa78xTiclVtn4px82UM4Ly98Bnty5g+4KMfewWV590FclY+320P/d/n+uhsO
ys0XYls+NhRH5M4TsSMfvBPH5RkP4qS88Qmckpdvwa48qICz8sMf4Zy8rIrz
8pD2r9O/vLkrLsq7z8dledFgXJWPjMTW3b67TcJBed+d2JYT9+OIPHsFduQr
1+O43LYNJ+W9r+KUvOo97Mqnf46zcq8jOCd/0v4N+pedM7Anx3rionziT3BZ
jg/CVY4fhq171O9oHJTPn4Rt+fRpOCKvuAs78s6FOC5nH8FJ+dhKnJIDa7Er
b/8dzspnbcM5ub1n+pebCqZ/uVA0/ctrSqZ/uXvZ9C/3+9z0P1f9f2n6l53D
pn+57pjpX97zX9O/fLz9m/Qvv3UKTsqx03FKntgZu3KXs3BWbjgb5+TJ3XFe
PtIDe3LfnrjIfHrhspzujavy0h9i61e+6/vgoNxsbP8/R+Q6Y0dOmfHi8nBT
LymXzHxSsmXm68oF009Wjph+c4xnzkdePn6m6V8eZc5nEZvzXZZr5u9RZbyT
TP/zfA+yTP9y0fy9bTn4rfn7y8cPmb+/nK2Y619+41Nz/cv5D8z1L9e/ba5/
2X3dXP/y5JfN9S/XtpvrX3Y2m+tfDm8w17/sPWmuf7njCnP946Xm+p/vu/Ve
c/3LHeeY/uXujaZ/uZw0/cujzOdHXJ5rPl+S8vjBpn/5YMT0Lw/qbfqXo11N
/3LlFNM/9Y6bz3/ZPWQ+/+WGsvn8lwN/NZ//8uSXzee/q/xWHJBHPWO+D+TK
4zgkD1pkvh/kxGwclqNTzfeFXBiNo3LvK8z3hxyvwzE5cp75PpFLnXBCHm6Z
7xc5e4jvxwZ5qfn+TDG/PTgtV7eZ71/qP4MzcmaZ+T7m/MzHTXKs0Xw/y+UE
biF+hfm+luea3xNtjH8O9uSGU3BBDn1rfs/IzWVckmuvm983nM9tuML5W2N+
78gdze+nmuzNwda9vkdNwgF593AclAMX45Ac6Yltue+pOCx3+db8XpSLn+Co
3Fgwvx/lyu9xTI6tNr8n5exCnJB33I6TcuF63ED8cpySl4dxWo53wq58+Ci/
rzNyep/5PS5X9+AmztcW8/tcbl6NW+S995nf67LViNvkjqOwx/kdhAvM5zxc
ZLxTcUnOfWPub+TxH+IK9V429ztyy0ZckwctM/c/C3S+7sIBue4mcz8kt8Rw
SO4SMfdHcmNXHGa8/5j7Rbn9Zzgqx98094/y0lYcI3+1uZ+UKwtwQg7eZu4v
5eg1uEEeFTX3m3LqXJyWsx3M/Sf9fcn9bkb29uKsXDL3y03ywTU4J1uLcYsc
usPcr8uRBG6Th19u7t/lxj64IGdOM/fzzO8wzwNKcuFDXJYP78IVufdGXOV8
PIprjJ/G1n0a70YckDsOxUE50Q+H5E1nYltu/x3PS8Jycj+OMP7rOCrXbcGO
nF+FY3J0AY7Lu2/FCeZzDU7KgUtwg9xq45TcGMBpOVw1z5/k4+/hjLz3JfM8
St6xHjcx30fM8ynOxxzcwnxvxHm5PBS3cb4i2JMHhXBBXmrhovxRhedxJbn+
L+b5HvN7AVeIrzHP++hnMa7JM+80z/8W6nyMxwH5YL15Hih7P8Yhee1Z5vmg
HP0vzyvDcuwvOCJnnsZRuTQLOxwfwzG50g3H5ZYveZ6akLOeeb4q5x7DDXJh
Kk7R36U4LTedhl3Z+ZTnvRk5+DzOyrUMbiI+AeeY/4W4hfmchPNy5AOeR7fR
37Pm+TXzvQ8XcMI8z6b+j3FJTrQzz7c5nyWel1fk9LO4KscX4hoei637fTdE
cEBu64CDsvMxz/dDcmmLWR+Q84txWM7dZNYLGC+Ko3IgaNYPZPdz1h9isv0S
jsuVFThB/Wk4KVtDcYMcPxenZO/frI+k6beIXfrbgDNyZIFZj5Fj43AT/dWZ
9Rk5dDpuod/PWd/Jy9k/4ja8GnvMbxYuyME4LsqZ83GJ+Z6Ey5zPMutRFfwi
rsrRFbhG/TvM+tci3/bVOCAXwmY9TG5pj0O4zHqaLZd24LAcWYkji+gPR+XU
tdiREz/BMTnZEcflpgOs/yVk61WclLNrcYMcW2DWG5mPWW9My85lZv1RznTH
Gbl2jPXLLPH3cRPjt+Ec4/8atzD+bJyX3dG4Ta7WYY94F1xg/MOsxxYZfy8u
EW8167nkm/XdCuPPxFXmPxrXmN8As/77gMbvigNyrMZ6c1B238chubrdrFfL
6WYcZrx5OCKHkjgqh4eY9W053hvH5HwHHJejB1hPT1D/dbMeLxeexQ1y8VGc
kq1ZOC0nE2b9Xi4NxBn6ORdnmV873CRn32K/QE4ur8AtHH8jzstOGLfJqX+Z
/QnmfOKC7C3ERcYfgUuMdzYuE/87+yMq+Flc5fi7cY3zcwW2MnInHJCtT9if
EZRzG3FITs3BNvF6HJYjnXFErv6N/SFRufw77Mj2XByT81fhuOyejRPUq7A/
JSkH23AD4z+AU9Qfg9OyE8Yux9fYL5PhfLyGs+Svwk1yvNHst6H+ZbiFfoI4
z/nax/6dNsZrxZ6ZLy4QH4uL1OuHS8RPwGXy32E/UYX6G3FVzs7DNca7zuw/
elAO44Cc+479TUE58jYOydYGbMv2fByW3dE4QrwvjpLfzuyvkoMfsP8qJmd/
j+PkP4gT5N+Ek8z/EtzA/DvhFHGzPyxNPy+b/WPEm3FGzs/CWdkZiZvwj8x+
MznVDreQ/zH72fLE23Ab+Y+Z/XHUn4YL5A/DReI/xCXOxwlmPx3jfcz+uwrn
/wVc5XytwDXZm4Gtxao3Egfkcl8clJ3v4ZBsVcx+Qbn6Cg7jp3FEji/AUbmY
xI6cuxzH5GwPHJe979jvmJAjH5v9kIy3AzfgZpyin7k4LdsTzP5J6g/GGebb
A2dl9z/s12xiPA/niN+HW/CVOG/OH26j/ltmfyj1l+ECx4/BRdwdl6j/KftR
y3g9rnD8NFwlXodr1D9u9r8u8e29jANybgkOytYoHCLe3eyXxf9gv21Ydp7D
EbmcxlHiP8cO8dNwjPmU2N8bl+21OEF8utkPTP5A3MB8AziF32F/cRqvwS7z
n44zsjsYZ6l3Km6Skx+Y/c7MbwNuob80zjO/obgNdzX7o5lPhf3VBeqZ/dZF
/CAu4bG4zHz74grz+w/7u6v0V8Q15ve02f+9VJ6NA3JyGA7K1jk4JJer7De3
ib+Cw7K7EkcYbxqO4p9jR7bPxjHZO8T+9zg2++ETHL8aJ5nPHbiB42M4xfE2
ThP/N/vxXeb7Z5whvh5n6c/FTcQTZr8/9S/ELfh7OM/4+3lfoI3+d2IPN+EC
48/ARTl3NS4xXhiXcTtcYb5/4/2GKvnbcY3xV5j3Hx7S8XfigOyOxEHcF4fw
KdiWvc94HyMsJ3fjiOyswVHsmvc7yJ+IY+QPxnHZ7o4T5B/hfZEk+a24gX5m
4hTzvQinOf4b3kdxibfiDPFZOEu9i3ETxx8177tQbwduwfNxnvGG4DbG64A9
xtvD+zUF8h/BRfJH4RL53XAZl3mfp8Lx63AV345rjP9T8/7Pw8r/jveFArL7
Cg7K1sM4RHwMtmXPxmGO/yfvJ0Xw8zhKPRc7jPcLHOP4rjjO+OZ9qATx53AS
z8UNjD8Mpxj/LJxmvP3m/Szyt+AMXoCzjHcNbmI8G+cYr8r7Xy3Yw3n8KG7D
k8z7ZPgiXMAn4SLzeY/30UrU34DLzG8uruCRuEr+ebjG+EfM+2+PKP46DuAn
cVB2ZuIQHoZt2e2Bw/gw7+NFOH4PjjL+U+Z9PeJzcAyPuG/w/wBgm8nV
          "]]}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.004583333333333334], AbsoluteThickness[1.6]}, {}}, {
         {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
          0.004583333333333334], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
      AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
      Axes->{True, True},
      AxesLabel->{None, None},
      AxesOrigin->{0., 0},
      DisplayFunction->Identity,
      Frame->{{False, False}, {False, False}},
      FrameLabel->{{None, None}, {None, None}},
      FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
      GridLines->{None, None},
      GridLinesStyle->Directive[
        GrayLevel[0.5, 0.4]],
      ImagePadding->{{42.78545833333334, 6.487458333333393}, {
       13.159202025497564`, 0.4999999999999858}},
      Method->{
       "OptimizePlotMarkers" -> True, 
        "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), "CopiedValueFunction" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& )}},
      PlotRange->{{0., 1001.}, {0, 5.3025028812214714`*^-11}},
      PlotRangeClipping->True,
      PlotRangePadding->{{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.02], 
         Scaled[0.05]}},
      Ticks->{Automatic, Automatic}], {967.5, -116.80842387373012}, 
     ImageScaled[{0.5, 0.5}], {360., 222.49223594996212}]}, {}},
  ContentSelectable->True,
  ImageSize->{
    UpTo[600], 
    UpTo[360]},
  PlotRangePadding->{6, 5}]], "Output",
 CellChangeTimes->{{3.81305426624684*^9, 3.813054291955165*^9}},
 CellLabel->"Out[84]=",ExpressionUUID->"b28bdd49-4536-445e-a4fa-ee4c720a28f2"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1366, 685},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>1.4 Inherited,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 158, 3, 92, "Section",ExpressionUUID->"fd76de14-de24-4a76-939d-182b44b1f370"],
Cell[741, 27, 23468, 389, 371, "Text",ExpressionUUID->"87dcf015-e835-4d58-a9f0-e1cc62932ea9"],
Cell[24212, 418, 3362, 95, 551, "Input",ExpressionUUID->"1b56c63c-48a1-429d-a6c3-7377e0ae88f0"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27611, 518, 154, 3, 92, "Section",ExpressionUUID->"cd97b77b-6eda-4198-be3b-ddb0d0e4e616"],
Cell[27768, 523, 12942, 217, 455, "Text",ExpressionUUID->"f084202f-33da-4926-824d-86be44fbdf50"],
Cell[40713, 742, 3862, 105, 497, "Input",ExpressionUUID->"274b5909-0cb0-4260-b141-d58a5d1a171f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44612, 852, 207, 4, 92, "Section",ExpressionUUID->"995c2255-176b-4e1a-aab4-2ba8523595de"],
Cell[44822, 858, 2776, 52, 629, "Text",ExpressionUUID->"eb5c5bca-f5fa-4910-a0f0-d971673b8378"],
Cell[47601, 912, 349, 9, 40, "Input",ExpressionUUID->"f3597dcb-7738-47d2-82ab-7769297cd647"],
Cell[47953, 923, 428, 12, 68, "Input",ExpressionUUID->"8464c760-38b2-4e23-8e97-ed2c0348e6ae"],
Cell[CellGroupData[{
Cell[48406, 939, 2117, 59, 293, "Input",ExpressionUUID->"3a6a925f-56b8-4a77-a813-c11c80d05961"],
Cell[50526, 1000, 37372, 665, 205, "Output",ExpressionUUID->"65d4e4ef-a5ef-461b-96ad-2d035da53501"]
}, Open  ]],
Cell[CellGroupData[{
Cell[87935, 1670, 1204, 30, 125, "Input",ExpressionUUID->"60afd7b7-e7bc-4158-80c0-baf5b2b32e9e"],
Cell[89142, 1702, 15618, 300, 205, "Output",ExpressionUUID->"8afb6170-cee0-4697-ade9-386237506070"]
}, Open  ]],
Cell[CellGroupData[{
Cell[104797, 2007, 2224, 59, 293, "Input",ExpressionUUID->"f0a62991-e83a-4f24-8ae5-33f6ddfd941b"],
Cell[107024, 2068, 13655, 267, 205, "Output",ExpressionUUID->"b28bdd49-4536-445e-a4fa-ee4c720a28f2"]
}, Open  ]]
}, Open  ]]
}
]
*)