(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 29927, 577] NotebookOptionsPosition[ 29182, 547] NotebookOutlinePosition[ 29577, 564] CellTagsIndexPosition[ 29534, 561] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ "Lyapunov exponent approximation \n", StyleBox["(Homework 6 computer exercise)", FontWeight->"Plain"] }], "Section", CellChangeTimes->{{3.419682249757717*^9, 3.419682256008585*^9}, { 3.420287751523087*^9, 3.4202877688960752`*^9}}], Cell["\<\ Here's something for computing the Lyapunov exponent for the logistic family \ :\ \>", "Text", CellChangeTimes->{{3.4196823322059183`*^9, 3.419682341194043*^9}, { 3.420287696473742*^9, 3.420287716646884*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"simplelog", "[", RowBox[{"a_", ",", "x_", ",", "n_"}], "]"}], ":=", RowBox[{ FractionBox["1", RowBox[{"n", "+", "1"}]], RowBox[{"Plus", "@@", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Log", "[", RowBox[{"Abs", "[", RowBox[{"a", RowBox[{"(", RowBox[{"1", "-", RowBox[{"2", "#"}]}], ")"}]}], "]"}], "]"}], "&"}], "/@", RowBox[{"NestList", "[", RowBox[{ RowBox[{ RowBox[{"a", "#", RowBox[{"(", RowBox[{"1", "-", "#"}], ")"}]}], "&"}], ",", "x", ",", "n"}], "]"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.419099206149913*^9, 3.419099241399033*^9}, { 3.4191009351170673`*^9, 3.419100956473486*^9}, {3.419101015208037*^9, 3.419101015558031*^9}, 3.419333753596129*^9}], Cell["This allows a more general choice of f:", "Text", CellChangeTimes->{{3.419682267473666*^9, 3.419682330361566*^9}, { 3.4202877250255833`*^9, 3.420287732621842*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"simple", "[", RowBox[{"f_", ",", "x_", ",", "n_"}], "]"}], ":=", RowBox[{ FractionBox["1", RowBox[{"n", "+", "1"}]], RowBox[{"Plus", "@@", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Log", "[", RowBox[{"Abs", "[", RowBox[{ RowBox[{"f", "'"}], "[", "#", "]"}], "]"}], "]"}], "&"}], "/@", RowBox[{"NestList", "[", RowBox[{"f", ",", "x", ",", "n"}], "]"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.419099206149913*^9, 3.419099241399033*^9}, 3.419333746215425*^9}], Cell["The graphs should all look more or less like this :", "Text", CellChangeTimes->{{3.419682436711038*^9, 3.419682471594778*^9}, { 3.420287741081884*^9, 3.420287748457333*^9}, {3.420287790564229*^9, 3.420287791391385*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Random", "[", "]"}]], "Input", CellChangeTimes->{{3.4196823597964067`*^9, 3.419682360633583*^9}}], Cell[BoxData["0.8475510875738919`"], "Output", CellChangeTimes->{{3.4196823623037767`*^9, 3.419682366929373*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"simplelog", "[", RowBox[{"a", ",", "%", ",", "1000"}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "2", ",", "4"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.4196823471563597`*^9, 3.4196823775793877`*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJw1l2c8l+37x6+vdV2XNCgkkawoCS2UOu5uGsRtpEEh0dCgKGVFQlaI9hB9 k5ER2eq4ygolI2RlryizbH7f/4P/g/N1Pvgc53ke7/freHIa21uvtiUIgouz uv5vB8JwVo5/3oyhC8J+KZ65DYTRx1JpnZUM/Wbdxuub64AwWWX/8r0oQ4d4 Lxg6aAqE6fec1ykCDK1fnWCm8xcIcz3DnPBxhqq44axulQXESRW3WvPPDPlz 4uXA8CcgTgVLH/lWxJAVfBN7yoyAON1X3Kz/gSHTHaSkY+yBOPtCqHNnBkN6 7auO0lgMxEXh2BGZKIZc2TKYI+MJhOtUzcL+KwyfNSszankGEO6HXMJsHBg+ vdreucFDQHi8lVrWfJbh29wflM/9HQiv82fFv1oxfPzlnZa31YHwbWUppukx vJm2r2Nk84EIK1yv4yrF8K6U8BaU4/T/KsTHnb+M4Z7ZtCHR9gMQsTVLjGNa Ge6BmERlygaI+BWPFXaNMdzN6t3mk4+AeP0ypdpVkuF+n934Vj0XiDe5TYr9 Dgy3t8TR2BW3gMjtVav9LMJwL5k1Z288DUQFtKoEWzNcO1x1w3dPAFHpe5ZP 0ZnhUuPbkRToAkRV2Z/GwiCGS46dICt3EYiaQwt8Z9MZroUrxuuUPgPRcH5L 8wWKYTV1D3T09gPR+TA4wCCJYV3/rFg8KAzExIhm5+JJhqh0eSOYLQnElKr1 SK4BQxRrf5SNKAJi2iGQdTKaId6Jty79lQnE7O9GyXe6DBE79O7dkmJgcfW7 HrF7whBeqi7bQ3YAS6Azt7wQGGLjmF6CZQWwpOq2Z7v44vy/Ht/V5+8BS1rE 5pNcA86rdnPVeHcCS8Y0qLZyA85LWYVsdKsBllx109ia7zg3GyLz/J+HwFpb 4a5co4hzWV8yDjtmAWtTyTv2hs84p5ZE6fkvAdbevB23u5birNr8mxOqpsDS FXFkO4/grPQ281dt1cDSc4jN4a/CWcFbz+UWDQDLQE6oSyUUZ4YO2apPxQLL JLhnm9sinEn68DeyYA+wLCxDepdSOKOcaLDD9hCwnLja/oE5nNZk3f+iyAAr cp/36JNenLxaWfosKRRYz4OemRu34+Q5ex8JA46PqK9Z+VQjTlqJtEZk+wOL bfrrjlM5Tu65sVOLfQZYcdaHVPen4aRIKc+jGRFgpbmuvTDjgROZgjbGAQ7A Kn5d8fOoMI7PyVgffzUGrKElEm3iu/Dvi7JVwqEGwLXn4eQtt0M4tvCNbNv6 3cD1NFDGMfEIDre0Ct/5LQ5cv9KjetitOLTc5mnjgmfArf3CWvi6NA483zxr /EgBuO+WXvm8qg/7Xi1/ZBTeB9xdh+uLQp9g9+zixC1iT4FnW7yHSyc/dh5I 8QtzPwU8QZudM2dnsZ22Hm651wo8Hd6pSVMJ2HLpomjpSCnwamhs6yvejY2b GxlZTRp4/RXEEs1sse66u79WYyjwtjRorK0ArE7L2KXNkw18KoWJqgKRWGGm ULR953ngu+VkPbwuB0uqHs7pEQTw/Qg4+WfFB8w/czkg8DQXkMr16FWhj7nn r/I+3PQWSJ+BI59OR2CyflXt/gQXIP1U/b1jFmPShOKn0pjjQN66kr2jIwCT ou+YKvcbARnIWp5u4YWJ4w/cZYsdgAwVros+cAFfR7tyPajeCeSjnaZusAfj phO+y2X+BTLpjpGK6BSy8x5vW6sxAWRynWe/6WVknzHPVdneCuSblSkxd4aQ LaJ2j8nxBDLt1RKJhT34wkGW3ZFeA2RW7ld+7m8YLScrmN3Eqf/Ytb9jMBEj 1/SOSST1Almnvvd+sRXeXXZJQ/usIJDfnYWDloljxDWte0pbs4GsT2/3tKrF 8B9XiMzVnUA2bXQ7PbUf78Tvka97+hXINuUUzfVbMXSvzk2dpelAtp93V742 giHW0eUHfMaB7HitK130Gm9vUZQOO8g5363YucBKGgPbialv3u1A9suJNoUv Qj+jYL1l/c5ADpzorGgpQd+N78sfvuQG8lf0m0Klm+gjJmaX/vsikENSeomF U3hjUEZburoDyGHL5dFC6ejVRt4IEFoE5MjTrnuWDuj5PXytm6YZkH/Er3tM dKF7znuzgKY9QP412++o8xxdm0TWn+D2BXL8odipO2boIqrjPzdwDsiJum7z FmG86rRJbdYkDsgpkTTDdRV4Zc6vWFlmK5DTpp7aVwPxcnF4XDJ/A5AzEfrq hbvRqTfMUeR5E5Cz1SvWC7HQyftGlo7TMJDzgj2rLfLQCVcYhZoVAkUYvhVO cMbLOfc7n8hvB4oV4kVPqKGzRwjZZTINFFe5/pz2L7w2IXUu7VoLUDwLV4yE xaI7d75emPsZoHj1erp/nEDPHep+l3eoA8UX8LZxnSTekF61K0EtECiaMsgv CEdfjSMNhvdrgOLfI54laIC3rILLeFevB2qBT+9rCwoD7irZbY76B6hF3Dci xj0wRKJNV/5bAlBCcxkn19phhA6fCteBF0At0/I2c5bDu9fKBU7w2AMl7GZo kN+K997wxEQ7TgIlOvlz67GD+FBe4pmW+UKgVv6RpEL/wWdr6t3Ors8DSr7f J/7vcmRrmVv2Bl4FSkHHMmatBbLHNtvyLioGSjFSI8qCjS8T2ouSc28CpWT8 637xBnwlkWgg3B4HlGqO6c0HuzGev/TSkK0IUNv95Sw0L2MKr5AxK2UbUMby hYIeVZihP2sPfQJAmXhGCqQtx4z0hrqb9nNAHWi4RvZYYKakxN8+e0ugDt1W mvuvDzMH62SOSk4BdfRPRP9qArMfSIes8XsC1KkCm4J8Zcybifj01JjD627N c5UMxI8Cry4sZHF8erDzCy1+4sd/w8yfqXB8Xe/xFsrQxY8uKc/qmzj33TjP SrTlx499oQGaDZFA+TnPtX8MwPySltG9ZiZAhQVP6rv7Y+FdQYGXJ48B9TL7 t/SoH5Ycfv5uUZASUDEziQ663Vhyb7KsaPIPULE7z72L3o0lNXfLQz2OAxVf 0H/YiA9LTR4Ui0aKAZX8tTc4yRfLjOu9i1Jlgcruap8444NfLArEFHSzgCoX rPvS4o0V33V7beqSgfpqz5aWvIAV8zt3FWrnAFXxxeHKscNYKX+KTJXg+K72 pyUblbHysnF985sooOoJTfvaBqwSjtfXmncDquP3U8FyFay2+vb7+JQoUJNl 1qZ5zVgr517lFLcGqGlF5bjpT1h7gHT7UaAH1Izv1JxmKtbeqA+4FBoM1Pyu 8JisW1jbZrZj6eYFQPNkF4y/3YR1Ud3DyWsMgV4cu+ZhQgDWq3ErrPPnA1ru 5kDzo63YdHe+aW+rFtDyTQfteq2xqTwweafnfqAVNuLEliBsJr1ULNkpQK9t DxeqasNml2RTL3VPoFVAY++CQPxxUjvvgsUKoDWnfVLdWrBVsX9onc0w0NuN B3eW8WPrvs4zR/ZoAK0Vf/iL2CZstZNaHWcaBTQcWdeT7o+tCbZpWkWSQOtk VK38vRHbVPc9/3dyM9CGlyT9LP2wXS/0Q2WeCdBGJbeWJaZiu4NhTcKhi0Cb SA1HTTdh+90ztksFi4E2/VqQ90AV21vjn6wJSgfafL3dcEUjdrg+i7vFkwC0 bW/m0V0bsLOQ13zzmUigT4FUX6gZdg6pDH4SqQT69P2AKz9uYpd4rZJ++Hug z+4+FuZSj11OD6wIJUGgHdjcRWne2K20X6NONBloV0tDVdk67Mk0DPz7OBZo t8jIJfdWYk8veWmkYBxo95ZfQ6QV9oodTlVIaQDa0yowqe8n9rrXv5rf4AO0 z/FipaQZ7DMoGRk5zunfN1pEQAqw72axxkG7T0D7tdv0h/lgX65X56vBKaAD TnDFOy3Cn+u0L9+OTAM6xGaHgoYU9i97lnOPrQR0KDuYjLfBfiP5xfKndIEO 62zsEY/D/tuP9ZYEcXxH2F57OaeGAwLnt8b2OgH98GS6TL4O/loWXiTkKQf0 izNKq/TtcNAyXEcmkzMv7DiXufdJOJhs4dunweF/2VfyY8MoDhEn3j7V+gl0 rN2pp0JuOHTA3FdQi5Mnnn2xou42DpMPxn6XuACdlDA8tbcKh7fybFksdRXo 5AGozxHF4VMRYmzVZqBTz/148OQ5Dn8+asnrwOHLvLBC5HgajsScu2P5+DLQ WUmn/1RN4EhrDV8JIwF09u/MGu3tOLpymveVPocvz940Qr4IRx8NqjWQI0C/ S2Y7PuDH0abp4XLdK0C/Hxw1of/DsdXaogLhJNAfHMIE+7/j2Ju1ihJ5j4H+ mNIyfGwljs3KGbjGcN7PH1auKD+Of/aLZrhWcuaj6OLnkJSf+Ge0lSx1ewp0 caq4vbQy/tXfx45x+g/oTyN2BuGO+Df+27RpDMdXmSO18PIMjjtOmlj2fAP6 c/mJ0fLPOP6Ta/exeB6gyxWxYc1jnDh1+LG17Q+gv/qIM552ODEwco7PmcNT 0eocU6+Jk66TLTFunLxqW3WwGj9OCTdNNCQUAV19f4NjQD1O5Z1oU+Xh5N9G Ao90xOL0xUqVshojoGv1e2D7VZzREBJNcrkOdF3cv/J39+KsCFVZGlkCdD1P pMBvUZxbxPOhvWQ50A2WU6O7u3Fe4XbQuRpHoBtzD9ZHZjCE+vl1OaPVQDeL pOKEL0PMGBl7TGwA+selRTFGpgxrcHtt5uJWoFu+2AXFyzLcG2ROc4UFA92m UOTINcbw8mRa99cfBLr9pvQR8wKGuoW2fhViQHdwjYiMx/////F/us/3xw== "]], LineBox[CompressedData[" 1:eJwVV3k4VV8Xlspw7zn3Ukjxk0KUMUpS2QspMlQaUKGQ0CzRSMlQppCQqIyh MpQo4SDzUGbXzDWECplChm9/f93nfe7ea6/hXe9ax/iy1YazHBwcnBwcy5v+ /wu8fZbBP+KY+bxFQb832wcAb7966ybSKp9LWMWam4oB3oFQVcJWOX95U9G3 2IMIeH/8CRqnOPM57fVD+1mmwDtk8LtlTW3+MjP27VklXeAdTtSjrrzK5/C4 pWdTXgC8v5bHx5dfoZYSfQdiD98D3t+Wy/w2ArV45+h/icl2wDuSbe54m0kt OFcWbFzBAt4xwc+m9Z3U/OO1gWWbsb0/VwWRXAr1j1pvcYwhCLzjVVelPF2p f5xOBZ1xGE9KV9M7Dam5Y7WpXunCwDvd6cF6/IualZDS3P27GHhndnbnDeZQ M8+ejh/6eAp4Z5/ujtf0pWaEC03mKpWAd15/8urEFuqv/LEDPjH4/sLrQ6b6 c9R0wabbqncXgXeJ861GXDk1feLGyqSgIKAt+2xDP25HTT1/FiG34zzQlgtQ f1LVqCmtkJEbPnxAW3FFhMXDTU3+PjO+WWwz0Lg31cd9TqAmDU+5M7WVgUaE aGuI9lETKqbSs8o9QCNHX0pez6DGZ1gvpl/VAo154B/t2wNqPPVL1XXT00Bb tex9872N1LhSQ7zuvTygrbm88UrvaeqPASE0byQFNOEKV5PdStQfobrUK/6h QFsn1bonZIkaY899sIk+BDTR9mDavhfUmNuG13JrI4G28cCy2KR2atRH3EvK +xLQJAK02qVKqNEdU9vnupYBTbLOQyA6lRr5YbeRmv4ENOmTXJ7P3KmRA4OK Rtp9QJO7QD/3SIb6vTG/QiVjP9BUAwRl7a5TP+nxeUy5k0DbUXfcus+CGq4+ u/bbsZdA2ykUFnlalxoO5H8k6YD92/VyHWkqSg0LC/6xK/ABGqSJje4vpIaU RWyuvNoBtAN1Mh+kGdSPR2YOxl+LgGYgZP8zZob6YaIi0ocKgGZ44o2kWA/1 QyqhXnSLDdAO9so/FcqgBooMW27+LgXa0UllF66T1ADPrgm/zmNAsxTard6f SPUlcsQOduB6XTtxsDBOk2KrthwxXt4KNKdwhW2KDhR79e1Cl5+4ftebiPjs EKpn7MJZI3kvoN04XP6wZoDqeaPyIFRWEGh39bSNFvyoHkn1DWs/6wDNW021 7XgT1S1xzElL6iPQngv9N8nrQLW/PvMj81EV0CKPzJ8NeUK1X5MU/OArDLSo oNbm9TlUO1yUy7DpBtqL7591kweodr69ouq9T4D2igzLVmVSbQOG4/xvHIEW re8sW6BGtZVLpDCl3IEW8+holIEV1XYtNVW77TDQ4lfy37PKoFr/jqs93iUL tJTJlcSXRapFK7zlzjcSaIW1v6QuBVNNsuqk2O0bQCva0hzVs49qnDlzZqRS HmjFD74KHZmjGotbiZId4UAr2x7BvdOaajyTLXiUdgto1c/2D63YRjVErZrK Et4KNNbp2LfPGqn6Dbusp/KYQPs9elKlWJiqMbuUJHgqB2ijuvve7qiiahQ6 am4p6AFtLFpZMvkeVbPigXDv+xagTRzhFXw8RH0/3RZg/Q6fn8nM/Hsim/om aHfr90Y7oC93XZU9dpKqDJIEyd9rgb6WqNASeUGVVH1K2yOzFejrjJ4T6Reo kodSdn6/fgFdJPB88351qkS72cKydBDoYgLkBcdmqjiHh9pBiANdQuRQaAUf VZQa1JuzWRjo8pubft5wpwpjSstFmjuArrmvN6TRlsq9UNq+52M40LUeZlie 30bljC9I5a6sB7p2heeWZZxUzs0XKYJxGkDfZ7SJkntJfXkocbbfDdvTP2Y/ 9IBFfU4WM1NC94F+zHpMQ/kAlblMLdr6Fz/Q7e/N/3isSL0byf20MagI6A4d rZ+Jh9S7ZdInvvONAv2C+iefR2zqrYDFyP016UC/NOkk7xFCJe9RnrzjzA10 p3M/r92coRLCjcT0nzwBuptRy5JNIfXiB9dX03e8QL/3JquGLUJFSaYW5e6P BLo7z9MYSyfquXU3kXLoGtA9Cg/rnJCmwgZ65FcINQLdZ3uZ70E/KohDZE1T pTvQfYMTTlX1UwFPnbLh5Veg+416KBzQoHwVwr4v27QD6I+TUO3eMcrj0qPd wzuNgR7E9V9MoR51f80Oy622iUAPtvrnhGKo24EeR91TPwP9CdWikzNPXSek HvheYwP9qWjWGvVj1MULVfnLxy4CPfRmyGBWCnXmpeanW1c2AD2syTF7Ow91 +OrXEc9dAUB/pnLI9/0Z6sjaVwoP3W8BPSJQwVwxm7LTQaOczTi/z38TCu9W U86jXldOeQcDPUpvmGPLBeq+975SRasMoL9aHh8rJUYFRF/fdsLVB+jRpx84 xbhQwYVMq1TJfKDH5J7Zt76GCj1gf1LV2hXo8S6iQ+seUFHPLBTkt7GAnjwc osgcpBICiZ6yAyZAf7vfkcNfk0p8xvNi1fxLoL+LO1hHi6CS34xuFdJgAj3N gn59pSGVMnZmXfxZXK+Pte7Zc++pDwuLFsu6BICepXDa7wad+qiwhvvd1Tqg f/LVMJ+yoTJtrmsJmccB/cveOY4xIepT3/T3kdsKQM//dHVf/20qbze13tYx BOgFi9scPEcoKuDMakslfP6r9t+ATVZU/rWyK4n3cDzF3+82ntOlClizm69m M4BeMeBjPSREFbHKoXX5HqBXyRl6P3pEFcfLRWjxGgK92pHvzeZ5qsSFrJrY XQz07wuh4w69VJmc6+REx36gNwrE3vuVRlUOip9X8cZ8bzphG+cvQVUNm3xp yMH5bn61uUw+lKqeMmyI6YsFeqtsGvPSXapm9cyFyXfZQG+7ek2FMUbVJD64 L9GmBfT2LNXjqTZUrVbLUFLlBaB3aeVEjR2g6nzqhYJGrwK9+6FbQWAuVb8z +NeibCvQe75p9W9VourHLHgCVmP+9JmVyzmuoRpv7rPX/4nr2f/S79AqH6rJ 0GCbt0gn0Af6D157v0g1y++OaDu/DOhDVxqzJ/qoFqHCuWfMaqAPZ4Z3hphQ rZI8KzM0sP8/509xbqug2vbm95g7bQL6b631Uo0aVPvtsbzX0jhfI95svevp VEe1OfoBD4A+Wh1/UUiS6trTOPiygRPof1bbBWaGU93239uibnoCfdxMNsOE TvXoB6S9U0sA+sSLkea/rhRb788ermac38m+9H9hf6jei/t910/bA316y/X1 ajZUX6F4Ey0mEOh/r6hpsVjUwKlKX86vFUCf+fjP9qY+NXjG9kFgkSzQZ//l +azNo36u06sIvoDj/ad5/122MjVqaLmftywJ6PPee2tPxFOTN54VOHVJAX2h mnvynzA1ITo63emL+3dpVaXwc19qxJXqU83hB4LDNGD3rkVq+IXtMVj0B2LZ i8OW7Y7UIK/+WLa0KRCcfQIP7vRTA2ujVyT+//8Vm5tf/2dK9fkL2L2MKAJi 5eWIitwqqldrUUS2iACC66P5iAWi2JtF+j4mdgPB/W8D/+J7qmf7eS4ZGRMg eKFv+0spqtuWtp9nxycgaF6vzVA41SXcFcf48wwIepXDnW6C6oiRjOLQFgaC XCX/6p4b1W74Trg47wMQDJOxIvFxqk0g+NaPw2eAYEZ9+FFgS7WMX/m8q2IN EKs2qytyGlDNA98khgp0gVhtF+GbvkQ1LaDeSZN4IARezw2eeU81iSfl/XO5 CsSaTV+iC4Sphkfm5dvHnYEQthVZvFJJ1VclRJ9pVQZibfydk+KuVL0wS/b6 0S9AiEruEXDro2rL+wzVhLB//1lHXVUMo2oVHNpODOoBIRazUN11gKqJ5JQf ElwAYsOGPG+NNOo7ITu6FlUCsWk9ml24Q1VEb6m12jIFhLTFq+PvFKnyeLnj T49IAiHzguPDKTZVln60c8XUJBCyogUXcnSpkq5L6mON7kAordPqui1AffX8 Wbev3wyIrWaxu7eUUoW5HBkFHTh+5WfLw1tvUQUzw1+frNkKxPY1RcY7u6l8 12Np4qE3gVAX1CmaeUN9cdJK9bLdAsSuowkbEi2obMnY0fd0YyB2h3C5mvJR n9paqh+N/wACrSrdkeVMZZpXzWj/eA/EXj7dZGctKj2Zf/O369+B0DmYxC01 RaX5bl/7qfMEEPse81o3JlIpvbYxy6JbgdAjK/7bzqCSI9+ql2QFA2FE1388 0UrFPHJI99ITA+LggTe/Yv2pVy802wbZX4E45EPoHQEq6naJmULHcSCO8FRz vk+gQuFH6rGoO0Ac1ZW3tDKjgqsOyWb7jQFxzNs/ZxVBBUxaSBa7RwJxvGRE uICiHr71te9QHwDClMvo+lVHyv2vV6bMexkgzHRS6jZIUXf/6187oIXzc8KT qVDDotzMd7Fut2G+nCy67HvPl/I9vXHxozUXEObLvw8qaVCRm8a1zb5gflto K+3t+kNFRvoo1F7RBMLSPTD6cRwVc2Bw/JKHAhCnC8YWkQmV0NDvH3ywEQir ZYdPjNKoN2xN2cJ7mN/WJ/PyvHSptKXMwN0R2J5NptzG/7yojKuAfvlgftry PfPMKKI+nXo0aCaP4z13nmvIAO8bvUVzp39j/+yKnQx6gSr4b6jr2wk/IBzW 96TdcqWKzNklY5K4HudvHRTgz6HK7pW4H43hBuJCQ45L4iz1Tb6hvyYf1+OS 4uY2UKXqo6JP8kngfF1+FKbR7ER1vSuZjTurBsSVvhUxl95TrPURfE+y3IBw 1HDkWjlG1SrLvM3NsAXiWniXfZQCVSl+Ad7cwPx0mjCo3ob3r02xmVFTM0A4 G2YrVSZThSeaDv+SxvxySZQOsRqk8h672lPPngJxkzPk79wm6rMVK1lE3hCI W+acJ4NsqIxhmb0s3k1A3M66TMnEUGmzN4vCxL2BuMvfsZHqppIHIhgC5oeA cL1wwMtEjIrz7Mi5ciMUCLeSrKGRU1TU2TQudTkc//0NUoaeEVRwmsrQ8XOJ QLjfDk4XbabuJ1varj+/G4gHTRwCGQKUUdLAK0txTyA8lS7e0DemAkR7ExQ5 lgPh5dPaxg6k4ieNeY1NcP29+3XRzW/UF/kXV/R5sP1H6GMsP0F9nlItmz6E /fOJkOBKPEDFHX26pFV0HgjfyUAH9JC6v7d70jnuMRD+RgvVTSVUAI+SWcVH VSACks5vvbSCcjyp5O0WjvUzcDnr6UotFCL/H0N3x1sggix0ZiLvoVhvXn3x BUUggj99OKWSh5L2lxtk7T0CRMhqcapiHr1tT5w0Xo7tPb0YIGGljt7Ur88e 0rYDIrT0n9fsDZS2r/FNkb88EOEb7YYDP6I0nZlPvxM6gXh2p8lIehxlvDbc JcXJC0REs3Y6pYSyth7o7dbD9iK3pgsev4SyvHXPvKnIASLKT+zG77fo02kU bX8f+/tiwLfdYxil9/d+CxNYBOIVzIKoDMqyDp7WM6kHIvq5bewHW/QlPqUo 0TMCiJipBu4DcShHLtB1gXMJiLhDmg49bPSlba9Km7M1EPHJKd9uiqPc7avu W82nAvF6pagynwXK3dekWbXtERCJlo+evo6kvi8UqrTRcX6SPk/ParRSdxzv 538LxHr5RsDmVJMwei2A+CyGcH7eXqrNv3gcff7nEhao+B8Q78qR5IoQlBf+ JNjAXgiIVIm33s/rEMUIKzESwv2U5rr2pwofKqhSMKud+wdEOsvbqMIIFQhr /H1HpAPxQXny/Rk/VBAwzbV2xTsgMvytBGcqUMHeWZ6q65j/H398vxnIgwpC kyKr1f4CkaW1u0N6HxU67fyx4nMPEJ8ikyHvAco/53Ng6VQzEJ+75qdvtKKv swnmvqNYX79IHHyrooyKLvIMSMzsAiLHNubMbx9UZHz0awL7NRC5SVNrXrNR ycqunYcnmoDI+72/Gte7mLy47v70USDylSIeiAajkvMfqlfsmgWi4NpvtaZh VNTqV+waJgpEYRYaCdJCxec2u50r6Qfi61xwnH4EKhEVvqfAwvOseE//iZXj qESt7rnEXhxfyX01Zv4BVPLsiVLiIOZTaZFP8a0Y9LVmVmavcBYQ5dwdt7fN 4fMtfz/k+AJRoa+kNGqMip9d9yMsNgJRGeA+kJSMit0UuqW6XwBRVdsYacOJ SgdkDL4hCyC+CcoYi51AmXalr7RbxIH4bnqbm/UelfdXVTvX7AGi5nl17hMa Krvgm1vgg+dFbZf4NUMrVB6UZPHzMZ4v9RuvbebORqV/rdJWjUgB0XC2uLNw FSrLXb2NVtMCRGOScMgdB1Sm9uj5rjC8rzT9On9AtRCVKs8+iboeBwRLMXfp zzoqllcy4Pv3bCBarvF9fONIRQcszah3FADRmml9/mwFyvun4vc65gIQbbOZ G8Q3opJGC6sDVRh37OFpar2FSm7ol5WeugZE572TfiF1qDThooJAiyUQXV9T tA5uQeWyqZOv2vF+1cPNMcPrjsrb18dbhn0Dgn3gyLuvrVR01NQbzzrsX69/ grWrMiqpFtJuC6cB0Vczu1bNB5Xdi97kt7kDiAEB/W/jbFShf6H/synW3x8m LzzeqaOKZfzX/X7h+TAY8Uf9XDCVkVioPzCF+TbUqT22YRiVpnMnZ6fhfP/c EBrfroUq9uz+bF53DohfNkMnQyNQ+duuTcLiskD8TtzNf2gclRtZJa5PwPZH fgaU0g+gyq9e+lGPsX6OKfTcKY5BVff3VkmE4H7547hN2W0OVUVvzFA+mQvE +EevQXVjVCn3qlRr+QMgJmZaoiaTUUWv49Xi09NATO2WO5KK63+387/W1s9A TLu58dqfQBUmujIqibgefwtrKYn3qGrH5oZ9/RuAmOWSdOqkocq4IyuTNBqA mNNz2RJuhSrGjehmndj+P7/ybuNsVDm7/qJvTBAQ8zWioeQqVJF3LuEdDevB 4urL+qUOqEJdxmuh/C4QS8cLl90vRBWa29LNM88ByREhkLV7HSovXRdq+V4S yGUd5y5MO+L6XTYuZMoBuVw8e2NaBUq9qEbzW/4JyBU2BOv8Rqpq8roN/5E5 IFe+tgyQukXVRX4Q99B2ApJr+L12Vx1VqFO84lfPGyB5FFbMRmyhMvNWlAv0 pADJe9Uk9eh9KnRu39QaQ3yf9j3r/LQNMpnTXLymPwokobBGJnwflXKwUvel uRCQpJ9zv/pmKu2ggWfmQjqQjOGm6HaCClXRY4gLDALJp6dqcXcUJa5WXjyl 8RNI/tehIuvr0McVPDq6h4yAXL1yqrngI8ok5rjY/QNACtgcC7EOR8U0f4lJ YxUgBQs/Hl55G5VZhHQcpmN/12wQJF9boEqO6EZP/itACrs5VehpokoL27GG HceBXNvR4PVTElV0zuaPrjEFUmT3Nm1/blQu1sp/PMMVSNGIEA6FYVQ1tm/P /t77QP43M5FTU42qekZ2KCbuAnK9yZGbjmmo+vnhZN/sMiDFP35QFXiCqsvm a4vssP2NAqvGM53RNyUPWdqlJSAlHB1TTc1Q1YPuVWfbsP+SNXXn53ajbyKb BubGDYHcpKgsE7kefa9RrX82gfMj7R/Ur8GJvgXOzL4R2QikzM8/0V19qGZN cqG6zAkgtxw4bHG/FNUYNWfFCTwCUjYxXUQiGdWo/WcjxfYHUp6br7nYD1X9 R1vfJyIFpMLZKyHnLqNq3T3Gi/bngVT8WnOY1xh931I/5FJ6A8itG5XIN9tR 9R6RmPzm3UAq33tcYSiMviexl8W5fABSpXPUa+Qfqnl14xi/NH5/+56D2oGd qKb15Z8dT+KBVH2eyqFcgGom4x1cKk4DuWOOkdMQh2r3Ho2JicR83Gl66aaz N6qtF3tq8sIeSPXMb6rCDqjGtjaOcSwcyN2C8uPZhqj2piy/xVn8/p5r/qmn lFDt5mfGSGsHkBq1v88vrkI1fzNCu0/+BRKUDGVeTqG6sHXP3Wm3gNQMeNuv yUK1fp6KGVVPgdT6TUT3ZqO6UImHad8pIPfqX7DwjEJ1ZqM/A1RwfXWSqkSk 76Fa4cLpH59DgdzPI9tcbo3qjjwLtGQ7A6lr6xtyfh+qG9L5Lty7AKRe0c/D 5GZUu2glOOemCKS+hD6ZSqBab17T3ARpIA3uJ1ccGqUeb5hpn/2vHEjDbprX eB36vu3402O7o4E8qOGgHfIR1W6xGX9RHATkocgKDtVwVJun55bu4g3k4X+b c1i3UV32Y2pvWD2QR8we3bxlgeqOxa3hK1YH8mjWkKoooDoPdRHfAfzecSHd 8TwJVHd58JMkJQukiVNi6mkuVP/06fTpczFAmtbzXFg2iGo4bd6+U7AF8sRW O5nYSlTf05vVqX8JyJOPS/t1UlBdKUPjisJKIE+NSEf/CEL1O2SG7V4LAmlh 4G3xyAnV/3bervQd59cy+YeIrAlqCOZ70ZPkBeQZ3n3N1eqo/s+lgYXfMkBa ndGudj2FnIvm8gYW1gNp/Rm+Krqh+lnRbfz/bQHyLP+eT90xqGEtTXzTAaxX tvbqKUFFqHFzrMgqAaw35wp2xGkNoAbBuG0HYhGQ9mu3PZvkRQ3+N0MTi7F9 h6tbH8fLoQYxl8zESlyf8+UKnseNUKO2iupGR8y/ixtkb3NfQY3JbT8C1x0E 8tJNmSufnqDG1BSlsTIOIC/XStnaZ6JGK6P5B5y4/65u3nhqHQs1FtQUunhj fjneX3+48h9qkuxVfbQuE8hrLaL774qhJiQt8uvTdyCvb127RwFQ453lV3bm XwPS+ZGQcpc1alyIdyqQDAHSpWe1TKAXapLePNLqhOt1cyefmGYiaugVCm/7 /Q7IW0Hk6vFK1LTX4pGyyiEgbw/ReGNHUNMGRwmDvMtA3tXkWjrGh5pSc0Zb jayBdH22fIpLBTVavxd/bBkMpNsfjuGs46jus3v9oQvdQN7XXei2u4GabCWm VNavBtL91VzT2gjU9IGWf8gD238w87eyIhc1erhsuSb3BUjPQ5MFt7tRs8D1 lx06OP9eiX+y5DlRs0bZCLvZA8iHHCNvOyVxPHmzXSIBQD4y/RnzeD9qYlQ5 33LjBdInbTAcHFCjcezhjH5uIP14+v3/+KGGDWPn8qffAul/mv0gJgU1C4pN B/1f3wI+dd08UovtFy5VWGF+BfK1X14xiZp5+hXu1WwAMsiuxSZTCDVNfhJ6 FNQLZHB+04lzO1FT/tONHu2YDyHC9YeET6GmUhdtt9G7QD69UqNT7ooa5x5e 2r49A8jQsupdt6JRU/DT0RyJSiDDxSu2ymI+5T4dcjzLD+SzG6Wb2gdQ42yU i4esOZARNUWiAbyo2VVunNajD2SkTMEqJIcabSVucxr+AjLqXh73mBGVuvzt cpH9+P4LVvZC9FX05oPbyyIpFpCvlLImjENwPY6eHQmaBjL6YcbQ8kzUNMWr q8yZB2RMd3pnBubT9JTexDzma5xaSsNZzCeJ/J7dy3C/xQe+qViD+SS2nfWt 9gWQCYOJVBmgJrnWmWQVPSATIf7jTSvUVDDulVl/BMik8Jg3WzxR09CHRc4W rHfJYy9ftSWipsqx1mu38P9vdSND/SpQkzX1X/hPYyDfvXrmt+cXYsm3iDUs eQKZMhN6f5SJmrYN3w0px/VMO/TE5ZUyapwsr3lozwVkemLgxcNHEUv7zKUc uatAfuDwt+Z0Qaw1Tx/tVMV6mGHqY/ohArFWPfaO18B8+pjmbWSTg5oPL9u8 u9MPyCwej72CXajZjbs55OBhID/p2foOi6Lmb/mxa0bdgPzso1ufZ4qah0YH L/pj/mRXyooEP0Wsi4bm24QwziFIa9s61NTlu3eHjgWQuQajb3YxUH1skkvu Dex/nn/tBFMf863Bbb39HSCpbxnqvQ9RYz2310qdHCALmKHuWUWo+WPCI3kO rJeFh25U+C1DLJ4VLl+V8Lz7GnRy1WkN1NTz3a90FPOpqG7PiW23EOtumL7Z Vx0gS1avj+HJQqxNYTJZ7/mALD26bKgd83NC1nLtbTwfykL6tqYroaa0pkyJ lCggyxtLbnpeRM2qlwa93XE/VgolFpglo2ZdN9UNjS5AVpn48soPoKZX63ce E8b9VR1+8TCnBGo2Ymxy5ygB8lvLwfCm06j5ga6VuSb+v2bd1u7kSMR6rEMT +Yj5UXtSQMa1BbV4DJ2IurIHyLrn01eOCCGWZZ3z5ThhIOvbWZ+kjRGLlC6q YBwFsvG/7KV/j1ELR6vYn2PfgGyyiNr/vQq1hC1//MFDHMjml26P43hQc+5f tz1+YkCyus8039BBrKw0LnFGP5CtG/auN3BHLeYDy1KFJoBss9p0TjwPNQaP RSiU4/nSHsuTMvkPNX/Kj58+/gDIjt7h6XI11Kxg4mTrNQlkl2S1RpQTYhVn dLv/h/Wo+2yq19U0xJou7neNPQZkT0LQN53fqDn5mXrobtwP7B9OQuu2oOag lzu7VuL51Cdz3HzEFjUlnlmTsOsHkP32avGFsajp6R6DSw1DQA4kr/sd2o0a DyhbBkxoAPljeGH7eVEqhnfk4NVNp4Acku26i8yoTqsBAUPOYiCHLxQUrw6l MlfQqFGXr0D+fBdH/KijPn3TSjDJw/H/GvE6msOkXt786SPerADkiKJdZKAB FaGdeKLtKp4no1cO9Nk8pKpl159dCsN6M5YuL6dWTN3mMCva2NMI5J9xphPJ SZ1Zs9WuwP0lkBPKf770aKCP/7Sv24zj+kxea1j+8Tb60nJ79u8IngdTGZn6 jz6hbz/GiKeKI0BOT4U/MZ9EdQeddSqam4GcUb3VprwV1dVO3rnvgu3NuphL cF1CDR+MEsvq7ICc+4QcWpNRQ4U59SfaAch/sxvep/xAzaZ/tfdtwfN4QX35 3ANJ1KSmOH/oI9aHxVsDWiaYX2hVgnYqjm/pS5mPbBRqtmWtmRusBwbHfHLt UitqFmps4RZsBwbnHv+1DVh/qTeHnjw9DIzlrpfPJB5Bzf6TLWcvHgPGCupw 0p1A1Hxx9fBvHgdgcHGo/DlUhVq5jlnGb1oFDG5NoZ1SvIj16F0YD9c/YPBE rfTurkEs1/4xfx5OYPDOTDZFhiNWiK+4wGFBYNCP9kmZnkYs90bmtq+OwCDS 6p0EZFDTzJI00eAFDAZRWPh9DLWYrtRpblcFBtPuPb8v7s+HcGR8tygw+Iqi T+93RazRwxbzXvuBsUo8MHX5PsQyW11kfNUCGKvvuC1QDMRaapQujvoIDAHW ZYNbjahFiW5w8mU4MIS2WTxXxf3Gv261+vUHwFgTaDg0boOa32N1lEbAEP61 Ry1FDrEMY2Kix/YAY52unLf9BGJt3z/fZy0HDJE4kSapbNTKs0XzV7QUMP7j oEv2uCOWtdFwoxY3MMROzl2LOoBaxdavtgujAWN91lChGT9qlTreqvSeAYwN q1v4BVio5aarRTDPE2BsvFxmWfMKtTxLNtu5PgsYEpVZKX52qCU3Z7nSTXlg SEknLOgqoRauipIWMXFgbHoQarB8BrH6P3ZxzZwChnSXZwRFoRZR57kF8evA 2LzLaei2F2r1Pt4v6cUPjC1hNmo7jFDL1hhDBeEDwJCdOOI1IYRaCjtjLnBj LH9QqzGlE7XqSwUdMUgHhsIbZUmHBNRS1qWY/SEJGErcG65tuohaD49eXO3f CYyt1nwF7G2o2emsRdTK1cBQpjj4ouap6D3f7ys5mQBj27pRS7Mi1LBGmZad mACM7c6dKYK+qHXovIW2TTEwVOuq52uxfr1ZwdWxwxoYagq5+v7rECuGd/N2 MV1g7PR5G6HLRqxCrvuO/cHAUB94PrQiCbXKzWxpVToBjN1avjvyr6DWuyGx us6BwNjz4pbXHTXUkjT3Y9wZ81Fjzr5RjQO1tGWffGf8DBhw3ExiohS1WZuV fcnjA4bme13H1ADUwqLU1vjoAUObsaPg/HHUmv95l+/IVWDsdZDmkxbD/Ft7 NJDrJDB0SoQs2AOoNTWU9/b2FmDs37jy3YsUfP7M5Jc1l4Ghe3dy/sR11Ppu YOLywSJg6LX06QvtRk0eq0RWa1YBQ397/bO6FajlTJPNFhbuH4OgwkF/3D+q R1v/ZcUAw/D3+x16Iajlx9anNJ0fwDioF+218hRq1728I1c1DxiH4gMbCiRQ 6+IrhZUH04BhvOyexJ1h1HbxY4RcJebHEfPLjmrvUds5r6N3WE+BcfSzRf7k LdSWZe8j7lkOjOOCRsw0TdSmWaV7+lcsMEyu7rG4wINaF6o+Wuk8BoZptdw7 6RrUlivax/JcA4wTm0X+9YajdnPJGb5dOP6TnvQDLy1R24V+viRfbO9U99yz k5tQ2z1Zng682zMsdg8NCo2gtpMDpbY+uH8s7+c9nML2vg3xZfipAON0yZPN Demofc0R2KF/GxhW9HMVH4JRW93xmg3HcT2sD+0+H+yE2q7MvFXK2wsMm6f8 xNXjqJ2VU9I+j+tztqX/3SE11HbXAHHoHAXGObFsI8V1qHVUoyywCvPNzjpg lPyH3sSeWhyJqwCGfaL141/tqC3H1EdceyUwHH6rKVXmoTY15szyj9j+BWWi NvkVats4MWsWg/vvonOP4yN31PbUbkHcA9f30pfM1XY2qK32Z9OMSB0wrnD4 ZOzXwef50tzKhoBxVcfy2Kb/x99rInpcABiOPirTK3lQG63VVltPERjXvvOE 9uL6JIf1vCiwAcb11R07CqtQh0jOaT19aWA4m6azolNQm+GubInTPcBwifK6 ee8xaj8scaD+4BgwbrBPrrPE/N6//Pg9XwoYt6QVv2gYo1aZjbOc53cC4/aF Faf+24baXnvUa+phft9JY83PC6E2U64z559oAePu1Luo1hnUFnbLaLQQ67ub urvG5xbUrib29cnQBmDcczveFZ6N2nbtfN7KFwaM+0Wy91wiUYe6UdhXFxYw HvBybDjuiv2HitFszGcPw4bC7adRm9dk7ZGVWJ89g5OsBbRQS3CGE01MAhhe za7LJyRQm8gC55DRT2A8FDWOw3xvi1m4OLsD68+jM9I6aQOobf7DDcZxf2D4 xP8beFyG2r96/MlZwQSG73CN96Vk1Dao87M8AmN/xXgZQz/U5rnx4uYWzM8A p5vlchdRm1jNBdbzKGA8/mzkQDdCbSPm3GJup4ERuChBG1ZCHSsuG3CVlgIj WOvvm/JVqHXFk6s0Lay3T7yrDF5Poraz3xKuG2wCRkjVq99eTajN4ad32K14 YITyXw84m4Va26L+S2jF8yTs+AHFveGolWO4vU0Pz6fw52I1EnhftPjiU7XH HBjPuiaucp5CrA+7VLbdYQPjuWQZP57/rDr/gkOe+P1I+8j3+eKIVfPO7s1F zO+olKtHXi5DjRpNU/tZQcB4MaEzeRfrW/v5vycHXwHjldq6p6eKUE2gxsCN nEZgRN8ZVd2VQBU6yMNtfl9gxBR8bV73kKrpPW8a3YP1Po4r3GXWgWqhreZx ZHUBI17/ojBLn+qh1I45cykDIyFQ83OWHDVsPXUydPk9YLxuFDoRyqD+jusL /1K5CIyktcP/nMaowezt2bUZrcBItqAij9RRbP3CMGGBZcB4ExuyRzmD6uC2 SWfFygDj7aBdJ/9TqjnL88jlPjzfUuT3uI45U/UJexU8b/wHjFTHVeu/m1I1 Fr7fX5BvgJGW9mA29xZVeUrpv1rBAmCkj0zVv3lOleRtDxjOIIHxQf7cu2c5 VGHGkG6wJ94vMs63eHt3UnmaVupFJ/D9j8n6Z64vUVlN+g+DR7HeZQ7m7rIR p97zbMq/9rISGJ+kFQWNtagkt+7kZ9zY/89nX42CNfWykD+P/xPmS3bcqnIF D8p/L7+mSQHmyxe2R6xoAhXskREkapsDjFzx6bu0UuptZ9yvBcD5y7O0Oz4z SH2pKj39/CDej6ioVqUfNKq01eFRTVA0MPLbDWiNslTjR77c90lYLwrX5fZ9 NaDao+skWgDz76uZYl76JarqwPXLB9g430Vh0eEvA6n88yNqek6GwChuWu3o n06lp4d+1w/PAEapgKfB7ToqwuzzHV13nL8y47+b7CepyE/qEtcisL6UB9lz mApSpfmRf+oicP4qvre27lOlLr08fFHWaRwYVQyDjG2m1KMtbItNp3A+qw3y AjbeROnb+xY8djQD45uvkh1fBPo0/1X5+kQ3ML6Xx2gu5qA8vVJW2As8b2p5 BER+daD83aJcqttxP9Tt85xqXUQUy8hG1vktMOo9/n4vF0eFi69TBV+dAUbD V/ukLE0qTY9xP/8Tfq9pWduDeCv0VdJDd3bpFjCawdD8iQcq9gl53UbHfGe5 UjvuJ6AyeYdM03u/gNGSu5X/cikqH30aem2lMzBa/8X8NB9EtYNG0jnyeB9o VxcoNqCh+tyQNwdHtgGj44bXS3VZVN++VMrY9RoYnZkzN2UMUBPsq74etQSM rimHI0KXULPhxHNUg/PXo9Iuv+Ix/t4alri2cysw2I5GXONpqOmimUydshMw etOo7m78fTmfnRB+Ae+PfaNbs7/hfa6l2kGP1geMAfnYkFwB1GqubaZ45gow flwQvPRGFbW4pHZItGG9Hkz21n1mgloyHsmu0Mf7z9Dg7EbvG3j+EP7V6q7A +CntMH89As9npyV4hOfLr7PtTdZfUBv1WnGm1QMYv+OM0g7jeVTjsmhzCOvp CDvfBy2itghnk2Ms7M/YBmUbhfWoXeCBmEgq7r8/lnEaooDaadrZvzcnA2P8 hZAwzQq1xY5tj+vAejHR7v1n5gFqHTONiN6E+Ti1brZyIA61WTzZOi2G58e0 2fn4hhLUevTilyPX8X73N6zDrfAHapfd9KnyCd6XZpoOmqXzoJbFZzP+xgeB MSdQoPJyC2pTEPX6W9UEjH9HVEh/fdSx5qEz3W8BGPNBcT9uX0DtD1JbKuzf A2OhRqjAPgC1+67VPDeJ5+0S4+Fzk1TUnjY4dETcCpgcdtx35E+hDunHVe/G jgBzWYGX+XIe1JE5NLnzIQcwl4us3MPKQJ2cZ84FzC8D5gonD7GU06jjqZfe dSVZYK78xrnkQaCOv2p3HircBSa39P1us8+oc8T+0pNFKWDy3OcoULRFHd1M 2XWELzB5W91iVq5CHSZG2Tu/RQKTrrLo3pqH7xdpMs7vAibhf9c6zQF1Rn0R rvvuDExyYH6vlxDqVD1ifXLBE5hMuC118ivq8h/pGuFRASbfszmurVdQF6dg tAvDDZj8Ezd+cIuizkZjE6FWO2CuNvhb2l6GOgu/qqefUgemQIJz4nsn1C2Q cPy7TBQwBZemHj0URx1idj59OtzAXGPm5GBehdr63GMtz2B/hN9P6KvcQJ2r PERqUySBuY5wlOOVQh3yudIt1tieyNk/RGcd6iKVOZLJaWCKUpd/Z7iiTqnj dcLuW4ApJjzyzWcL6nxfFy46dwGY6x0vplo2oS6ravnYuBvAFK/8FbjdHXVO /Snm47gIzI1S56/SFVCXUMKhAof9wJRwHT7c3Yo66zRct5ycAKYky1450xt1 ztKkJl3Dgblp6+Bqv22o65Xjsiv3goAp7Ws7eaYbda9+0ptVfxOYMn39jTv8 ULdTl4SGA35vi4ZNJqmGusmth2UIjGXDekPZfaibT8Z554A4MOX+WLl8CkTd 65sHAqbHgKlwoMc0YDfqKo37t3SZE5iKcad3Wg+hbln3t59tQoGptNC1bmco 6j7RW/Zz1XNgKpuY/2NqoS7FnfzLTuB6q6S1t/WNoG7DPWY6v84DczvtZE52 BOp8LbNktPk4MFWtW6MC96HumzKfuYsUgLkj18z17Djq9j97+tpbnP+dQizL XS+wP7nua+9YAFP9ignw66NuVcG8D6dxfLvKm8QHZlDXt+fNl9hXgLlH4tiy nDjUbR/LXrqB86Nxp54dfAh1tDMkb+V1AxM1GX89N0+lKeTZOkjoAFNTsTZu TyJqnrm0OnYuGphajw55rj6KOg4Ma9ZElwBTm/397BAH6pJa5n06eyMwdXYb 7ct7hzqZFY5WoWrA3Pe0WjrkBOpOcqwy5XEE5v5RAx4HLtT9ojSBz74HmHq6 FUPoPepee8tXfQb7eyBGr0LQAnU7FMef3oXj0f9XlvyThjp21OrKq/wCpuGx /b75WajzbIC3fNEIMI1SSi6E2qDOjJSm1xGWwDzEo2N4gQ91B0SGfln8DszD Z4oUNHNQd4ty3WSvNjCNv2gz19ih7nSDdfM5DsA8KlAw+lsAn9fmE9NdAcxj l6CmsAB179w4lN+F+Xw8LX75l4eou7B0aypfGzBNJnhVM4xQz7ZbQjpT9sA0 237R/h2+39R45mWBBjBPuNRGxrei7l54oXIJ1+tk9vbvL16h7kmtM8z8R8A8 NR/BGXYO9axN2u8njuO10FjcHqiAur+tXtt61hyYlvet7B5Oou6CQ4W6fj7A PP215Pn9bFQluuugrOorYFpxbfl26x7qLh31jVt7H5jWuo+XXcP8+EkL+NAg A0wbn4ltF0jUs+7l5GLIPWCerTY5Z9OAupdK2ImvsR6cY36JMI9APQfUf37w FAamnfH66uOnUY+LvO3+6Wxg2od4cBzchO31bylbhuvt0DSosv8X6rlUqK9x RRSYF9Ya2MIH1ONx9lmYZTEwL55Me7bzFupx29Rp+Ry/d+mFQJUyoJ6J8ouT 7nzAvNx9Y0mWC/U0r1rTZJ8JzKsbO5QlqlDP0VPvNjTEAtPxrOZZ0WDUc7Da 9PNr3G/XXseHC5qhHl/eFT8FVwHTaZi3kiGOeqJO/Ft5czMwneUuLnINoJ74 lef9YACYLpfqlJbeop5XpoPGZzCfbqSrWs84op6/EtuPaFYA8+ZEROgfNdRz g9Zbwf8UmLe3L5YPLaEe58CrvHvWAPPODat5dgli89qdN1P/Bsy72aWKbX6o h9X9ozDjATBdF2StGoxRzyjXhet7q4B5Dz1+Wr2Wcml+ydm52RaY9+9PlBV3 ox7VklMfJLC/7kUm//ISUDd7xYdnj7OA6cGVo5B1AbGFKs7tcKAB01NP/Eya MupJ7b4qKUsC08vXIyRxBvW89P9LWysPTO/qwdJoCvX0e77msFcC5iM+g7kI L8TedfC2ygPcDz7G6fJPDFDPs3M7Brm/AtP3qeBpv1WoZ8EtUEbvMTD9mm8+ 8WAhNie9ujzoJzAD1naU3H2BenLtA/bXYX8fn9KcdT6LemKn/d66Yf4GvkiQ uyyH2Kve8hqvx/UM6qFZnhtH7CeMYo+8M8B8InEp+PQnxLYlHFbIYH0LOVtX bOaG2NqR3WopesB8mqg6Y7wPsR/M9D5WjQBm6PBzWX0CsSPG7+iJXwJmuNyS hXYdYidU10R/KAPms8vWQbvDUc+k4dbIuCJgRqSXFm23QGzniJaoHnz/+aTs XwUp1A0zyrI8mH9RqoFbNv1CbCn6Gte3mK8vbkyar3+P2LtnDPe/fwfMl19M A4VvIHbuado0aQbMVws5X/k1EDvpSFur00lgxoD4NG0FYuvf8tR4iPs71t1z M2clYhcdVdk/g/szrmjo5L9gxDZ9O1+5E/d3ArdhwKQpYi+31t1Vsh6Yr23y yhqSELu7IPCNPtaTxEIlzow5HF8A5/1U7H/y+ujdTw0Q20Bxm8VHXP83d1c5 X49C7PDt3gdNMb/etj5IOzaC2DfIbp/gY8BM2TE5rAqoM7hEbF0vji/1qa2k UBB1I+nQP8HAQGCmjTebT7NRTo5efYYV1qf3h3TDmlVQ96HIF0Q4AuaHd59r szxQz24+udW0JGB+pMvSw5pQT1XuFsW6MGBm2kXq3JBG7GWCW9ZNvgRmVgnp Zorzk36ggphsAOZnCdfPahWIbXXkb3LOIDCz741OrBVBbP41rTp3PYD5pfOM 3Czmq+dO050lHcDM3VVn25KL2Ow2nfmleWDmPdN+lc1AbNbKFcJPLgOTms5o jTiN2C9Lah2McL4Kjm4SuJWO2FUbtawisN4VpocZneREvVzBx/fU432niMnz cNdRxM7X3L7jPp4HxRduForEo17NktGf2XgfKCkf/vdvGrHrlCw6/dcCs0z6 1PZ2XdS74fP7ymJeYJZ7VF/OeYZ6c/V/0a4wgVnB1kiKHEa9aZLPrfmxvlSh 1N67u1HvQ2nNN7S/wKyO2vCfuT/qVZYc43uL7X+bCzbZ04l6y/gVE9bifqgx XR4kpoR61X2XcQVhftR+dKpYvI96z1yP7A7ZC8z6VQMr8P7SG/tr81FZrM8N V0w0KEnUu3yj4yIyBGZjddmNl9dRb0acitd5XP9m2Z3v3UpQ7wOBY8YNuJ9Y D5N/nRZGvcUi7qNszJeWAdFNYI/6NB1OcalgPWrT9j8tno16oxweClfifml/ tRixjEB9Rhfqt4ji+nQsXq7vMUd9nHH+aw2wvned7CELUlCvJ5enkyKeB92f jffHcKA+i3iXB5HNwGQLfb3vfhj1hXsWG37iB2av07YvVjGoT1hMOf7mIWD2 1cZPaU2iPtLO92wA3i8HFNcoSuigvjMnv7+8jPXkh99D++WhqO+EckzsXDow B4dmY3p/oL5b0uV7ivG+Obzfob1IDdvPGngehuv7M65dKM4H9cUDLaH6BDB/ cxoe8mhDfdeCzRqbsH6MWOb5nJVDfTyzGwQ/4v1xNFexSMcV9X2ZjpUpbgfm n3WvFqW+o74gKcuFk6XAHL+xSo1LHPUtemurHsH5nWh8cHXgKur78KLDbAb3 35Ty5JuSQtQ3dGzsrzXeX6YDz/a/FkD9e7ZnF21IBubf383rvc+i/kM5uZ/s 8byf1dc1O5eJ+hL3v1GV/wTMucTPT3R5UL/LBWsxCfz+PNeWahkz1G84xbp8 HvN3wTqSmycZ9d+P9ZRas4jxJfI+ZwnqT9ifm7UG6/fCDde5+V7UX9n7t+Xs DMbuY05/OVD/l0fmRel4n17wPzPyRxTfPz4Wv+oHxmH1dj/VUH+FxCGFp58x jtnL7j+G+l/88sv3xvqw8DbzVJcj6q+6W7NUgvephSzpppYA1J+aXLpktBrj gmeH6rE/74p/aPRrYVxFq6guRf1vxjkEpkUwbrqzt7QP9X/U8nJ1xd8bC92/ 8wqxP5/zNpNv9mD801ItB/vzorKqfNk+jKdq3mfuRP1tHYelnfuAucihJZd2 HPWHJdnf/X//LdI+JCRjf+5u+8njjONbFJQSj8P+sP+zYOqMY7w+NOLFG/ye RNPIxCzGW3gEwstQf+t47If/829x282AYOyP+fbpS8ceYqzxk8cf+/PH2qU0 Gud/Ue+Uu/d/qD/FrPDPAp4/i0eq5++ro/509U5GSDzGFsj5Nvbn0oN7g7/+ f94ubew69id0p1PXWjzvF69tdLj8GPWPSm+sM7mF8d0nffZvUX+nTc/uU7jf Fh+utLDG/rRP9Z/VqMM42Jll3o/6M89P2z86iHHkoLEpJ+pPUm9Xl8ffJ4uv zaqMxVC/z2h6LHc5xumV+wx2of4ysVKZL94Y5+zO32eK+hdSLnN5+GFckqIO Tqi/92vxVgbm22Lt+gz1IFzf/Xd+Ox7GuC1IYVsKGuC9lfR1DH/PLA5wJipU oP78daY9/HifWBxz2ijzA/U3u8zLp1RjPNcfKbEC9XdFDGxIrQXm0koTof/E UX/3O+5/xriflphlgWt2o/7xd90RIljPl9ap0/gxX6cCIvM2amIs+caD7oz6 Z1eeI8rw/rmkKLq4MhjnIylRwFEa450BLksp2N7lcMtO/N6S9tKf2UrU/88w IGnkKsZGV89PDKL+vvUKPlGYr0um7P6RFWiAv6/r3WIjxtZHLQexP4UsudVs nO+li8Ut7D2ov8a9jRDHerd0Y8eR9hOof/gB9zXV/593T6xuwv5M2gqOe3Fh 7L9uf00wGiDKpbaxcb8thfkWVKRif7klNm/B++1S9MKuoirUPxLCPfIRf08v vb30MW8QDXCbfbg0iuu/lNmt+Bn7w6ExSxsNxrjgcNKHDah/Tjjk10M8H5Yq v0q800ADy9bK7duJ96Olpm1Rr0+ggRWS5v6P32LcnbAm2hkNLOf7b2UW1o+l n2uCnj9BA4IgZ8JBx3jqEf1pGhoQ/kW2r8X9sLQ05/m4Cg3813hhU//U/wCs IttL "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{2., 0}, PlotRange->{{2, 4}, {-3.2209814575446165`, 0.6911418068258955}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.419682387174488*^9}] }, Open ]] }, Open ]] }, WindowSize->{640, 579}, WindowMargins->{{60, Automatic}, {Automatic, 24}}, PrintingCopies->1, PrintingPageRange->{1, Automatic}, FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (April 20, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 253, 6, 89, "Section"], Cell[846, 31, 221, 5, 26, "Text"], Cell[1070, 38, 892, 26, 62, "Input"], Cell[1965, 66, 172, 2, 26, "Text"], Cell[2140, 70, 611, 18, 62, "Input"], Cell[2754, 90, 231, 3, 26, "Text"], Cell[CellGroupData[{ Cell[3010, 97, 122, 2, 27, "Input"], Cell[3135, 101, 114, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3286, 107, 284, 7, 27, "Input"], Cell[3573, 116, 25581, 427, 233, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)