(* Content-type: application/mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 6.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[
NotebookDataLength[
NotebookOptionsPosition[
NotebookOutlinePosition[
CellTagsIndexPosition[
WindowFrame->Normal
ContainsDynamic->True *)
(* Beginning of Notebook Content *)
Notebook[{
Cell["445: Solution of Differential Equations", "Title",
Evaluatable->False,
CellHorizontalScrolling->False,
CellChangeTimes->{{3.2*^9, 3.2*^9}},
TextAlignment->Center],
Cell["David Gurarie", "Subtitle",
Evaluatable->False,
CellHorizontalScrolling->False,
TextAlignment->Center],
Cell[TextData[{
"There are two basic commands that solve Differential Equations in \
Mathematica. One of them, called ",
ButtonBox["DSolve",
BaseStyle->"MainBookLink"],
" does it symbolically, another one ",
ButtonBox["NDSolve",
BaseStyle->"MainBookLink"],
" does it numerically."
}], "Text",
CellFrame->True,
Evaluatable->False,
FontSize->14,
Background->GrayLevel[0.849989]],
Cell[BoxData[
RowBox[{"Off", "[",
RowBox[{"General", "::", "\"\\""}], "]"}], ";"}]], "Input",
InitializationCell->True,
CellChangeTimes->{{3.626*^9, 3.626*^9}}],
Cell[CellGroupData[{
Cell[TextData[{
"Solvers in ",
StyleBox["Mathematica",
FontSlant->"Italic"]
}], "Section",
Evaluatable->False],
Cell[TextData[{
"Output of any solver is replacement rule (",
ButtonBox["Replace",
BaseStyle->"RefGuideLink"],
}], "Text"],
Cell[BoxData[
RowBox[{"Solve", "[",
SuperscriptBox["x", "2"], "+",
RowBox[{"a", " ", "x"}], "+", "b"}], "\[Equal]", "0"}], ",", "x"}],
"]"}]], "Input"],
Cell[BoxData[
RowBox[{"Solve", "[",
SuperscriptBox["x", "3"], "+",
RowBox[{"a", " ", "x"}], "+", "b"}], "\[Equal]", "0"}], ",", "x"}],
"]"}]], "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell["DSolve (Analytic solutions)", "Section",
Evaluatable->False],
Cell[CellGroupData[{
Cell["Logistic equation and IVP-solution", "Subsection"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"General", " ", "Solution"}], "*)"}], " ", "\[IndentingNewLine]",
RowBox[{"Clear", "[",
RowBox[{"a", ",", "n"}], "]"}], "\[IndentingNewLine]",
RowBox[{"gso", "=",
RowBox[{"DSolve", "[",
RowBox[{"{",
SuperscriptBox["y", "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{"a", " ",
RowBox[{"y", "[", "t", "]"}], " ",
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"y", "[", "t", "]"}], "/", "n"}]}], ")"}]}]}], "}"}],
RowBox[{"y", "[", "t", "]"}], ",", "t"}], "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.6252`*^9, 3.6255`*^9}, {
3.626*^9, 3.6254`*^9}, 3.6257`*^9,
3.6256`*^9, {3.6253`*^9, 3.626*^9},
3.6255`*^9}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"y", "[", "t", "]"}], "\[Rule]",
FractionBox[
SuperscriptBox["\[ExponentialE]",
RowBox[{"a", " ", "t"}], "+",
RowBox[{"n", " ",
RowBox[{"C", "[", "1", "]"}]}]}]], " ", "n"}],
RowBox[{"-", "1"}], "+",
SuperscriptBox["\[ExponentialE]",
RowBox[{"a", " ", "t"}], "+",
RowBox[{"n", " ",
RowBox[{"C", "[", "1", "]"}]}]}]]}]]}], "}"}]], "Output",
CellChangeTimes->{
3.6254`*^9, 3.626*^9, {3.6257`*^9,
3.6255`*^9}, 3.626*^9, {3.626*^9,
3.6255`*^9}, 3.6255`*^9}]
}, Closed]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"IVP", " ", "Solution", " ", "with", " ",
RowBox[{"y", "[", "0", "]"}]}], "=", "b"}], "*)"}], " ",
"\[IndentingNewLine]",
RowBox[{"Clear", "[",
RowBox[{"a", ",", "n", ",", "b"}], "]"}], "\[IndentingNewLine]",
RowBox[{"iso", "=",
RowBox[{"DSolve", "[",
RowBox[{"{",
SuperscriptBox["y", "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{"a", " ",
RowBox[{"y", "[", "t", "]"}], " ",
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"y", "[", "t", "]"}], "/", "n"}]}], ")"}]}]}], ",",
RowBox[{"y", "[", "0", "]"}], "\[Equal]", "b"}]}], "}"}], ",",
RowBox[{"y", "[", "t", "]"}], ",", "t"}], "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.6252`*^9, 3.6255`*^9}, {
3.626*^9, 3.6254`*^9}, 3.6257`*^9,
3.6256`*^9, {3.6253`*^9, 3.6253`*^9}, {
3.6256`*^9, 3.6256`*^9}}],
Cell[BoxData[
RowBox[{"Solve", "::", "\"}],
":", " "}], "\\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", ButtonNote -> \
\\\"Solve::ifun\\\"]\\)\"\>"}]], "Message", "MSG",
CellChangeTimes->{{3.6256`*^9, 3.626*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"y", "[", "t", "]"}], "\[Rule]",
FractionBox[
RowBox[{"b", " ",
SuperscriptBox["\[ExponentialE]",
RowBox[{"a", " ", "t"}]], " ", "n"}],
RowBox[{"-", "b"}], "+",
RowBox[{"b", " ",
SuperscriptBox["\[ExponentialE]",
RowBox[{"a", " ", "t"}]]}], "+", "n"}]]}], "}"}]], "Output",
CellChangeTimes->{{3.6255`*^9, 3.6252`*^9}}]
}, Closed]],
Cell[TextData[{
"Outputs are ",
ButtonBox["Replace",
BaseStyle->"RefGuideLink"],
"ment Rules"
}], "Text"],
Cell[CellGroupData[{
Cell["Specific coefficients", "Subsubsection"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"{",
RowBox[{"a", ",", "n"}], "}"}], "=",
RowBox[{"{",
RowBox[{".2", ",", "5"}], "}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Solution", " ", "f"}], "-", "n"}], "*)"}]}], "\n",
RowBox[{"sf", "=",
RowBox[{"y", "[", "t", "]"}], "/.", "iso"}]}]}], "Input",
CellChangeTimes->{{3.6256`*^9, 3.626*^9}, {
3.6255`*^9, 3.6253`*^9}, 3.6256`*^9, {
3.6257`*^9, 3.626*^9}, 3.6257`*^9}],
Cell[BoxData[
FractionBox[
RowBox[{"5", " ", "b", " ",
SuperscriptBox["\[ExponentialE]",
RowBox[{"0.2`", " ", "t"}]]}],
RowBox[{"5", "-", "b", "+",
RowBox[{"b", " ",
SuperscriptBox["\[ExponentialE]",
RowBox[{"0.2`", " ", "t"}]]}]}]]], "Output",
CellChangeTimes->{
3.6255`*^9, {3.6256`*^9, 3.6258`*^9},
3.6256`*^9, {3.6255`*^9, 3.6257`*^9},
3.6257`*^9, 3.6257`*^9}]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell["Plot solution", "Subsubsection"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{"sf", "/.", "\[InvisibleSpace]",
RowBox[{"b", "\[Rule]", "1"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "50"}], "}"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",", "Pop"}], "}"}]}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
CellChangeTimes->{
3.626*^9, {3.626*^9, 3.6255`*^9}, {
3.6256`*^9, 3.6253`*^9}, {3.6254`*^9,
3.6254`*^9}, {3.6257`*^9,
3.6255`*^9}, {3.6255`*^9, 3.6255`*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwd13k0VusXB3BkqK7EeyQkhHuToWSe3+/rHUvRrbhISW7IPJMhQ0UaiAa3
hEqaM1VKRaKQyJCh0miKKEPJkOH3+P1x1lmf9eyzz1rPsPc5K3b5bt4twMfH
J8BPLnJXSp7OzpQvMG+MfaXOxzdEvzd/lYK3ezJd6WXrxsULh+hPwnIvqbpn
0ReylmooSw7RK7X83YTd8+jHk43l7VWH6NVHm3ijbo/ptTc+fJjeNER/bry/
eMCtnt4Tqzemlj1Er0lp9+91+0jnSx6Uk7Qcpq8vEbHqcRukR6vU/bP82gi9
bpH+AsHtfPB4Lmej9mmUbpHMLrtiPR/pu75sfKU+Qd+9ZnLR2kXimH7Ly/XO
maJvXLV7WIoriVMX55VF/82Hwuao4OeG0rC7K3ozCwIwoanlCUbLgZ9x3v7U
LkG4/U7NyLFTxC6VJVGRYcLwrfyik7dRCbvOOfcG7piPww7Cf057qaDFbZc1
03YhXKoiXNe6/wVPg4hI5mJRPBGR+dvTXhW245/T32QuwlqJmNN9QWrIrn7Z
p6C6GE4/zDSsvDTQsE/7fr+BOM5nPpMbO7Ya+lyBVMtmcRQmGbc3bdbCF7XY
3yb2Egh8YLYpY8NatNqrua94KoHSmpyRxG3ayG8W1P9zJQ3uXKOSPXY6UJtY
q+gcTMM3VZlnBY46cJ/wOvJfCA0BCydqJ3bqYJZpuaQ+lIaIuvvvEvfoQM96
Ot4onIakLYZT18J1UN8WNzQvmoa7TgYmX8/pQEdR1sr3EA0CYXrFHp91kJb8
sST7LA3pV9fe9fTQBWctxyHrEQ2Vb1yqs3x0IXbhw82MEhpGFp5ufxWgi8sO
b1nppTSs957kN43QhRDfC87JMhom11ZYiR7VRUmXTlbUUxrsHm7tvXlLF7PO
xcKra2mQehm67PuQLs7Z/Kqd307D8R8lsX6herAy10gu/UWDd9NxZ8cIPXy3
Cjiyf4zkL3Bh8KL1oLvbeZo7ToOQ7wJ+xQQ9/PiRzKyfoCH869aY+tN6yF8R
uattiobdn/v3rS7Sg4Boou5LAQrG9TKRAz/18ELmq7+dGIXu68Eh7gH6qApZ
pjCsTGHr8pkPV0P0ERNq+y1bhcLT5HhuX7g+7gxXaNn+SeFiUJq0x359uIg8
lC7+i8JO8+IHnif1oRKXFxq2isK7xqkZnyJ9mOfcf/Z+NYXm8f0JwZP60PAS
NWQbUajgpP4XF2eACrPMyPQNFLpcH9aFxxtAsTnqvspGCoIJXQKBhw0w8lo9
+RYxp0rf2yXVAN3/CIuVWlGo4bYzWBcNcEhQx/7tJgpNPJV+oXIDTD9+v5PP
hkLH+iKzRH5D1O+xIeeGAp/1247kGEOUdB7sHPOkUF1qoxx0wBCapgoZbl4U
jq9udLE7ZIiOhi7BNmJFseouxeMk3qjl4x1vCqi721OQZQjlHLfnXr4UYixT
vr4qNURjYG16XQDJz+ONLJ02RPH2hw+27qUwiyL+83uNEPQx3/nAQQoOrC81
r6OMUCBkr95HfIcrfVIizgg9D6QubYyn4G4V/tf+w0YQNHmevSSBQsM2c0vX
dCM8dyrtuHCIwoWQypMaJUaIjeY7eeMIBebN1pXFfMY41ML892AKhYSlYxub
Eoyxs2y1qdM5CqnV1Z3lR4xRnntT+iZx5t6ze28nG+PE2IKH48RF7WY5J9OM
wS//x+aUDDL/WQembK8Yw3BLSd3jTAoMVdqN9kpjpJ7QUaRdoPDbUHN+t5AJ
dqe4WB3LoeDr4FI+tt8EiEvd2XOLgmjW10+mh0xwzs6aoZ5L4Xqn/2zsURMo
F2j3+RL3eMWYip4ywVJtpeRx4h3RGUWKl02QopRqKJJPYcOltuvrqk3A1C5S
lS2koPrd8kS6qClqdBzzlhZR+BSr+y/9lClcHKzDwkso/CklluR3xhSje8R9
Cog9rn+5dyHDFBly/Aa9xD+azv4heNkUCsmzFjalFIRV+G9XF5kirnUjV/Mx
Bc3KOv7Nr01xYqZcoaWMQvhCt6x/l5lhPr/drOBTCtSJM28TL5jhaIzkhYHn
FBa2UI0BOWa4UbXeXaqGrPfSpKpt18zgU5qyGMQDZ+PuaBaYkfpQMJNKXHnR
M6mxzAzl+S+dDV5QCCs0s5D5ZIY/O8U3hdfOnY/P167KmwMX8r37XlLIEVfb
W5VujhHD2Bj/VxRW99iPC543h2bwj+7TxEUPE0MtLpkjrpan8Yi40rUvuOSm
OVS8BRyFmyn0ProacLvEHClRQevPEGvsWemV9dEcqdrqe4tbKBQ+UXEKVabD
NflCdWcbhbIABfbKG3RciZDhdbdTmBwdPC6XR4fCk82Z/O8o6O4teydxm47S
myEflhNfjXYOmnpARwerXM6W+PjR7OzGGjpytuHEM2KnyysFIr/SMa7rE3L+
Pdn/A2lOvZJAd+PDBM5Hsv8Uz8uZbQVGZSPwTwcFscNrqzVtgffL+mR8iKN/
lAfK2wF/Nbd+PkC8s7K7ZnYb0LnD1ayQWMlLfW+5C5DlWqki2knh8r2iFm4g
0FD7/fEj4lyruqTNqYDOrESsRDeFkshJPvcGoGBy5s2RLxR+5fUMxTUB6nkd
1ueJtTqbPmU0AyP3H925Q3xp3Y2yV6+Bg6NHLN8TH5ZyjKF/BqTvL8/T7KVg
m18yKzUCqJbJO7wgHuyMmXlGMXCxtS1jso/CCkuRKeV/GLgb4nJDfYDMXwuz
SdaegVOBPDUTYo5TzFWJbQx8jlqesZ7YI3Bi6+wOBnrGHu3YQ3wnvf/WW1cG
FC82b748Fz9Q75QSwoDQB3dh+W8k/th/FTOnSPzh2jeC30l8/aojb5oZENSc
EC8cJOuh/CRnspUBfbnXHY+JPUPtypa9YeBb04JrdcRKCgmj298zMP30KK2P
ONmny6mjm4FLhiu2KQyR/IuydPt/MaA9/q0mkVjRUvLDlLQFygcqXW2HKRyt
nNFS3G6BhFcrDzSMkHro37N2uZMFXk/lj74l/rnspbasswX2CU3u6CZ29M/Q
ldxtgdqsbulJYk05U8P53hb4btmmpvyDwkv/vfShSAuEVRbVBhKLL/+5seyc
BdJzx3QW/6RwIrDPY+c7C5yxPfV+zSgFWU64/5YPFngnv5JlSHxR+o8wzicL
5F7UuwTiwhKNgxpdFrgy+m7T38QNIn5Z4/1kPPOTYwCx2Llfr47/tsC5tWKs
28SHns4zeyLLxOH5DX5avyhELpFfrGTPhJ32UY7MGIW6ggUO9G1MZLOtqxSJ
5a1GLzluZ2J4eTlDlfhxfK1xmjMTxpxaZQNigfFwV1EPJqhzQhlbiRPetpX8
CmeiwzLNPYk4NSvF68U5JuS25SZME19bJVQT9ImJqfOjP2rGSf+SY/z27mBi
h5fU+QbiqcVRGm5dTKgPPuG2Elv9+plk38uECq0q8TPxz/KuLeZDTASMfX09
RgzHp+9F+FiofXjmsvIEhddJ+4fPKrAwHCO6NIxY+Ce/zJMdLDy5UTwrOUlh
PP+/pK6dLFwWnScqS9zvvUZwvgsLX5xnKQXihi+Og1ZuLGjEKyxdRXy2/d6z
d74sTP7t8tWUeHWFd8B4DAsWLcuqXYhtU9+8WJPNQvjBmrBc4itaBfsyv7Cg
qR+3x/A3BZ8lwxP3+ljoVQz2MCXWndQKbuxnQWCdmDuIyyryPQWHWFgYa27P
I35tm2/nMc7CDscuyo54flSetv4CNpJP3RwIIXZ/fqunTo0N18hB80LiVbuu
W017s/GvmuqowhTp7xnG+SO+bFjdSppSInZ9UyPR68+G74k0/r+IL23qb24K
ZmNxX5qgBrEiXWPb1Sg2KkXivxsQS8nlum05xobn+1tbrIkFWgtirt1kQ5qR
4htJ3L6uuNCmn43AbK8NDcQW1pfsgr+xYXC1dLaJ+PrW5JmTg2xMr6jKbyEO
d9q9vvkHG2yteaLtxLLBEh2bp9i4tW/7lW5iu/N7JDYt4kCokRk4Sdz6S8Zv
/RoOGjtfeClNk/2RHa5JD+RAJDh4nRcxja+o2DiYA4XwlC0+xFsch9n6oRzQ
osUd/YibJd13aEZwwH/B1D2IuO2gTfKy/RxSz/y8I4nfuWkNj6VyIH54p8ZR
4m71njv5hRz4bfl+4Drx+O2/TZVGOFCJVu/oIBY11H9J/eRAb82i+C5ixUey
O4V+cRCgfU21h5j3tHN/7wQHXGVJjz7iM83BL3L5uTgm0/h+kNh49IyDiQQX
iUZliVPEkfodYVu1uDDYR01KzpB+cT/gbrwPF1e9jXwtiN8JFx656MfFgmzr
Hibxg63DzqUBXKhoqzmyiYOGfMV+hXARFKfF4hH3rvR22x3NRWBJppAVcdNp
NxnWcS6GHLqt7Ilzghyj+Au5KB8RFfAlttTisCNHuejyuamQRszX0sIfNMZF
QqiOyX/ERXt3l3pOcPGlrt3mDLFSxQGDbdNcZMifOJROPGlTscpYiAfdjXmd
WcTXIyA2LsnDw0W6/leJF1QZtwXq8rDCRuj7feKq7Wv2eAbyMKjQktJGzKmV
8TgWzCP/V/Sk18TPjAU980J5sNN7lfiGuGLpG68fETxcSQvb1078uDHOL+IA
D4Vj9ds/Et9jt4YeOcVD/PTG0R7iK5pR8Tfu8SBCWR8dJU6Yrsnu/80D48OL
feKz5Hv+pMzob2LJip4Dc1ZUd+eITvEgEcN/ZM4P7IT6NIhT4/TT5jxwm77a
h7jyWmbenDd53Lk/RPwu1u3TnGXbMupGp3lodWpnSBD3eA/IC8/wILRhmDvn
24ImflLE0eIiVnPeoP2aZkDsoKTtMOfoYzS7MOLqPfEBc+5iJnRMEu86tDJ7
zvlvW3T+mOWhYsL42pyj/FUOLiOOYVrlzVkq68lKU2KJnUEP57xu8rdXFPGt
ktJX/x9PWV96jHjBSMObOXesPLM4k7hqpuPjnCNt9AtLiUMuCA/MmTdwYF49
sZuB9Micl+x/tfUj8eOcVeP/f15G6fIgceyA8cycc/P9xmbn3vfHBkEa8f8A
AspectRatio->NCache[GoldenRatio^(-1), 0.8948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox["Pop", TraditionalForm]},
AxesOrigin->{0, 1.},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.626*^9, {3.626*^9, 3.626*^9}, {
3.6255`*^9, 3.626*^9}, {3.626*^9,
3.6255`*^9}, {3.6254`*^9, 3.6255`*^9},
3.6254`*^9}]
}, Closed]],
Cell["Several solutions", "Text"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{"Evaluate", "[",
RowBox[{"sf", "/.", "\[InvisibleSpace]",
RowBox[{"b", "\[Rule]",
RowBox[{"Range", "[",
RowBox[{".5", ",", "5.5"}], "]"}]}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "50"}], "}"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",", "Pop"}], "}"}]}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
CellChangeTimes->{{3.626*^9, 3.6257`*^9}, {
3.6256`*^9, 3.6255`*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwd1nk41csbAHBbUllyvlJEnRQiQrniOMf3dfZsqZBQ6spO9mSLULRYkook
KiJ1ESJFKTdkzS5U9iW6yL7lN+f3xzzzfJ6ZeWeemXnemV1/ux2z5eHi4hpA
hVPLxKw+ebjjpY70b0M1Lq4evEhAYaerQwxONjDbfmxTD/7hYnbaXocUvEo5
Tua5eA9eoephz++Qgz9Z9b7jqdyDV91sYs/av8d3WofbnbDswT+TworH7Rvw
oJwjlwNe9+DVt7o8Rux/4Fu+y3bRLvbieqXrjYbsJ3DrzTdnw4X78TohjQ18
p7ig5Mr114m/hnBqDKMs44gATNoqFe+VHsNtVZaE1IQ2gzVbut/m8wRuqGA7
Jc4SA+GG5c/xH37jeS1BPp81t8G6ijb1ibezuDZBMYcvWAr+op4uq2Yu4vbL
ccnp5kToqFaOkduwirtVDB/MMZQByp8Tf6lYc8F1C37ZVZc9QAwtj+315gGb
ygA7NQc5sCD5fGG58cGH9RJHnU/uhazKs40vTvGDmmjI3VFvRTA5Omqx11EA
rKcpSkYuStB7eUzuPHkjpD78JDUftR9M1hKNfbYKQl40qavpmCqYGulWlZcL
gdcbinGygRqsbiPYdCuLwLvq9N/XLA/AVFI4UZW5GRxYWqWO5gfBRU143a62
zZCUqfbK2Ukdqk9oiKWdEoXY6dLL7r5/wV67HGe/BlEYzPK54OCpAVd1JhWk
FQlQzoxLCA09BG+bQ0MNgwjAdaSzLyZEE+7TXFtvvCNA1TvT3d7hmhAvsp/n
wnsCxO5vtDGP1IRbv/JjrMsIQBSuGiDGasLunKB61Y8EgLpXQy9TNGE1a+Ng
zScChOjf+tn8ThMWhAo/DdSi+Gz2762rmlDy/Hbtk04CrEEhd6qfFhDUh+lO
MwSwoA9XdwRpQf25K1tglgAFrG3xoqFaED70lL5ljgAORv5yYde1QP/jzrLS
eQJ8sdTRt0vSgtettw6tXybAowsV8UqlWiAsvOdKIDcGtBdt8sVcJIh68+DC
tBAGEVvnDZsiSJCkNj25VxaDuKqq/o83SNAeFlxRh/zQ775ffgwJDpbayXrK
YVDYRUmPv0eCo0r+IcXyGAykhK+YZZBALV+ygKqIge5ewvOuChIEP585jqtg
sKypLDC4Thu2GwyH9Wti4GZh83E+TBt2m1snDOthIJjys4ccqQ0hxSk3A/Ux
yOr3WLt8Uxs6jzZtFDXAYMglhCx4RxuGe+yeahpicDo4uZD4VBtKHHJsQo9g
YJDWnnW4Shs8xOV1uUww2Puf/u0kQTKU3Un+kmGFQc9l9XP4HTLskjs8WOSC
gay4cLR7IhmW5whpKq4YOGUNFz1KJkNc22m+DOTppvub+J6S4ah2fFb8eQz4
93DnVxWSIYu5+Y6zOwbKFXXcxzrIsJgdt0PAGwP/jfYp57ZTQB77073qjwF2
O7Hz2iMKfPktx/gZgcHGVqzRM50Cj00Z7YxIDLi2RldaPqMAVz5J9hHy+P3Q
AuWXFPBY7zxrdg2DisfO0Y1lFCjeQ3R5fx2Di3kUqkQPBcTUnLrDojDobux9
lrlDByQd2oW74zBI36zoV5mkA5uZSdHNDzDYP3RygS8VWSJ0SDYZndfba77U
NB3Q6z677ItcYTfqU/pCBwTtvlltf4jBSEmmZ36pDvTVePNZp2Cg5CjvkvJD
B67J7cc6H2GQ92GPte9uHDaWpsakPsWgzHMnQ/45DiS29D2/HAyWZidipXJw
eFMbS3mNrO5X1i2aj0M7WbBgDjkz+Kz3yhscHHozRb1yMYi9+eRJYzUOFZOm
tx1fYmD9VJ4n8CcOX622NBjmY7A2fs96RAxg1WF260whun/EVCmKCQAkOu16
X4qB8HW1KmUzgECeU/VTyMHTH712mANYbRI23fMOgzMVg9VrlgANZlKTkcgy
Lvv8PtqgeCIL1sbvMXhaVNjK8gLofvxdp6MMg2yjuuhjcQACDluLassxKA1c
4nL4gqbjl6UafcZgLmdoMrQJYL+0fFkAsmp/U09yC4BxyJttz5DTDj8va+4A
uHxXyIK3GoPr4lYheC9Ay/csYiGyWW7pmvhvALw4vV2sFoOJ/pA/nzBdaP+0
o6ioHoNd+utXdp/QhRzxVnZiM9q/VlqT5EldmCr758tbZKZ1SKaopS7cCxim
fEd28lo0WTutC/5F9VUyLRgUJI3902mnC6KKyY+zkJnjDda3LujCFuYt31et
qH9UQvmfO7oQGfRo/+t21L9B4cbXFl3Y/JPd59OFzmP3h/SlNl1QF7v/Nh7Z
2de8bPtXXbDwFgnIR5bZGTF76htq1yx5MYkcc37Aum9QF2ZL3L45dqP4Qinq
Y3O6kP3q7tUT3zAg6ot9X9lGhdt+5AjiDwxuVvxRJZ6igqTCjkjXPgzAY0hN
2poKKoJbe8ORZ7bXH5A8SwVfvRL5B8hWHsnqYrZUmHhBv1SNrCxF1hRwpYLF
pY6jcv0Y1Hv44ZOBKJ4oI6UTebP0jGHZAyqcWmg6oDGIwW2vUacz3VTIvpGx
rWQYA0mmv8fx71QYWuQfqkd+vG3TRWYPFXQth9J6kfNKla4oDVAhl3x0nn8E
gy/r3VMWxqhgELtp0hhZ+MFcc+wyFdxoCaQ+5Mh/eSkfJGngoph4aWkUg8At
O0RkTtLAU29shXccg7qXGyxwSxqMSU3ViCHvMJpNszpFA7y296Ys8vurtaR7
Z2nAfSNvgInMs+BvJ+hEA5FZs1+RyBGd7aVz/jT4p6PIZeMvlI9TbrnUPKBB
1qypJ+9/GDxTWFft3UOD3dpJEw0TGDRJ6S679tEgrdwjoxt5RSRIyX6ABpGn
j5mOIhvNzUSfHKHBvkrj2zyTaL8/DhzXmaTBuUWJRHVksPr323ouOiRBRNNd
5I7osKn7O+kQJTCqajKF8t8Mt8SH03TYeOu0YsFvDBZyE6IHztDh23xa0Hvk
MVcVPgEbOjTgBpXVyF+GrSaM7OlAnP5G6UW+31X0qduNDsfGJ14LTaP8VO7q
uRBCh2z9ql/nkM3ivtaoPKGDnd/HSMEZDDJUX156OEwHmTedMtRZDM5vmVos
GqXDyK7UnfrI6kuqPo1jdHC5NiBmglxWnuvMN0kHExbvqC1yh1muudMCHcR+
MTSvIQsE5RzQ2MAA8xeki/XIDp//GapTZMCiYH7G8TkMFP7OMlp1ZYDy2Xxe
vXkMqpJJub/dGPDBUaTRGNnua7XoiAcDTk5LJ55ATjMea2nyYcA6bhsJW2Qi
rmSZGcSAjYlrncHI4lLZ9sejGBDIJ2RXgMzT9jLk2QsGmOooFEksYNB1uDjP
dIwBAer7JNqQqUfSzH1+MYB3TkuxGznLJOZP/AQDjh6X1+hD9re21WuZZoC4
Uo7Of8iSPqJ9x1YYEH1vjMC/iIF5qqOosRATVh/V1agjt81JuOupMKEnX/VB
NHLTE39l3IsJijXq8epLGBC4CotJPkxYs/JX0EI+bjXF0PBlwhOeXW8pyC1i
DqeVA5hQNPWgmYncfsU0ZnsYE3xLHn41R+62V52aj2NCdDaFHIg8uG+oIDeP
CdccxdzKkBfyj5JlfjMhJDVNGl9G77umRj02wwT7s/Z1VGRiieSZdXNM2PXD
yJ+JzP63P2xkkQkJg/E1hsiJLT412dwscIohGVghk2YTLbRFWbBi3lR5ETlQ
o++iiSoL3A8KFeUgr732fHX1PAsUPWnDW1bQ+vjzbjx2Z8FhrUsHJJDfmEyd
fefJAuE/0/5SyN6TbsJzF1iQUJPIuxt5RN7V3jaYBSZbRWZUkJvu2kvQY1kg
oZhy/DByurdVEHceC3aRRqT8kfVVmYzAWRYculxp1ILM1drK7T3PgjtyrNo2
5EI/23fOiyz47NTE+oosUx5+yHIVzTchqvUdecm0XIG0jg2ECGu+EeSsABBe
EGND+a8r+DLyhkpSu5c6G2wy9Y4RVzGoPKXi6OzFhqn+M5q2yMxaCacoHzZk
eP1xs0f+ROJzzvFlw6xwcbojcvnWry7TAWzga7yx6Tzy+8ZQ94BwNqRjLWUX
kIsYbb437rCh++WR5gjkDOWgq8+L2OgfmvM4EzlitfrJ2DIb7lyoIQ8im8RL
zC4jnzx9VY9j4j4HpuAKah8Hc47fmK8bVUL2XSnw4ng8H99/Hnnvz4Qsjo2d
Cl5PIpPXjmwbQpZsT66bXWWDNDFgmuMh1/Ed/H/YINSlyD2MnM+n7S6OXEru
FObY4EAH4RDy3zOHFDkOjiKYX0Tm/mviDMcDtIi+JeSEVYN6jnM7Ww9uWmND
hOx8F8dBHnuubEe2mEkd5Vg85YM8Gfnl1Wm+EeTDS8suQchOl+K1ORa/pfcu
Cjnvm/ZhjvvkE0UeIk/095lxHGiqkfcO2WRQxZNj9ng4bwNySE9rMMdbwppN
fiA/Cw2I+v94CZmnE8jNNcQkjrNz3efXkLlLPmVy/D9Tz3hK
{Hue[0., 0.6], LineBox[CompressedData["
1:eJwd1Hk4VW0XB2BEEilnl0xJGgmZkvn8juGckwpJJfHilUzJXBEllKFSlEih
0FxmSkWGQgoRRbwlQ0iFZBbfc76/9nVfz1rr2s/ae61V/3paOPFwcXEJc3Nx
cZ6yF/6mp0jn6Jslndbi4pqgP14gt9LD5QLdaChRSG3RBL3sWGbGBpdU+rVV
QjxuEhP0SmVv5/kuWfR4gdBt3eoT9OpzjexR5xf0aY1OyYVuE/TX2mFFP5zr
6e01bi/EmifoNbFt3n3OX+i6mf2PVDIn6SbF/KbfnAfpmo2r3Mrdp+m1izQE
eG25cGWCVz7Obo5ucMG49I7ZAoTThidLkrnhtGlqkcqiJQiYFNgaeX4edsg5
DYuylmJQxj028B4fcpuC/V9rikGz3u6UwlV+6NDks3hPSsHH1vbUvRQBOE/H
Jd+ykkHD6sHW8GBBeFb2qmXtkEVd5NNHlZKLEG09f+3fQ2sgk+tfTRUIw7Hq
+EEVl3V46Xk7tJ5vCcr4xXe679uAx5JN84WuLYGKSMiVfj95bJ0zG8jfIAK7
ET0F00MK8NDwUppNEsGNlFdS4+eVEJirs4iXj4bcGO22RgtlbNyZo91rQ4Pv
Uz3z5O0qKD0x4iOTRkO+dRVPq7kKNDOv86un0zA2ZV6wdLcKdpgueM/MoOGY
zgGJc7YqEBsRXuJ2m4ag59HdgYdVIFAgr3znPg1nSj4es7qogtSsTJ+ZXBoS
y71vUk0qMHf6tVW0goaSmlu/o/ar4mR4+YJVX2nIFnvnHWynihTp38/GiNMP
Tg15Oarip+Xu1jedNETymA3udVdF41lZS79uGnZpTwysPa4KCbOXu0p6aei7
Z9JTlqQKWuzSGo1fNFBRgx8nW1XRqqAw7/kUDS4srWJXKzV4fXDLk6dR+LlB
/FWOjRqEVjQl/iD2WTj5dtJeDTkV73ofURSO1z5pj3JVw7OhNS+VllGI2aU5
cy9QDX8fxJ2SE6NQYLdF5/t1NcgtuOs6T5oCz7HNRW5f1dDxVq5gmxyFa3dV
Ctzd1PHprUKZoT6FylbH6tTD6oiedtHqIP698Erbex91HFynoRlEp2DiMcWt
e1wdp63ba3JBYUqlwlTonDocBJ/vkDCkYPXMsu/hI3XwBbCXvmdREK07Kvlr
SB3L8vJ//TSncHGk+JTX0c1YYSN3/pMDBY/Giw42xzfjoQn/ZYt/Sf0cRwb7
5Ga0J2cZ1BDzeQpwy0RsRpfUzsQnjhQCv1uG1F/ZjPIY942xThScvg6cUCrc
jBolxbWqrhS068WDfvzZjCMdO7g3e1Houe9/xMVHA+NKeqUmwRQsV8x+vntE
AxuP/uksJH554QyrP1ADg+VRmbInKKT5JYi5hWnA7cFexiSxvX7RU/fLGqDP
f2KTFkKhvWFm9nChBrxWFqR3h1FomgiL8J/SgHmuoNrGaAoVzLjE0NAtsI9Z
W18bT6H74LPawDNbkOAZTZe5QoE3opvHN3oLtGhbvX2ImVUaHo5xW/D6UpS4
aAKFGlYbwyhtCzLeLgy1TqTQyF4zwFe+BUmSrE8fkyh0mhTqRXFrYkO/Z9ut
VApcZp86L4RoQlZipUPiXQrVJbtX+4VrwnT9mnVtxBeVGhytIjVhtVigaMU9
CjLC1d0yFzVxPXje1jRi1BZ8y0nVRE3fw2X37lMI2Rb7/X2JJry6qVOPHpL6
bPbv5X81YTHiFB+XTWEOhdw3ArTwNGkerfkxBWuj3pqWYC2k3KQWiTyhkM8S
uywSqgWhztkP24ldTAPXhUVroXnHvJYK4nf79bcdvKaFGLl5l7KKKNw8UnlZ
oVgL8kuPC594RsHw4Yf1RVzaeD65T3aihELE8vEdjRHaUGSOXnN5RSGuurqr
/Kw2/MMXKF0hTglICsi7oA2hdFpyBXFhm96tywnasMnwo0lXkv6nhs/suaON
R/HbbBuIGRtoD9oqtTGg0RulUk1hWlNxQQ+fDngvrR1qqaHgae1YPh6mA5Fb
NMbregpCqd87dCN1UM7vPTlMfL/Le+7UOR0oOpdckHhH4duhEF2heB2wm+P3
uBP/czK5UOa2DvKyYS7YQGF7xsf7W6t1ELU0StC4kcKGX9suXRPSRWR2gu2N
Jgodp9QP0ON1MeMbITDTQmGtqHCM11VdFIlb+ki1UnC73/v4ZrIucve6luoS
jzQmCfLe1sUZ5kuhYOL5a7jzqgt18ejLFZcZYsXKWm6LFl3YsszoY5/IvCx0
Tj0gqYd9qx5MfWinQF26+inqph7Y5zO4d3ZQWNhMNfjc0kOHSkGCKzHX8piq
/ff0sPlHvUQo8Y+k0HzFHD2s22fbl0tcmeYe01CqB7xwTaO+UjiWq2cg3qEH
rj15KQ1fOfPx9d5daX2EftERRxeFW0vkA6qu6cPHN3O2v4eC0rd9E7w39JFn
0zI+R1z4LOqoQYY+TEQSPy/7Ruof7PcvfqiP0b5l/gzivud3ffKK9WExmleV
QKzguv5Q6hd99LlfWcLopZBbtsbu6Go6TBWcEs70USj1WWm8/gEdov05Zzq+
k/00OnhRKouOsSDD9cPE6gGl7SJ5dOzWf/uUe4DC3ZMOfjNP6TCrGXgqS3zx
XHp6Qw0dYu1BMgeI7W6v5wn6TkeKz/WobuK5Hwl2fUuBeKWR5e0/yP8nc0NK
zxJ4cxpKd39REI5WqVbcA4gJ/IwrID45Uu4rbQXcyjT5WU5sX9lTM7cfsPii
GfYfseyhjQHljgBj3nN12iCF248Lm1m+gLX/6oIA4kzT2hiLOOA1K8CeMUSh
OGiKy+UdIJKyeDx/mMJY1reh0EZgY+jD36XEyl2NHclNQNhLo+63xBlbH5S+
bwGyH0Xe6yaOFrUJoX8l+YnGxdRvCnuyi+dEfwPyjUVLvIgHu0JmX1EMvFht
n7pqhMKqbfwzq/cy4Dq1n8/tD+lfs2GjxD4GhOTbRH2JmXYhd0X2MyB9wHhV
ELGb76Tl3D8McE8OycQQ518bePTpIAMz5XpvczjxP+rtYo8w0Mvj+nGcE38+
sWI2nkHmlL/nxCiJr5c729rEgPkzj4mAMfI9VpfdmvpA6hvcXxRK7H7UqlSy
lQH/ukVS0cSyKyNGbf9jQEbEYU0S8YXD3XadPQyEf23pKiJ2W5SqPjBG7iM5
b/UEscy2pZ9nxAzAd3G82nOcwrnKWWUZWwPkPt/6afsE2Yfe31RW2BngdLLq
qAXxH8k6VQkHA2TJMAX3Edt4J6svdTJA3Y6ZDU7EilK6mgs8DBBjvXhjMHGd
dwB9KMgAu5nLdzwkXrLiz47S6wbwtzeV5Z+kcMm3382+3QBHmErNWcQSzEDv
XZ8NMEDrKs4nThMTPMbsMIBbR1VaEXFuscJphW4DvNqu6FBB/I7fK3ViwAD2
F4qyPxALXx97f3HaAKMnHJJniCNfztMrkzBEuNbrUuMpCkHLpBfL7jMEZVdX
VUtcmyNgTd9vCNXgL3KNxNKmoxk2tobwo+iRH4hfnHmrneBgCDfamNYXYp6J
wINCboZQy7rqPUQc8elj8VigIQRuOqjTpsk+To099Oa6Ic7rl0vtIr4nx1fj
12GIdjXW8jfEjVKMaY9OQ5RFWznXEc8sDlZw7jbEubi6/AZi07E/Mfv6DCHZ
PcRqIf5T3r1Lf8gQJ0NYZj3EsHn5Hz+XEd6Nz3bNErfEhA0nrTSCsZSVmvIM
2X9/uMXL/jFCiWxT41niiezEmG57I4iYO76KIR7w2MS7wNEIjlOHC2KJ3/Xa
DJo6G8HPkjcmgTip7fGrdk8jrEj5d1U6sVKFh89EiBFk0isai4j3xLW+2ZRu
hOHFhw73EN9RzjmR0muEJlmfsc1/KRxeNjz5uN8I8SvebtckVp9S9m8YMIL2
qPJNbeLSimx33iEjfKlcZUQnbtmTbeU2YYTFAbU+bOIFwVmqGgLG+PMx38+a
2OX1o2+18sbYJPJ8WTCx3L/3Tf96GGNUmSlTTFydrJ3929MYp1gd3SXEB1tr
RPq8jbHkttjdUuIM84GmRn9jKJ51kX9JLENX2H832BgLx9Kk3xCLSmU67zpv
DLkkz8oWYp4POSH3HhpjRd4diRHitq1FubsHjEGPmBFaM0vBwCzDyv+nMYIV
j95YS3zf8sLs5UFSv++a6nriQDsnk6YRki/TZyFPLOEv0mkxYwxHgQdhysRW
N1xFzBcxcUL2/W1d4g9j4l4mm5jgGdaf3kXcmB6oSPdl4mnCJtGTxDSuwiJt
fyYuX+1ghhDvshk21jjKhEjwjSOniJuWuvyjeJwJ+S6z92HEH0/vviAZxoRr
56lTkcTtzsrD43FM9AZ8yIkl7tn4LT87lwnlCAPndOKJvJ26sr+ZWMhju+0l
sZCmRh31h4m6z/Z4RSzzXMKeb4yc37itXknMftkV1jfJhPqKKYlq4qtN/m8y
uVnwm9/w+Q2x9uhVax0RFo4bFhu+Jw7S6DxmqczCY/nDZ74Szz3xKThzmIU1
Da+FZjjvNz/3bJoXCx3Jt3s5fmo57FDiw8JfpTdlf4n9hjyFx46w0C/wy3eO
uG+9h7PTSRYUjj6o55kj/briLG50kYVIpWgnAeJbfjbB3LksCKr3CIsSb1Nm
GgeNsmDydPSIIjFXczO33zgLipGOUkrEhQFOJe6TLLT9Y1bGsWxF+Jb9f1mI
K0kTUCae2l0hp83HRo2t4WVV4vvHITyxlI3smorYLcQCVdoffdXZqPt1WJlB
XGW7ydXdl42BRz1LdhEz34q7nfdno6BLPoPjV9q87llH2XDJ+61hSVyxvPXQ
yHE2FN4I7t9N/KIh1Ot4OBv1Pp7Je4kfG384ejaejVcqVoI2xHcUg888eMzG
XH9AzAHiiL816QPTbFwv02P5EVteFh+dJt4eX2vCscxGF6bQDBvFO21NOX5q
xdevQDwQHrSb4x95dKXDxEeaihw5NnfLfzJE/J+YcgjHEh+Ta0f/ssGvI1jE
8TePH9LzZ9mQ3pD4nOM8Xh0vUeKFs2tKOd6u2kLbQvzvAf0qjk+ep1kdI5aT
9GrmuNswonOKOEalfpjj7E/NaoJzbOg+tB7lONh7zWlJ4qrF3yY4Fk0tW69L
HBA7M8fx1qnpQ8HELic3CPtzzmNNSs4TizPyRDjuXH91cQrxwwG9ZRwH7dbI
LSEOmbaQ4pj9I3xePXGp/X8rOV4W9t7yC3FPpvPq/+eLy94eJB7uGV7HcWa2
1/gccRdfkDzH/wOwVKCq
{Hue[0.57954`, 0.6, 0.6], LineBox[CompressedData["
1:eJwd1Xk8Vev3B/BDUprk7FJk+qEMmafM53MOzjlX0iQpSq4yD5kyJVNxK9Go
0kAoFUWIFN0oQ24yRCmN5qIo8/x7zvev/Xq/9tp7r72etZ7n//723XaAn0aj
qfDRaLyrbNJMxnWpB6ZcExNZGk0AxQuVpL3dkhjNsV/FPywQQHnI/UxFt1TG
UYPVIoXCAqjS8HMVdMtlTP9I7neVFkBNQhN3xPVfhp6iiES9qQBeGsaW9LvW
M6pqtBmXIwRQe6bNr9f1C6Nf/Eru4ikBWJYtsO52HWB8zE0faR+ej7qlekIC
e2iQpVUXdSxcCFaSxbOszQuR7xrILukSwgH1yaWaS5cj4V/ai8nri7FJ6cBv
Uc4KLJfPC+rZtRT5zRFBL/VXw1U77dDbxmUwoivnCkRKQPtuW7G69HK4Tp29
dtNOBtaem3q9S5fDt6pHO3eTLIzrhc19jEVwYrfg2hkveWw7pNyjVSgC5+pw
F023dVDQOrs+j6KjfIHYVs9dimh3arx4yZmO7gw1y/dOinA+v7b23X46lsDc
jOOhiOC3v8RFXeiwC/HVlQtXhHZrq9hpNzoGeivF2q4qgvW2ISXUmw7JWv92
y6+K2KTCWCUXTEdowqsAZVclvIzhP/DmOB2aIlHJ3wOVkS5bb/sih47HC98p
Vh9WRoOSmELZPTpYNLUnmUeVsT7zxaaH9+mwGWj74nheGXs7KlIz8ug4VLdB
6W2BMuol+vaFFNJRenzgScUfZYTrzsvhL6WDy+/47crB9ahSUKddq6XDcchE
xdpLBRqbuqiZLvI/dFpDc6AKslvtB790k+9rVATYR6ggJrlLpqKHjhU+7Cdu
iSrIyRuLPPqdjppea8vYPBXw6fmy+H7SofHZ0e3RsApknrYYtg3RQXsZnSkb
oYorlX20FXwU0q5XSoydUsOD0B/heWIUrjiovNuRrIZfIps3u4lTSBY/d6bw
uhqi8lqTpddQSEjeJ+ifqwb1tNbYBAkKIYlTA/0NarAyGPziKE1h6xHN5+2U
Ovw3CA31y1Pg33vN/fVldZQEewpnqFPITzRsa9qmgSgj2dazZhScriyb/9hW
A/VS41tkzCksv92udmO3BkJmrjvfI/YpPxHj+7cG5iX67KyyoKA8/EFpib8G
DKyDooY4FNJ3hYdYnNFAx6n7dCMrCmflS1c+qteAGLdfINqGQsBjky3XrDSx
iZa+ctSZQuHuav73WzTx80PjWNB+CqOTWx6u2KEJ8SchgSPEIUb7xRP2aEJ6
v53p0AEKh0tPdIb5aEJuOlHihyuFuKfvQuxOa+LScKlZpSeFSxV+N6hmTRju
KjLUDqDwtPbmn+P2WkjctEl2JopC3uoGvwhHLUwOBn/eG00hw2Vy8KCzFlQP
2Gx9RvwP/+aBnZ5awMc/9jExFLYbjvetDdfCtwLfRfOOUui9Y9lVnqKFyltW
E6NxFKjjA+8m3mtBSlg+8m4CBTeOQZm7nTZ6utVeWl6k8FNRrPKBgzaWLN8b
kk7sv2ji1cQ+bcR43fwyQRxe9+jjcXdtOJvvfZd1iULidv3pO2HaSNqZ4kdL
ofDQcYPRj6vaSFlHGuUqWb8Q3RKPb9oYWul8quQG6Y/bmg89PXSwcl50sWw2
har3zjWpPjqwXOr0wo34z6Lktjf+OvhyxOXSfWJL70k+43Ad7L19Ltkgh8Kk
5nPrJQk6+BlmPWh1j4LdE5venHs6+HDwEsM9l4Lo6+A1vwZ18NQwJTS0gMLp
obLog8G6UFV9s5PvCQXvptNODuG60BDTv6lHbPnAmcmN1MWSrbMvPYnn+wrx
ycTrYsGvdSktxGE/bKLqk3XRsknNLquUwoFvfUfUinSRLzyz3+wpBcN6scP9
w7rYMZ6RsLecQtfdoENu/nqQStNer1BNwUZy9vPtQ3o4lxDdaUv8IimO8z1M
Dxs5mofiiNMDL672iNVD0TU9uS7ifaYljz3P66E11TYrrYbCx8bpWZ8iPSyn
nQhYXkuheTw2PmhSDzcPPxl684rCc/bZSzExGzCYmi74u5FCp8uTurC4DSga
2ay7somCQHwnf8CJDcj3buIYELOr9bydzxK/zlkdTVzLaWOap2/ABeVlXsJv
KDRx5fvmV2yA1w3B/nXNFNoti0yO8+lDSLf8G/MtBdrmD+1JUfq4slFkav0H
CjVPd8gFHtVH+nSw0l/Ep9Uane3+0YfFP3FMF2KZZTWdMqf1oZLwWS2NGHUP
ux+k6kNVZjONaqMQtfHMjzdP9dG95UHCMDGNy/2zakYfH/OmSrM+UZhDEV9a
qAG+Hcv71viVwm7zntrWCAO4ZCXy9RMXclafF4kxwBqv7MWC30g/W4etiz1h
APX9vZ8NiBvsTTe6XDFA6Bmn2DTiG4eqzquUGWBvyP037u0UzHLeKpTQDPHf
0Auhnx0U4leNbWqKN8RJybB9Jd1kP6ip6ag4aYhbwYrKdcTXQ1NCC5IM8e2C
V+dX4qI2k5vnLxpiVWaX8sIeUv/Uo9O2WYZ4WtxFsyVmKtKz26oMsVbStu0X
8ZS+6sKu+UZgjA5eXvmdgu9u54qxWCO4++qkG/RRWJL646vxP0YYPKfG3kh8
t8NvLjrBCFqXTNociLu9ooyXXDBCnNmr1kjivZHXimRuGeHh7AK5SmKrzHd3
/6oxgo5z0axVPwXFXxvPXVliDHHxunXbflL4Gq2zn3HBGJut2CU6AxTWii5L
PHjZGDo6yqUsYo+7PcU3rhlj6S7Xgi3EQ00piwVuGWNNExXpRSwoz1dQU2QM
qVXMhAxi1ao6vm2txjitrhEnPEjmZZFr6v41Jhhs7dP/Qkydu/zh+A0T1EA3
zvEPhUUtVKP/TRN0OTTIehDTViVW298xQUOSanEgcX9KTKHqAxNI3uNUHyeu
SvdMbHxmgm1t2ccKiEPyTVhiX02Q9isnav4Qbz6+3bktZQr7v51SMohvLlcO
rb5iilkfqcyXwxTUuneNC6SZQjx1UuINcdGT48GsTFNIV3449ZG4yuV7UFmO
KdhPyncOEPeW3vYvKDPFzXOdCStGKKi4K3ilfjGF4XzTmD3E+eXyjsFyDPza
vrbqO/Ezf2kLhWwGTir+WvhzlOxPIwOnJXIZ6Gqr1Rwm1gl99lGkgAHhoXk7
pohvRzoFTj9mwG3m6gmhMTIvCRkZjbUMXL7Qlb+W2PGWAv/hHwzks65PORDP
9V907F0BzMW6n6wmvi6TJmFiA3C2PIw/NU5h2QnNGlVbYEjsreU54sihigAp
O6BTKVboMvG+qq7aOXtAvCkpNJNY1mt9aIUzcDWrSf4J8a3iohZOADDJ9/tu
D/F967rEbWeB7UmpCSYTFMoOT9LcGoCSKlXbT8Sjud2DMU2A92aRiHZijY6m
r9eaga689NQe4sy/sp+9aQW0/D69+018QtQhivENuLF/1UrBSQq2eWVzon+A
jWbFlarEAx1Rs5UUEy0eX1aFEf/fxgXTcjuZmApz1V8wRerXYtYkvosJ+dc9
kouJ2Y5Rt0XsmZiQtqcJE3sETNjM7WVCiCNfKkpceKXv3gcXJuYPxAiv48X3
1zueOcTEH4EjCua8+FOXns9eYGK8wVflCC++Xunk+2Ym/k4rEvhBfFqu/Obk
WybezLlH/CT2DLZ7tuY9E/XdykODxLLS8SN7PjHB5nxpGSNO8ul0bO9iIuXf
nSHzp8n7l6bq9I0ysUY+R1uGWGbjis/Tq1kY33H1/HbihKpZDZk9LAQcnTiS
Twy/bk1JRxYGs39LPiQeXvNaS9yJBWZhy+NiYge/azorDrAwMOX6q5RYVcJY
f6E3C81qHcbVxK/9QhmDh1koW3c9pI14ueTwpmdXWZjZrXOIf4bCuYDvHvs+
spBda/vZklicHea3/TMLlmm5E1bE6asXh7C/spA7t5/aTJxfpnJMpZMFK/8f
zO3EDQsOpo73sXDiVOxJe+JlV0ffnJ5iYey6YJMn8T8v5pmUi5vh+NLPOieJ
D6+UEpbdZYZhPqNzlcR1D4R2M+zNcNY/J7CaWMp6JNNhjxmyH8XbvCT+N+6V
4UUnM9yRjBauI+YfD3NZ4mEG9qFM/2bi+A/vykbDzOAj9uZXO/HZ1DNe/101
g2PNi+pZ4jtK82sDv5pBdKsmtGbJ+SXBnPJuN8NY+Uy2NvG0cISKa6cZkJMt
qktsPTqcuKvXDPb+ub0biIcrOrebDpoh0Mw30oQYDi8+LaCZ43vemRAucWti
7O8UaXMsYFql7iEWHOYTK99rjtypur6jxON5lxI795nDK9tXPY64z1tdYKGz
OcTlZ/ziiRt6HAasXc0xIPll6DhxSltx5Udfc9hWt/5IJFZ77u0/HmUOVUHv
W5eIbc++/089wxwyq1i7s4mzNB4cud5jjsyqmI46Yp+VvyeKv5tjnfWaztfE
OpMaQY195gj472NHPfGz53meAoPmsAv+9LWRl79tnp3HuDkWyn1vaCFeGJGr
pSdkga07Xxz/ROz28l53nbIFZp/AuZ9Y6e+71jPeFqhxN3y/YI6c79cM8/74
WqB+p9WVhcQu72tFev0sUBFX5SBEnLmlr7kpyAKHo0c+LSKWYajY346wwKlX
lu+XEotK3HfdfsoCPqn6hRQx/9sHUXdyLGDDbzorRdz2V0n+jj4LvHd/PKlN
zNqcaRf00wJMtUdZOsR3bZJmzw9YYKWlnI0ucZjjAcvmIZJfiFm2HrF4kEj7
tmkL1GmUbDMgtktzF9mylA3du51xpsRvR8UOWqqzsfJTeDyXuCkjTJURwMad
VE9+e2I6rajEMIiNrubFx3je7vDbQi+YDbN08wUOxM0r3PaqhrMRYRW9aA/x
u2M7ktbEshFeob7Ekfijq8bvsbNseEZmTf9N3LW+uzAvn41VjxTOexCPF2w1
lv3DhsRIw70Q4iX6eq+pYTZKaR9EQ3n1KhXfN3+UjQPxZyJ55r7oiO2dYCOD
f+mWMOLLzUH/3efjYCLJ+Gc4seHI5d1GIhx0OBaujiQ+rNceYqPBwUxHru4x
4rlH/g/jfDhQvvfP+dO8/ATzT6Yf5EApcPcfnh/b/HZ66s+BkMmnzWeIAwd9
l40e4kDn86jQWeJeBW/XA5EceCVHh53j1SvZVcz8NAfff6huSia+GegQwZfP
gZxf6uMrxBs12BaHRziQHlj7/CYxraWFL3CMAxdV4UW3iItCDzz1nOBgb8Gz
LTzLPj+6wX6Gg10WEZ94ntzxXMlwPhfNLKfhLN56h2PZ+AouBBmTq+4SC1Ub
vgvQ4cKJmbzhPnH1HnV3zwAuBnaZfHtIzH4l5nEqiAt+/z0yRcSVhgKeucFc
9FVzHHl+vuq911A4Fwqd4R95/rcx5mD4US7i3ro0FxMXW7wNPnmBC7lun8cl
xFmqEXHZxVx8pYftKyOOn6nN6JviIk/qsuQLYpvzYiNTxFuGpKV5llnvxl4y
zcVU2S0Znh/bzf+uQuypWyjHc38BQ82HuDrntRLPWzwKHw0S3y3j0+NZ/N21
upEZLlrc91nz3O3dLyU4y4UMrXszzwUCRgdFicMTPLfybKXVSt9A7BF/yIbn
yFN0uxBiWcGE3Tx3msW3TxKzFB668Jz3oUV78RwXFVKGbjxH+MkfW0O8a8m/
7jyLppYrGBPXvKzx4vmvySmvCOLCujb//90/Y/n0FHFv0L5AntsVLgtfJ9Zd
0RXE8+EdevlPiVdo/wrhmdt/dF498Z0C/zCeV8a+sflCvFN5LPx/z4vJ3hog
lrwYHsHz/byDY3PEUxNzR3j+f85Hilo=
{Hue[0.36934`, 0.6, 0.6], LineBox[CompressedData["
1:eJwd1Hc8VW8YAHDKKGm4R2VVkp2d1TXuY9yBkgbJiFKuvfe8oaRESSikbBpG
JVnZpAiR1bBT/JC9/d7bX+fz/bzPOZ/nPc84fNn5zNUtDAwM2owMDPSnQPR6
2uODBeqbn6r5GBjY4e02sUOONtEEn/4V93gWdqjyeZkuapNCGK4deSuxkx3q
ZVypLDZ5BHHhLMbTPOzQGNlOmae+J2xdYaPR5NnhAz703QT1MyHMUIfiZ8MO
Tff6XMeoPwnxW9vs21vYQbecVX+UOkUwwrd6jcbuhOadituZzBng0OhaUpfA
btCMJlZmndoG8xE9SdJKe+Cq9MpO2Z17IHo8mml7wx44KXb17z4yJ3ipiD1i
OcUBhR2Bnh+UuWAQ3/fOsowDVHDieUzBfMDlMJIYcBgH1NWY5Axjfoh/9CqO
wxMHsfoTNb6W/PD+gMIdKS8cVD4l/jlpww+nJ6UbdL1xwEVZUlzw5oe5ssOs
13xx0Bhr3kqK5we1Q8NVw4E4EJMS3TrWyQ9dz7j8Im7gYNyy3Eb87GGoSP2V
qxuPA+f6X8fyTgpAYV2Be1UR+j48t5cwEoCYHe0P69/iYE+pS1rORQEgj89t
NBXjID9vGZfuLABbDx652VaCg5mEHTMJ9wSAm6CW1lmBAy976fyQTgGw+VF7
paQeBwF7vCWMzI9AjkviqeVOHNwyYRFadxCEoT/KO6xncKAuw77w20MQ6m7c
lJCZxcFfZlzD1wBBEA1/XrGMbFxwwDY/Epm1U/H2PA6Etim8sHouCHN652lZ
SyjfN1YKTeOCwE44Ule3gYOFPVXEeDsheLph6pXIhoFVg7+1rI0wJN0qYKjg
x+CKqrD+grMw/KdzO8fqMAbWBa0Kpd7CIJ2x4/s2AQxsEoWYieHC8CDs6Mjp
Ixg4OX9OM84UhhNeu0N+CGHgt//IQPCwMBhzyEgPimMQQ20ybbEUAWWFD4YP
5DGoYuU+bX9BFIbA7nwJCYPRNCndnkui4AycaQZkDNhBW4tsJwpKTy/6jiIb
+zgrHPEXhZprQg17dDCYGqvj7ksSBYPNrGkLPQwONLkN6vaLgjKjVu3YKQx8
Iz+5i1PFQMicJTfTGANZDlrcbw9xqN/28/WwDQYl27pEGwLEQZicOXjeFgNN
BqnS9DBxeNuhU9qEfG6q76dFrDgw6nM45dlh4NWsJPb1lTiIEdj6PRwwKIuY
Kq2eEQd9W2/RaWcMKFssBhJdjkICsac3xwsDi1k1CX0HCTCyvdKhEIrug2No
7fCQgMTn8d7JyCUy1e6mgRJwu7yvnikMA04nUqlNlASw9cTbtSM3junrhuZL
QKr1QLPtDQxkfljYFM9JAMlOSSUiAgOGD9fSBQIlwSRi75XIaAyePK7jW7wj
BV4dKRv5jzBINJPoMoyTAt1s6urWRAzieO7fe/1YCpLq770xQo6Ms2Rxy5OC
W8U9l1eRfaJWpyZapSDWvylTIxmD00GyNYOYNHB46TFUp2Cw5WKybctDaZAP
Czx7Px2Dwih8X/sZGQgxM+8OeYHBpcRdzCVGMuB1rY+nGnlP9qDUUxMZWHB+
Isn4EvVH1a0Q58soviKiNghZfK5XjN1NBiAmtdQ7D4PUC/4+xHsywCA4PGFZ
gPpFsGxv8WcZ8EvJJnG8wcC9RM0g+YQsrLT0j0yXYvDapGFLj4EsWDt65YiV
YbCwYvCG01AWgs8YqV9G9lG5whNpLgt9Jre+tSMHlN0a9nOSBUPhL5MF5Rjc
qOjyMb4rC72n2G5T32OQUO36FOuQhYfXw5ILqjGoaMqYiTCVg3G8ZF5jIwb5
XK2ugRZy0CKu3DiLnGa9Mu1iJQfPmtbeH/yAwc0tp6bO28uBz/X9JzyQz+KX
xoX85UB56GTPgSYMxnJ0R6oeycF1W9Pvdh8xwCKmupZ75OC8Mknhv2Y0T+Tj
5bbGxyD7dLliQjsG/4ly1xWYHYNybae8CmQ3tuVPy5bH4HtRK+sIsn9z8bcI
22PAFMNpLPsFg6izyms5fsegti/AtxH5jYWSyp+kY1A8m/91ugPVz0fhnd3A
MeiI8X8h3oX6I1v2jb2dPGzaEQQofRjU91g1pjjJw8Ad+VJr5Bm2uL4vbvJQ
odODv46s67jCqOovD+PbwyarkFdka/TZI+VB53aSGP4bms/Sc2PPX8hDaPp9
BuHvGOxr8eadnJaHA8ZHy8Z/YHB3tvyai7cCkIRtIo8NYuDYfveSmb8CmPva
vNdF1i2w0qAEK8D5Dff+S8jMztsZ+cMVQCRsaCIa2e/POdrnOAXYG5UqM458
dWA8SKpIAXTbnZ0fD2GA/8wdMDGnAFH9b4SXhjEYyfX0snFTBFbRNoL/LzTP
BzZ+ZHspgvbluoVo5NroG+TffopQ2qSflI6c6hHPZReqCHIHTV81I1uqvyux
j1UEa8b/cIfGMPjWtrbhVKQIyVM26RXIHUuh4Z4rilC7ldHt728MakgxCSEh
SpA17ZYsOoHBsHVps98NJVAsU6EoIzOFD29xv6UEgQ0Fo2RkUoOio1WMEhjs
+rRhjdxE7tPQTlUCtu7apHTkdorgOHO1EuzW3JFw4D8MBnWL1CIYlYHNKOYT
yySa/1O9g9E0ZSgn7OErnkL7osLwiEeYMvzl6SiuRb4r1WZlfFMZRDJ69dqQ
+Xc1DvPfVQZqRvO5P8jQ/Ga0IEUZhmV4Z3mnMaDp3fvzpUIZFhcrqYHIDBTK
zP51dM7kb3X8LwabUMT4xPc4xBpcFoqfwcBE+1dTd+BxqOdJtU9Ffk3miuUI
OQ6Ntx2zXyDb6PsJh946Dp2/6plrkVtN1fWsE4/Dl5qLBtPIT73qYyXKj0NF
56+j5FkMtJ5/FXnHgIehFbzsJHL4/sWT7eF4kPjUz3pkHu2Dxsah6tt4CLwu
Hn0U+bHvI99X0Xj42VK3Wx65qE8tIzYeDxy/mNe1kYdTwtaMsvBQy6920xpZ
QxT3rK8eD/eJW6yzkVeVJbeNMKvAFdrkJZEFDJxNrKoXQ1Uga1Kcii2ifZ3y
p1/1pgqoj9/f5EHOHXLdvBapAsUjN+8JII860FTZH6hARK5EqizyxeDkIv5M
FcjRn/U6hXwivStXp1EFPq9IzEYgi07q3U9kVwWx0DGTFeT+a/JXCA9UoSJW
fuPDEgZC+3ZFuTxUBVqmInMrsl3ur7dPk1WhoVGG9SvybPujHUyZquC3emBt
AJlFkPFVY5Eq7KYZ5C4jS9Y3M57pVgWPLS95RZfRvLBRU67wqoH2JxOna8jY
/Ye9EU/VIEmUmiO8ggFbJ9bmlqEGN3/JHz+KzLA/qsE0Rw3chj3qpJEnHoW8
lixQg/upRl+UketT7aPaKtVAiHGyVhfZp1BNk7tfDaSGPjI6rdDnYyAn+6A6
MBt2zBUgZ+wR921IVAe/wh36MqsYSI1eWGJ6og4+fj8t5JGLSiO8NdPVYbJX
00kZud76t2f5c3WQX6+hAfJYWbbbq3J1GLNmCDRAlrAVcUj5qQ5cv390OCMX
VglaeB8hQHt1yt7nyJVuh4gizwggIijhzLuG9tP81F2+PAJ4x5NCDiLL+1Z+
43hFAP/u2ZjDyNnBlzzWSgggZrPrhQjy3ci0tLYmAuj+8ak9hmyRKbIl4A8B
qtYarukhb07EW4xxAvxkE9byRX7M/4RP7RyA4CizVgvyrluyjZJGAM/9iWWt
yMGz1e4HjQGC1x8d+4JsWT/StGkKgOXf4etGFnA46lttBdCzWNI1gJz5tqiT
7A6gEec3N4f8Ur856kwMgIgArY17HYPygBUGm1aAxdY3Vy8iL+SNToe0A9xp
KU60RJYZau9P7gDQde7/fBk5XedZ5ZduAJWONDkq8q19ZjTCAACnr/CgM7JR
fvnmvhmAwYuuf2jIU0O0jTpMAz6MReY+RT6sx7p25LwG6PeYKv9Alu/Uaue5
oAFUU8Phn8gkC1o2h6kGCNhsRg0g27kvn9u8qAFJKrb9w8ivE8df9FprgOu5
fPdxevzEZ4t7XhrAvnXSbIkefyehZuOBBthIsWtwbKD4z2K3ezo0QO3+9TAC
8t0jVRkrXzWgIqacUQPZ3tu4krdHAy6kFQRqIgscCp83/64BeK12VyJytNOw
xeCIBpiZYHp6yHY7U+THFzTgezet2giZX4/zxxqXJkS74z85IEfWb8jwm2sC
19LEtQfI4Doqe8BCE1p5C5/GIc/xtsjxXNIEzW8xlfHIZq7J8pxXNcGBNXf1
IbIkn6ryNkdN2PVak/oYucXVlzAdoAmJm8WcWch7DsydrEzShJUM3+a3yPfd
f9tZftOEjRNxvV+ReUh+rmd/aEIPV39VF3Iq1w4fUr8mEMY8s7uRC8slrksM
a4J8Vbh7L3Irq0vK0rgm8N87y/gDeVfSwpe7q5owN7N1bRj5Zu1WtSoeLTB2
9nCdQQ7Ye3C3wAUtSLz33GfHJgbNBdtNCKZa0LPBS2ZHPqg/n25mrgVZkt85
dyK/v/EJH39JC27xnn+5C3nLkp81u50W3JQ63sOBHN7bVb7gpwUjBvK7uZBj
Uu45fEzSgvCZH5tHkHPEmJs8+rUgIvoUBx65nU9j1XFQC+6G6ryge213oAR1
WAssYhopKsj6C3NRF8a0AB96PkgVea56+Kz6tBZkKNQMqSODWe13VgZtuMQk
lKCF3B0V+vfRIW0QtUjKOYHMMsfIXXVRGyzDq0LMkZfyE6KGLbWhXIO28yLy
uKM00zYrbbgc3hhPd+svsyl9qjY8pizmWiA/6ntb981ZG2qSuj9eQpaqcXRb
omlDgRbj0lVko5iej9Jp2sCuEL3miJwlUxD0+Jc27A4uawpAdtr7d/ntb21Y
PF6mHIgsvyLj2TauDebvLTLprqzJt2ea1oawi/dpQfT8jfKN7Za0AX/mjQwN
eVtgnpzidiJcNxemhSLbfHgx2ixOBLkn1d8jkMUu5+qvOxKBGqSX9gC5MRmf
P+NMhKcXAhbotu5p4hhzJcKEoINOHHK6wXhHuycRdp+PmKSbnyBhmh1IBJ38
9woJyPv4XlLP3iFC3zi5+BG9vl8LaDnPibBFjy86BblP512h4TgRKEzbo7KR
NU+lG3v+R4QnDDLddOeei96InSJCMR+XQA6yn8VV3Y5ZIqhHSLyhm8eTY/DM
GhGE1u27c5GNn9hyGOwkgS6bO/YC+esCt4uuNAk8H7CYF9D7I81PkuCOXNdP
eYeMYyh6h/ckAeuFHzS6z5r9JSp6k0ArIa+Y7g5Om4uS/iRYymgWKUHuum4Y
zRtKAuGW1a2lyN+oMn8XY0gw56n+ugx55Ojo6/xCEhhXU+bf0/vj1WlVgRkS
UA18GeqR2ZUVW7A5ElhHBknSzV/GY8m8QAKPr5YmdFNqh0LHlknAz1j+iu6H
HZ4fXzKSYTn165UGZPz8QxMVDjI8ss6sbEQOUBz0OSdDhgJJJ4uPyJvFbm9u
OJGBk7iH0ErPj6XwdqoLGRYekqh0l5z7e6nCjQy3cKbRdHtMO+9a8CJDRf/B
n3SPiThSrwaT4d32vYFt9P8VR+XWvksGh60dhe3IGR5mgYyFZKgKNmXpRNaT
IRED5slQ0qF+tgeZobOT0WORDGeoaa50F/lerbBfJkPe/ulougVqwpRM18lQ
xKf5ie4VwxoxPDMF+quatHrp9faHXUucFOiJD5LqQ97egO9yl6eAUtDkzDfk
BnNpW3t3Coy6Far0I5M+cdvd8aSAb2alEd11eCb7PG8U//mtK901+3scZv0p
8Fj/Uhbd79tCXPzDKBAknMYxgPyW+NX79gMKqO2eHaA7SzLwxrO3FEjmJbsN
0ffPelPa+CoFujayDEeRz8Vyz68i19NOGdHNf9SGxL5GAZeNhX8uMWb+LYE8
P0E0pnviFUHKCbkhbsiEbgO718XTyA8PHLSkm6cruXl+nQJVS9F2dI86Thxk
2aAAYUbRnu5XTCou+5D9f33/5xNy3TglZOm6o450B9/BGfsgh59scKZ7WCt8
cAV5Z8y6B935vZ3HdmxSIICc7kl3oKvgdV7kZ0u6XnTvS6kSUUXW0U/wpltn
ZdUhEPl2vZzfv/N7uhV3kNcu9PzzoMjD3Y+RpX8H+9MdYKhYWIHcvfIpgG7K
RNjWz8gnA90D6d4b+uXcT+Trq9xB/97nFsicQg7yrPznl/kui5vIyuPWwXT/
{Hue[0., 0.6], LineBox[CompressedData["
1:eJwd1Xc4Vv8bB/CHKJX66jmyooTs8diyntt+IpL1tUey9yrzi4iUiCSJCAkp
I1RERlllZlQqsjIqW7bf5/z+Otfruu7rvs65P/f7c05e8DJwoCYQCPZUBAL+
5Erczn1wvExZsk5gL4FAhBd0Aic8nBPJyaLcti5URGgIfJbH75xFNrRTG+yh
JkIzycdpr3MJuepV8vFsWiK0xvdSVpzekHPYitglDhKhTT7q1S+nLnLMizck
YSYitCcN+Uw5DZMjxrfjQoSIoF27T2/SaY4sxlZ0O9OQCOK7F1aY5ufIltxD
SrNGRIjntOnkVpwnl5XUbcqaEMH/QnWA0dt58t7exxWdpkRQn/R+W9m3QK4J
pggsWBFh9Pc3u8CVJTK7QH3iAWcicG6/yNiWWSOPGuV8Sw4hQschmf00VgSI
appOPppNhF4dYlpnKAEcyi6Wsz4kwqe4P7z3MggQk1N5giOHCOO0j9VFvxKA
X8Zz4mQeEba2mSNMLaiAcUydn7uACIJ/1leLzajhBvHqGE0pEa511o4bm9DA
mwXPGa03RFBN1Kh/fI4O5JvCCn4MESFCu6Kw15MOiIXZuy++EqGOlvv29k06
kNRZbrz5jQjyoQRHww90MDFyM1F2mAiSLjUHdyn74TTGox09SgQeNfF/TVQP
gI0TJFFPE2HvX/Y/NNL08CkyoyllhQgOYhuHxA8xQFpxowjXIQzODsU+YjnO
AM6xmR4TyFKxR5UIYgywuqnOUnAYA5phknuXPgMoFySWCDFgkJPg1O6RwgCq
9NwLQhgGw7/7YoqOHYGHOSwsR1gxMCsu2eUWJMIpp75Cbx4MdAUcFpi0GOHC
jP33e6cx2LgRUPPNhBHSqNuYSfIYFPy5ejXPkRH0G7fH3iHTVOSzSsQwwuTC
X+l5BQyqladUdN8xwtlpYitZGQNeI7fkKLWjIJIZM9auigEh3EdynswE/YNd
48k6GJT3hQW0ybHA1ht66nUL1P/qxQujZ1jgoUFiTIQlBg9kdM5tmrNA8WDb
g31WGFxPYxEQCWOBXN2CG4zWGNhbPv96q4EFhNjf2wjYYsA4PqVmqs0Khf77
ZMkXMbi0ZEj8acEGXTnQyOGOgQJRsIQmnB0WZ45VVAVj0Ell3taRyA4LLImh
3CEY2C7EjaVms4Mc+50/icjR3dPMAo3sMFZv2e0UiuoTCiN0aTkgm+lnEvYf
qj8oYJB6gwOibzWqW0Siehr+Fb6045DWPHziaSwGTpvJmY9MOWEql/LJ6TYG
KXq/moJsOSGl23euGrn+ocaMrjMneO1zbj2UggELZU1m9TIndNk51Jcjt6ZY
dWve5QQhkSsd63cwEBDl3zPVzwm+7OlXQtMwmLWtdRY0PAneZhJXjTMx8Gr+
KVmiywVeLcnr1vmoPxS7CZtwQfxnk4xCZIYa79xCay5w2TlEvYxcWrJOzPPi
gpgPR5mvPcZgMe3gYloSF/gtvXN/WoDm5yZWeqWfC1JGD8z/KsIglOGysIkV
NwSZkNOhBJ2H+d5T2+480P7uirj6CwyUSfSr0/48kDxs+icUeYGW2DIQygPP
OxLCK5FNyzhcSuN5oOtMxgnelxicopN+al/MA/21LWa0r9D7VtpLt8/yABbj
2FtTjcEqQ4PGXddT8HDjYun+OnTeLSGO4s68MOacXyn2FoOLirx6q168EL+Y
nWSF7FjWLV1zmRfIQl+1byA73z9FqxHLC8U2IcaTyJ5eXbmm+byQ852vM/0d
BsHM3D/Cx3mBd8nEfrsZg2SndotOWz4gn7FtL2zDoGEf63k3M36QGXldm9eJ
wWSuqPZnO37QVTz4owOZHtTVtFz5obbt4ugqsmmglzR3CD+c0XkcTunCYG7q
HetQBj9Q0TryzSBztPuOao/wgyGDoxpfDwZB8R/8BJ0EIODAP45xHzEQPxKR
Ou0vCMMsup8XB1Ge6Ab5W0IF4cShxsPETxioEkRr8qIFoZDMJCyObDQ3NGyT
Iggdbx4d90K+1CErMPBcEPiHbCamkV/HzdU0LgrCae4D+4c+Y0Chtvlx31sI
emJ7eYqHMLBZUhLWcxeGLRcru0PD6HuIhO4+f2GgLSpQFkCuJjX6WYQJQ+dW
HI0GMqOnZo1zgjAsXrLUDkVundLTjioVhtI0y/ppZNJ3G+eXy8IwkhPgWD+C
8t4WmccVJgL6ObocVqMYZD94x/73pijoP3zAbDqBwX1L4UHjVFGIDDMGT+RU
tttJFQ9EIZ7nlUk0cnyq7V7fElEIamw8W4ocmLA596tbFMiCETb7JjE4/594
0ygmBvy80cLlyNTWmS6d98TAI+t828ZPdL8kyA/1GpBAZFP+sMUMBnb3D9NW
m5DAp+XZFzdkhoJR0YfmJIjc5E0PRfZsuH7F6wIJ3lifXcxEFlz+IkDvSwL/
tnC6EeQcs5BAjSQSxOYy9dvNon3heX30ZRcJGoUXxS1/YeBXraSfeVYc3LKa
NkT+YFBh3kL9WV8cbo3abysir27oVzIaiwNDheuSDnKgwkW2eCtxuLLZ/8oF
OfT19fFgT3Fgptcqy0OOqRsMNL0lDmH1qYYscxikNfo8xPrEgUvJoHgVua79
0WKchQQckN7XkrSA8szS7RNmIwH2xx6ezETOddyY97aXANpYA78C5GvU5+b+
dZMA40rmrTpkQ/m12VMhEkCu/JU7izxVqD3RkC4By8eqDdQWMcDi5gbXP0vA
rkOaxS9kZ63TtS6mkrCPOyucdxmD3/ys78osJYGZ1JVOQvY9sP5h3VYS6lod
nsojh3S8/BrnIgm7L8de6iEnGMptFQZLQiDtb5dLyJU2sgozGZLwbG+A/Ftk
6kDpV64/JKGEk6vZbAXtR4F4pZurFGyOXNv0WMWg+bN9a5anFHxu3UO4hLx4
IHXoo68UnPFp3w5D1vbYoFIMkQJlk/MTN5E3xJv06OOl4HTJyMWnyKY1RlPF
T6VAI7tbcBaZqfPysT/zUmC83/6v3V8Mbi3VRnpflgZ/hvxSpTUMPHpv2VmG
SMMpo5oGNWTtMnsVSrg0aDAydp5BpvXaT8UZKw0MzQ8+GSMHzxhFdKVKA/2l
8EoPZIcfs/+JVklD72ZaUgayfBdr6K9laShjXppdRZ4oCrjk7CsDFn0vJx+s
ozxz7HwvuCQDcVTm5/KQ3ybGaE0Hy4BtsePzQuQc/7ssrlEyELjQ41GBbKv8
qtotRQa+1FsUtyF/7dna8aySAWPBBMMl5L61qNiADRlweW9yUW0DgybN5LQr
V2QBXD4aDSKPO9Z0BMfIwgzfI9chZJrYcWq/67JwL1EhbBhZs0XGwz5ZFt7G
16f9RG7XGlJRz5GFJHHfwr/IvRSeWdpGWeAlP5Jk3sRgVLtKKY5KDg6LemwZ
IhPOfRlNjJCDwqmA0bfIrXXG3P7RckD+EC7UinxLtMfe9JocmPBv+7xH5jzc
Os55Sw5URM3/9iBDR+VkWZYcJJv1jX1HjtBJmvlYJwdMV0I01vD+FMoi87Yc
+PN9kRPYwmAXqqiyg06D2EFjjWhkc/Wf7Z/CTsNgSYxULHKFFkvKkSunIXqJ
fPI6srNeMG/U9dPwIu/8UiJyt4WyjuP900BwEgq+j/zwUnOKcO1poI6PYCxH
Vise4HtFkIcyQnrcN+RY5r+6vbHycCgrYlhkG90Hra1jjTfkQeCYATMJ+UFQ
etDzRHkQNHfQlUCuGlJ6lHJXHqjvmVTIII9nRW+ZPJYHE77rHmRkFX7ik6Fm
eSi6fCFBH3lTToRuglYBpGbpbvsge5nbN/6NUgAjVo8/z5Dps2ZGFK8pAHN+
/89S5KIxn93IeAVYKkkcLkeedI9QpL+jAKpL221VyNbhmVWc+QpQuDNzvQ75
bN5g0ZlWBXB0V6nrQOb/o3P7Pr0iCB/1X55BHomUuki+owiTb0d+cu6g/zXT
4QTve4oQ7lcfyIXsWvTzxcNMRajs06fjQV7qTT9Ik68I//xzjZsPeS8P1fPW
KkWQ7vHUFkEWae6gMvikCBl7DlucRg4+4JR18ZgS/IqaStBHxm7f+xL3UAnM
RecDQ5AP9GM9vo+U4P2oTEEoMoE5ocWiUAn2/Lk1EIb8K/1KhUiZEoTMRJIi
kJtz3BJ66pVA/prb52jkwHIlVdYRJbCb2dy9uYPn40dhwXFl2DWYzstCfsQg
GNRyXxk8y3Z865FFJ83WaLKVYcEiKLIBuaom7rJqnjJIcv+T2Ij3d5wOqC1W
htJDiY/fIk+9LvB9XqsMEv+udLQiC7vwuWcNK0N9F9VCN3J5A4/NZW4yNMcI
3hlBrvc9ocH3hAxl/tqh28gbK3O32EvI8C3tDmUHWSqo/uuR52QgG5IZd5EL
wu38t6rJIHattJCwi/ISn5vb006GjPfMHXuQbfL5qENnyBBz2ml1P/Lur7s2
U4wA/a/tBo4iP+DMZlcyAtiaOj8hhHz4uniriAmAyeNHccLI4UuNfsdNAY4c
CxYRQbZtnmjftQAoKs73E0XmchcKarQHKBBh3SAh57+o6tfyAxg0ezknjfxM
ryPBIBlgoy0ggoxcG7pBcO4GiH1fpXMeebVkcv5KL4DZu6TvuEljvSOZfQBK
5xe8DZDzzjyp//gJ4FUtOcUQ+TqTZQT5B8BJ4VP9xsgmpbW7TIsAvmZuYI48
Nxax8w5Tge1r999cQD6ps2+L+18VeHxVJsAXWapfrZfNTAXCviqO49a0iSg4
YqECLid+Gfghu/qtG+1aq4BvsKaYP3LF/dmnXxxVYGtuYSwAr//VZZN0SQUi
Iq2VgvD6m2lNO3dUYHrEtDwcr+8SuPG5TwWqP0aV38DPg7vh0caACiyqOR6N
R3a7bFp/7LMKfMtZvoyb60TsitU3FbDSZlK8iZzoOW4zOqECNNQB9Ql4/0NZ
UrOrKjBk/e1VEjKnDuP3LRZVaDsu53MXOb55h8RppQrHfc4RcpDBZ1Kcw0YV
dgsXtXEvH+uUYLNTBd1s8zu4LX0ypRgdVEExMkggF1mEXVGOzkMVvElrennI
nT5B5PlQVeh5cjIxH5mBY1m3PkMVbk3bThQh3/abdrX9qgqPZ7YWy5HZNIN9
DL+rwkL9muhz5ByWg4GaI6oQNPHOFXd5rfBV4XFVuF7U9wN39z7vrLVZVXBj
bu2owPctY/XjrU1VUKk6lVmFfO3tHqUGNjX4T76UrRo59Ojxf7jM1KDjvG5a
PXJH2X5zsoUacJwZ6MB9XG8lz9JKDc1JYk8D8puYD/J37dQg+s4FD9zUa8GO
9K5qINr6hNyIHPtlsHY1WA3U/XaHmpCTs5Lc32eowXWrJ8vNyIUCtO3+I2oA
liVPPyD3sqtseoyqwfHRkwO4t/4JE3YaV4OaWP8d3HqrywlmU2pwL7VarwOf
d+O4ofK8GliHffiNGyzffttHUIcmag7eLuRPCVEL6SfUQTWDNaIHee8yFWuD
tTq8WIn73o+8VpqWMG6rDjzS0dQDyLMeYjR09urwssaGF3f3T8s5PSd1kNhX
44E7fejFu69e6uDuTbuJW7TJw3ctQh0eHyQwfMLzkfz5vViuOrQvKvB8QX5M
KvvvwU91EO3au/kN2fPowvqLaXV4lBjP/B3PxwYpoGdWHSwcNyRx1zeVutHM
q0NufoI77k8mpaaua+owkHX/K266sBIJmf0awFpg+nIY2bnt6WSHoAYkG09Y
/0AWuFCkt+2hAf+eoLiPI7dmypcuemmAuviXSNyOn9uPTPlowLK9xV3cefqz
fb0BGmBuK9SIm5MsbFEQpgEyPHuYJ5CZ2J85Gd7UgEDRzFrc1ANlEYXFGnBn
8ePOJPLQmVflxrMaYPZhwXAaWfVcnmnAbw24/EHICXeRUeJOypwGFP1nHIw7
2MZBu29JAxJIgdm42QKOjBpsIWtvzeI2zXY5on9IE4IkD0bMIA+ssnpri2nC
U8fsjFl8P3KDRch+mjBwxKXoNzKRUPVKPkAT4oJ/v8RtaLmgIXNZE+5fdGzB
3cfobC0SogmKP6THcQ9eNU48FqUJwUFnOf4gf3UiLfxN1gSOvzY3cU8ITVaU
lmvCjoqVwxy+H8/PK3ItaoLOeNTqPDK9nEwntqwJ3uwk6gV8Xq/ZbGlXNWFF
/eMh3JS3Y1FT65pAjtrhwX2vL+D9MyotcAycNsAtv3LPXOGIFrDpMz/FHSoz
GmhE0gKJNSqLRfx/8dK3MsZTC5avhz5Ywt9vb/mNHG8tsE6dLsBdbbRgV+er
BdJfdJ/j9p/3Orx6SQvK+bdacE/xeTg5hGuBw7P5edy9qU6s6re0wKqCrLKM
/MjfMoyqXAsUhB8P4dYhaWqErmhByxjj+goyob+fyv+vFnQKrBJWkauCHOrc
1rXgqnIPHW6upmhZi20tCPrpz4J7w7hJQJ6WAh/qfGVxF4XA4TVGCkDZWX/c
+1vkB/2kKJBtlD6Nu8VKzMXNjwIi5qlNf/H7/QOr680ACsiyprfjfidP41Zy
mQJd3+/24G5i/uy+FEIBuuTIYdxveq54h0RTYOvU4U3cLzQGLt+4QwFROwbx
NTxvImExT15QYE1iKA137HZ77uwmBWJjbpqu4/MzpZ0W3qKAsayCGe5fz8mi
nsjuA1P/t75rxct55MhhNQvcbIOZHSvbFAjgWLfCfVbiE1F2hwKtkG+NO/wm
0TQQ2cfA0Ab3uFrs6AbyIDy1xc2U1cCnuEuBVA5zO9xnNjbdw5Arp/ZewB1q
LFNehyyd8/z/flbq/XcXmU3b1h73/wAElT9/
{Hue[0., 0.6], LineBox[CompressedData["
1:eJwV1Xc41W0YB3B7JaHfEUJC1it7ZZ37cFZEvRmZUcreI/sllJIVhZIoZCQj
QhmJQspIQmlIiCgUQuF9fn+d63Nd93Wf5zz3/X3O3pO+R08zMTAwnGJkYMA/
JVI3Cm6KVRt0FoWsbm0RoJ5Dbo+3WyoxlHEkk7pJgCehFYWybnnEmgVsMOUv
ATqU/V3Z3CqJqfUSUlJrBOhKGqAvuz4mVnAzfjvxkwDPdeIezrn2Ee2OPkqV
nSRA9+VR/2nXT0TKu9YXQy8IcCE783mGyhjRVOYgl1k3ARTciLZqlWNEGcbm
mo4uAgSwX44IKPtMPFngR2p6RoBNinrLQv4X4uLFarGGxwQQbA83+p78lTgr
tpizWEMA42Z2synXeWI57H17IpcAKlsnlwUW5onKFFezqRwCJIk79krqLRAV
Pm8e8LpOgKCTj4Itni4Qt6/d9o7MIgB5yu/pg8FFYh5XvVDFZQKMf/9wInT5
F7H5VLWk5XkCiG/U39jQXCUmDA07FPsRoGe7JieLAwNQLQPjtx8kwIAJf3Zv
JAOc0F+6+5tGgJGLP6Sv3WAAu6uSL8epBJhgLSYrvmeAj7EPfjaSCfB3Y1eM
tR0jdOX/aQ4DAsj/WFspt2GCUZpxhIQ2up/e5glLKxbwsek7ICFLAMNUSmvx
YQ7YzByj3uUiQIxxbemADwfMP9rJ+YqTAC2skhkbyRwg+vDu+d8cBNCJZHAx
f8kBO3+bVlPZCaDm3rhti84JgrLHq78zE0DKSOWYlSEXePGIdjpsYMD2W+QH
iwY3VE9gT54sYHBaaX27ynZeWAzge3PmLQaHRhOKBMV4IVrIM7V+BAP1BII+
gxIv/Bvdt7o6jAHLJ2WvviO84NXAMRQ1hMHtFNdu7yu80JSdVZz6GoNP3wfP
l+3mgyOCTtjrHgxsyiu3JOX5YcLau2CwDQNTudOLAjQMrK08dW9VYLB+Kbjx
gxUGlxR/MO9ELvlx7lyhCwapt3jI5+6h76u9I6R6HgP3oig5z3IMHhlMk0yf
YbCqFhOlV4aBtIVnepwRAUK0X/xYKsKAIdpfbYEoAOmND5Uf5WJwfzAq+Lm2
IPyR6NM7m4L6nzt1cvygIAjuSDdgQb6paXL4j60gSH2fqUpIxiAxW1Buf5Qg
6JdVTKYmYeBsX/M+7Ykg5FC24ooSMcAmpo2sjYXgZE299iQ615lf5vxf7YQh
/tC7xzXRGOjyy1eyRIvAbO1MUGgABr2Mts97UkVAYSzEVQDZafHil8x8Eai1
G6qs9ccgvn9ml1ybCKSXv4Cffqg+pTTGlFUUGigidv6+qH6b3NHMS6LQN+in
GO2F6llkl2WyxUDR6SVtGN2T65/03CJrcbB81Dvtbo/BFbO59jAncdCaJjmK
I7feonwzdRMHDRqJb9gOA0H6quZKiDh8sBKVoSB3XXHop2aJg2LOZ0zKFgM5
RVnm6Tfi0DU61jN3DINZp2Y3efO9QF4qcSsxx8C346tapakE3MowOd5igvpD
uaeClQTc6TmnE4XM2+hXUHpcArinrxzWQ66qXOMv9JWAH2wuko3GGPzM3vYz
+7IEkFwPeTUeRPfnqVQV+0YCWlovHOtAexHJG6Jg5SAJnu5KwEJG87Bl27fh
JQX2o5OmAvoYGChzr8wESUFOhPepST0MFln5O4cipeCO15XbtcjW1aLuVUlS
sMilf9UceR+Hxj3ncilwmXM7kqGLzvvAWaN7Vgq8MlmYhHUwWOF9Qsny2Adx
o1VnjLTQvDsjXFTcpCHItVeMQRWDU3rSZiu+0uArn7U0ooKBS3W/RmOINJSO
Ws9UI7vl7GOlJEgD2yFnpVPIPr59BdZ3pOHC04/y3coYhO+S/Bw9IQ0pfmzd
OUoYpLt22/U6yYCByzdGy/0YPGEX+tfTRhZUkp84CMthMFWgaPz2hCzEMWVE
/5bFgBvIRjQPWeiPUX88iGwd6qshGSELWZIxWanI89PPhEZvyALLe49SNmTR
7oBx4zFZeFHrqLcujUFY0stAeVc5YK6Z9F+RwkCFLyZzJkge2K/pruvvRXni
GJbtjJSH1lN1NZLIhgyKjYXx8iAsLBjLiWwxP/rJ8Yo8TCwfcngjjubVoyU3
VCMP7S4l23yQmy7ON7b9lAfy8/sHCvZgQGdy/Jzj9w+EftLxEBHDwPGXvoKZ
lwLY/3e2lbob/R5+hv7BIAXwuu3iroT8SLkt0C5KAQRTHKUFkTEfaqNbigJ8
2aS+nxFG+zltZhxXpQAMoeyaKcjKHx3dGpYUwNDGNv2dEMr787OFElH7gX/r
fthZQQzybz4T+Z2sCI6TVsksAhjk2CsMW2YqgtMa186fBAwyhTMu195UhIyq
f0s+ISdlOrEFVCqCluRJxkbk0JQ/83P9iuB6dI9XAPK//6m0j+9Ugn3Z1YIT
GAZMx3Pde68pQcTDiuCBneh9SdEZHTiqDGUHo3NH+DA4kcPD+shKGd6dO1zS
jcxbMq54y1YZ6NuxliZknyeJsb4nlUH09DGuW8jyS+/kuAOUgZwyJOWBfNsm
IpRyWRmMJmi2m7xoX6SaCA19ytCY9ClJETnwkf6R3EMq0Ot+xLKcB4Na206m
t0dUoLZPfzIPeWX9yAPMUgWsPzRHZyCH6p4STnJQQe/Mv4PhyJFNiRPhPiqw
9KNn/SDy+ZbhUOs0FcjHjml8245Bdpv/rZ2DKlC2oS2litzSXfTzop0qxA8d
Fv24DeVZsN8/ylEVFD6VPBpELnBZX/BzVoXHTgGOL5AvMB2eP+apCkF7iY8b
kM11Vmf3RahC6pTwcAbydKnx5JPrqvDC/NlJE+SdF+eH196qgqKAdNtjLpQn
2oFmd2s1eFmf79/EicF3WaFn1fZq8K5sQbgWOYBr7eWakxoUWLu9uIsc0dPw
/qK7GvAb1ZFykFPMtf+WhqtBE7Y3Lhz5gaOW7rcbamDE+M1cB5kpVOOhx2c1
oGXwszzmQPtRovLA00MdFDipp0fZMeh469yV56MO5+PiuweRf3Jljr4OUIfU
k0c1epGNvdcZ9SLUYS7x4u4nyOsq7WbcSerAxPlFrhjZutFiuvyeOiS16xoH
IQv0huz+saAOS9UC4fzIab+az/qFaMAEm0yxAxsG3gNpJ+wjNKDX07vUGtm4
2plEj9YA4tublebIrL6cjOIJGvAyvK7zIHL4N4uYvkwNiBZxltZEPv159j/F
Og0Y+bgtZweyTp9Q5NySBkgFfD30jBWDybLgM24BmtCdkaOui2whuvmx5Iwm
2FpmtmkiP009T5sJ1wTmvTYWqsi3g7IEPeI0YXFJLk4O2cng4SPPK6h+4+Z2
QeT3r/5u+tRpQv93/tfLLBgMrsYlBK9rAvv2y1y1yO3U9OzYWC0IKmxJ0EGe
cGnsCT+vBQJh8wpayCwJE0yBiVqQqdQ9qIZM7dT0dk7XgobuIIX9yN20URL5
thZkv3b+swd5gC41y9qmBa5hwwQ25HHjOv2LjNqwRfmgOMiM8n/43XhqjDbI
JPfGhyB3tVhKBsVrQ0BQ7kAQcpriK2frC9qwEXROIgBZnKdrQjxNG5a1r3d5
IUPPg6nqPG3gtSrdfxI5xuTyt9ct2vDGYsXGFO9Pp//ctaENN0cOK+9D3oI6
xvywA/DFbnvMOyYMbMlfu0eiDsCxbbO7RpBraYJX+GIPQEVVf/UbZDezcOm4
xANgynhvth+5387AxCXnANRwhod0It8603FFofkAmEzOfa9FNiofknnIoAMS
3xqMLyMn7PptOpCgA9zsvl6HkdO7ur60XdIBeuHcdVPkm2HXw2pSdSBxYKTb
BLluVL/oSpYOODiuKtKRJ/Li/1oV60D1iDoHCZkky393tEMH4sXvbqgi/9He
zzHJqguxrGtzu5B9bZ3bfsfpwvh7tYJJRvRe530b07ugC57ljKwTyGVf/LfO
JunCpzFlt3HkKa8YPe6runC9rUL1E/Lx6Nw68Tu6MFnMPjKCfKhwuOxgly6k
qfzj/BJZ9odJRg63HrAKWQXWIo+dVT9FvKoHpumqxAvI+wR4Uvyu6UGKRsfK
eWSPsq/1t3L1oHvGq/Ic8q+B69tY7ujB7Oy8ZBwymxRjTVedHrw/XCP0H/L+
jh7GoyN60NHvrhWEHM7lmndqtz7EfkzWOIG8M+Pau4u39CHAvZXHAJnrzc5X
AUX68NLdRkAfmWFXSqddqT4kn6kV00Oeux5bu79aH/Yeq1DRQe647ZnyqlUf
RBrznTSRQ+/rGwqN6UNx+8ZnRUY8H59LS8QMYIT/po84chGvfFhnjgGcrOK+
zIqsOGWzypJvAFZeyTUsyHWNF0MMCw2gsu35EDPe32UmuLncAIjtAXuYkKeb
SgJqmg1g10DQw00GDBTcZbzyPhnA/ZdVO1aR7z+RcgyRJILG7jGYRW4N2EOR
uUuE0KjgzNfI68vzaSKVRODNV54fQFYPa33PV0OETmY/Ou6S6BNBfx8RIXou
7G8/clpSQcGrbiKcn6rw6EV2vCPDFPmNCAmFRo7PkbfmshynUWywFZs7j5Fv
iueL6FsAZKSSs+4h8ySqdO23AuAjN/Pijv7VFihmDVCsfzOxHNmpY7J7yw5g
qSkr+i6yhNc/YW3OAPxSTj6lyHfq697QAgE6xRgdi5ArzHpSjqYDvEs9lZiL
3By5zuDWD2DXFLKWgrxSObUQOwDAdpR0Frfyl4Gx3EGASxu+nLgLD95tfT0C
YG74RzgZOVHAPob4GeAjrdXwErJVVfOWwE8Al6ONhQnI819iNp/tJEEnj15V
DPJeE/a/ksdIcI0eMhiI398bowFhGxLcOpcVhpvqGFPCZ0eC7h1bYrg9Atcs
to6T4M0/BR4ByLU5s/feuZBgyucAqz9eP9fnePkMCV7u8KH54PXJ2e2bV0mg
8D160w2v75O79HaQBOnsPAqO+DwknxStD5FgM3HP6HFkzxDr1t1vSbBrPDgR
t8SehGWHDySw2Lsw44Cc6jPhOD5JAp6LWLk93n97nvrsCgnond/0bJHFTbCP
fwUN4a08JdUSOaljU1ncwRD6wrhmDyGD/5SKqKMhMB9Yvod7aXevqvAJQ1gK
5/XHbe+fq46dNgSe78O/TZD3i+hpc3gbQo/tHXbcvf5hxIVIQ/hAFlc+iMwr
umTaesMQwsrOFVCQMwJnPJzeG0LTDkduIrIwNdzf/KMhGMXRRw2QbwtuC6WO
GYKpknUZ7vvNCucUJgzh98MZOu5+dr+81VlD+OswkKCP79uNlddpfwzhnlk+
tx7yhafM+k+EjeDtexeVA8iRBLEdEjZGYLla9FkNuaea05ZoZwR1I0MPcIuZ
LRfaOxjBtXqFRNyPz7/UyTphBOxLVDXcTKvhLtweRqDBrJqgipzwbrh5JdwI
Zu91aqggp+dd9npxwwhS2MXvKSKXyrF2B40ZwdEGk0055AER0h/vcSPgv636
AfffHVEKrhNGYHxbrQm32cpSis20ESxkxIXjXmqbMDdYMAISg82qLD4P+6cf
2BnIIGTetSqDPJISt3h9DxmuditySyOzLTEKPTlOBvHxP6ESyKtV2SkTTmRo
sb5kh3vWW4mFw5kMj36pGeDu/2o/b+ZKRv9j/cy4r4/WP3vvS4alL6ppe5EV
270DVmPI8Gss+Z44no/0ty+UCsiQyrJtXQy5WLn6v5tfyVAo/f7NbmQfwuJa
/QwZPLzPteBWX1cOfjVLhmGKYTHu1vYqT5YFMvB6rITiHrGqsvZYJcOPeqoo
bo6oSlVNTgrYWIa5CyO7Pb831SNPAcbqhh1CyHIny8w2vCmgd8c2TQC5K1en
6qcvBTSJtf/hdnnbzTftT4GiO7u9cRcemR0cCKZA04OdxrjFiQp2JVEUOFNS
zIJbQKTC1TyZAj3ZlpEEfL5D1TGl5RTQEjwbjCGPHnx433KWAv9wziXxIxse
LrQO/k6BHdXKUbjLLFI3r8xTwLUq1ht3uONp48FfFNC5YGaGWziYb/zoXwps
nojnxW2d7853ZDsVjnA6ZvEhD60I+RkrUeGhW3wFL74fBeH7iYFUCCo3Y9iB
zM9Q91AnmApyA95LPMjm9osUzRAqNDhmTuMexNyO74+gQv4V5le4h89Zpu6O
o8J/maa3cb93VV78nU6FzivXqLgn/5mqrbpPBVL4p6vb8f2o+VdP4icVYLrQ
lBuZW1uzd+cSFXindpFxizcJO7GuUGHKNF0HN/3pl7jpNSrcJefK4L42GPyi
gpEG0YrbmHDrLF+z1eWjgRmjYcM2PG+a46EWyjSoXlmSw73VEPDgvA8N4tzp
e7nw87Hdv3Tbjwav0lsFcT+yWDzREkADxVkiL+6gBV+elTM0GOe12OJEnpbx
dj0dTQPswqsPuAcyXYXIaTRgXN59A3dRkH0U430akKgdorhNlKmUyGUanFHX
U+dAZnjzhjHoNw1yS5/vx10XdrrFc40GM0r2Mrgl2uO17DZowOyUJox73bJd
ToeVDimf1Blxl0UAzypGB90Fjz52ZM5OneFAdToEXPPzxd3poOTuGUiHtDJC
Mxv+vr8U8kgOpkPHSet63M90WDwrQ+jgsplXjbt911uvXxF0IPbp38H9+FWs
X0Q8HYqda1Jx11OGQi5dpcOuZ7HOuIv3R52/W0+HlqRMHtwJG90Fs3/ocGVb
iz8r8lwNUdHnLx2kxm754T7iUduwgFzw45wvbuHh3J7lDTrQHpp5445O5rcO
3aTD0WA1L9wTRgnj68hinoKeuA+u//GK2qJD3vUNd9wVVX6/t5BfMH9xw/0/
AspectRatio->NCache[GoldenRatio^(-1), 0.8948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox["Pop", TraditionalForm]},
AxesOrigin->{0, 1.},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.626*^9, 3.6255`*^9},
3.6257`*^9}]
}, Closed]]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell["DSolve knows many ordinary and special functions ", "Subsection",
CellChangeTimes->{{3.*^9, 3.*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"sol", "=",
RowBox[{"DSolve", "[",
SuperscriptBox["y", "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "==", " ",
SuperscriptBox[
RowBox[{"y", "[", "t", "]"}], "2"], " ", "-", "t"}]}], ",",
RowBox[{"y", "[", "t", "]"}], ",", "t"}], "]"}]}]], "Input"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"y", "[", "t", "]"}], "\[Rule]",
RowBox[{"-",
FractionBox[
RowBox[{"AiryBiPrime", "[", "t", "]"}], "+",
RowBox[{"AiryAiPrime", "[", "t", "]"}], " ",
RowBox[{"C", "[", "1", "]"}]}]}],
RowBox[{"AiryBi", "[", "t", "]"}], "+",
RowBox[{"AiryAi", "[", "t", "]"}], " ",
RowBox[{"C", "[", "1", "]"}]}]}]]}]}], "}"}], "}"}]], "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell["\", "Subsection",
CellChangeTimes->{{3.6*^9, 3.6*^9}, {
3.6*^9, 3.6*^9}, 3.*^9,
3.*^9, 3.*^9}],
Cell[BoxData[
RowBox[{"(*",
RowBox[{"DE", " ", "as", " ", "f"}], "-",
"n", " ", "of", " ", "friction", " ", "a", " ", "and", " ", "forcing",
" ", "fequency", " ", "w"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{"Clear", "[",
RowBox[{"a", ",", "w"}], "]"}], "\[IndentingNewLine]",
RowBox[{"DE", "[",
RowBox[{"a_", ",", "w_"}], "]"}], ":=",
RowBox[{"{",
RowBox[{"x", "''"}], "[", "t", "]"}], "+",
RowBox[{"a", " ",
RowBox[{"x", "'"}], "[", "t", "]"}]}], "+",
RowBox[{"5",
RowBox[{"x", "[", "t", "]"}]}]}], "\[Equal]",
RowBox[{"Cos", "[",
RowBox[{"w", " ", "t"}], "]"}]}], ",",
RowBox[{"x", "[", "0", "]"}], "\[Equal]", "0"}], ",",
RowBox[{"x", "'"}], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}]}],
RowBox[{"(*",
RowBox[{"Solution", " ", "as", " ", "f"}], "-",
"n", " ", "of", " ", "friction", " ", "a", " ", "and", " ", "forcing",
" ", "fequency", " ", "w"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{"sof", "[",
RowBox[{"a_", ",", "w_"}], "]"}], "=",
RowBox[{"DSolve", "[",
RowBox[{"DE", "[",
RowBox[{"a", ",", "w"}], "]"}], ",", "x", ",", "t"}], "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}], ";"}]}]}]], "Input",
CellChangeTimes->{{3.*^9, 3.*^9}, {3.*^9,
3.*^9}, {3.*^9, 3.*^9},
3.0003`*^9, 3.0003`*^9, {3.*^9,
3.9997`*^9}, 3.*^9, 3.*^9,
3.2*^9, {3.1996`*^9, 3.2003`*^9}, {
3.9997`*^9, 3.*^9}, 3.*^9}],
Cell[BoxData[
RowBox[{"(*",
RowBox[{"See", " ",
RowBox[{"solution", ":", " ",
RowBox[{"Re", "[",
FractionBox[
SuperscriptBox["e",
RowBox[{"i", " ", "w", " ", "t"}]],
RowBox[{"p",
RowBox[{"(", "w", ")"}]}]], "]"}]}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{"x", "[", "t", "]"}], "/.",
RowBox[{"sof", "[",
RowBox[{"a", ",", "w"}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.*^9, 3.0003`*^9}, {
3.*^9, 3.*^9}, {3.*^9,
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Module", "[",
RowBox[{"{",
RowBox[{"tf", "=", "30"}], "}"}], ",",
RowBox[{"Manipulate", "[",
RowBox[{"(",
RowBox[{"cso", "=",
RowBox[{"Evaluate", "[",
RowBox[{"{",
RowBox[{"x", "[", "t", "]"}], ",",
RowBox[{"x", "'"}], "[", "t", "]"}]}], "}"}], "/.",
RowBox[{"sof", "[",
RowBox[{"a", ",", "w"}], "]"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"GraphicsRow", "[",
RowBox[{"{",
RowBox[{"Plot", "[",
RowBox[{"cso", ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{"-", "5"}], ",", "5"}], "}"}]}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"ParametricPlot", "[",
RowBox[{"cso", ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{"{",
RowBox[{"-", "1.5"}], ",", "1.5"}], "}"}], ",",
RowBox[{"{",
RowBox[{"-", "1.5"}], ",", "1.5"}], "}"}]}], "}"}]}], ",",
RowBox[{"Frame", "\[Rule]", "True"}], ",",
RowBox[{"PlotPoints", "\[Rule]", "200"}]}], "]"}]}], "}"}], ",",
RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}]}], ")"}], ",",
RowBox[{"{",
RowBox[{"a", ",", "0", ",", "5"}], "}"}], ",",
RowBox[{"{",
RowBox[{"w", ",", "0", ",", "3"}], "}"}]}], "]"}]}], "]"}]], "Input",
CellChangeTimes->{{3.*^9, 3.*^9}, {
3.0003`*^9, 3.*^9}, {3.*^9,
3.9997`*^9}, 3.*^9, {3.*^9,
3.0002`*^9}, 3.*^9, {3.6*^9,
3.6*^9}, {3.6003`*^9, 3.6*^9}}],
Cell[BoxData[
DynamicModuleBox[{$CellContext`a$$ = 0, $CellContext`w$$ = 0,
Typeset`show$$ = True, Typeset`bookmarkList$$ = {},
Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ =
1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
Hold[$CellContext`a$$], 0, 5}, {
Hold[$CellContext`w$$], 0, 3}}, Typeset`size$$ = {576., {120., 125.}},
Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ =
True, $CellContext`a$1450$$ = 0, $CellContext`w$1451$$ = 0},
DynamicBox[Manipulate`ManipulateBoxes[
1, StandardForm,
"Variables" :> {$CellContext`a$$ = 0, $CellContext`w$$ = 0},
"ControllerVariables" :> {
Hold[$CellContext`a$$, $CellContext`a$1450$$, 0],
Hold[$CellContext`w$$, $CellContext`w$1451$$, 0]},
"OtherVariables" :> {
Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$,
Typeset`animator$$, Typeset`animvar$$, Typeset`name$$,
Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
Typeset`skipInitDone$$}, "Body" :> ($CellContext`cso = Evaluate[
ReplaceAll[{
$CellContext`x[$CellContext`t],
Derivative[1][$CellContext`x][$CellContext`t]},
$CellContext`sof[$CellContext`a$$, $CellContext`w$$]]];
GraphicsRow[{
Plot[$CellContext`cso, {$CellContext`t, 0, $CellContext`tf$1449},
AxesLabel -> {"t", "x"}, PlotRange -> {-5, 5}],
ParametricPlot[$CellContext`cso, {$CellContext`t,
0, $CellContext`tf$1449}, AxesLabel -> {"x", "v"},
PlotRange -> {{-1.5, 1.5}, {-1.5, 1.5}}, Frame -> True, PlotPoints ->
200]}, ImageSize -> Large]),
"Specifications" :> {{$CellContext`a$$, 0, 5}, {$CellContext`w$$, 0,
3}}, "Options" :> {}, "DefaultOptions" :> {}],
ImageSizeCache->{624., {178., 183.}},
SingleEvaluation->True],
Deinitialization:>None,
DynamicModuleValues:>{},
SynchronousInitialization->True,
UnsavedVariables:>{Typeset`initDone$$},
UntrackedVariables:>{Typeset`size$$}], "Manipulate",
Deployed->True,
StripOnInput->False],
Manipulate`InterpretManipulate[1]]], "Output",
CellChangeTimes->{
3.*^9, {3.*^9, 3.*^9}, {
3.0003`*^9, 3.*^9}, {3.9997`*^9,
3.*^9}, {3.*^9, 3.*^9}, {
3.9997`*^9, 3.*^9}, {3.9997`*^9,
3.*^9}, {3.*^9, 3.*^9},
3.*^9, 3.*^9, 3.*^9, {
3.*^9, 3.0003`*^9}, 3.6*^9, {
3.6*^9, 3.6*^9}, 3.6*^9,
3.6*^9, 3.2*^9, 3.2*^9,
3.*^9, 3.9997`*^9, 3.0003`*^9}]
}, Closed]]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell["Direct Integration", "Section",
CellChangeTimes->{{3.6004`*^9, 3.6*^9}, {
3.2*^9, 3.2*^9}}],
Cell[CellGroupData[{
Cell[TextData[{
"Variational problem: Fastest slope (cycloid): ",
Cell[BoxData[
FractionBox[
RowBox[{"d", " ", "y"}],
RowBox[{"d", " ", "x"}]], "=",
FractionBox[
RowBox[{"h", "-", "y"}], "y"]]}], TraditionalForm]]],
}], "Subsection",
CellChangeTimes->{{3.6*^9, 3.6*^9}, {
3.6*^9, 3.6*^9}, {3.2*^9,
3.2*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"Direct", " ", "intergration"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"int", "=",
RowBox[{"\[Integral]",
FractionBox["y",
RowBox[{"h", "-", "y"}]]],
RowBox[{"\[DifferentialD]", "y"}]}]}]}], "\[IndentingNewLine]",
RowBox[{"Apart", "@",
RowBox[{"Simplify", "[",
RowBox[{"int", ",",
RowBox[{"0", "<", "y", "{{3.6*^9, 3.6*^9}, {
3.6*^9, 3.5996`*^9}, {3.6*^9,
3.6*^9}, {3.6*^9, 3.6*^9}, {
3.6*^9, 3.6003`*^9}, {3.6*^9,
3.6*^9}, {3.6002`*^9, 3.6*^9}, {
3.5997`*^9, 3.6*^9}, {3.6*^9,
3.6*^9}, {3.*^9, 3.*^9}}],
Cell[BoxData[
FractionBox[
FractionBox["y",
RowBox[{"h", "-", "y"}]]], " ",
RowBox[{"(",
SqrtBox["y"], " ",
RowBox[{"(",
RowBox[{"-", "h"}], "+", "y"}], ")"}]}], "+",
RowBox[{"h", " ",
RowBox[{"h", "-", "y"}]], " ",
RowBox[{"ArcTan", "[",
FractionBox[
SqrtBox["y"],
RowBox[{"h", "-", "y"}]]], "]"}]}]}], ")"}]}],
SqrtBox["y"]]], "Output",
CellChangeTimes->{{3.*^9, 3.*^9}}],
Cell[BoxData[
RowBox[{"-",
RowBox[{"h", "-", "y"}]]}], " ",
SqrtBox["y"]}], "+",
RowBox[{"h", " ",
RowBox[{"ArcTan", "[",
FractionBox["1",
RowBox[{"-", "1"}], "+",
FractionBox["h", "y"]}]]], "]"}]}]}]], "Output",
CellChangeTimes->{{3.*^9, 3.*^9}}]
}, Closed]],
Cell[CellGroupData[{
Cell["Cycloid", "Subsubsection",
CellChangeTimes->{{3.*^9, 3.*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"x", " ", "as", " ", "f"}], "-",
RowBox[{"n", " ", "of", " ", "y", " ", "and", " ", "height", " ", "h"}]}],
"*)"}], "\[IndentingNewLine]",
RowBox[{"fs", "[", "h_", "]"}], "=",
RowBox[{"h", " ",
RowBox[{"ArcTan", "[",
FractionBox["y",
RowBox[{"h", "-", "y"}]]], "]"}]}], "-",
RowBox[{"y",
RowBox[{"(",
RowBox[{"h", "-", "y"}], ")"}]}]]}]}], ";"}], "\n",
RowBox[{"cycd", "=",
RowBox[{"ParametricPlot", "[",
RowBox[{"{",
RowBox[{"fs", "[", "2", "]"}], ",",
RowBox[{"-", "y"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"y", ",", "0", ",", "2"}], "}"}]}], "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.6*^9, 3.6*^9}, {
3.6*^9, 3.6*^9}, {3.6*^9,
3.6*^9}, {3.6004`*^9, 3.6*^9},
3.*^9, {3.*^9, 3.0003`*^9}, {
3.*^9, 3.9997`*^9}, 3.*^9}],
Cell[BoxData[
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxl0Xs0VWkfwPF9HMfZ0hjJyCVKlFQiSpjG8+gtCk0pxIxMpCIxkku5RMqt
cosQKrll5FaRS/idkFxKB4VxCwlJyG0k5+x3z3r99+619trrs/bz+z5r70fB
7s/DJwUIgnhC3/8+/+/qS3p7cFmvnt13yfoQ6T1ApIQY/qE3u2QzILhKPmNX
mGjlxdRc4w02QHCEM/9cLor+994eCKlVkUBJIN0y62ef1jgCUeLkoZ8ku7Te
BQj3ksXj+uuX1rsDEeU643NMBTV16lSVJHsC4dvwebZIdWn+IhBWUh9WhGou
zV8CQvSVyYyE3tL8VSAU1ke4XMRIyilr79rfg4CotVh1elp/qRcMxK8df29j
713qhQFhqMVgPjBe6kUAMVci5XXAfKkXC4Qd6k8+aIF+efIH26yMNnO9dKX5
0aX+LSAYZz6XOFst9eOAkPERW914bKl/G4jiGO2XCSeW+nchoMtfzGany1I/
HQjjbU43eC5IRHdW5ZErbeHQWqcXf6LCC63XOl+mA+feMRmpXFfEmos0UfXM
gIAHo/nPjN1Q1iTJbWnJBLxMZ1Q33wONf1xok7vxF2DRkILzQn7Im9s7+ISf
BwE8dTX+ZT+07seKvT3m+cBJJ8/q8v1Q44GkTKHcfMArDXJq5i+hiI+bC3NU
CyDl5X43QZ8A1Gc34kF1FQDHVoo4FBqIrhyz/Zau/Rg4M36BEX3BqP7wEWLy
ayEQBtc0xxRD0PF3Qr/flykCbLSlJvx0CJo7WlZ0+D+0N+umJk2EIKVjCk5F
sUVAaEXs6xYIQwGnJ95573wKxO1WqmrHdaTle/0hy68YUmSDHJobI1FaRrWF
LFkGXNWMk0PtschvXjPXQB4gwNI9Qa71Huo9YBhdfb0GotZwcqIhE6k33hyO
n6mHgp7lSpsHcxG34Ewwg/EGHGQsxH/d8QSJ27qmvNnVApwC/a/DL4qR/Kf6
Ccuxt6D9Ii9/o3U5Ms1w/JDj3A7JTy4fTW/nIBuvjkGZ9Z1Q55/JueBXjW76
hV87ktcNJn0uKGZXLdrItE2t3fQeuuMic4Lt6tGBN4zBv8b74NBx+ymmzysk
vXOu7p36AJSUTgZpMLioRaWlUDTiAzSc6r23p7UZ7bKu+HtwcBAS5LTZpxJb
UVlW46EUpSF421FVWHr3HVIKeN7k4D8Myd5mnb/HtaPUc513JatHINF/MdyQ
+zf6YbzTo1lhFCzbC85+lepGKtmjQSrenyFVtkH5gmkvsgqQGzn1bAz8E8Pf
jTv2IY/+naaPVo3D+KOCSvfUfpRh7bG833oCTIO+1Cn0D6BqEVui/cAkRPo6
D2zZP4iyt37nVS9MwqFghxatBx9RUm9mumLyV5BlbuVlawwjAbPaJmHjKZhT
FXcTqRtBux2db1XOTYFX2K7Ro4dH0Q2RDayhuGk4f6n+7lOxMWRy/vB8zJ4Z
UHZhaqjXfEGSDlGxhUMzIOjj2Tt7eQKVLJq/TguahTwhjrt25SRqtogMk9aY
g81iaVoST7+i6V86xxmtcxC36YFNPHcK3RhNaHe9/A+Mm7w2KBiaRmpKjo6H
leahSVP8Jx2hWZTN48QL1MzDlGuQcuuhOZScfGO/zrlvkHz9t5WWt/9B8XWO
6QPiCyAoH+0ZI/kNCZeVBYpXLUC29bo59sMFZGYcwxNz+Q5HJtKeT+xZRFmO
fgqBoouQG3Z/vE+Vj84f6Kq7WL4IeoaWLuaSBDasNJIetOFBnmQM95IMga8l
8heP2/LgrYTNb9nyBH7l+biv5wQPHjmmTbKVCXxYTSar3YEHex56+fXsJLB1
yietRjceHEqx75mxIrBbYKjZ42AeOO1mRozdI/Bdg9rIgDwebLT6Oa9Xk4Fn
mvSF5Pg8WFCxIqV0GNjEmynIJfhwi+hWNtdj4Iz1LxhXmHww3773Xec+Brbw
3ccbIfngYfuqSsKGgUtUDs4UruTDGEv4B6MwBva7aj1gosIHjTy96W8DDCyk
61Xpa8aH3Wdkm8pHGNjmo3a52lE+JBarJV8dZ+DiqIXSASs+iMpR+ooLDOww
7Fe07w8+aD7V+fHuCgHcEHs1R+IMH5q63ujc0hfAERPRiTn+fLgTKktOZwjg
VRk5nt3ZfDDVLrpD5Ajg9DhxZJvLB4FhPR2JxwJYPfQCeyifnmd9rj1QKYCN
nPYkTBTywaS8SH15uwD23dZdKgB8SNkSKvpSmIkHKkR4G1v5EGm/ameLOxPn
tZ254vGdD77xme71Pkz8cx3XeIHHh3Ma9wvrA5m4rlRLwp+goKNqtvpjFBP3
JwtkhLIo+Lb3hXBULhNLnEisTfqRgpmdxU/fjDDxxYl64SpFCtRSz940dRDE
e0iVKFETClKKBc96nhPEQ3XJbvW/UvBY0b4211sQh4WJmV81pSB8l0iFebgg
blo2L71gQYEwp2Im/okgtvzhZdpHWwp6Va5p2QmysLP4yeJyLwparl7vHhdl
YbHWjtte3hSA2srSeGkWfhxj4qvhR4FDl2S2vhoLz0ts188KpKBKeMxC3JqF
A1cxG2PCKcjvN5OsLGfhuNX3e8+kUbB8BeeXpHoW1u6ReL4+k4I+hbnjD9tY
uOtOaFpfFgVXQl47u3xl4bVrXB0s8igIEn3awFURwtkKaEq/lILxkTCROG0h
bDLw+O3iMwpE60TuZxkK4YnUDcXFlRTcREGD9aeE8HYlUV/VGgrOKXJ35P4l
hGFDD0v6DQXtr4PvXH8mhO2GD420NlMgtui0aea1EBbMqmmIeEvBbt2c9sQZ
IWykkhMp2ElB82vt6gWSjcc+rXGDbgraouXrVOTZODI7xsz7PQWXFJQZp4zZ
+O1mH+nJQXp/l6ikihNs7Dk2/j17mAKO947eS35sLJ1r13tylIKTMuyuyQQ2
Lndu46z9QkG9s2OIx1M2ttlqlNY1QYGb7fdbiW1sTExUBMVNUdBYZz6tu8DG
afnbHExnKbjI/dDdsI7EBq4ZRsvn6f/P/Jq55QiJR9SlVV8u0Pvt6GhqNCOx
0it7sW7aGk5bDFwsSGx7umB6kraH0X7PCisSd94xLJP5ToEva8gs5DiJG5d5
GrrQ7jG0d+hzIXHeYIvtT4sUZAu/8OC4knjUX37vJtqVK1BmhhuJlWXPbES0
26WPrwnwJHGKKWPCgfa6942vvC+ROLpSzbec9ha5Jq2KSBK7J9yIO8Gjz7NY
vWbZTRIXaHZcvEDbrXlwxDaWxF+aFI+F076QbDG19TaJT7GeKRbTPnla/ZpZ
GomPun0qWManvz9t2QYqg8SxP+yIXUPbaPZFTUkWiblZAV7baV/5En7wWB6J
971fpWdDO24wsD+6jMS6JgaNj2iPXWL5dFaQ2HM4Oq+WdtvDkET0nMRPAnui
u2ivZlzzP1tH4i2l7pYsigLV6JmgXa9I7GDG+VmGdqhqHkuXS+L0CZE1arTP
eiUEd3WQWG5D2kdL2ge9WCGxPST+7fl4nTNt5Xc7zj8YIHGctW5OIO1a6+jp
zSMkbvknKDKe9id5u6MG4yQWjWl2y6HdECVaqTlLYqOtchbPad8IjDJ34pE4
pN5Bp43697xOvF5LCuMa+8LVn2kLHvzJVnO1ME67mzrDp/1fN466fg==
Axes->True,
AxesOrigin->{0, 0},
PlotRange->{{0., 3.755}, {-1.6277`, 0.}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output",
CellChangeTimes->{{3.*^9, 3.*^9}, 3.*^9}]
}, Closed]],
Cell[BoxData[
RowBox[{"(*", " ",
RowBox[{"Rollling", " ", "wheel"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"RW", "[", "t_", "]"}], ":=",
RowBox[{"Circle", "[",
RowBox[{"{",
RowBox[{"t", ",",
RowBox[{"-", "1"}]}], "}"}], ",", "1"}], "]"}]}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Marker", " ", "position"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"MP", "[", "t_", "]"}], ":=",
RowBox[{"{",
RowBox[{"PointSize", "[", ".015", "]"}], ",",
RowBox[{"Point", "[",
RowBox[{"{",
RowBox[{"t", "-",
RowBox[{"Sin", "[", "t", "]"}]}], ",",
RowBox[{"-", "1"}], "+",
RowBox[{"Cos", "[", "t", "]"}]}]}], "}"}], "]"}]}], "}"}]}],
";"}]}]}]], "Input",
CellChangeTimes->{{3.*^9, 3.*^9}, {3.*^9,
3.*^9}, {3.9997`*^9, 3.0003`*^9},
3.*^9, {3.*^9, 3.*^9}, {
3.0003`*^9, 3.*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Block", "[",
RowBox[{"{", "}"}], ",",
RowBox[{"Manipulate", "[",
RowBox[{"Show", "[",
RowBox[{"Graphics", "[",
RowBox[{"{",
RowBox[{"RW", "[", "t", "]"}], ",",
RowBox[{"MP", "[", "t", "]"}]}], "}"}], "]"}], ",", "cycd", ",",
RowBox[{"Axes", "\[Rule]", "True"}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{"{",
RowBox[{"-", "1"}], ",", "4.2"}], "}"}], ",",
RowBox[{"{",
RowBox[{"-", "2"}], ",", "0"}], "}"}]}], "}"}]}], ",",
RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "\[Pi]"}], "}"}]}], "]"}]}], "]"}]], "Input",\
CellChangeTimes->{{3.*^9, 3.*^9},
3.*^9, {3.*^9, 3.*^9}, {
3.*^9, 3.*^9}, {3.*^9,
Cell[BoxData[
DynamicModuleBox[{$CellContext`t$$ = 0, Typeset`show$$ = True,
Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu",
Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ =
"\"untitled\"", Typeset`specs$$ = {{
Hold[$CellContext`t$$], 0, Pi}}, Typeset`size$$ = {576., {115., 119.}},
Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ =
True, $CellContext`t$476$$ = 0},
DynamicBox[Manipulate`ManipulateBoxes[
1, StandardForm, "Variables" :> {$CellContext`t$$ = 0},
"ControllerVariables" :> {
Hold[$CellContext`t$$, $CellContext`t$476$$, 0]},
"OtherVariables" :> {
Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$,
Typeset`animator$$, Typeset`animvar$$, Typeset`name$$,
Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
Typeset`skipInitDone$$}, "Body" :> Show[
Graphics[{
$CellContext`RW[$CellContext`t$$],
$CellContext`MP[$CellContext`t$$]}], $CellContext`cycd, Axes ->
True, PlotRange -> {{-1, 4.2}, {-2, 0}}, ImageSize -> Large],
"Specifications" :> {{$CellContext`t$$, 0, Pi}}, "Options" :> {},
"DefaultOptions" :> {}],
ImageSizeCache->{624., {159., 164.}},
SingleEvaluation->True],
Deinitialization:>None,
DynamicModuleValues:>{},
SynchronousInitialization->True,
UnsavedVariables:>{Typeset`initDone$$},
UntrackedVariables:>{Typeset`size$$}], "Manipulate",
Deployed->True,
StripOnInput->False],
Manipulate`InterpretManipulate[1]]], "Output",
CellChangeTimes->{{3.*^9, 3.*^9}, {3.*^9,
3.*^9}, {3.*^9, 3.*^9}}]
}, Closed]]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell["Pendulum", "Subsection"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"(*",
RowBox[{"Analytic", " ", "solver"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"DSolve", "[",
SuperscriptBox["y", "\[DoublePrime]",
MultilineFunction->None], "[", "t", "]"}], "+",
RowBox[{"Sin", "[",
RowBox[{"y", "[", "t", "]"}], "]"}]}], "\[Equal]", "0"}], ",",
RowBox[{"y", "[", "t", "]"}], ",", "t"}], "]"}]}]], "Input"],
Cell[BoxData[
RowBox[{"Solve", "::", "\"}],
":", " "}], "\\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", ButtonNote -> \
\\\"Solve::ifun\\\"]\\)\"\>"}]], "Message",
CellChangeTimes->{3.5997`*^9}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"y", "[", "t", "]"}], "\[Rule]",
RowBox[{"-", "2"}], " ",
RowBox[{"JacobiAmplitude", "[",
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{"2", "+",
RowBox[{"C", "[", "1", "]"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{"t", "+",
RowBox[{"C", "[", "2", "]"}]}], ")"}], "2"]}]]}], ",",
FractionBox["4",
RowBox[{"2", "+",
Row