(* 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[ 28507, 561] NotebookOptionsPosition[ 27271, 516] NotebookOutlinePosition[ 28002, 542] CellTagsIndexPosition[ 27959, 539] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Assignment 1", "Title", CellChangeTimes->{{3.433567408527817*^9, 3.433567415725386*^9}}], Cell["Due: Friday 7 November, 2008", "Subtitle", CellChangeTimes->{{3.43356747830616*^9, 3.433567490982802*^9}, 3.433583481267744*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ "SetDirectory", "[", "\"\<~/Teaching/MathModelBioSciences_Winter08/Lecture2\>\"", "]"}], ";"}]], "Input", CellOpen->False, CellChangeTimes->{{3.433567727678491*^9, 3.433567749540262*^9}}], Cell["Consider the 2 species, Lengyel-Epstein model:", "Item1", CellChangeTimes->{{3.433567500149698*^9, 3.433567548592442*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"GraphicsArray", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"Import", "[", "\"\\"", "]"}], ",", " ", RowBox[{"Import", "[", "\"\\"", "]"}]}], "}"}], "}"}], "//", "Transpose"}], "]"}]], "Input", CellOpen->False, CellChangeTimes->{{3.433567717141712*^9, 3.43356772056262*^9}, { 3.433567770112527*^9, 3.433567843270501*^9}}], Cell[BoxData[ GraphicsBox[{{}, {InsetBox[ GraphicsBox[RasterBox[CompressedData[" 1:eJztmotfFFUcxTmApoEZIGaI7xQtUySlsrAwtUYrENFSU8oUS0uIQElNy7DS ytVIyVc+krQiClHxUfYQU8BKRQldHgJz/49m9jHssPfymcGZ3UvM9/Nh7s69 d+Z3OJ+ZnZmzM2RRRsqiwICAgGnSXzwCAoiFhYWFhYWFhYWFhf/JmepvBZ2R 7Cn+VsAzv70/Cpi1+Wrbfi/X7Au7IyaviZSvGgCMyLrsI328UoxQ0buXcqyV 9F7kaBfjritmi+Kecwin9NLO0OKIk3ITizSTJXUCyjW7Rvb0ryCkDPjZbE38 o8M18t7Ia2QuHjJZUWegjWvp8xzERjvbeerrRNqjl0PwqU/18Yku15qnDkLo DZ/q4xM9ZyghN0LxorIyHQ6+NUcY1+hz7UcgV1lpmLLaZrN98Io5wrhGn2sp GBpWqKzdXHhNWp6m3O/939HlWvUd2HQ47LSyfn3uvybJ4hxdrq1Hr1qyJapK 6bg03W6SLr7R45o4HPOlZsXo1iNsp9AFz8+qljKqa/SkqAgolRrxuYQ6d5c4 aptJ0jim4u3BVNeoND4GnJc/HEF8pbtz2fAueLCRi3Pv0TizaIEgCKlfErJ/ jiAklxNyXO5dg66ZGv3T0Q3TbkpnaCKqKUP+zILLi6QHmlvHeH2EKUksKn4B T9GGjMmCOxSFvgn0yRE3o2+pSQVulz39gMnUUkzX1q6k9TLV649Ca7uhRNy3 ciu2LP9Fw3R/ZK3N5/6mD7Bcq+sbfpM6wFLPikJPZKv5yj3wB6KlZWkm0sWD Wv4FnrJWlmsbgHfpIyz1jCh0+xA16e6BaoS2SM26wPEapXKUtTJca4hCUr86 6hBTvd4odBROSUt7ZHf6Qe0NP1krw7WPEVY7JY+xDUu9zih0HZbJzQp809pX Pk/FYQ0FXlJt4aH5Aoxiu7LPdlPNW4OQQ0qiGxn/L8MenVGofWCEfJSdxdLW vvZdoxZgumYG7bpmw901hCR8wtiWZY8qCnXxXaoaj30WBWdISzE8VrNoWgF/ QD1Dm4chU2q+H9zE2IqhXhWFujiQqGatx9jnPYqlZUwg7fabCq2AP6C6th29 5BRTjM9nbMVQr45CNVAY80XTqTjs1jpfdwGToLnWMhJvOD4cGtFM34quvm0U qgF7/usfXs45oHF2BwqYA821XQhxfsuJY3fKzbGejivJw8oEhnqvKNRoTC+g FYpr4oPKVW3vaEe6VPikzWbbknLWPYOunhKFGovpBTRDcW0ferrTlZb7nQ9A u+VEs+WMq5ehnhKFGovpBTTjnRSJ47BYWSkY54wyt+3wmEFXT41CjcT0ArdD IZDxkZu8Hq4b98y9ygS6eq8o1GhML3BbpI1V4XpGb4o77hrnWz1n7LqP9XzV iTiKVLkpHA28eoyQ/BQgaLbzq/vGO7FAyGs72t1eN/agAmN36A9crpHGmaiR 21q47kxlGobhCQNrnamXFmLIAg1TbRE2Awsbjts10vR0snyZO4BuHmltNqIM /Bnv1AZHwZc1TJ2D2cbVNR7FNVI3Xj7IpuEZj9ETwK8GFhM2luVHBR7VMPOK 91vqfHBnsEwQEBzsDO6vDN1EKoNwyGNOSyTWGFiyZkYgwjr3T/f9Zi+ReBYj lyxxRcfnIg9mYIDqQXsWJhpa9Orvtwzdn8/ZP0P+yipWzlB5JSwCb6kmFaBb jW9l8c58+VxRuUayEKh+xigF1vlUFPfY4y60dW0SIq+r5mRjRNQPvpXFOycS mtWulQNjHlfduk9A5k/RFb4Wxje5uWrXlqJ//aRUjzu06iCUkvwxXfP1SBbN k056utbQB1mk5gGP60EBwqVLanpSV3xjjU3FhMMerhUg+BIhlfduVXpmIUla Nj2y2g/aOGZHnIdrEzFZbk72cr9lIt3kfia3F0MPeW/alUlude0IMNPxIb77 Jucp+TVcTwqJvcv8oI1faga6XKtfNUMQhOVVpDpDah3vAW17XhBSNxLyZ5Yg JP3lR5H8UXXe3wosLCwsLCwsLHzPfxV4xXA= "], {{0, 0}, {310, 70}}, ColorFunction->(Apply[RGBColor, Part[CompressedData[" 1:eJztmMdunVUURg3iAei9iDegiV5eAULvEhJjQo9EHTJMcbpThBBCCKHEdood 0puSEOLey+3V9970CgH23voGR0sIJKZZkzVYcqJ7/f/n7O1733lv3rtXt/3N VW1t17T9G8efudaZD6f9a1j9VucbOP35hbD6jU47rH6n0wGn//56WP1651s4 /f+/g9Vvdr6H08//A6x+u/MjnH6+n2D1u52f4fTzb4DTz98Jq9/gdMHp99MN q9/kbILT728zrH6LswVOv9+tsPptTg+cfv+9sPodzjY4/f38Aqvf5WyH09/f Dlj9HmcnnP5+d8Hq1zm7YfXFi4zdcPp87IHVlyw29sDp87MXVm9fYuyF0+dr H6y+tN3YB6fP335YfdlSYz+cPp8HYPXly4wDcPr8HoTVVyw3DsLp83sIVl+1 0jgEp8/3YVh99SrjMJw+/0dg9Y7VxhFYfU2H8Sucvh9HYfW1a4yjcPr+/Aan 788xWH39OuMYnL5ffbD6s04frL5oodEHq/f3OXD6/vbD6s85/XD6fvfD6gP9 Dpy+/wOw+jxnAFYfHHBg9eedQVg9zoNBWD3OgyFYPc6DIVh9eMiB1eO8GIbV X3SGYfWRYQdWj/NkBFZ/yRmB1eO8GYHVR0ccWD3Oo1FY/WVnFFaP82oUVh8b dWD1uK/HYPVXnDFYfeUKYwxWf9UZh9XjvByH1SfGHVg9ztMJWP01ZwJWj/N2 AlafnHBg9dedSVg9zuNJWH1q0oHVY56ZgtXfcKZg9Tjvp2D16SkHVn/TmYbV 4z6YhtVnph1YPe6LGVj9LWcGVl+31piB1WdnHFg97qNZWP1tZxZWj/tqFlbP zDqwetxnGVj9PicDq8d9l4HV5zsZWD3uwwysvnGDkYHV477MwOonjhsZWP1+ Jwurx32ahdXfd7Kwety3WVi9c6ORhdXjPs7C6idPGFlYPe7rHKz+gJOD1eM+ z8HqHzg5WD32gRys3tVp5GD1mBdysPqpk0YOVo99Iw+rP+jkYfWYR/Kw+odO HlaPeSUPq3d3GXlYfWjQyMPqp08ZeVg95qECrP6QU4DVX3AKsPpHTgFWj3mr AKtv6jYKsHrMYwVY/cxpowCrx7xWhNVjXivC6h87RVh98yajCKvHvFeE1c+e MYqw+sNOCVaPebAEq3/ilGD1mBdLsPqWzUYJVo95sgSrnztrlGD1mDfLsPoj ThlWj3m0DKt/6pRh9ZhXy7B6zKtlWP38OaMMq8c8W4HVH3UqsPoCpwKrx7xb gdV7thoVWH18zKjA6hfOGxVYPf4eUYXVH3OqsHrM41VYPebxKqze22NUYfWY 16uw+sULRhVWj3m+Bqs/7tRg9Zj3a7B6zPs1WH1br1GD1WMfqMHqly4aNVg9 /p5Th9WfcOqweuwbdVj9M6cOq8c+UofVYx+pw+q/XzLqsHrsK3Ow+pPOHKwe +8wcrP65Mwerx74zB6vHvjMHq8ffwxqw+lNOA1aPfakBq3/hNGD12KcasPqO 7UYDVo99qwGrX/7DaMDqsY81YfWnnSasHvtaE1b/0mnC6rHPNWH1nTuMJqwe +14TVv/zstGE1WMfbMHqQQtWj32xBat/5bRg9dgnW7D6rp1GC1aPfbMFq/83 /xNX+pV+pf+f/hc5ExZJ "], #]]& ), ColorFunctionScaling->False], Background->None, ImageSize->{310, 70}, PlotRange->{{0, 310}, {0, 70}}], {193.725, -41.6}, ImageScaled[{0.5, 0.5}], {310, 70}], InsetBox[ GraphicsBox[RasterBox[CompressedData[" 1:eJztnHl4FEUaxnkDhAA52CSgSQQhAcKxZBFFYMOxHBJARGACq1wBdI1yKocI ghrBABEBFxRChOAqyAKLIOdymyyiAhtR5BQ2QRAQAjlIQjJTz7NzJDNV3TWT 7ume7n/69zza6W+q6yveqamu+r7qbjZ+6tDxfjVq1Bhk/e811KhBDAwMDAwM DAwMDAwMDKRy7u1KLlVX8v4xLdrjmbPXtPJ0vUqXkypXvBOVHKim4N2+2Sq7 9oLcHhc18vTfKl0yVK54J4ZusfOb53IF8btV9kxIVojsS862+0X1ZnDJd6iS 5APB35RS7MFTi1V2bCUrUP41e2Nvqd8Q9yzTS/CX+1u895GfPhCIX/mT0C4W fHlDBC/OJcVLuwKhM07xKnuj2wPvWyIbvQTf0OC6Ii9l9fC92Mrp4dfj2tuP X7m9qZS1m6GoKfLQSfDckJUK3UTiR7GRN6T8FmP/F85AO3c/qeyaRxQ2Rgb6 CG7p36pcoZsoqYKTCxF7CCkNx9/d1pUYfV9ha6Sjj+A78E+lbqQLTr5peJJs RN07buv62e9dpc2RjC6Cl8dGVyh1wwqekWxnUC3HMZldUm2NvNwdozxUNijw htL2SEUXwT/DEsVu5AhO0poCnsbpA5imuEES0UNwS1ytahZFEpAxpFg99kQL 5y0z0bHgO0h9bm4SlK+4RdLQQ/DD6KPcjSzBi0KQ4BS89Om0zMzMFVPpAtPx gfImSUIPwUdihXI3sgRfizpt5znPiiba7p859CwxG20UrMPkoIPgRfVwTrkb WYI/iefzotY5T/NfuicoUB4MtYN4btBB8K2IUKE3yRH8FHCU5IS5VppXE4sF RQZIC0goRwfBkzBEBTdyBE9Ga+tXvC/stNOy2ST4zlPQToVGSUB7wS1RWKDU x8375gZcwbnh2YIgLLIdV0U5I7GWtp+yZXYBymdOUtBe8IvAl0p9FH7UDTzB uVjmAp/Y/iipFelMecyOZbv4ZWCz0lZJQnvBNwJnlHspSA2TKPjF2SaTafgi M/lhismUuJeQ72zW93GTKWUOgDYxQ+0FnwEUqeHnnrcBp8lW95YE3GatLdGT U3bBM146cYtywYvPsnEIWnDLuUPnhDOSAQhV6FEhxxKOZI1Hf4G1Fx7mlJ3/ tMhUlJ2tZFXKCn5uXiQQPf9rGRUsbjGzY3s6tE8Jvr/lgB1pvTLMzAUxaO06 KZ7oj9Yry8mF1CZAqxT2Z+4rtkUA/YSuRgDCyTnhCb4xOHz0kAl3KUtxqwKu m8yhAMYuL2Gtwh5+vhbek9LoKrbjODG/4bfFZaEEn4KV1p9+n6fowKi5DrrR FRxr8Kr9OA0hmuUWK86LZyRTAFHGjiP41Tr1r1jvC+OoAWkZ3GVnx8BftOIQ DSm9MVJ85b/ms+Q4P/lrmLVKc5cw15dMCb4TtoBFYZsOVBe4AQxk6s4Kt2ca n8BkN83WhnnAPrFVJPhaDLYd7h13Wkoi0Ei4jKokC7ggtIkEX4KHzMJCZEoM y1bnJ/3tI996fO60UIIX+tvziccxwlXVz8DzbOXbHvkfISeB00RP0oANYqtI 8IV4WWD5COF9l/PrLP8DPhbaRIKfAS8765ZZsA0EJ6h4Mn3T7O74hgdRsdD/ AOMFdSxvc4e8hE4yvPqAFQAnzSoSPAOJrKGsMVK/iyrlVzoMonuASHBLE7wj o505+ND+/0lOCy34Akx3mFyjyAEgWVjJ5B43Ax1rE4pXkmloMaAWdJdOB9Ko 05kOt080rvT/e6X9FzRmh+UMhN4jCav56qxDPeHkVTwtTMaT/Kv5JD1svSVu ojK0tODfItC2G6K4ZoBz68ceYIKwjopnmiKoUGB0L7hP+ARMxMGd4ORZMBv0 yqPxrnWwbsrf2nIN2CMwiQVfBL8Twgs/TmLJcn1U2DWhlPSDa+cmLXhFaOh0 W4cIRV6VaRdHcFIYhL9xW6wda4EUsVU8LbzcoJ/1X2T59JDj9B8Isc3Ke2Ry Ky0PxLOCeYpY8N7o+qgwpZqawLKf+qxkWp8MmFznzEozERkpZnIa/s77+D7g FVHDsoGF3BZrRwa3CZyFz9Eg6/iRD0dcsiIW82zHAy24WfEs1O4iqFUkeKE/ Di/pVianraf96lMbIhnB12D5ru5jwjHTaTkETmcejeZhnEmZlqwCL7PNEZyc iZ97x9LGsVLdhCD7jNzSaSOv0jnocbM5u2dVJPh2BJdZXhR3Qg+8ADrWyQh+ 2Xq7vPb+tC9dv6sTwBhhBbcDkL4/LEdo1pRl3H3EPMGJJXvBm6vtN0PzHzHb YdvZ1j6ZFiSo22MROdOYyW+JBE+2pQfK4tOlN/V7dsnCBq9aBAluJ5cgnFjZ GhFcQNZF5QntWpIKbBVbuYK72Ib6lSECy5+22Q5sgtp6z/zBKkgbOhAgFNw6 K1xjPfwaJX2z/HMdH5BL/Hk4IRMh2CJSBPQVVGCJta8m5sVxghmaMQvIEls9 C27p4IrpbnnM/itmEtTrEGn7Y85AaikpFNy67rFPODY/JLm/xawl5IMmzlNW 8O2i9E4oOgosRwDbStnyXG/ttviJSAauiK2ew7O7UdcZlDG3dozVdIJ6GEbb Dg/iKEWEgr8OP/uEoqJ++9+JNHrOJ/nNljpPWcFvicLMHfAoayjrCdif+ziA +FyJPtVnMGrL3SZu6YwOq5wM6ezo1q4E9Y1AzLL/sZBa1bGCP3jHZDJN3EnI 19NMpklEGjktu0VNds2KGMH3D0WAoNuORAAzNT0yweozaZP1wnEm00hRsEcr OtNRY2n8+DhLZQqqKkF99AWTaXjqLVKaNsJkOlp1kRoZH8sNWlO2h5uv7BGE 6xdqla6VR4T4Xu4logQ1g/YptoMQ3kdVYF8c8Nq3hHw2Gqg5dofdJiubUsyG UhQhTFAzaC/43ZpY57GAV5SNgn3yVQDKvYxsyilATp7LHdwENQNX8MxwfmhA EtVuBOqEV72v3S3lg0fa+tUu1KYGLG42hct6BMpaXbuBn6Cm4Qo+DmO991mt 4HMR733t7rnfxeZ3EIZSNm42hctEceTaG/gJahqu4LfWSJ0ScqhW8JOoW+K5 hHfcap5OcmtiL2WSnk15DGtVaQQ3QU2jw95CSzQOqezTwflGu99CEzpsJzmb ctvPX6UUNi9BTaPH/vD5vtrklBUWLhBYajblc2Yk8iV6CP5r7Rgf7X5PgR+7 VuVlU3iYeCl7n6DLQ1Wj8I3KTivpg0bsMouXTeFwu05bjR6A0Efwn/yEiXt1 uAB06EXP7qRmUz7EFz5pEAd9HowdFeCT1f3reKSkTxLVVyVmUyqi4yTOHpWj j+B59WdWW0Y+pQ2RQu62o26bErMpG3DYB83ho9PLDZbU9UF65wvUukpIbuT6 KoPEbEpZjKenlFVGJ8Er/jxcZbdW/oJ+tsOp4Kq0rcRsynsRHpbiaqPX+1Ku hH2lsl/bni5HT433T3eM49KyKWeDtBtQdHwj0L8j1H2fWknqMJPJNCeP3H7L erQ9Vywxm1L6uFYPIdvRTXCSEe9mB6S2WF6cpNUU3I4vBG/m2J7FfcEU7Xq0 ZnMxDywZr1ErLjpUifWB4BLfW0hWL62uhO85OFWrb91n7y28n1dJ9TFYpa9h UoEKzcaTsipdVHmGz8DAwMDAwMDAwMDAwMBAwP8B1IjiKA== "], {{0, 0}, {369, 78}}, ColorFunction->(Apply[RGBColor, Part[CompressedData[" 1:eJyll8dunVUURg3iAei9iDeg9/IKEHqXkBgTakioQ4YpTuIkThFCCCGEYsd2 4u70nrj3cnsv6T0Bzt76BkdLCCTWZA2WnNj3/v85ez/44ceLPrq+4W+ua2i4 oeHfOPHijcZiOO4/wOp3Gj/C8c8vh9VvNRrh+N9vgtXvNZrh+P/fAqvfbPwE x7/fz7D67cYvcPz7/wqr3238Bsd/3++w+v3GH3D892+F1W8yWuD482mF48+n DVa/zWiH48+vA1a/w9gOx5/vDlj9LqMTjj//Llj9HqMbjr+fHlj9PqMXjr+/ Plj9AaMfjr/fAVh9xfLAABx//zth9ZUrAjvh+PnYBauvWhnYBavfYuyG1RtX BXbD8fO3B1Zf3RjYA8fP515Yfc3qwF44fn73wepr1wT2wfHzvR9Wb1ob2A/H z/8BWH1dU+AAHL8fB2H19esCB+H4/TkEq29YHzgEx+/XYVi9eUPgMBy/f0dg 9Y3NgSNw/H4ehdU3bQwcheP39xisvnlT4Bgcv9/HYfUtmwPH4fj9H4TVXzIG YXU/DwZh9aFBA1b382IIVn/ZGILV/TwZgtWHhwxY3c+bYVh9kTEMq/t5NAyr jwwbsLqfVyOw+ivGCKzu59kIrD46YsDqft6NwuqvGqOwup+Ho7D62KgBq/t5 OQarv2aMwep+no7B6uNjBqzu5+04rP66MQ6r+3k8DqtPjBuwup/XE7D6G8YE rO7n+QSsPjlhwOp+3k/C6m8ak7C63weTsPrUpAGr+30xBau/ZUzB6n6fTMHq 01MGrO73zTSs/rYxDav7fTQNq89MG7C631czsPo7xgys7vfZDKw+O2PA6n7f zcLq7xqzsLrfh7Ow+tysAav7fTkHq79nzMHqfp/Owerzcwas7vftPKz+vjEP q/t9PA+rL8wbsLrf1wuw+gfGAqzu9/kCrJ5YMGB1v+8TsPpDRgJW93kgAasv NhKwus8LCVi9ZWsgAav7PJGA1U+eCCRgdZ83krD6w0YSVvd5JAmrf2IkYXWf V5KwemtLIAmr+zyThNVPnQwkYXWfd1Kw+iNGClb3eSgFq39qpGB1n5dSsPq2 1kAKVvd5KgWrnz4VSMHqPm+lYfVHjTSs7vNYGlb/zEjD6j6vpWH1tm2BNKzu 81waVj9zOpCG1X3ey8DqjxkZWN3nwQys/rmRgdV9XszA6u1tgQys7vNkBlY/ eyaQgdV93szC6o8bWVjd59EsrP6FkYXVfV7Nwuod7YEsrO7zbBZWP3c2kIXV fd7NwepPGDlY3efhHKz+pZGD1X1ezsHq2zsCOVjd5+kcrH7+XCAHq/u8nYfV nzTysLrP43lYfYmRh9V9Xs/D6j6v52H1C+cDeVjd5/kCrP6UUYDVfd4vwOpf GQVY3feBAqzeuSNQgNV9XyjA6hcvBAqwuu8TRVj9aaMIq/u+UYTVlxpFWN33 kSKs3tUZKMLqvq8UYfVLFwNFWN33mRKs/oxRgtV93ynB6suMEqzu+1AJVu/u CpRgdd+XSrD65UuBEqz+rFGG1X2fKsPqXxtlWN33rTKs3tMdKMPqvo+VYfUr lwNlWN33tQqs/pxRgdV9n6vA6t8YFVjd970KrN7bE6jA6r4PVmD1q1cCFVjd 98UqrP68UYXVfZ+swurfGlVY3ffNKqze1xuowuq+j1Zh9WtXA1VY3ffVGqz+ glGD1X2frcHq3xk1WN333Rqs3t8XqMHqvg/XYPU/rwVqsLrvy3VY3anD6r5P 12H17406rO77dh1WH+gP1GF138frsPp/8z/x//tfT6Qb9w== "], #]]& ), ColorFunctionScaling->False], Background->None, ImageSize->{369, 78}, PlotRange->{{0, 369}, {0, 78}}], {193.725, -124.80000000000001`}, ImageScaled[{0.5, 0.5}], {369, 78}]}, {}}, ContentSelectable->True, ImageSize->{399.45, 178.4}, PlotRangePadding->{6, 5}]], "Output", CellChangeTimes->{{3.433567795983915*^9, 3.433567843671119*^9}}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJzd3XecXVW1B/BLZtIJSUjovRcLUhW7ICAq0u2oSBFURBG7ICCKFbCLvfde sVcUe29YwN57r+et7z5r3TkTJ5OQ/PPem8+HO5N7ztl79fVba+99OOLEs087 5YwTzz79pBO3P+ysEx942uknPXj7Q888K76a2GA02uCS0Wj0tu1H/u7iz/xo P5f6yH/8P/z7Lu3XRPePf/yj++1vf+u/0bz22X/3n//8p92xTd7n3//+97+7 v/3tb6Pz8zv/dm/8Ht1n8N0f//jHaeP9+c9/bt+7Y8X06b/2ta91Z511Vnfm mWfG7T4nuq9//evj27dun5NtelN973vf66688srz82uUf/WrX+2uueaa7l// +tepOeqvf/3r7tnPfva0UV/3ute17/+LiMnuXe96Vzdv3rxu0aJF3SabbDK3 23TTTeOJD3zgA2M6tmqfcxsd//znP7svfOEL3V3vetfu85///OjInDUI6H73 u991P/vZz0Z75f1/+tOfute//vXdYx7zmO5b3/qWe0an5P3f//73u+OPP77N NpmfS5cu7SYmJrrTTz+9u/rqq1fOTOhxxx3XvfWtb53bvfOd74yBfv7zn4/1 teWAUAR9+ctf7m5/+9t3H/3oR0dH5cSk9rSnPa27wx3uMLpB3v/73/++e+EL X9idffbZ3be//e2m05Pzfvr71Kc+1WabzM9zzz2323777ddA6EMf+lBq2zG/ RtB3v/vd7vDDD+9e9KIXHdi+XjhWLrp+9KMftVso+qc//Wn3hz/8wbXR0lVo P+KII/K7xe3ZX/7yl+0ZzxqDDZrO2H//+9+7X/3qV90PfvCD0Rx3xDg//OEP m7KMvXHK4Mc//nF30UUXNYM0Rjw7um+S/va3v727/vWvv3b87jBQAuXc7W53 6+5zn/u0AYtpVsRsn/Oc52BG5OmWLVvW3ete9+re9KY3sY5ZmF7Wnn30ox/d 7bTTTt3ChQu7m970pmyi2bixv/nNb3bnn39+t/fee7exlyxZ0t3mNrfpXvrS l2J+RVLOkPF03nnnsYF1Ynf7Abts5YlPfGJ3netch05GN81rHOkTn/hEt2DB gu66170uJx8tn4G9ZXl/xJk2zc1vfvPu/ve/v3Hjzv5a6LF74AMf2M2ZM6c7 7LDDsDraJK/ROfEdeeSR3Ytf/GK6H52+Lkxtt4rHhwN1W2yxRfeqV71qdPN8 5HOf+1y3cuXKNuh3vvOdmfjpv9uw8fO2t72tua6A8Ne//rWFY8N/+tOf7m59 61t3K1as6M455xzfbZqz/+IXv2iuxoU/85nPNAs6bX3YoQcmebOb3azNGayN 2YlgJv7Nxk7vKotaXHv1q1/d2Hnf+97XXBCLL3vZy7rFixd3BxxwQPehD31o 8xz5wx/+cDO+BzzgAYZbJw62HSiERz/rWc9qRn355ZePbjEzB4uSA/r7zW9+ g/G4yedEd9VVV3UXXnhhd+ihh0aO6L8TJAzd89j7Ke4vu+yy7k53ulNTFL86 +eSTG2X3uMc9SGK0eZImoAgg++23X/fa1752kBjXkVfaEol22GGH7iY3uUlL H2vg9Sc/+Ul38cUXizdxk8+JFlMOPPDAbttttw1e++8uuOCClmn7EOpzvszV nFe8esYzntHtuOOOLb684AUvyMQy2bLijW984+6EE04QoNeJvW0GqmQ55hQY IieMbjk7e2yXX4dT9FkovnvIQx5Cjd0uu+wS7PXfsY54LL3P57yWG+S3zTff vNtggw3aM5Hr5ufQQpeoKpxMts+JFmrYdW8LPdXSkin33HPPFk1p+dT1EQMt y18777xzt9lmmzHD1YlhJp/cJNkzDF8MMTTHECJDwvNzJBlc4t1oo42AjTnE ECOZULCMWJY23+sl4ly39dZbdyeeeGIjaD04nGwOKFyANne+851vde2Z+8tf /sI8uvnz56O+e+QjH8nd5qevCvsf+9jHGs9Pf/rTuze84Q0tDfBCpuDv4KoY 5Mi3ve1tu3333bdFZWH1lHVhcOuBxMREg7occXB0rbmc34KLoDJ37txmDcce e2wLPrLf/DRCskT/S17ykoZYEqK28fwbL8Unb3n4wx/eYCWJyYbrxafJuJFL DCkevDZ8rkweqOsd73iHgNSCjGAjKx5zzDFNdZEdk9+5LSBfccUVglDxO1oy IAdq3XDDDZsRB8ReP/5EI4hMeABhCPNa87eoyV3aLoQCURfCljIF0ZiqWGTc CBK3RRezLslZ3/3ud3c3uMENuutd73pqmHXibli9AB5ExRUjL2Q6mT/GOUQK 1wAjyemMQG1FhlbPAXZywsYbb9w94hGPqGKkoZenPvWpLXdCLzF+RSNZ7ZJL LmmhPxynFErGqGOwBCjhnbw+LGMH2WAar4o5xyxTNk2pBd785jc3GL0qy7T2 ile8AgZN8DOvsae+eO5zn9umpJqA3gVP5RD5AR9f/OIXW8k8b6ABhk8iNEvT G+Y1ueaMM85oaF0GDumtH++8TKQwoLg5hZTmt2t0oEwGLqmOmeF9l3bPinaP +gv0FYk4/I1udKPujW98Y5NL1ZigFVmIy6z7k5/8ZDNpP0xK1JkYONFd7nKX 7kEPepBrxTtxPu5xj2ukHnzwwYZYP97pFm8GpNfwwHTkXi6UZEKCjr+T5x4E oD+uxb8inYTq1Ebvf//7my3hW5IN45zTLm/YLFzg9iNAKbkrKBMRlqPUbjqP x4tll8R2mlERRapKC+upMFRINtsLqxPD4Xk1ezPVu9lwZuksGMztbmL3mxlu NIi3ilWSE1flnoOTLtqU5mj0gIzlKkixOWhq1xSgPGaPPfZo4MHzFcyACujL f095ylOK3Y9//OONrUhp4wKwGDKk38L+jWcVxSZJEAbAKf0WuvjSl77UkC6B Up+Ysme6gftcl2ConR6h+Kkk+l8y7CPYguYZ3/jGN5Q5zSLgK4gT/wMNiLe3 vOUtmzi32267hclyhdu+JxEkBDEIzqfHoVo6UCv390y0e4RVddtJSUlJXN8C IAXtFWDiUmCmUqzuwO67796QtvhEossHGjjttNOaBpbmsFxY4sK/5CwM+Df5 xLMlHopSOCie4u+VOSK9kfgU8ZPtaQwZORhcslbqVNC8/OUvby7EHgmSbfFe JLGvPfJe6oVzDjnkkFaI3f3ud28dBJBQpMfywtnUurDBD42tO97xjt0HP/jB Vmww4yc/+cmNfoIf4mS8bLXVVg0x7L///jl8n08FRlkMlDMVUiJ4hjyMsqDV oPIprVSgAoFAu0haJ+dQJCBX3vCGN2y3MqfnPe95vHfstSyQzg1z73vfu12b Sbn9/T1gZrmqcghMM1Nt9LCHPayJIKVVwsKJy+G09F84DlB4whOe0ApeDJqd nkDuanfOruXNkn48yi0Q00knncR0xvTzNPRorC5pn715VXO2igAVBuCg25B9 6Qbm5H7W4plQy4LZDGBR96hHPaopk8IrjIuTin6SB2SCr6oG6ZlHTE5OFtzL erpPQqoSZcZ973tftjJuYgElko9nXI8xxr2CSjQSdDBfIdFQt7rVrZqUmWfg lll1LE6YVt7WCiRHP6bm0DQS3/Xy6GGNryjTNMQY31Vn7ZWvfGWzPgTDsUx4 6ay6rbbIgx/84KYXPScuJLguS7rhLjYXeCB063Oi2WUV3aUVNq1UxVPBS2KK arJhFwsJ8e9ZdTvZRka/0eLuOfm1dAumCysxUKnW+CI6eYPLxx9/fIK1vuf8 kY98pMFeFWS2jRptGvuKcI2wcJjRo3O8qsCNd7vb3a77yle+Mov+pip2hUj4 5ETeSjSm0LMq3VHpmWee2ZYwiAMmCQsq3cnN1WRmCQr6tdMdrKdKY0MSJI8u 3T3/+c9vCxiBH0N3Pnt9JvYbldTvec97towvS4R8SuqGJqJLL70UT7Oqbm7T 0fLly7tdd9218seoRpLBEaj6Hht0j/N4a+kvwnhGo0VNXxIim+RX7iVYMcT9 eI5Ct1QHwx111FFjfBqqmVV1EtNBBx3UPDvwfj/rgmYd+BPJWTOjA/kNEddW DginTfmAgAPFrKW6lDLV8Rf6wgYS60y2TMM+qDGct4QHvxEeFkLIq1FDfV1W BDANdIkxXUaXAuQNNQA5+1qAjDlSRD28YEHEpEwCxVQiJo3AuFFKl6OpjKyH LG6f/eJNreUVoGIftGRF7tD22fd4Bf2p++Y2REimSKL8CAmjai4KkMY2y2TO xUrF6fCo8LP+O9mwks/c/GSDRx99dHN8MHeqslvQYjlwxtoFrPe+971Nv2BJ pJbyUigaVGDixD+FytekduKXEdQSiJUZ9kxLF8FkeBVXVW1+eIY2oQyxZr0j VVbUnQrfmTOYWDODNCVHJJfHUj6yS/lh6GnF85tWwFQK54YIF7fCV0aleaGe UepXL86utQUE0h1q3gKw0jX+C837nGgQgkaGmnef0hM9uuGhhbHm5TpjD7vm 6ncZ4Ra3uMW4a/6a17ymWdpQ84QKT6MfGFLvL07hk8973vOeNp/ekYJP6AFo WQRfHDi+iB/TNXHBUADD7BawRT75pCc9qU1A01AZpFvLh1xFwtCsW9I++75G mfCCFNPjH//4Zr6yJCIRzyaFIdRgbFYzWdy8C/EsXWApY6v+qGt0PdXQnFqy KztRVZedeE5uO+WUU9rqm+S9zz77DATTeyN4d+GFFwYFPieaOPRX3HGPnEcN COaaB1SO62P9f/azn23PeHoyxzBnIHpeHfrvv7v88svH6W1e6lj6dx/bMseQ Z7aC/hDd3LxdkCE1vicQpZm2wMMEn/nMZxLVypyDonQOUH3qqac2d15j7bRF PkxrMi8NYi/RU0nOah0oFFIJw/A5YXm/FifHhoELaUoGZKISuRJLjxDwAgpm NYxeyVH+Nq61COH7OS1cTDbDUNlRbtQBQ7sQfcsu7ONYNrhmZrUs5mLI0bl5 zfdDPQdFs9S5i5qE+a17tYwyY7Q5+C4sFbY8mTKtkCd5EIBHxhFxstWLQm2u 1q7Ip4hQ7IRD6Atj8d3satw8tcBRPSx7W7QFOu1nYUpjW5hodZugcr/73S/U 6bPfxDJ9AWt++zeqBUm6IECmZSlzNdXxwpQVaylsoUBi79YDUOP7bbbZhn8M Vcg9SoVh36nCec3ICKs2Pqi7gqFzZojWCkPReiYtLkvW9BrshyAnxKvJeKrI rPiIBDaazLGVc9oEYI96u0cLPeOgsg0HOAwFrUiZCWaQClmJ0LImkKcTtuFa KZHvgWasnOMRWTUZQhbXSdpkInQJxNpLfGtWpSxsSUluhuxhWUCMPuEIy7cB /EZDP2TUpZGgddzwYdm6OExMr5yZsXLrbiHKwsJiPMINYddJmNiMmlmew0o/ 0ogyj4jYiz6CqBZanRj4FeZJVyoVB8NSq152WdODgcN5QkdVRESLqvC5kJLP fuW21vLnJiWaY7vttltDCyh57GMf2xwgCjGCWbwWelzYropksm1m2kaYdMFK 1MV75LS8Xijacsstm5ZYy2raiqVMdiyqEXwtDDMXDSCJGd+lTNchiCpwwjtT mf0yoXUZSlY3ybCyF6VCXEh5VI6DbNjVOIxc930mjW6czkAX6JQrIBr08Ves l75EUf7IteUN68tzk3VNgAhPTWxhWOV6YiJXVxmFhOfow8V4hK2FYONWJVHc +U46UV7q22oR8XaxnPYHPcrZek+LGztWSYkJIBUC/JhEmHCN0XLF3VMEfAid NvYJM2KqDA53rSYVFtbQXxFGbI4cLmRqmJUmo0pIfNOTJ9IDB4V/GJYAqikz KGFHj8zxLPtoMRhPVaz/OpNGtxo4nzFonxjA/Ri/IK7LsASRwytTxXevCpCc pPSFGJeYLxspOflI3D+Z3LCbWn+UQWRfPzTIJtiTIGSKinuztw7JQchi6SBs xJtR1SdijyxiJqiFenZLFXJMsqmWOy+w8CENrkGFeCVaK2pLBjIyFMao6fDD Dx93Z0VduxM1iiXxiy66aEwiRwEx7EZRByPxEXmNqqjEXFbWYPfVq7GPN7KA +KlSFwqrCyOZyFp0FCKbyDghRpKv/8TRanmSVq2EVN4UP1gvSKnRM+zBqGn1 ye24CJbWTnVgtyyFP8uiUyuFc1tgY+roBfFNtmuqTjllsZznsT9e6DvtjjWo rtoRESnGpZVhLF5KohGAlgycgkRrHZFjhIbmDC7D9lIt++MYwcDDk37ISzTz pIyswVu6cys8bi/tloPh8AxYSwB2Al5zzTVtV1buOvIVHC7WTuW/3gPEdzpJ iD9ecSYeEVH1JFJcccUV4+dcU4UzLQHNtosFsyqub1UsHe9s5W8gBtGY3XdG 9cOwlZ1mpUgKVreVEnk7bgQIUMoGN8qepcysUp5h1so2gxzuCSGMkjvSCKPP 7MvGuwFMqV2BbIs9fImN+XEPYQOAjMJ6hrmCrIflFGRj91dWqy3yXj+v8Sfu rIFo6BNOOCG8LJwyuOEf1oBE7eCm3EoIlApF0QxG6XLzm2lqY1nJso4hP1FV JIRwPXzPa5SDNrqKpsDdoOk9myb7JV2aCxeKASPVjVp4atUznOxH6IDg5HPJ kD2FnZUilaRgt7nR53KutaxJkapJuBAwtlwp9i4ZSLJ6ihC84LB4QPWVV15p +0hSPb8hOfaoLe9HSubtYeF5z0TbbmKpySgPTdJELDvdiDdI3yunkITsagh4 MWogI9s/bK9/sP9UAmEBno9oMZkeYoOBdTw+S2aok1W5PspqT7HMVCDjLW95 S8tes2vudnm1VhgZ/Hhv22riXr+mhHOYMIynhKy9CowTsqijwdmfE+jXOHgC e/KYaG1JqvZ2UJAho+AOAccd2UjBOvuczBRTByPAPp5kLMoTxoyvGbRz0oNM XWFZElS3ke7QvIZhCpT1qgTYYlZRrcUGGZi6NsjoBR+VoYJTgjAWSFmfxqxy tcAhXG6xQLu1npeCwKSoJ4p3+UMrQIZCX1W7ih94N1DBmHdyxZ97hAlnNg5K a6Rk+vCvSrnMxWIE+B2szC6KOkkjwEA6wGoMeODsEqo4rz2jJAyHGB3drvX7 TYUmbQTrhf4eLwn0cdiuhFqvD1w4LlxIRcTSJpT/475MC/0pDA0Uz5GSHRIH pTQVeRQw3F4jeCnoBF3pZe3EIASzL54boWw1Yqi78QKHwXaBoY4ZGCvcql9g IJm62CdfbJCa2go0GUY9tqKAUIMy5F3ymlpUW5ADGFNNVezLoIxGVuVQ5nCP CkmokjnXjn13S5xqinDhWdnvI7GckPiopZ1eBv1GArBVZ4jxsi9ykVItXIqs UQbUlgmxQeInFk0Evr9Lzi5Ya8KxiAid/aajBS3kUK5tWKAYjqjO0migmS3X imHuDJEZLES6BoYxwEk4NW+1Fe/YVAIvxJyyghiNKLDD/CCMaFjcyoKiGaBA 4THmtEjHc/URsOLardu1fieQZAs6UZNSJTt5o9n5XYt9ypxPvxpzEUJGx6Wc GbnHbPVjlSKfThR1hZaHm17QwuisNjgcUDzBVuoP5S3jp91DkiL8yqkKi5hq 63Vho3bMM0h7hlmRbMCBjsuwpFQWdwRsTiT2xCS17U4IURmLaB4VMnbK2eRm WJjdRqpf7039/FStpJfH/4/Pa9RLyvxJWLEcX4cUJEXxgslQP1GHwe4wGFYI wBULlD/Wic5t82tYRE9N8o8QfVxaILcGT5DOUBTKSOP6i9JgyhKgbJWtNo9k IITFfXUsDEbhu9AwpVHWepDcy0CTX0lDjIDhsUmTyQQLomWhZGRw1XdMvDBv 4wNqDQG99rqArfY4RcLZbmDr0Ar/NFMA4XU70LPdQNoitsIjYNUxGVPMwzCV kbUmS4pKS8AxUNKoGoQiqITDxwTdOuwHRGYAHpWkxB+xSX5Uz6zzgaTtBqKX gG1YBfcFpKMzE4gfQLCkzjx1BwunOy2k4A319OB28fjkm4akYolp+DGF74Ud AFoMGYYvCBMuwDtssE6HybbPrwVzk1i3CHM/OmkjY7LnhPooifnHJYUErWLj 1FVoG0oRLmQKr5EjxmDAc+gDFhVWTp1ulSQYW3PXXpC4Z50OKw5PYNI52aBb YDlyQJ9dgmK2Uifyy6h6756xJK2ClJOyAVcFqp6CsA6SIdPBzuohVBfA97nH f7wobUp9IeqifbF2nfgbRj+erHbPbUOj27dry8apTd6WO+VVMMHfgJKc4Xn5 AXKmI2CG+oOJXlVLxyIkKnHYMAKbmF29SkMpXJiAAIkHG+eH2zEUxVwXwmKm 63SydsccsI43SRa5Ez3LwP7YuXa3PDwh0IY+gTV9WPmYC4LmAZ3mtMvzW81p uPopr45bRv09i9qKOoheO9uZqaUqbTfR3PWgdrMBy6xGlOG54aWj+10rlrXx w+/n5mrcYKvKaLipn4vZQyNcyiQX5DVqlRitABCRa3WAnfspPo08mZ/Mp3S3 eqI0gPbZZ5+5UFgMZNfTqq8mYAwKQb3xSPjnJz2iIRnzDS5PZXXyX4lDSAad zE9JVA6eiZ6J1hiUpwIRhHp99t8VMfWaBsRovsN19ZIE3osf2ov76yUJKneJ fzikTkC1R1Z5UwMZGgkUntc+J+roQLtj27yvDoX413n5HSrzMP74RRH+zTyH 49FhPZvTN9v7X/TejP+Lf1/dfvUuRbvavG+cfkk5LDlZpctLn2y/+nrHk2Ko IKgZCNjydkY+3GlkGLWKYfrm+vIx7ueAnhUYGVK9ZcRvqE0AEmzqHnVCnU25 eDpJLkEXwGb17cVoa0eGAAhnI8mw0J+SDIyyrblWiJBESFp9PI5PVt3rvDg2 pO8k6VPtVw/p6/Y8bDru9g9JqA0sIrL0Aknk3uVxShW9pR3tASHEdJcOp5ts FKrPVDER6DeeYaZ+9sVNnoQzOAPRjIMQzI5igow6dvS06bPIcWuYpV4IALHp lEsZ3JoYtTlkYPWD0Bzpf/YJqrvLHHBok9Bk+5xo2wGwYcltXn4n3Vf8qaJE pNOEpgRmKwdD2/QXtvH0tSKAcQjQw+2Acq8gfcwxx4y3A+KwunJ1iB6eQKi2 sra2EoNG4fK1mx3p0mu/aOBzolUfQLd55+V3UuCqswufqgD3qo7ZkQrGEkUd dCdatZxtxZPts98jWKH3/GttBbWDleLgeQ6DYuuVC/J+Dgbn9y/XiXoxPu2L 0DC1NnbBOs46rzX89Yzw65qOfPBSRweISD9DC9bjGoYAIA+XndPrrz3XK9Pk yE321ykid1q3rJb7AepkkWAJHlHHwHIboM0e/jqIvo7imEyrSBuK8aho9TmU VbmuO95xaV3H0rw2VlIyuvBaz7wiDQ8f+KJMIQtIIpJ66YeColaSiwIJwYYv qiveHzsM8vMabjUMsKS4WCMlfVtZ2VNw0mJGlWAqU2U8zdfiVK2zoNgClj5m TLsKITr9mnxWCmSMVQgB7rCiNzA8Iazw0YNSNNlBFUIfrhzrNihjoLXcjDs8 JGy9z6i1nHvBdJo4Glit0DOJjYhBU99DX9EEIUfbk6QzwzL1IbP72cyOEaAd Mtd1tdgx2NjeViOmslK/A3yvvfaqLtAq7upJXTc93VD8rikcyUfP1urbHN2V cFpsmUya8GN5CfBU8+i/iNi+qyM1rLRSih+/ZQNFrBaVirOUIcHZqiNlntFT eEo6CFsiF3xlpGtDmaZwzAk5DCBauybrnCs2TMccXFdzyZ8irOw3JIOsbCwE a5OMfllx+bhMBQDUyzgFHorDOnQNDillh/fwGtPsnJ5Hha4RvGjGfnma71UH bN1zyuQdcn7qpSygw9jmuDqP7g7hGPqMXfeQap2erYOnxtbmFUuYFwSpKafp Fm67YiANWIc0Ns3vBF15iNRrj2ydy9bhMa0OWKvM+70OlCQY2Ms0J3f5ucUw 9L33dEEbjfnZYqATCAQoHIUc19gBBvUHXOenQpgtQBZanKDsG8NLmtAUldyf LznWw/XU2IKHBgdzhBmQBIBun3QQqGfdw1bkQGaeqK8p3nXtNGPInZbCLQhr nUHrdRIU3hAZ9PP4rWfwyChmknblZfNwO7Ko4Nfb6rIm1npjF6yiByRYihMM gkYoF021a9T4NgSEoaTU986Z6n0JdnEKcagvXMs87W7RlUc5CW+YVJst32AS 8u7hHRupwxbzMs7QGiSBigj/cxps7eGlmEDIoipQ4Tudl2qisBOr8bKAa2ym SAPWMDjY+zZulDIWz+ax601mkbOYQL8WjHNTaSboyca9BK0VGjoY7sGpXdXA sl2B+/Yy3SfdnEtLomKrJJSYrz0K7orwHo0QMKq94IiVIu24XtI++1c+CQtT 8lzSxtbaoxMWV0FWqBE+CAumEV7y3XzpFYvHW66FgOS2PSuiKpXsQpO57Qw4 fwaprUyrZp1aStI2BWiZFSwZggxr/iSIlozO4+O3+Ad9HNPYrxdfL8VFjRXB W3QSWpDJtJzMZvCWNkK/JTlcAU2OfSxpn/3rGDUEpyTX93Q1kdRP+vVMpsAG nwRAFKcgpii3S5KjtuX/JCeo4t798gUJICv8ajaBIdFqCicAPDFfe8XrFTS1 OKPNCYrUu11qO74wFqSkrPbLiG5kBLNII8KqKieMBIbePW0ZhNY8U/FFUVVs m8MeCWEU9ADr65SWa/mSxnB6nxMNZNZZk8VJgVwjjMKwMC3oBOFYUYjCqI6u 1PseDbQ4hxvsTRlXHXxFSyAi6eiw9tlrOFHQuCbyXY03meMpGOzKlGrm5XfD lle9Io5YkUt7VAGEWKEBZOqoatUKnuWXMklG91Ytivqh5v17deyftLEVDU05 SQ1pBvGZR+2eQpeEqYuX6uGztdohUsdWeXH8XbqQ3xGhmNwwS1mmPP2FafMb gtca4MywNnuzIyuwwfj0lU1ryBue0lNtsZ6hIsyp3RhQMxThc6KRm3vMxorA Zp36q9JbSSlFB0Yfl96DbSBjRbAnZolOmEQgEphE9nqpab1qy7PWJGyg4R+a FMSrc7qKFujMMJKLmMbj5S8eX2IwFb1bBljSPvsXvOaSwbjwYXfO3iPJGwUc 4i5tyX9sDsqOGFfaMoSdt5FOQzM++41ovGOoLRwJ2RoVahQeLg/Kc0WmrYPS slEW51g8LHc1je6e92njAAXumMz7VFa80MLRvPwuIVzKfypUS9EkShvuqbQi Kghw1dhbORAxV3UqVvQQoXi9xb/UxgEpwdpXRFmC4GWXXTZu9jEAezICgIYW fE607CRrThnj4pYdGRqUImAN38ZpUUcOZ8QC6m45b23BlzLs8wjMliXr9CBd b5NlDcTFGC2aAOHDM9e0LOKE+Ao+mF4Qq56jE28HDLmf14QpT8t1eh4UHC40 PDPLIz25JN98zHsr7dc5LqAOFlL1Eg7mmLDGlx0L+da/3XJWnSaxHi92psHH IbMLZmB/eT6Cl0qPlhzz1XTDg8MmFzenG0O/sQwgq1eiMIrwjWle2R+Kg71h cV7pKJKek9wUtO05iAm6E1bdwWx7PKABXgJKqyZ2zSEhA+r3uk3/iTtSv0Ad 1jlTSl6e/MhHfLqONdB27RCF2iF2Qg+q69QSmRvZjuzJ3Jedq1opiKmUQUSg VJ1alz7yoOh+Q7EsbKNyAYQo6/XJ8aRUTfCYR8MXtGAk5hO2tiVCISRRJExz 1xySJXNcQ0I1xCgrCeeFcs+bQTZ10q9Oj9AlsebhnFpdrsPr4k0eOx47ilYe BCOW7bTTThOCQvZiJY6pje59FtabYKM6gSA++OLxRIwFrPcbRCplGLNBnNxW L12ixFwRaQzecfAMmRqW2WWEbePALQNYN3rMaoXSn0zRKESuUCSZVJFT74zJ BmId1aAe8uekPBAZfsPv9EAacwdkMmd6YjQCscrPd9QLRycszDojpTJ3/GoM xajHbM4st5UnFJJSF1iFy9rVBOqI7Taw2Mgy3ABBEZAETuvV2ufOIJ2lA5NR dNKKYCyVDjpkdO0QkQCT60hZyi1ooEX6qsUcP+7PN88MjzSpDOz8EGuteg9x rJghpBA6/Wedu++AQHZcO1Xz3TbDDp9Iw4KqSA9Cjx1c5mj1Nhh5JegdPm0R QMaFxtlB0HROCgx+oASl4tIkGbxikvIC7AmtLpAtwpxkS92CfHtmhSWYQhGf Lcb5g8lZSZ3N5fXhBEPaWIhWgsogj0b075FPU1rWlEIB+nyikP4kKO3u2kxj FF6j7ycsiLTiX3x/dA7DpKgeZJFLDSMbeCw7A82qtQhrhcS9zDQEdnEKBxUq GKhIIIvMF8Jx4mGy+atr4l88UkcUgSKGmY3ruUlQ7blHQL2y29PW1snRj3v8 bUQNH/fx6GBmr6GUJtrggiksNafhqsXNg8QU7QI/bJQ3giNodZ1529pzdI5D sXhU/NY4ggJ8Len4ATU0EXxvHECzDpJcnOOweVg+vDtQpc+JBl1qkakiNv+w TxKHta4/mWNUoe59KROJFdmwWtRPbWWJ7+bkHiSXBQUD9MdoRt/I+WpDD1+K gD+R0U2SowJ618/zdzZH2yRUJLGpk+a4I1OrwC4BDk8mMFPXEGUsv4W2aise kaTUnldeH0wePJ1QDg4VKidCdRNpL74mbaeFAFaVjSZynazR5rT5C5DlbNZP uK/QHu5ddEpK3FXGtaFawGJ6hrYMoNC+TSpg8H8+GKdRz4PjDCsiQNJ+o8FV 0qK2kNzzZ2ap3n9AxiKKI12ARnpIxU5gAW7CDVfVn6u3VLEj/idDAEuZtA5P yikZHcMK1ug2FAoO/025RMtgRb5wktVQXiGfVAQ8sM/Kh+KiyJZH1e6ARu3Q q+36XEuA0tWjPOng8NSMrAJgi1KskjWZQxzWb8PLKkQr2BGNukBez0ruYRcA zSGU9g60XEuuFUXp0PU6xtqb2KLxgTMpGXW06bt681tuXhtuwJf8mRD12cJ4 23atP3KIGXFA80Iv396IRKWr8GESCpDbIwgXH5r/QqKXzbfD28GHgZg/G/Q3 P8q3po19kd2yZ+YvN7IGhMJyCGLjxYR8rhpXXOinWujrmZjfxhEs1Sl6JTYD 6JKSWb+/fvSl5MB00gae2/a4mCdGIAb6Vl/ZVWt2eqb3epWAAC8FKitQXKIC F0JU9QZAvl07wLnX1Ma93rQFkwQwh0ynTQiS5MXoyBzlDyId8J+mvTA9ix50 aAhewgZLwvNq0ztMhARWma8wGj2zn+/LKX6gUCCm62rkrjLZWPzqJ2EIds3/ mcV4/7qSjCNRDVPg4xEynjF9NrZjzc+OnXzj9X/NVq8rYdCihuwtn5q9jknW tgdZSuzkKJYYyiXJOzddfCU1J9PQMPuDHWaU6IKWf4U45qf0tZaO13znSu3G Fg2tKmFXdMY+jD48xp6bXr6atEoyDFTsN8uQ53rjDgOnDTwxKz4meJNoha6q 5fFZG00oF2iqhSIoMgulS6dTUbWuKlJqgeWmU7GohRqorfrMcIbRSReqlX7H b+5b1FaNJDKQOk/9NzNChNIqF/uTkK+lweh6EYdoKdYlEdWqUZvqBAm+Ul2t CnsMTWquXC2e9pZtHuAZOAUYrkCPbxBRBIYsBPpLeoq+nreIPAxOhi2RVP9O Mw4vjJHIXKtd1CzLYoJcBejmxoZCTWxWw4kTQl/9nsf+GlNht0yoQk3uMCuE UUtUbhsntf5IM70QspxXrxVkJWSE33zRcOOZLKUGnYfxuzeWtnjGvEXceq+F UFDrLHiROgwBzKi8KhM9dTqhbtXWUXpeddVVczMvCHMwVLhq35QKiUlASnxK 5YvidBTec/JtW/BqoupR7c8mWMszbNB/ergMr45va8Rw7XybV6s4ZVyFbhL6 zaRIAM4liNE7+kttiP9FWzfX5+/RBv8DRDAT4g==\ \>"]], Cell[CellGroupData[{ Cell["\<\ Derive equations for linear stability analysis about a uniform equilibrium \ solution (10 pts)\ \>", "Item2", CellChangeTimes->{{3.433567896430321*^9, 3.433567897316579*^9}, { 3.433567939123716*^9, 3.433568021148382*^9}}], Cell[TextData[{ "What is the condition on the parameter ", StyleBox["b", FontSlant->"Italic"], " for the system to undergo Hopf bifurcation? (10pts)" }], "Item2", CellChangeTimes->{{3.433568074212306*^9, 3.433568114113836*^9}}], Cell[TextData[{ "What is the condition on the parameter ", StyleBox["b", FontSlant->"Italic"], " for the system to undergo Turing instability? (10pts)" }], "Item2", CellChangeTimes->{{3.433568120592667*^9, 3.433568130332626*^9}, { 3.433568233579821*^9, 3.433568233620543*^9}}], Cell["\<\ Choose a parameter setting such that Turing instability could be observed \ (5pts)\ \>", "Item2", CellChangeTimes->{{3.433568234481335*^9, 3.433568263304813*^9}}], Cell["\<\ Simulate the reaction-diffusion system on a square domain with periodic \ boundary conditions and try to obtain aTuring pattern (15pts)\ \>", "Item2", CellChangeTimes->{{3.433568273517758*^9, 3.433568347083156*^9}, { 3.433568391411892*^9, 3.433568407894048*^9}}] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{973, 839}, WindowMargins->{{Automatic, 139}, {Automatic, 14}}, FrontEndVersion->"6.0 for Linux x86 (64-bit) (April 20, 2007)", StyleDefinitions->Notebook[{ Cell[ StyleData[ StyleDefinitions -> FrontEnd`FileName[{"Article"}, "JournalArticle.nb", CharacterEncoding -> "iso8859-1"]]], Cell[ StyleData["Item1"], FontSize -> 18], Cell[ StyleData["Item2"], FontSize -> 18]}, Visible -> False, FrontEndVersion -> "6.0 for Linux x86 (64-bit) (April 20, 2007)", StyleDefinitions -> "PrivateStylesheetFormatting.nb"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 95, 1, 80, "Title"], Cell[688, 26, 139, 2, 32, "Subtitle"], Cell[CellGroupData[{ Cell[852, 32, 233, 7, 20, "Input", CellOpen->False], Cell[1088, 41, 129, 1, 26, "Item1"] }, Open ]], Cell[CellGroupData[{ Cell[1254, 47, 457, 12, 20, "Input", CellOpen->False], Cell[1714, 61, 24281, 411, 193, 10147, 178, "CachedBoxData", "BoxData", \ "Output"], Cell[CellGroupData[{ Cell[26020, 476, 234, 5, 26, "Item2"], Cell[26257, 483, 233, 6, 26, "Item2"], Cell[26493, 491, 284, 7, 26, "Item2"], Cell[26780, 500, 173, 4, 26, "Item2"], Cell[26956, 506, 275, 5, 50, "Item2"] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)