(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 14715, 378] NotebookOptionsPosition[ 13988, 348] NotebookOutlinePosition[ 14355, 364] CellTagsIndexPosition[ 14312, 361] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\n", RowBox[{ RowBox[{ RowBox[{"a11", "=", "1"}], ";"}], "\n", RowBox[{ RowBox[{"a15", "=", RowBox[{"-", "2"}]}], ";"}], "\n", RowBox[{ RowBox[{"a21", "=", RowBox[{"-", "3"}]}], ";"}], "\n", RowBox[{ RowBox[{"a25", "=", "5"}], ";"}], "\n", RowBox[{ RowBox[{"w0", "=", "0.5"}], ";"}], "\n", RowBox[{ RowBox[{"alpha0", "=", "0.2"}], ";"}], "\n", RowBox[{ RowBox[{"h", "=", "0.05"}], ";"}], "\n", RowBox[{ RowBox[{"tt", "=", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{"qq", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"wz", "'"}], "[", "t", "]"}], "==", RowBox[{ RowBox[{"a11", " ", "*", RowBox[{"alpha", "[", "t", "]"}]}], "+", RowBox[{"a15", "*", " ", RowBox[{"wz", "[", "t", "]"}]}]}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"alpha", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{"a21", " ", RowBox[{"alpha", "[", "t", "]"}]}], "+", RowBox[{"a25", " ", RowBox[{"wz", "[", "t", "]"}]}]}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{"ww", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"wz", "[", "0", "]"}], "\[Equal]", "w0"}], ",", RowBox[{ RowBox[{"alpha", "[", "0", "]"}], "\[Equal]", "alpha0"}]}], "}"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{"alpha", "[", "t", "]"}], "-", RowBox[{"2", " ", RowBox[{"wz", "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["alpha", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "3"}], " ", RowBox[{"alpha", "[", "t", "]"}]}], "+", RowBox[{"5", " ", RowBox[{"wz", "[", "t", "]"}]}]}]}]}], "}"}]], "Output"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"wz", "[", "0", "]"}], "\[Equal]", "0.5`"}], ",", RowBox[{ RowBox[{"alpha", "[", "0", "]"}], "\[Equal]", "0.2`"}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"aa", "=", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"{", RowBox[{"qq", ",", "ww"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wz", "[", "t", "]"}], ",", RowBox[{"alpha", "[", "t", "]"}]}], "}"}], ",", "t"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"wz", "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{"0.15180194939380343`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "5"}], "-", SqrtBox["21"]}], ")"}], " ", "t"}]]}], "+", RowBox[{"0.3481980506061966`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "5"}], "+", SqrtBox["21"]}], ")"}], " ", "t"}]]}]}]}], ",", RowBox[{ RowBox[{"alpha", "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"-", "0.4237229365663818`"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "5"}], "-", SqrtBox["21"]}], ")"}], " ", "t"}]]}], "+", RowBox[{"0.6237229365663818`", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "5"}], "+", SqrtBox["21"]}], ")"}], " ", "t"}]]}]}]}]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"wz", "[", "t", "]"}], "/.", "aa"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tt"}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVknk81PkDh5mRzZHwbW1JjlJTiNqt31aL96dQqGzlXIrI0Wrp2FU5tvql HButSoVIFEkiLbUoM4xzxj0nY6LG1cSWlZzl1++P5/X88/z5GPkd3RdAU1BQ +OML/3fGGT/l9f6R1oqfPfLm5npQwzPc5W74E84HGP008bYHc3WV0YaGPyPq V4UftcQ9yGAlGswzDMeNX/NUHz7uwS9WYSETBvEwhUVv6KEe8GYOJ783SMXF b16G0Vpf4ivJa8aIwQP8peGVxX4uhZj1Imx46AGKnJ7lJhVIkZ+bwn5bnI/7 Z3M5p9OlcDrh5CO3KYD3Jd97qVFS3FApvzEYVASRpsh8sbUU2ukNtFHREyjc u+zwsrYbWZzM428NyvBYx25g96AEH1JKTXoSynDSzHuRSqcE24O4so6pMuhs PRwu5UgwTJ9wLeeVQ3nFhmXthRJssnLaEh/7HNvs3V1enpKgrWiWxnjHROp2 dY1ADQnOcswc1/uwka6cmS7f2YW84merHl5nw//U+9PHSRc6UrbRVzaxoWJs +0ZtYxeMgzyeL95cA1p8zFCEfhc49AvmClQthMP9xlb/dkLHSqLdWleHo2/S CtdndKKwKF7yizkHp3Y+Gn09J4ZcW57g6MtBkyTl47cTYjBOOlqvTubg27hJ JP0jxh1LtSzZFAfnb+ZZnpCKcbUxIdCzlou6ct117RVinJT9ObpjfzNyE0YY AxFiZKpxo7Xs2uDRP32zRVEMa6vI/ee82xDMUN6TNi1Cd6jpxnen2mDvXXHh 2JgIizsuDTTlt0GNUfO1Zb8IV1J2OcRptoO/tr5sd6MI0StbNOYk7fio8uSW 4zURgtCeOnyCh/+k+f5sbi6CQvbHcsklHu5Lrvxpu1qEFCW9bs49Hl60Ly3y Wy5CQ2Og/gMhDwcMFO881RFhtctsdsAWPkSq+55VzAkxdJhR8FJRgOwQjwvm PCGCr/5e2ZYkRKLdCprnWSFo49k9zAdC+OYtUnKNECLNvUGhqFoIW8vLB9zC hODoUTaJH4T4O3B66MgRIUzv59U7eIigahq94pW7EMMVvLZqAzGupalGyNcL Edq/RlZS2Inc4xXuFnIBer0751c1dqLbVWFjdL8AzuI48+a+Tkw2/UST9gqw iTsY3q/bhYEjpdbZIgHoxTmaOrFdMLPV+XS8VoCUKCPrUz4S+P3hwtPLFqCK WpK6WVOK+DWtIW4+AmxIbKi0M5XC6PTtg85eAuQqn+7bu12KHAu7SRd3ARKm hBbBX773CHw7GvCjAG49yXVpQ1LYt+/7uhICyPM1P8ywXuL+QcsEGAmganSB 0bOpFyR8ysRpkI8GyysWZ+x6caJ2q6ZaHx8xHre/X7avFzb5bkbcXj5oSX/v 8DrSi72pKo9duviY/jwcJMroxejzV68Smvl42+2a10p7hVyWUlNCCR9NKYw1 rKZXsC7uCs69wMdlTa5J1kEZWn0WEUszPsq+X1bqFCrDrXkstYA1fPR5H8Vs pAzfdKlvSlrFxw+PKBf3mzIYRwTcGTHkY8hh/+8LWmT4e96i72q+5sMmeqQ1 /Ic+mN2uu949x8Pk+MKwfTr94Fo6l2cIePCTuDDpTQPQ0k43vhHDw/HaC35H GXJE7XMwrZzpQPExdefFDiPw+XbIaCy2A1mRR+Tr3N/jzg9by/cv6sD5zNm9 H4ZHkW7rtme2sB0KvF7yPGIM7WL9B1p72rHqt8DTLibjUGXvXvluuA194zfd 1F58xD9PrjgfuNYGx/7g368en4TNPw3XZza0wakvyVRffxpLBzfmnytthZzt d67h0QyWyMji1VqtMNcSLrPx+YTG2TfnR/1aoIpx3+6xz/AJiAkdZTdjx4jn xHJXBeKl15AYsrQZ4mp/RU8zRRJXNGtdEN6E0rgab8NxRUIZPvaZ6ucicF5Y b0sLjdTQVzyzcuDi0InIUvU0OtFgDJxULufAIsi0ZZevEjkWN6YUuY6DVm0z tbDN84ju7b+2BtxuREXh+VtmdGVydEMb/d/FjYhUNNnTJFUmbz9NRh681ADV 7gj70JKvSPZAFD+E3oD+XKPkzKj5ZLmppHtHbD2SVy9RKnNVIRUvuBMR9Hok pQXW+xqqEh3rj5/UL9XBM/lAi/qgKlGGG93vmzpM2uZMZVaqEUGF6uoz6bWQ 9W87XR2rTp65MgKszWqRbbV+abPvArJy/6Fgo9IaPLoYphW0VoOwaZEWdvY1 iHU9tJM+rkHk8rq+pz1sbOzOehDbvJBYD9pcKzzGRsqaW4+ad2qSkkUTPVZq bJxlJabsLNMkFblL2TmZ1di91DWHY6xFntbu1XPeVI2tofVeWvFaZPiV7cKP nVXI1PMzbhrXIj33qMimsCoYs+/fcPLSJjlHf9az061C0tpqk45KbbLLO279 /moWAseU2htNKBKgFOM27c+CX/G/hvF/UiTTjH2/YCELv5x3nnW8QpF4jRVl 32mwEOZcIlK/RhHBmOZIuToLceNhl5NuUGT+sj+eN6qwULh5auZmBkVyHVKn hugsTFXNCXMfUmTDTHOH13smkjrUE2vqKeKof4C6W8VE2t2QwzGNFOmM9Wsu q2Ti7m8tNvZcimR4lnPbKph4qpM0zW2hSJa5/fW5Uia6PanDPAFFDriu9TqY zwRDtsTmtYwiMfdyCsyvMbGuJEL/Xj9FFjrI8hySmNhyUTIVMEiROf+Fb/wT mdjFyHj8Rk6RYl3djemxTJw4Yqg/OkoRh207ZxZEMRFl+d+pJ2MU0T3kpGMa zsTFBa/5v41TZM9mjTD7k0ykFN29NDlJEZbFw0/Rx5jIPqcUVD5NEaeBSEZ2 CBMP9wZsi5qlyLsz75NZwUyULK9bZv2ZIt/X0bb3BDFRObZqam7uSx9fsOWz PxP/A/SFPDg= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0.15}, PlotRange->{{0, 5}, {0.12263478214712117`, 0.4999999183673652}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"alpha", "[", "t", "]"}], "/.", "aa"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tt"}], "}"}]}], "]"}], "\[IndentingNewLine]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV1nk8FesfB3BbKkuW0+1mSbSQJakoa58pZS2VKEm6qKgk/KRNl1LpVLKk QkmEIkt1uSLJVraDYzsODmecc6wRkkLFb+4f85rX+/WamWfmeb7fzzxq7mfs j4kICQn9po7/zgl/u4uvP3ppiyzvq9GjuGpUNKvuPKB6ELfWWEjmTVVj7lNx qKrqCdwMlfsse6AGCSXhy+epXsC12hWeMXK18DY7e/rHcjrs9p7KH7vMQPNP r5ix5XH4ltgep2nQgMBKFb3R1/EIDFssJU5vgOK9ltoR88cICPPMNOI0wF2L EP3smYjfQVp5jgQT5zMlC762JSOqPyJ4hGRifidPY2R5OvwNtq00pjWBXfL+ 7PBAOuSNH7/WNGxCRlps+efXGbhT3Wdc79IEO3+7I0PmmTDJDPMzTG3Cg4WF D/o9c8CWW/346YZmyD+uEhlvewMHm17xTUQLkmoS/T4vL8CcYfrYhCYL32Lz tLh3CpBnNWPTa8WChWctv2m6APrq+dJRniwMi/5wLGwuxEorRVHFFBYMzeyM 6WFF+Ms8aERIsQ3MnF8iGqMfIOpmsjZuug3yWOgqrVuCAw8L9fwWs2GiHmMc drQELmKB05Lr2Lj7LX3iYmMJHpgH1t32YEM/qvWY+8tSaL/bP1xazUZwjY7N +iPlECmt3OF7ux0vXuerv7xfjns5QfJZye1oit0muppRji47ceHcgnas8nQq WmpUgc2xb2fU+ttRI3pNV4j2Ed71QhtuoQNLzDrlGz59QvKH5YH9fR3IzqF3 euvWQCQ/3o4uz8GQ/NAdG7caZPxaXPByBQcagTZb1sTUIG3jwcr0DRw8NZVM 4k/X4JYSXWBkz0F09Z3jzh9rkf11OO51BAeB/IhxS5c6eK2pkTUS68Ibi7Hk 1RF1YMVd0TGQ78Jo+h4H0bI6jETenC+n2gUvX7n8Yo16nDZ9Ne5l0gXn39FB BhP1CNi41vCpbxcSJWtD5XYw4ZT7ReOv1i5sMbvkEuLKhMY/V1cH9HSB46Nt MHqOiVmFwbBTI11Y2nS7j5HBREFCvuO8ed2Iit1pfVO2EW7SD6eGNnQjdHX9 ornORoy4srUy6d3wRGPcsH8zMpWURd1WcSGU/L2w83YzWFkt1X+u5SJWTJlT k9IMy4VFY28NuKiqPq6SzmpGhYaP8DsLLtY4/Eo+ZtwCj6377pZ5cjHgpZHZ LdyK7o6xo8JpXITU7qyrU2wFu2z2kW8WF0t1/b8UbWzFinVBjxi5XFhNFOk9 OtaKazLG12zLuUi/bJ93oLoVAg/7Hi8uFyejLxczI1nQv64gpSpBQmQymfsh nYVR9W/LLKRJxB+oEsopY+HcGf00D1kSNco08/BvLKSLaDtfX0JC+/mLSmun NnzdaC95bAWJ4XfNzLLlbOSqbVv0zJDENZXp8deGbNxIlWmyNCGhfEWFlrSX jToTE6FeMxI7LU44BoeyMdbMaVlgTiKLOdtu0s8G1BsHf+wk4dOryc/NbsfB L7zS4CMkSNf2BaXV7ejeKD9DupHYx76pWydoR2vKNgvjoyQMa/sv9Cp24IwG /VCrFwnR16myS8I6YNxTZJjjRyI2SG3LuSOdyAleZjd3hYTEd6ZH6MVORFmV qQyGkrh8JoQecb8TTpefTtVeJ+Huzm15XtsJ690bhS/TSehYJZxib+ZgUeK9 0fBIEqU0hTgj2S4Em17XP5dAQj+8qniHdhcybHzlpp6QSBM/L9hr0QWD/KcN /k9J3JlmrTsZ1IW5It0V+5+R2M+N+RQ/0IWBQdFLrS9IVDltH04T7camY83v 1TNImDRNyP2j0o2b+/7c7v+ShOpH+8O1Dt1wst1mM5VFYihD9tvPkm4o5HCu lr4h4bKqRGEBpxtx715r9P1DoiHhDBb/6Mbs45Xe8/JI5EbW39Kh6m6KZSqk n08iOPCO2uE4Lngqczfs31HzoXZNg0utm0amo+OTMup9TKPW/b2DBNuYO+dR TuKG05PNy+xJhPL+CF9dQdVL5FvLQ6dIiJ13OvP4I4mZ2WHPNmoe6rujdJ2r SLxVnDkTSH3HHz931cxRDtw0//wf1Li7u5puJFWT+Hpa7eY+Jgn++fUh7TUk PnMcXzSI9CCVqgblOhIZU+6vfGR6UHN6U+4zyl6Lfd9KK/cg82F+iEY9CYHt rSpbgx58uvTdclUDiWTPh8yhrT2I5f2zJ4HyX6EpbLpdD3yPVTDlqXE4hcUD lZ49ONq5y36Ccjyrdux4QA8G5VzdXBtJOH1lT8270gOO0pjIJ8otmhMLtsf3 wFrddOntJhKMWA3NEkYPki21zEeaSdzK1V9/pL0H2584PTFrIWHF3Go029uD 8xzJkluUK+a7WJvO9cD4l5KvUiuJqytP7O2U5GHzpwGJI5QJBB68uJQHtbOK 1xMpFwVGnni7ngcHDSm9P1kkLkYn+B3YwkOWLdfbjrJhdsaF7zY8bC8zS7pK +Xt1/pX7B3hoqFnc/A/l3N4Kuv5RHirPWSzgUfYXbopq9uUhpDzdQrqNhN4y bpz/ZR7EUiViDChnOUynv3rAg49fsf9lyt6+4m92P+OBO+Evk0BZ6w6t8EsO D4MeA40FlAeeq5aFF/GwnpP7oZlyWvnaGp1qHky/eHQPUT7KNW6qbeXB6t/A jXOUV/y07DjJ42H6i0eFLJvq3yWOvIWjPPDEYh8sp/xkg/vQi588TGqezdGm 7GJ35qvlAj5KvPYtNqCseDJopm8xH3fyQ2pMKLOv00VuqPFxUHJrGyg/SHog sVqXjzcn67GVssP7Z/IVxnwcEzq4kKAs3/5K0cOSDynd+UamlJnf3q8QceAj 9Mlk83/Pvytbq5X0Fx+aTpu5OpQLNi/Ls/PhI61NaqcqZYHrGfy6xIfYuEL5 f+8ve6O0Op3OR79pts0s9X0mWTSHAw/5MPGRGhygfLzlWLdYKh9h3idTmJSj f+Z7vXlD3e8zfjWP8vsVEhNHSvi4Lld86+F/82ntclm6no+V58nSQMqL/bLn v+vko8DVRXcfZcQKRXsN8qGzWZOlQ/lBX0pauZgA+Sn7R9uo9S2V/qHnJy/A 5bAT3umUh/Wt36moCpC2ZoPpecrmoSMNF0wEKH/5aUiask8GnDWsBSg93NfS QtVbfGOUoGW/AOwEjzVxlMeWb5pe5y9AxGiimiJlZcuboV3BAgTuOtDEourZ 0qdD+na4ACuFDotEUX5S9PeKvhcCbAuunvxN9UM1n5l5718B/vzeUP6K8jeJ lZu3VlDXmx9Vd6Nse7DS9jFXAE3a3PtCqp+mJmXO2i/phZOS9qe9VL+tXOY+ N7eyF4VXFKrGqf60255Lz1rfi48Kf4ZFUE6NdnqyYFcvtk69311B9bf9uqTK 4tBeyASIRNOoPMjy2qCoPd6LtcwfwzupPGFHXEthz/bCueLqtxoqb8TyWbo3 pPpAjy4UsaTsPO+SeY9GH+w7HbcbUfm0ILnc+6FrH1SH9zeJf6Lyv9Phgyij DwlWV7KtS6l8fNRz2q+9D/vvp81mlpB44+yjzO3rg0IBLW0RZdGOsAuFwv3I OS25sK6YRAq7cKPf5n4crKu5sqmIRF+r6ovu5H680Ao3/vAviZPM4ciCCwP4 7EvfVUzlpHzkBWiEDWDj4PGxr+kkCneLf4mJGYCouPwedcqSDcttfXMGMB41 sIf+nERmnf08jd4BbLEXSyVSSHypeXshZs8gPMtSBgKpHPb7eM39jMYQOlUj XGgRJFQEFuLfNw5B8QXzhf5dErUiCzOCiCFILnVVdQwnsZoIH6cfHEJZgUvs vdsk2gvvh6TcHoKI3gaO+E0SW1+lPWkfHUKcqLR2VQgJuUdVHeZvP+OHcIFo KfV/fe0rtW+p9Qg6f9w40bmXRNKlU0N6B8ZwSPf6/1LEqfxL/LX32/A41K0C wmfOUfuhZpIoujiBDg86/V1PN9QDjp930JpExN098e5m3RBMPtwv+f47hEXP 7bia0gWb3pOXo/2msPNVv80FoS7YCSK1VVRm8FbmzYdGN2p/We4eUpX1E+k/ WW0OzE7oyrGWmR/5De9zV5V8DTohgUk3zsQshGTMbLOfd8ByxPnHCkch4vbn CwW6ch3UPuqosLOOMFHjs2swIqQdeTcrXFUnhYlEGT9N719sHJ93lqyvFyG2 qp8OXuXLhof/pTypeFEi087Sd2S0Des8tet3uokRI31d6z4GtKFBXkfyrNE8 4uPaQ30T4yy8y776SEdUnIhQtlH2u8jCJWGtPYwuceKPJ0Td85+tkOBctPLJ nU+EJxXt6bvSit40tZjEoAXEwMf5RXnyrYhZoyBW4LiQ2Jx/xDPjcQsi449X uqlKELK/s4Xr9VrgHHO4Xqpfghj909vaoLoZU9tTpxOLJYmAG+3Rqw41g9+7 7XxZmBRRNWKf9Oh7E5LN1ivVuUkTSbR7TcO3m5B1/ayc59pFxKHYzSxvrSaE OXrYik4uIiQWKVVXMhphwElKD6uTIXqP7Cl0PNGIWM1HWXW2ssRXIw/xS4sa EVwSHmtbIEv8HVy1e1k2E7uUHFNrVskRdvs7vkbvZWKrT+UhOboc0ag9/+Wy mQYkKruvYkzKER/qnhq+sWjAqvLnD+wOyRNbTb8I3ifXI3JtmVZTsTxBM5RX 3jGvHscnxBqrtWjEHvGAttnDdXB//VWVHkEjzCS0t6//yID31X2/bKJoRGPX Nfm7ZQyc3ZfbJnWPRlzwr1D5/IGBm5Nn70Y+oBE/dLrHUgsZyDaa/vkwgUYk ru66qPKKgenSOVbaSxrxqHGxmtIjBiKbpMIrKmnEH+n8u9v8GIh/dtrrRjWN 6NATlnnuw8CzgHpzq1oaEdC3myHpzcC/SyJnautphOmpk7OtxxngONO8mltp xL7B7zjtwoAGX8Gcx6cR99vCNj63ZEAv96JKSi+NkN2w9JrMDgaMr3dOH+un Ec82hcid38bATo2EV4NDNEJhtIRnbcaA/ylVlfFxGtEje8J1bAMDQaZXpt9M 0AiTTMX7znoMXJfmtQRM0oiZv84t+LiWgdicZ7enpmiEMFspNX4NA8khYp6F MzTCMMipfp46Ay/3HtsW9ItGvPutqe+3koHcFZ+WbZmlEVX1t1gcVQaKJ9Sn 5+ZohIr72TwrFQb+D7F0ANM= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0.2}, PlotRange->{{0, 5}, {0.20000019387750045`, 0.5263983007638073}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output"] }, Open ]] }, WindowSize->{1202, 705}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, CellContext->Notebook, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 1440, 48, 232, "Input"], Cell[2010, 72, 617, 20, 30, "Output"], Cell[2630, 94, 206, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2873, 105, 328, 10, 31, "Input"], Cell[3204, 117, 1365, 43, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4606, 165, 197, 6, 31, "Input"], Cell[4806, 173, 3307, 60, 231, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8150, 238, 244, 8, 52, "Input"], Cell[8397, 248, 5575, 97, 240, "Output"] }, Open ]] } ] *) (* End of internal cache information *)