Skip to content

Instantly share code, notes, and snippets.

@ketch
Last active August 29, 2015 14:07
Show Gist options
  • Save ketch/707ec2208c8130deac1d to your computer and use it in GitHub Desktop.
Save ketch/707ec2208c8130deac1d to your computer and use it in GitHub Desktop.
Reproduction of IMEX RK stability regions from Ascher-Ruuth-Wetton paper.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "",
"signature": "sha256:92812c5315577f461d522827e28a855b9cba55cee0a072012cf59e4c740ab771"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"%matplotlib inline\n",
"import numpy as np\n",
"import sympy\n",
"one = sympy.Rational(1)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from nodepy import lm"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"CNAB = lm.arw2(one/2,0)\n",
"MCNAB = lm.arw2(one/2,one/8)\n",
"CNLF = lm.arw2(0,one)\n",
"SBDF = lm.arw2(one,0)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print CNAB.order()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"2\n"
]
}
],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print CNAB.stiff_damping_factor()\n",
"print MCNAB.stiff_damping_factor()\n",
"print CNLF.stiff_damping_factor()\n",
"print SBDF.stiff_damping_factor()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"0.999999996\n",
"0.33333333037\n",
"0.99999999\n",
"2.2360679775e-08\n"
]
}
],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"CNAB.plot_imex_stability_region(N=200);\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmYJEWZ/z/vXMwww/Q93T3Tc8HAcCiiAq6oMCu4IIiu\n7uKxrnif67Eu64HKT0QFBDkUVFYRFxcVYfFAd0UQbEF0wQO5b+QaYIa576O74/fHm9EZFZ1Z3T1d\nVV1Z/X6eJ56MjMzKyqrK+mbkG+/7hjjnMAzDMIrDpPE+AcMwDGN0mHAbhmEUDBNuwzCMgmHCbRiG\nUTBMuA3DMAqGCbdhGEbBMOE2ShCRXhF5Z862BSKyUUQk3ldE3iwiv6zlue4KInKyiHyrSsf+gog8\nKyJPVeP4442I7C8if6iD8zhQRG4e7/MYT0y4x4CIPCoiL0/qbxORARE5N9rnNUn7d5L1Rcn6xqic\nkGy/SkS+GR3jxyLy1ZxzaBaRS0TkaRHZICL3i8gngu0DIrLnKD6WS8rQDc497pzbw6XO/4P7Oue+\n55w7egzvO4iILAu+ow0i8oCIvGdXjhXjnDvDOffuShwrREQWAP8G7Oucm1vB44qIPCIid2ds6xWR\nrcn3tEpEfioiPcm2/xCRrwf7ThWRzTlth4rIPsnrV4rIahG5RkT2id7y88DZwes/KCJ/FJFt/vrO\nOMeTReSLIvIiEbkuOfZKEblCRLrKfO7W5LrflPzP3uS3OefuANaJyKtG8BU2JCbcY8NF9YeBE0Rk\nctD+VuABhophUyKCvlyZtH8AeJ2ILAMQkTcABwGfzDmH84DdUcGYDbwaeCjaR0b3sSrGWN53efK9\nzAY+AnxdRA6o0HlVgwXAaufc6tG+UESmlNl8OLAb0CEiB0fbHPAvzrk9gL2A6YDvOPwmea3nYOAx\n4GVRmwP+BDQBPwH2ATqBW4GfBufYDSxL9vEsR8X8kjLnfyzwP0ALcBGwMCkbgUyxT/gasA2YA7wZ\n+IaI7B9s/x7w3jKvb2hMuCuHAM8AdwJHg/YagBcDVzNCEXPOrQBOAr6V9OK+ArzXObcl5yUHAz9w\nzq1PXn+/c+6q5P1vTPa53ffqkx76z5NezxoR+ZmIzIuOuUREbhGR9SLyExFpSY7nnxaGXDfJE8dN\nGe+7QUReLyJ3hj2kpLe3SkSeN4Lv5BfAamC/5LUiIp8UkYeSY/zQn2Oy/UQReSzZ9hkpfTI6VUT+\nK9j31SJyt4isFZFfi8i+wbZHReQkEbldRNaJyOUislvGZz8KuBaYm3zPl4zw2B8XkTuAjVnfacJb\ngatQEX1rme9ofbKPv7ndBOyXXIMALwUuB2aKSFvS9jLgd865fufcH5xz33HOrXPO9QHnA0uD7/UV\nwJ+cczuC9/yxc+6n6G8zhOS1+wC/d85d45y7yjm3yTm3FRXml+S8bibwOuAU59wW59zNyWd7S7Db\nb4AjRWRq3nfSyJhwV57/Ak5M6m9EL7jtGfvlCrlz7lK09/5n4BfOuWvLvN//AV9MhHPv6Di+x3Vg\n0KufBHwb7SEuALYCF0bndSLwdqAb6AMyzTRlzj9839nOuSuA7wL/HOx2LNqrvr3csURkkoi8Gu0R\n3pY0fxh9sjg8Oce1qBCQ9Mq+Brwp2dYEhKaLwSefxBTw/eR47cD/Aj8LesAOOAG9ES8GDgTelvF5\nfwW8Engq+Z7fMYJjg14frwSanXMDGZ99d+AfgB8CVwBvzBAqP97QhordLck5PUFpD/twVMx/F7Xd\nSDaHA08759Ym688F7s/ZN+9aPhr4VWBai49/V87r9gH6nHPhk+PtpDclnHPLgZ3A0pxjNDQm3JXn\nx8AyEZmN9hAuzdlvVdIT8yW+AH8LtAKXDfN+H0IfGz8I3C0iD4rIMXk7O+fWJD2lbc65TcDpwBHh\nLsB3nXP3JL38U4DXi8hYzS3fA44TkVnJ+lvQm1wec0VkLbAF/U7f4px7ONn2XuAzzrmnnHM7gc8B\n/5iYqP4RuNo597tk2/+j1EwVfo43AD93zl3vnOsHvgzMAA4L9vmqc+6ZRMB+hpqtsoi/n+GO7ZJj\nL3fOZd3YQYV4Q9LjvCFpOy56z6+KyDrgWWAW8C/B9t8ARyS/3aHA71HxPjxpOyzZp/SDqJ38QtRm\n72kCNuWcZ17Co+PQG1Z8/APR6+pjOa+bBWyI2jYCe2S0Necco6Ex4a4wzrltqE3vFKDVOfd7snsk\nbc65lqAM9maSnvNJaM/x3HI20ESAz3DOHQy0oT2zK0Uk84IWkd1FB64eFZH16B+3KRLmJ4L648BU\ntNe4yzjnngJuRgW2GTgGFfM8nnLOtQCzUXPRp4JzXAT82N/0gHvQJ4NOtJf9ZPC+W8l5lEd74o8H\n+zr0s4emo2eC+lZUVEZC9wiO/UT8ooi3Aj9KXt+P2pdDc4kDPuSca0afBhaiTzKeG9Ge7XOBR5Jr\n8+agbQZJD90jIh2o2edrzrkfBpvWMlQ4B182pEFNP0cB10TtS1Ax/3ByQ8piE/q7hzShQh2yB7Au\n5xgNjQl3dfgu2lsZrrc8hEScLkYHHT8MbAY+UfZFCc65jcAZwEz00T6Lk9BH0UOdc01ob1so/fMt\niOo7gVUj/xS5XIqaS05AbatPD/eCxKb6CfSP601QjwPHRDe+3ZObw9NAj3+9iMxAb2hZLEfFzu8r\nwPykPfN0hjvfgKdGcOzc4yW93pcDbxX1GHoaeD1wbGC3huR3c87dhXYWzgzs5TcBz0N7vjclbXcn\n53EccGtos05s0tcCP3HOnRGd0h3odZNF1uc4BHgsHKwVkYXAdcBpzrlyN+0HgCmJyHueR2BaScZl\nppFvvmloTLirgHPuN2hv44Iyu+WZHt6PmkhOT3pp7wQ+nmFK0YOInCIiB4vINBGZjnpgrCW9oFeg\nHgeeWWjPcX0iAJ/NOK9/FpH9EhvracCVOXbKcsTvC2ryeAF6Q/ruSA+UmDzOAT6eNF0EnC46eIuI\ndCR2cID/Bo4XkReLyDTgVPK/6ytR883LE9vxSagnw+9y9h+NueiKUR475i3AfahYPi8p+6BPE/+U\n85pLUQ+jEwASG/FK9Jq4MWlzaC97sA0gMe39Evitc+5TGcf+FfCC5Dv1r5mcXHNTgMkispukHlXH\nAj8P9p2HmnsudM6VuLsm298mIn9NznEz+qRxWvKE+FLgeEpNa0cA1yfXxoTDhLtylPg/O+d+7Zxb\nl7UtYZ2U+nH/ayJEXwTemYzs45y7FxWtIRd7wgDqVvUs2ps7Ejgu8EI5Fbg0MSv8I+otMAPtQf8O\n+EV0bg4V1f9Ee6/TUKENtw/7+TPe15uRfoSaOn6Uc5y897kEmJMI9FdQT51rRWQDars9NHmPe1C7\n/+Vor3cjKl7bg+N63/P70SeAC9Dv7zjgeP/dj+Az5p6zc+6BUR475kTg6865lUFZgd60Tgz2C99z\nJ/rdfDzY/hvUzBWaJW4COigdmHwt6qH09uCa3JD0/L230w3A3wevOQUdg/hE8lm3Ap9Oth1LqX37\nXehT4Knh8YPt89FxHc8H0Ot0Jfrk+r7kv+B5c/JdTEhk9B0pw9h1ROQUYG/n3InD7lyZ95uFPoEs\ncc49Vov3bFREZD/gUufcocPs1wn82TkXu5mWe80vUbv3sKaPZHDzG865THfCiYAJt1EzEtPMn1AP\nkd8Ot/8Y3ud44HrUtHEOcIhz7oXVej+jlGRw/QXR4KZRQcxUYtQEEXk3Oqj4i2qKdsKrUbPRctTO\n/sYqv58R4Jx70ES7uliP2zAMo2CUy5EwZkTE7gqGYRi7iHMu05PJTCWGYRgFoybC7ZybUOWzn/3s\nuJ+DfWb7zPaZx17GS8OGw3rchmEYBcOE2zAMo2CYcFeBZcuWjfcp1Bz7zBODifiZ65GqugN6r5Jq\nvodhGEa18Akpa61hwfuaV4lhGEYjUFU/bsMwjCJTr9YC63EbhmEUDBNuwzCMgmHCbRiGUTBMuA3D\nMAqGCbdhGEbBMOE2DMPIQUQGfarrCRNuwzCMgmHCbRiGUTBMuA3DMAqGCbdhGEbBqIhwi8hkEblN\nRH5WieMZhmEY+VSqx/0R4B6gPgP7DcMwdoGRzkhTa8Ys3CLSAxwLXAzUn9+MYRhGg1GJHvd5wMeA\ngQocyzAMwxiGMaV1FZFXASudc7eJyLK8/U499dTB+rJly2wWDcMwjIje3l56e3tHtO+YZsARkdOB\ntwB9wHRgNnCVc+7EZLvNgGMYhjFKhpsBp2JTl4nIEcC/O+eOD9pMuA3DMEZJracuM4U2DKNhqNdc\nJTZZsGEYRg42WbBhGIZREUy4DcMwCoYJt2EYRsEw4TYMwygYYwrAMQzDaGTq1bHCetyGYRgFw4Tb\nMAyjYJhwG4ZhFAwTbsMwjIJhwm0YhlEwTLgNwzByqNdcJSbchmEYBcOE2zAMo2CYcBuGYRQME27D\nMIyCYcJtGIZRMCxXiWEYRg6Wq8QwDMOoCCbchmEYBcOE2zAMo2CYcBuGYRQME27DMIyCYcJtGIaR\ng+UqMQzDMCqCCbdhGEbBMOE2DMMoGCbchmEYBcOE2zAMo2BYrhLDMIwcLFeJYRiGURFMuA3DMAqG\nCbdhGEbBMOE2DMMoGCbchmEYBcOE2zAMI4eGzVUiIvNF5NcicreI3CUiH67EiRmGYRjZyFj9FEWk\nC+hyzv1FRGYBfwL+3jl3r4g4qF9fSMMwjHL43natNSx438zu/ph73M65Z5xzf0nqm4B7gbljPa5h\nGIaRTUUjJ0VkEfB84Jaw/aabbqK5uZnm5maampqYNWsWkyaZed0wDGNXGLOpZPBAaibpBb7gnPtJ\n0uYADgPWA+uS5RZgNtCclKagPljOO29Q7OMye/ZsE37DMKpOvZpKKtLjFpGpwFXAZV60Q14R1JcB\nL0UFPBTztUHbWuDRj36UtRnb1wKbgVlkiD2lwt/S0jIo9r7e0tLCrFmz6nKk2DCM+qKWgt3b20tv\nb++I9q3E4KQAlwKrnXMfjbbp4OSY3mEo/cAGskXdi/26qC2sb0N7+S2kQt8SLFtOP32I2If1qVOn\nVvgTGYZhpAzX466EcL8UuBG4g1SjT3bOXVMt4R4rOxkq7KHgrw2WWdunkwh8VjnttEGRb21tHRR9\nX6ZNm1bDT2oYRhGpunAP8+Z1KdxjwQGbgDUMFfg1wfoaSgXf13cjFfnWqN76xS8OCrwX/dbWVlpb\nW2lqajK7vmFMEEy46wgv+qGwhwIft60J6pvQAd3WqAyK/rnnDoq8F/y2tjbr5RtGATHhbhD60Z57\nKOarKRX61cF2v20tatrxQt9GJP5nn01bW9ug0IdLs+Ubxvhgwj3BcehAbij0q0lFPl76+jpgBqnQ\nt0Wl9fzzaWtrKynt7e3Mnj3bPHaMhqFe3QFNuI1MBlDBD4V9JGU7Q4W+3S+T3r0Xeb9saWkx+71R\nl5hwGxOC7aRCv4pSUV+VsXwWtd+3EAh8WM4+m/b29sHS0dFhPXujZphwG0YOfZQK/LOoqPvi28Jt\n21CR70AFPlx2XHjhoMj70tbWZjZ7Y9RMaOH+CjAVmBYswxK3Ze2zW7KcAlg/y9jOUDFfBawk7cmH\nZS0abduRlDlBveO885gzZw5z5syho6ODOXPm0N7ebkJvTGzh/iAa9LIjWPr69mHa/Lpv66NUyMN6\n2Ba375azPpIyPVrGbdOwGSnqnQFUvL2Qr6RU2FdGbauBPUgFfk5YLrxwUOg7OzuZM2cOLS0tZrpp\nQCa0cFfyHQYoFfMdpKLu13cE2+P9srbllW1l2rYF9R2kN4TpQfHrM6L16VF71npcnxG1zwiK3TQq\nzwBqq/di7pcrMuor0MRpscB3AZ1nnz0o8F1dXXR2dtLe3s6UKRVNzGk0GA0n3PWIQ8V7K6nA+/q2\noOS1xdvCtq05bVuD10ylVMizyu4jWIYlr21yhb6zRsObblYEy2dIhd2vr0B7/i1AZ1C6gK6zzhoU\n966uLrq6umhvbzePmwmICXeD40hvFHllS7QM27dk1DdHbeG2acBMSgV9ZlSfmVOflVMP1yeCRPWh\ndvgVlIq6F/ang7YN6KBrF9AdLLsvvJCuri66u7sHlzNmzKj5ZzGqgwm3UTEcpTeAzUndL319c7Qt\nrm8K1jcF7VtRM9CspMyMlqMpewTLIgf870B77V7Mn07KM0Hdr+9OIuroFFTdwNzzzmPu3LmDxQS+\nGJhwG4VhgFKB35hR30gq9puitk0Z+21Ae/F7ZBQv7LMzts0OlmF9JvXp1eRQm7wX8qeC+vKobQYw\nDxX3uUl93te+xrx58wbLnDlzmDzZDGPjhQm3MaHxpiQv6mEJ2zZk1OO2Deg4QyjovjTlrDfllOnV\n/NBl8AK/HBXycOnLk2jKg25U1Ht8Oecc5s+fT09PD/Pnz6erq8sGWauECbdhVJA+UiFfT6morw/a\nfD1sD2d9ElIRz5u+L2yLJ/3Yner2/Lejgv5kUJYDTyT1J1A7fSewAJifLBdccAELFixgwYIFLFy4\nkObm5kK7SZo7oGEYg2yjdOo+v4xndMqb5amPUiGPS5j2N1xvRU0llWAHanp5PChPJMvHkuUAsDAp\ni4BFZ53FwoULWbx4MYsXL6atra2uhd2E2zCMirGd0hmbwrzueXnefcKwSQxN8RtnfwwThLWj4r8r\nFu91qIg/BjyalMeAvyalD9gTWJws97zgAvbaay+WLFnCokWLxj161YS7IDj0YupDc2DvDOrxMq77\nMpCxPlCm7oJ6vO4ylkTrYVv8XY/muxfSx2+J2sIyKWc9XA5XJkf1yVE9LlOSktUWl4ngUrireM+g\nOMVvmOo3TA7m6+tR043PCeNLB9lpBOagAWfDsQ4V8EeAh4Plw6hppgfYOylLL7iApUuXsnTpUnp6\nemri3z6hhXsTQ6MO8yIYwwjHkZYwPD6uZxUvzFn1AVQUpjJUMOJlXA9LliD5pQTbsgRPgm1Zwuh/\nyVhAieoEbcMR3wDCtrgM5KyH7f1BW3/U1p9Rz7vpZd0g4xvozqgulAr51Gjp6359WrDuS5wzJyt/\nTpxaoVyqhawUCr4+lfr0VAnpQ4Xe54TxqQHivDBh2oAZqIB3BsvBYKOoZA3W7kBF/SHggaTcB9yP\njh0sBfYHDgAOuPpqDjzwQBYsWFBR08uEFu445Hu4fCFTKb34p0b1sC38U2X96cI/afhnzftDe8E0\niotPi+CF3At7Vlu5m3xWx2A7Q3PqjCTVQlbqBN+R6WdoOoRyaRDiFAjlomPzSrV9QRzam/aRo2EE\naRhk5H3QZ5L6n5e4KZJ6tXSSmms2APcm5a6k3IG6gj4PeD7wwksv5ZBDDmHp0qW73Duf0MJdJFOJ\nYdSafkrz4mylNBrW17PSIMQRsVlRsmEkrA+GmszQyNbdSQOewqCnOBBqj4ylL7tyQ4h90LNcFJ9I\n9ukm8V4hGewktZEvQG8WtwN/Bv4E3Jq87m+AlwF/e/PNHHLIIeNuOx8OE27DMErwvu2hkG/OKGFU\naxzYtImhvvCb0KfXrMClLD/3PDfIvCCn7aSuiI8FxdvFV6DeK/uQmlGeg5ppbgduBH6NDpAeCbz2\nsss4/vjjmT179i58i9XFhNswjJrgBz5D3/bYzz32a/eeMaHL4w5KXRvjeU/DgVGfibE1ed0jwIOo\nCeUe4O6kvhA4BDgMFfSHgB8DNwGvBj50660ccsghVfhWdg0TbmPCM8DIPYPKef/EbeUGZ/MGe8sh\nUT1r4DkcwB6JB0/eQHk8uB4ux9szZwepG+OaoO69XOIB0mfQp4FO1JTSQxoQtAgV7T7gTuB3wA3o\nb/U64FWobfx8tHd+/n33sXTp0up/yGEw4TaqygBDU9Rmlbz85qPxJhrOayj2FPLrjtF5BuV5/4TL\nLE+f2Msnz70yC5dRj90/YzfR0INnIKiHbVmuqVneOuHNK/TM8d9POc+cLG+ccl44w3nfDJe73g+8\n7hZ8n9tJE28tJw0CehTtXT+MCvrBwEuAJaj9+zuowJ+KmlO+AHweeH+NByNjhhNuSzTQ4PSR2iSz\nSpy9L07lGrZnpYbdjv6BQk+H8A8Wu7/NIPvPugdD/8wj8RzK8hjyrpx+3VIljR7/lBK7y8YutFme\nOPHNNfa+CesbGX6ykrBjEJY+StMJh0nDmlAzy3NRe3Y3ekNbC/we+CJqC78waXsT8BG0R348sFGE\nj4+zeJfDhLsO8bbCMBw6thPGJU6Q5AeN+tCLOvQAmBmVMDf23KAeupTNCNpikbaAl8ZjEulNs17p\nZ2h6YP8f8BGlq1Df71+hPuEPor3uL6E39g8AJwF/BI5FBf8G4EXAwTfcwJFHHgnU3h1wOMxUUkW2\nk9rm1kQlDEdeF6z7wZrJlCYhijPNZY3ch0sv0tMxv/RGIvRD38HQnm9fVA9Lf8Z6XvRvlnllJPb/\nrBJHBmdFBw+3T95Yw3CBWwPBd9cB/BPaE/8c8CngCOBQoBf9z70INa/0AqegdnEwP+5Cs5l0aqp4\nctlwpnFftlGaC6IlWIYJgOLMb82MX9pPI5sBhvpE5804NFI7fxxBXM7O78XYMdRUFEeCxm2hHT9e\nLxf9m2X3zxr0FIaOC+RFAudFB8fHmBS1lRt0jQdcJyXf1VRKnwgfRM0g56CeKFcCPwU+hA5mfhwd\nsHwNcCLqfXJf8tp6E+4JbypxaI/4KdJE82ES+meC4kjdj8LcDF2oLS10U2pHe7/W26092xnqjhaa\nj8JlOClDODNPaPvfkhwzjFDcnVKzUdaEzr742eLLRQtPy6jHNn4vuo1+Te1kqD+5L95sGLoT+qfV\njcn6A+jNDuBc4KNJ/WLgpcD70QHIA5L2ueiTL2iQzh+BtwJ/D5xZhc9XCRpeuDeQppiM008+gYqz\nn/LJh9p2owMXh5PO89eFmh8a/U9TD3gTU1bWu3WU+v6G6VD9H9qhYulNSt6EFEb7zUKfbnootfuH\nkYPxRMkT7bfvY6i3ULkninif7RmvzXo6iec9dQwdg4kjNH3wzr6kT6r+d34YtV2/A3hf8lnOB34O\n3IzeFH8AXJVs25C8FrTn/dukfuiufW01ofDCPYAK8EOkbj8PoQMRj6J3Xp8P2IfKHp0se9BcCDYD\nX/XwI/kroxImKAqz0a1GfzNvSvKmptCM1IM+4YTRd+EYQCObmfopP69n1gTQ5cw6eWaerehv558k\ndouW06P69Izts9Cnz7AtK89KnFtlLEm3voG69P0YDbYBHZg8C7gFvZauRa+l5yfbHwWOS+ozks8O\nmpGwXimMcG9Hs4L5hDL3oY9ED6M/wt6ob+YS1LVnMep838bE6ynVAocKrp8dJZwGKzQv+Sxxneif\n2C870N/sxahZKcwB3UhPNg4V1dAbyNezpk3LCikPTTrbSXOKeO+fMN/IjGDdi2Ez+iSZZdrJKl5o\nq53Nw08msR69drwZxH/ulej1NAc1bQzn1nkZcAbaY16ctK0D3gn8J9qbBrgEeBfpNXYf8G8Zx+tG\ntWVNnXmUQJ0OTq4D/pCUO9CR3UfQZDIHoBFO+6LmjL3RC9WoLAOoGcnnRv4rqcnpCfSPtjv6Z/BZ\n3EJTk0/fOYfG6AHvIB109l5CfhmadNZRatLZgApg+FQQ5uvImrw4XPeunN5kMF4mG4d+B+ViAsIb\nTlaJXVgdQ3OW+O9kJnoTnw98CzV5vKPM+d0F/C3wG3RQ0fNvyXtdnKxvRq/Rh9EOw0Dyfs8ky6uA\n7wE/Qm3tM4CdAwM1n6Wn7gcnHTra+2vUMf7/UFF4AZpb4Djgk8B+jCwxuzE6NqP5HO4izXn8IHph\nz0afYPZCn16OQP9IC1Fzxe61P92KspXSwWj/lODTkIZeQ1spzZURege1op2K0EvIB4DMZnx8oR1q\nIsnq1eeJa1YvPyxCaaBLGB8Qt3WjnSrvnhqONfgyGlfVm8kX7gFU2D9PqWivQCMj7wvabkC1pT1Z\nf5L0ZgH6nfnr2ptsdu7cybRp9eXRPi7CvQX4JXBNsuxDo5tego4AHzBeJ9bgrEdvjH8G/pKUJ9Cn\nl+egTzBvIjU7zRqf06wIDhXch1AbZjww/SQqRn4m8zCp/6GUmnQ60D93LftcDhXRrCnI1jF0kDYs\nXqSnMlQ4vdCG/v5d6G8ep2qdGbxmJuMXjDNcgNf/or/le6L2S4F/QH9LTy9wVLD+MNox8Wym9Al+\nKg0q3CJyDDpoOxm42Dn3paz9HBpO+h30ceRg4JXAh9G7ZKPYNOuJrWgqyxuSch/6vb8Q9Wc9BU1/\nWd+ZicvjUDPOXWg2uHtJZ0mZhArSIvQpYX/gGNLE/O3UNupzG2nvfgVp7z6cOSYcqJ1GaQxA6PPf\nggpubGYITTBj+V1Ds8hysiN04+x/4bp3rfQ+6NuSz/MS4LXAW6jcE9tFaPRj/FteAZwdtd0CnBas\nP4k6Knhi4Z4M9Pf3V+ZEK8iYhFtEJqPh/kehv+8fRORq59y94X63oT3pp9CBgrtRO5NRHW5D7YKX\no2J1FBp08CKKb24aQEW6F7Vn3oQKwoHoZz0ceC/69NCefYiq0Y+OA3hT0yOk3k2Po4LWhV77vnc/\nB33a8b7/HaRpTHf1t9qK3gjCKN24d54ltqHJpJ+09x1G6obLPZLz3ZPSKF6/vx/onIKaRTahv9d3\n0KjFD6I26KYRfKa8cbKt6LXw/ah9PXoDf2nUfg/6fXtWUNoj30bpmIxQf+HuMPYe96HAQ865RwFE\n5HI08KhEuI9BQ0zfVYE3NPLZjH7HNwPvRgd2e8b1jCrHDtQz4Ay0J7kMjXI7j9IeUy3pR4M1rkOf\naG5BxdebmvZE/yDeFbWTsfXwN6E3g0dJzT1Pkdrkn0V76gOktvfQldL3zOeiY0bx4GhoSplMaQBM\n7HY4QDpf5+7JsfygdN7T83S0t/3a5HOcht5sf4o+CeYh5Av3Xeh3HU+FcG/yGcOnDh/12hG0rUO/\nE49PmubZCDQ3N9edeI9VR+eh15DnSbRjV8JKNFrp/UFb3tcwkrSXtn/2/m9E/6APkPYainT+o9n/\ngTo7n5+jvbu491ip43uzwmLSuIR5wHdz9n9qlMd36FjHKlSod6ICviFn/9eRztm5BRXAp5LX5x0/\nZE/0RiwOYyx+AAAZDUlEQVSoE0K5/eM85SPhGfRGkrX/pOD42yj9zT6fLE8d4fvkeZvkCf1o989j\nrMJdX7ehCc6TwNtoDPe7onEs1R2n2YLm1PAePj6A7DMVfI/70SeEHtSEMwv4cs6+70GfOPwUaOtQ\nW+kZFTyfsdDPyMTNUT/ja729vfT29o5o37EK93JSv3aS+pPxTiegEU1tIzjgaO8Etn/Kl4E3ozbu\n96F/wCKdf7n9r0fdQrejA1t/BA5i6AVcy/P/KfAzNBKvB7WvH0BqKtmrQuezE/W+ehB90riO1FQy\nldL8Od5l8VSGJjJrQm3tTagox2abW0kHIv3A5CkMNZdsQXOATCbNt+5NMJehZqx5I/isN6PmjKOT\n441mkgnP9ZQONnqa0JtJuP8a9DdZG7RNJc1rAvBp9DOdkqwPJ+qj7SmX23/ZsmUsW7YMgM997nNl\njzNW4f4jsLeILEKflN6AepSVMAd9DPt34O3JulF5jkTd/c5CB+sOQ3+QIyn+YPCRqLBcC1yNPlk8\nic7efRgaAn8A+gheq3GU1yTFoYORN6O21StRkX0EFbdFaI9mPirwXaT2YO92WM7ZbCo6xVYW20lt\n3GGA0GpU5OOAIB8UtIVSt8A4HXC4XECpD3Zc34ORT1YxAPwC+Co6ePhF4J9H8Lo8Ad0fDdAboPRG\ntAR9gghpJg3/92ku9kDt2J5YyOuVMV3jzrk+Efkg2iGYDHw79igBdTv5AJq8fB80kOPNwN9ROjBg\njJ1FwNdRN6j/RnuF/4qKwzK0V/hCtFdYtAkQBO2dHZ2sP4u6mP4e+DbqMfA0Kt5LSSNrF6E24flU\nx6tGSNMthDi0lxdHnD5A6hK4AhXbmaSeJWH4f5izJSsN8G6kOXhGQz/Z4fZZQTrPMtT1L95/OqXB\nNT7K03uVzEB7679Fb14fQnt4I/k9yvVpfRK4W9GbuGchKuY+4hr0et8TvakemLS1o4O9nunk2+nr\niZqHvG9EeyRXoj2U56ICfhQqKGafrTz9qIfJr9Hv/Db0z/gc4HlJeQ4qdh3Uj81vV9iC/jHvRwUy\nDMBZjoqe7/mGYfpzKe0B13K6swG0JxxOgutD6uPgG1/8rOhTGDrRRiigsfteVmh96EmyKzhSoQ9F\nfSulftxTUT/u0d5kvo3eoL+ds/109Df/TtT+bvTm/bGg7c3Ay1G3ZFBT1zfQIB6AC9Br58JkfSaw\nYuNGZs2qbThaXU+ksA3167yWNEDkANQt5UXoQMkSitczLAJrUTG/PVnehV6wDu2lhnbavVBvhm6K\n/Vv0oz1cbx9enhQf9u57wGtRgQ9zrrcHJewN+9JE7ee2jKe42xAt85JYZfWuN5Fm9ItNJaMt4esq\n8Z0MJ9xr0U7H9WhH0HMzalK7LziPi5JjeW+cu9AxOG8m+A4atOZvAtOBtVu2MGNGbXOI1rVwx2wB\n/oT6w96alFWkiaWemywPQB+PitwzrEf8pBIPoj1VH0jyV/SRcy3aSw3ttWGSqXmo0BXdV78P/R7i\nGY5Wk51kag0qgLMoTT/rSzjl3GyG9ojDMl6B1QPo/y8vK+HGYLkpqsf7bSK1o/uJJOLPHZb4e/Ih\n9u3ok9DH0A5DHAUZ8h9oIqmbSb9Dh06M8C5UwEGv5RehN+op6I2vNTnvKcAPURPjlaS+6n39/Uya\nVNsuS6GEO4v16F3xzmB5N/rn8jbMfUhtjEvQH8JEvfJsQ3uqjydLb7MNe6+r0d6ojwz05oc5DM3/\n0UHj5ELvR8UqzhAY5xHZGNXjiMVJDE3c5LMEhpM8xGX3nBLmu671xM7+ZhBHZ8ZPBlmTYm8hTfTV\nij4Zlhtgd6hv+R6of7j/nLeiA8i3kzpFHIy6Lb4iWV+MPvXvjQ58fxP1y9+EXrObxyH4pvDCnUXY\nM/QuUr6H+BD6o4WBCmFZgN7JTdirw060d+rzcITZ9nxeDp+bYxVq98wakPOlJaM0o2LVaL+hQ71E\n4kx9mzOWWS56vh7OjRnOj7kNHQzMmnZtes56OAVbvJ5V/IBkvP+uPoX1MfLp2ragUdqLUfH1A5+f\nQE0hP0G14auo95UPk38l6jxxPGrrvgD1fHkUHcx/3IS7+nhR9/mjfXk8WG6h9DHfD07NQ+24vjRK\nb7Be8RnwVmcUPyAXDsqFPdntpI/Www3OZQ3Ehb3a+sr7Vj0GKJ1CLGs6sbgtnBFne1QPl3nTl/n6\nZErF3Yt6PPuNL+GTRPh7hb+tT50beqZsRs0ij6M97/3QAdIj0TGzc9Brak/UDNiGJrrbE/W++ilq\nS78azYHyGeC3dSjcRTdHDkFIB5GyQmlBf9xwgOoptLd+I+lA1TPohRX73PrHfj9w5RMD1TrtZyMg\npH/CxcPsG7OT9NE6Hpzz5oh1pMmd4hzUoU0W8uecjGeXiScKjnuv5Xqn3jVuvK6TSaTnX0sc2nMu\nd8MIp1gLp2RbiXbCQtt5OJPQWvSm0EpqS98HNZHsj2YIPAEV5JehHbaPoj3zK9FAtc7kfUCfAluT\n+u2ozVxEGi5XSSGZiQYE7VtmH4deFGEKTv+4/xClg1Yr0Z5HOMN71uN/6Jfry0Tp7VWaqaR28rHi\nU5huzijxvI2+bU3QFgtS1qS6vrc6QOns7tNJZ3XPWvoZ3uNluTIlKPH65Kg+OWibnFMmZdQnjbBI\nUPfCWkm8Z80a9Ddci3pH/U3SdlCyXytqt94TFfZjUbPI+9Df3yejupN0MobrKnyulaThTCXjxTZK\n8ynHy7WUmgD84/9ulOZYDpdZZoDYFDCTYrvoTTR8fg8v5ttR4YjrO4P1nRnLcqU/ox4u+5LSH5R4\nPSwD0bIf/U+Hy4GghG0uWkKpmMfiXq7Nl8kZ9bwbTni8R9D/09dQd8DZaJqIl6DpFLx9/Br0f3YQ\naXh8rXvcE87GXSR84EJow/V23PVk51COPRK2ouKdN9NJlk039Ebw67FdcUI+ihlVJxbz/mB9JG3x\njSO8mcRlB2mH6Vn0v7YaNZssRW3gp6O97N8AJ6NP0z9BU88eSJot0ITbqCg+dHm4AIvQFzf0u43N\nAn45hXy3stCuGw825dl246Uv06L6NMbXDmxUB2/n3kH6VBEvw6cQvx4OcIYmKW+mCvOEhwFH65P9\n2tCxqB7Uq2xf1KZ9JzoI+VI0SvKTaOKlXuBMNDT/xm3bmD5dY7nrTbitY1VwJpMGMFQKh/6ZvIhv\nzVj6euylsJZsj4OsP+aOaN2bCAZIRTws3oZbzs47JaOE7SO168aP4xK1CUMf6f0/TKLtMPRGFK87\n0g6Oi9qyStgbzTNTxD3SgWg9y1wSFt+2M2jbGSyHKzuiMpl8O35o84+XsTeKn/19BqWeJ774KNY1\nqEvfw6g74J1o9sK9ULPIFahDwkFoEq+r0SR4zwLXrFnDbrvV73xR1uM26g5vjw1tvF7U82y8WcKS\nZ+stZ9uNH70HopJnu/X221hcfVtI3v/BCz4MvQnExd9Msm4e5QYUw2XWjWxqtC288fltWTfLvBvr\nbkHbWMZiHNqR8ONGPqeLdxhYQeoRthztMCxAbdZ7oT3tvZPzvxW1Y9+NTib8NtRv+xvAvwCf2b59\n3CcHNlOJYRh1wU5Kx2niqEk/vhMX788/mdRLy3sU+VzknaTxFz3oDexeNNr6DjSVxt2obftwNFNm\nH+oSeA3weuDkRx5h8eLROqZWBxNuwzB2iZ2URmnGJU7/OtwM8DspnYnel3A+zNCzyqez9SUMiHNo\nL/tR0lw6Pnr6fvQJbT80v9FzgRegJpVb0BzUvcDzUR/vN61eTWtrK/WECbdhNBg+ND4OYolD3PN8\n0LPC5EOB9vU+1Ibskz7FGQC9N1OcJjb0cGoKljMY2aCzj372+W+eJM2N8wRpbvNZaK51bw7x2SyX\noiaa21CzyC1oRsAZaErXoy67jKOPPpr29vYRnM34UBfC/RGyPQrKeReEAxdZ696GZt4HRi3w7mXe\nzu7t71leEnkDr3EJg3NiL4qsEPIwqGcqpQN2WdGcoUfQTEoH80L3z9AdNHQX3Y3K/r+2k9qkn4mK\nt08/nZRZpGkofO50P8+mnzhiJvq7PExqErkD+EvyPgehYe4vuvxyDjvsMObPD2dZrG/qQrjPJdur\nILzIt5H9J9hJqcdBOEjVR3lPA982haGDKbEXgl8PI8viaLOs5XAlL2BgtN4KIylk1MNlFlnz+Y3U\nqyGvHns5xAN7sYdDOV/cvMHEeLkzagsHK+PBy7AtyxMiq+4oHWjz9XCZF/2Y12Hx7TOCtjBJUxw+\nH7pa1jr3dxb9qO3ZDxLG0cRepP1yM2qX7iTNHunLXMrnCdpEaga5Lyn3ognmutB0zwcCB15+OQcd\ndBBLlixh8uR6+JZ2jboQ7mq9g2NoVJkX9PjPl/WnzfI+CF2hYjHIWpYrWZFmsXCNxFuhnIDGbmRZ\n4jscklEPl5OCeuzqFt9YwhvSpKitXNRbXvRbnvtefAOdGtXDtqlRW95NPC+0fBr1IZTVpJ80WCUs\nYQRw6M2xCh1MbCZN9TAnYxmm9G2hvGfJerT37MtDSfFzZy5Bw9X3Bfa77DL23Xdf9t13X2bOnFmJ\nryCTQECr9h7DvG/jCbdhGKX4aMGs7IprghJPBrEBtUO3ZRQ/e3w4C1AHKsSjCQTZSpqt869BeSRZ\n7kBziQzaqy+6iCVLlrDPPvswb968mk9mACbchmGMgAHUGyNMdxCnQwjd5mIXuh2UemT4elvGMqw3\nM7anCp+U7XHSFMpxWYcGziwOyqLvf58999yTxYsX09HRMShY9YIJt2E0OAOoLTYrt0y8Hs+O44V5\nIzpQGLvIebe5ZrLd5rxQV2uCiQ2kqZDDZejpMRUV5oXo4OHCM89k4cKFg6Wrq2tces1jwYTbMOoQ\nhw6M5/khl/NLjnNDb0ZFN/RVjn2Ww9IcLX2ppT29H7VZL4/Kk1F9gHTiET/n6PxvfpP58+cPltmz\nK520dfwx4TaMCtDH8AEhcXBI1iQK4foUhvogz85oawrWw3o42069DGJ604V3sXuKNBw8rK9Abxyx\n6928iy9m3rx59PT00NPTQ1NTU92ZMWqBCbcxYRggDfgYSYnnUYy3hfMt7kB7tWG62qzgkHgZ18O2\nIk1m4XvI3vc5q4QzOHWTTs/XDcz7ylfo7u5m7ty59PT00NXVVdfJlCYqJtzGID7rX1a2vzjKbkuZ\nEs8ME6eF3U5psEfezORxmZVRnxXVRxp9VyQcam6Jg1LCQJWnk+Uq1J7dTTqtXhfQff75zJ07l+7u\n7sGy++61nqTMqBQm3HXKAGlgUTyxala03HD1vLDnOAXrZIZG2GXNpxi2xfWwLW9uxkYT19Hi0MFG\nn7kujBhcEZVn0IG9cH7TLr/81rfo7u6mq6uL7u5uOjo6mDp1avx2RoNhwh3Rx9CQ5G2UD1GO80ln\nhSnn1bPWfVsYLZcXJZc1SUHWJAZhiSc6CAW5XmywRWQHpVGBcXRgXJ/O0Emmu4DOiy6is7OzpFQz\niMQoHnUh3HeQHeGYtR6XcPv2nHrejBpZbTA0N8r0aD0rLHl6RttIwpb9ckbGfsVyjGo8+kintVqZ\nsYzrm9DgEx8J6MO355xxBnPmzKGzs7Nk6WdPMYzRUhfCfQBp6PC0jHpWDoi8bXGeiLBeLlGVb7cp\nfxqXbZSGZfuyiqG5NJ5FXfhaSEO0w3Dtzq9/nY6ODjo6Oujs7KSjo4OWlpbC+SEbxaQuhLueTCVG\nMfC5M1ZllGcz6s+iT1Q+HNsvB0sgxL60trYWOhGRUX3MHdCYsMQJjLzo5tV9AqPZpALsc2a0A+1n\nnUVHRwdtbW10dHTQ3t5OR0cHs2fPnpC+xkb1qFfhNsuBMSq2MTSDXFZZFdTXoyIcJizySYvaTj+d\nvRMR9kLc1tZGa2srU6bY5WkYWViPe4LSR3amuHL11cnrsjLItQFtZ51Fe3v7oAj7YiYJo6jUa4/b\nhLvgbCdN3bmmTFkd1beQJiZqI53Xrw1oPe20EtENRXjmzJlmjjAmDCbcRi4+ci4vf3JefQ1pGs8w\nRadP59l22mm0tLSUiLAvs2fPNg8JwxiGhhRuETkbeBWqHw8Db3fOrQ+2TxjhdmjCojBPcijCoRjH\n6+tJU3m2ZizbvvQlWlpaBktbWxstLS20trYya9Ys6wEbRoNRbeF+BXC9c25ARM5M3uiTwfZCCXc/\naW7kOFF9nLA+LuvQQJs4P/Lg+uc/T3NzM62trYOi6+vNzc02EGcYxiBV9Spxzl0XrN4C/MNYjjdW\nHGq7DWcMiYW3XNsmNGNcVqL6ZqDlC1+gJ+n1hiLc3NxMS0uL5ZAwDKMmVLKb9w7gB2M5gBfedWMo\nUxk6S0hYn/flL3NAc/Og2IbLpqYms/sahlH3DCvcInIdmhsn5lPOuZ8l+3wa2OGc+37WMS4lnaop\nXIbz6vm2KQzt8fqZQVqA9jPPZO+WFpqamgZFt6mpaXBpuYUNw2h0hhVu59wrym0XkbcBxwJH5u1z\nDmn2uwPe9z5e/MIX0pz0er0A++W0aUVKa28YhlEZent76e3tHdG+Yx2cPAbV5SOcc6sytuvgZI1d\naQzDMCpBo7oDPogm31uTNP3eOfeBYLsJt2EYhaUhhXsEb27CbRhGYalX4TYXCsMwjIJhwm0YhlEw\nTLgNwzAKhsVZG4Zh5FCv43PW4zYMwygYJtyGYRgFw4TbMAyjYJhwG4ZhFAwTbsMwjIJhwm0YhpGD\niNTlDFMm3IZhGAXDhNswDKNgmHAbhmEUDBNuwzCMgmHCbRiGUTAsV4lhGEYOlqvEMAzDqAgm3IZh\nGAXDhNswDKNgmHAbhmEUDBNuwzCMgmHCbRiGkYPlKjEMwzAqggm3YRhGwTDhNgzDKBgm3IZhGAXD\nhNswDKNgWK4SwzCMHCxXiWEYhlERTLgNwzAKhgm3YRhGwTDhNgzDKBgm3IZhGAXDhNswDCMHy1Vi\nGIZhVAQTbsMwjIIxZuEWkZNEZEBEWitxQoZhGEZ5xiTcIjIfeAXwWGVOxzAMwxiOsfa4zwU+XokT\nMQzDMEbGLgu3iLwGeNI5d0cFz8cwDKNucM7VZb6SskmmROQ6oCtj06eBk4G/C3fPO86pp546WF+2\nbBnLli0bzTkahmE0PL29vfT29o5oX9mVu4mIPAe4HtiSNPUAy4FDnXMrg/0c1G+GLcMwjHrE+447\n5zI7xLsk3Blv8lfghc65NVG7CbdhGMYoGU64K+XHbcpsGIZRIyrS4849uPW4DcMwRk2tetyGYRgN\nh+UqMQzDMCqCCbdhGEbBMOE2DMMoGCbchmEYBcOE2zAMo2CUDXk3DMOYyNSrK7P1uA3DMAqGCbdh\nGEbBMOE2DMMoGCbchmEYBcOE2zAMo2CYcBuGYeRguUoMwzCMimDCXQVGOv1QI2GfeWIwET9zPWLC\nXQUm4sVtn3liMBE/cz1iwm0YhlEwajIDjmEYhjF6bAYcwzCMBqGqPW7DMAyj8liP2zAMo2CYcBuG\nYRSMqgm3iJwgIneLSL+IvCDadrKIPCgi94nI31XrHMYTETlURG4VkdtE5A8icsh4n1O1EZEPici9\nInKXiHxpvM+nVojISSIyICKt430u1UZEzk5+49tF5Eci0jTe51QtROSYRKMeFJFPjPf5hFSzx30n\n8FrgxrBRRPYH3gDsDxwDfF1EGrHnfxZwinPu+cD/S9YbFhH5W+DVwIHOuecAXx7nU6oJIjIfeAXw\n2HifS424FjjAOfc84AHg5HE+n6ogIpOBC1GN2h94k4jsN75nlVI1wXTO3eeceyBj02uAHzjndjrn\nHgUeAg6t1nmMI08DvjfSDCwfx3OpBe8HznDO7QRwzj07zudTK84FPj7eJ1ErnHPXOecGktVbgJ7x\nPJ8qcijwkHPu0eSavhzVrrpgPHq6c4Eng/UngXnjcB7V5pPAOSLyOHA2DdozCdgbOFxE/k9EekXk\n4PE+oWojIq8BnnTO3THe5zJOvAP43/E+iSoxD3giWK8rnRrTnJMich3QlbHpU865n43iUIX0SSzz\n+T8NfBj4sHPuxyJyAnAJ+khdWIb5vFOAFufc3yT2/CuAPWt5ftVgmM98MhCO0dRfGrldYCT/axH5\nNLDDOff9mp5c7ahrTRqTcDvndkWIlgPzg/UeCmpGKPf5ReQy59xRyep/AxfX5qyqxzCf9/3Aj5L9\n/pAM1rU551bX7ASrQN5nFpHnAIuB25O0nz3An0TkUOfcyhqeYsUZ7n8tIm8DjgWOrMkJjQ+xTs2n\n1FIwrtTKVBL2RK4G3igi00RkMfqIfWuNzqOWPCQiRyT1l6MDOY3MT9DPiYjsA0wrumiXwzl3l3Ou\n0zm32Dm3GP1Tv6Dooj0cInIM8DHgNc65beN9PlXkj8DeIrJIRKahDhVXj/M5DTKmHnc5ROS1wFeB\nduB/ROQ259wrnXP3iMgVwD1AH/AB15jhm+8BviYiuwFbk/VG5hLgEhG5E9gBnDjO51NrGvEazuIC\nYBpwXfKk8Xvn3AfG95Qqj3OuT0Q+CPwSmAx82zl37zif1iAW8m4YhlEwGtF/2jAMo6Ex4TYMwygY\nJtyGYRgFw4TbMAyjYJhwG4ZhFAwTbsMwjIJhwm0YhlEwTLgNwzAKxv8HoXhTScB+lYwAAAAASUVO\nRK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x109d25310>"
]
}
],
"prompt_number": 6
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"MCNAB.plot_imex_stability_region(N=200);"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXecZUWZ97/PRCYxw2QmMjOdZkiCgOu6yoiKKGbFsKuo\nu+u6xt13eQ3o64ruGkBAFxEDQUFBRTDhmtC1VRSzC0icIUxgYBLD5GF6uuv94znVt271OTd033Ru\nP9/Ppz6nTp2659S5fft3nvPUU1XinMMwDMPID2Oa3QDDMAyjOky4DcMwcoYJt2EYRs4w4TYMw8gZ\nJtyGYRg5w4TbMAwjZ5hwj2JEpFdE/iHj2BIR2S0iEtcVkb8TkR81sq3DQUTOFZHL63Tu/xSRrSKy\nqR7nbzYiskpEft/sdlSLiFwoIv/c7HbUGxPuChGRh0TktCT/BhEZEJGLozovTsq/mOwflezvjtJZ\nyfEbReQL0Tm+JSKXZLRhhohcJSKPiMguEblXRN4THB8QkeVV3JZL0tADzq13zk1zhUD/wbrOuWud\nc88dwXUHEZHVwXe0S0TuE5F/Gs65YpxzH3POvakW5woRkSXAvwE9zrkFNTyviMgDInJnyrFeEdmf\nfE/bROQ7IrIoOfZ5EbksqDteRPZmlJ0iIl3J57eIyHYR+aGIdEWX/A/gE8Hn3y4ifxCRA/73ndLG\nc0XkIyLyFBG5OTn3FhG5XkTml7jvis+d5F8iIncmv5c7ReTFQdULgfeJyPis67UDJtyV46L8/cBZ\nIjI2KH89cB9DxXB6IoI+fSMpfyvwMhFZDSAirwKeBLw3ow2fBCajgnE48CJgbVRHqrutmjGS6z6c\nfC+HA/8CXCYiR9eoXfVgCbDdObe92g+KyLgSh58BTATmiMhJ0TEHvM05Nw1YARwGeMPh58lnPScB\n64CnR2UO+CMwHfg20AXMA34HfCdo45HA6qSO52FUzK8q0f7nA/8NHAF8DliapN1AqiBXe24RmQtc\nC/xb8nt5F3CdiMwBcM49CtyD/m+0LSbcw0OAR4E7gOcCiMhM4KnAd6lQxJxzm4FzgMsTK+6/gDc7\n5/ZlfOQk4KvOuZ3J5+91zt2YXP8XSZ3bvFWfWOjfS6yex0TkJhFZGJ2zQ0R+KyI7ReTbInJEcj7/\ntjDkN5K8cfwy5bq7ROSVInKHiLwgqD8+sRKPr+A7+QGwHViZfFZE5L0isjY5x9d9G5PjZ4vIuuTY\n/5PiN6PzROTLQd0XJRbaDhH5mYj0BMceEpFzROQ2EXlcRL4mIhNT7v3ZwI+BBcn3fFWF5363iNwO\n7E77ThNeD9yIiujrS3xHO5M6/uH2S2Bl8hsE+Bvga8AUEZmVlD0d+LVzrt8593vn3Bedc4875w4B\nnwK6g+/1OcAfnXMHg2t+yzn3HfRvM4Tks13Arc65HzrnbnTO7XHO7Qc+AzytxP1UfG6gA9jjnPtR\n8tnvA3uB8I2vFzgz63rtgAn3yPgycHaSfzX6z/RESr1MIXfOXY1a738CfuCc+3GJ6/0G+EginJ3R\nebzFdVxg1Y8BrkQtxCXAfuDSqF1nA28EjgQOAalumhLtD697uHPueuAa4LVBteejVvVtpc4lImNE\n5EWoRfjnpPidqPX0jKSNO1AhQERWJfnXJMemA6HrYvDNJ3EFXJecbzbwfeCmwAJ2wFnog3gZcBzw\nhpT7/QnwPGBT8j3/fQXnBv19PA+Y4ZwbSLn3ycDLga8D1wOvTnnd9/0Ns4CXAb9N2rSBYgv7GaiY\n/zoq+wXpPAN4xDm3I9k/Frg3o27Wb/m5wE8C11p8/r9kfK7ac98OHBKRF4jIWBF5CXAgKffcA5Q1\nEvKMCffI+BawWkQOB14HXJ1Rb1tiifnUHR2/BZgJfKXM9d6Bvia+HbhTRNaIyBlZlZ1zjyXWzAHn\n3B7go8CpYRXgGufcXYmV/wHglSIyUnfLtcCZIjI12X8d+pDLYoGI7AD2od/p65xz9yfH3gz8P+fc\nJudcH/Ah4BWJi+oVwHedc79Ojv07xW6q8D5eBXzPOfdT51w/6gudBPx1UOcS59yjiYDdhLqt0oi/\nn3Lndsm5H3bOpT3YQYV4l3PuV8D/JGWh1SjAJSLyOLAVmAq8LTj+c+DU5G93Cmqd/hJ4RlL210md\n4htRP/mlqM/eMx3Yk9HOrMmNzkQfWPH5j0N/V+/K+FxV505+x29GH3AH0N/amxPL3rMbmFHB9XKL\nCfcIcM4dQH16HwBmOuduJd1qmOWcOyJIg9ZMYjmfg1qOF5fygSYC/DHn3EnALNQy+4aIpP5IRWSy\naMfVQyKyE/3HnR4J84Ygvx4Yj1qNw8Y5twn4FSqwM4Az0H+wLDY5544ADkfdRe8L2ngU8C3/0APu\nQt8M5qFW9sbguvvJeN1GLfH1QV2H3nvoOno0yO9HxbESjqzg3BviD0W8Hvhm8vl+1L8cuksc8A7n\n3Az0bWAp+ibj+QVq2R4LPJD8Nn8VlE0isdA9iV/4x8BnnHNfDw7tAKZltHPI7ztx/Twb+GFU3oEK\n7juTB1I5yp5bRE4EvgA83Tk3HjVErozccNOAxyu4Xm4x4R4516DWSjlreQiJOF2Bdjq+E/XVvafk\nhxKcc7uBjwFT0Ff7NM5BfYOnOOemoz9yofgfZEmU7wO2VX4XmVyNukvOQn2rj5T7QOJTfQ9q8XkX\n1HrgjOjBNzl5ODwCLPKfF5FJ6AMtjYdRsfN1BViclKc2p1x7AzZVcO7M8yVW72nA60Ujhh4BXgk8\nP/BbQ/J3c879BTUWPh74y3+JugfOTPIAdybtOBP4XeizTvzGPwa+7Zz7WNSk29HfTRpp93EysC7s\nrBWRpcDNwIedc6Ue2tWe+1nAb5xzfwJwzv0BfSA9O/jMSuB/K7xmLjHhHiHOuZ+jP5pPl6iW5Xp4\nC+oi+Whipf0D8O4UV4qeROQDInKSiEwQkcPQCIwdFPyRm9GIA89U1HLcmQjAB1Pa9VoRWZn4WD8M\nfCPDT1mK+LqgLo8T0QfSNZWeKHF5XAS8Oyn6HPBR0c5bRGRO4gcHuAF4oYg8VUQmAOeR/V1/A3Xf\nnJb4js9BX7V/nVG/GnfR9VWeO+Z1qF+2CxXf45P8RuBvMz5zNRphdBaAc24tsAX9TfwiKXOoqA2W\nASSuvR8Btzjn3pdy7p8AJybfqf/M2OQ3Nw4YKyITpRBR9Xzge0Hdhai751LnXFG4a3L8DSLy4HDO\nDdwGPN1b2CJyAurHD/tPTgV+kHJfbYMJ9/Aoin92zv3MOfd42rGEx6U4jvtfEyH6CPAPSc8+zrm7\nUdEa8mNPGEDDqrai1tyzgDODKJTzgKsTt8Ir0GiBSagF/Wv0xxyHNV4DfAm1XiegQhseL3v/Kdf1\nbqRvoq6Ob2acJ+s6VwFzE4H+LzRS58cisgv13Z6SXOMu1O//NdTq3Y2K1xPBeX3s+b3oG8Cn0e/v\nTOCF/ruv4B4z2+ycu6/Kc8ecDVzmnNsSpM3oQ+vsoF54zT70u3l3cPznqJsrdEv8EphDccfkS9EI\npTcGv8ldieXvo53+B3hJ8JkPoH0Q70nudT/w/uTY8yn2b/8j+hZ4Xnj+4PhitF+n6nMnnfcXAN8U\nkd3ow/sjSaexD2VcSXEoY9sh1RtXhlEZIvIBoNM5d3bZyrW53lT0DaTDObeuEddsV0RkJXC1c+6U\nMvXmAX9yzsVhpqU+8yPU750VuTKSc18IrHXOfa7Sz+QRE26jLiSumT+iESK3lKs/guu8EPgp6tq4\nCDjZOffkel3PKCbpXD8x6txs+XPnHXOVGDVHRN6Edir+oJ6infAi1G30MOpnf3Wdr2cEOOfW1EtY\n63nuvGMWt2EYRs4oNW/CiBEReyoYhmEME+dcanSTuUoMwzByRkOE2zk3qtIHP/jBprfB7tnu2e55\n5KlZGlYOs7gNwzByhgm3YRhGzjDhrgOrV69udhMajt3z6GA03nMrUtdwQB9VUs9rGIZh1As/SWWj\nNSy4rkWVGIZhtAN1jeM2DMPIM63qLTCL2zAMI2eYcBuGYeQME27DMIycYcJtGIaRM0y4DcMwcoYJ\nt2EYRgYiMhhT3UqYcBuGYeQME27DMIycYcJtGIaRM0y4DcMwckZNhFtExorIn0XkplqczzAMw8im\nVhb3vwB3Aa05sN8wDGMYVLoiTaMZsXCLyCLg+cAVQOvFzRiGYbQZtbC4Pwm8CxiowbkMwzCMMoxo\nWlcReQGwxTn3ZxFZnVXvvPPOG8yvXr3aVtEwDMOI6O3tpbe3t6K6I1oBR0Q+CrwOOAQcBhwO3Oic\nOzs5bivgGIZhVEm5FXBqtnSZiJwK/F/n3AuDMhNuwzCMKmn00mWm0IZhtA2tOleJLRZsGIaRgS0W\nbBiGYdQEE27DMIycYcJtGIaRM0y4DcMwcsaIBuAYhmG0M60aWGEWt2EYRs4w4TYMw8gZJtyGYRg5\nw4TbMAwjZ5hwG4Zh5AwTbsMwjAxada4SE27DMIycYcJtGIaRM0y4DcMwcoYJt2EYRs4w4TYMw8gZ\nNleJYRhGBjZXiWEYhlETTLgNwzByhgm3YRhGzmiIj/scETouu4yOjg46OztZvHgxY8eObcSlDcMw\n2o6GrPJ+PrAWWJNstwFHAZ1AB9Dxmc/Q2dlpom4YhkH5Vd4bItzxFfYD91Ms5muStBUV9Q5U2Ds/\n8xk6Ojro6Ohg6dKlJuqGYTSMQECbdd3WEe5SpIm6z2+hIOodQOell9LZ2WmibhhGXTDhrgH7gQco\nttIrEfUlS5YwbpyFrBuGUR0m3HXGi3pooXth3wwspeB+6fj0pwd96ibqhmFkYcLdRA4w1FLPEvXO\nSy8djH4xUTeM0Y0Jd4sSinrsgkkTdXO/GMbowYQ7h4Tul9ha34KKeifF0S/eUreOUsMwhosJd52I\nRT202H1HqY9T77Q4dcMwqsCEuwnEIY2htR4PPuq87LJBUV+0aJGJumEYJtythhf12EpfA2wHllFw\nv3R89rODPvXFixczZoxNLWMYo4GWEO5FSVocbX06EpsYHGAfQ0Xdpx3AcgKf+uc/PyjqCxcuNFE3\njDaiJYT7QeBhYAOwMdj6/DZgLgVRX8hQkT8SGF+3lrY+eyiOfgnTTmAFgahffvmgqC9YsGDwR2AY\nRj5oCeEud4U+YBMq5FkCvwWYTbbVvhhYwOgU9z0M7ST1aS+RqF955aCoz58/30TdMErQtuGAIrIY\nuAY1mh3wBefcJcmxmvm4DwGPUizqscB7cS/lmlkATKhBe/LCLrJF/QBBjDrQ+cUvDnaUzpkzx0Td\nGPW0s3DPB+Y75/5XRKYCfwRe4py7u9Gdk17cvZh76z0U+M3ALLKF3Vvuo0HcHydd0NcA/QSCDnRe\ncw2dnZ10dXUxc+bMJrXYMBpL2wp3ygW/DXzaOffTVowq6QceodjHHrtlHmWouMdW/ELaW9y3U2yp\n3xfkxxKJ+rXXDor69OnTm9Riw6g9o0K4ReQo4OfA0c65Pa0o3JXQT7HlnrZ9FJiJCnlssfv9drTc\nHTpnepalPgnoIhD1r3990P0yderU5jTaMIZJ2wt34ibpBf7TOfftpCyXwl0JseWeFi3jLfdyHart\nIu4Ovec0K/1+YDqBqH/844OC3tHRwaRJk5rUasPIpq2FW0TGA98DfuCc+1RQ7gA+GNRdnaTRgLfc\n06x2n4997mHyZe0QCjmA9jmkifqDwBwKot510UV0dHTQ1dXF8uXLmTChXR5thpFNb28vvb29AHzo\nQx8C6ts5KcDVwHbn3P+JjrWtxV0r4miZNIHfggpb3KEaR8vkdRBTP7CeoaJ+H3r/CwlE/ZJLBi31\npUuX2gyNRlvSiKiSvwF+AdxOQaPPdc790IS7Nhyi4JbJCofcSvEgpthq94OY8jYTykHgIYpF3Qv7\nZgrzvnQCXckUAV1dXTaa1Mg1LTEA51/R+Oo5KdsjyJ+YtCJ9qLinuWXWJ/ntwHyGCnoo9POAvMhd\nOO/LfdF2J4UY9S6g66qrBkXdYtSNVqclhPsi1CLciorH1iDtQsU7TdR9PtyfAxxWtxa3NwfREaqx\ntR7uP45a5mmivghYgv4NWl32dpPuelmDvsEMul6Armuvpauri87OTgtnNFqClhDuUlc4BDxGQci3\nRdu0svEMFfM0gfdlM8iPFdlsDqCdiLE7Zn1QtofCfDJpAr8YDZVsVXHfzlC3i99OIRB0oPOGG+jq\n6rLIF6OhtLxwV4tDhaMaod+LRm6UE/hw3+IYstlHaX/7BtS6j2Pb49Rqtq1D3U1exO8N8g+ibiQv\n6F2f+hRdXV10dXVZJ2kb09bhgCUu3hKdkwdREY/FPdwPy7ahllcs5rEbJ9yfSutamM1gD8VumDRx\nF4qFfQlDLfgpjW54BoeAdRSHMd6bbB9BO0kHRf3znx8U9SOPPNL86TnGhDtHONTXW4lV71M/2eIe\ni/wc1JUwmt03Du1AXE+6v92XTaLYSl9CsYtmETCxwW2POUBheoDYWt9Hseul68tfpru726YHyAkm\n3G3OPkpb8lujtBvtlC1nyYf7zRaoRuPQ78772NMs+E3o9xgK+xIKVvsSmhsGuRMV8HuDrRf4KaiY\ndwNd558/aKWvWLGCiRNH21+7NTHhNoo4xNAIm1J++m1oNE0pcfdlc5P8FNrffROOTs1K29AwyNBa\nj1OjI2W8Pz0UdN9Zup7CoKNuoPuyy+jq6qK7u5uFCxea66WBmHAbI8K7FtIs+C0pZVuTz5QS+Lis\nXaNvDlI8xa+PbQ/391EId0xzyywBpjWwvQ9SLOY+vxt1vXSTCPtXvkJ3dzfd3d1Mm9aoFo4eTLiN\nhuPdN1lWfSz4YfRNlsDPDcpnkd9h9jF7GWqpr6dY4MeTLuxhp2q9o5G868WL+T3Jdi0apdPt08UX\n093dTU9PD0uXLmXsWBvmlidMuI2K8dE35ax6n9+BisVcsjtkQ8HPc5ilQ8cbhNZ6vH2EwvJ6sah7\n3/tc6uOS8ZN43UuxoN+L/r2Wo4LeA3R/6Uv09PTQ3d3NjBkz6tAaY6S0hHA/HX3NPDzZ+ny8H5cf\nTv5nxWtn+ikePFUupYVZzqW02OdplKyf6ncDGjoYu2XWo64O75KJRd3v13rW8n2oRX4PBUH32ymo\nmPcAPRdfTE9PDz09PSxZssSs9CbSEsLdi/5gdwVbn/dpZ5AP642jWMhLCX25B4I9BJpLGGaZ5pf3\nZVuCsgmUFvZ4f3LD7mZ4+MFLXszXBVtfNpmhor6U2kfJODQqx4v43RREfRs618ugqF977aCVPmVK\nq0TXty8tIdzDvYJDJxIKxbyc0KeV+fx4Krf2S70FTKN9fLutjEP/drGYpyV/fAzlBT4sa7XIGx8C\nGYt5KPCPoeIdinos8CPtptyD+tG9le5FfQ363a1EBX3l5z5HT08PK1eutMm7akiuhbuWZD0EdpL9\nFhCWhds9qCWYZunHAp9WFpZPxWZHrBXxdAixqKeVQWUC30pC/wRDrfbYcp+IirgX9VDcl6DhkcOJ\nIOpHp9m9O0n3BPkxwCpU1Fd98pOsXLmSVatWsWjRIhP0KjHhrgMOfeVNE/W0slIPib3o6MAsoa/k\nYeDzrSAqeWMvxUIei/yWKO8YKuaxwIf5ZrhuvNW+nqFWu9/fScHXvjQlLaa6jmSHzo9+N3BXtN1L\nIubAqgsuYNWqVRx99NEsWbKk5edMt3BAI5UB9IedJvppQp/mDvL1DpDu6y/3MEjLT8QeAmnsJV3c\nY+EPXTexwMfCH5Y1qjN2P8ViHor6OtT3PYuCkB/FUHGvtBN1Byrid1IQ8zvR3+1K4GjgmAsv5Jhj\njuGYY45hwYIFLWOhm3AbdecQpX38O6NtqYcDpIt9JaJ/OBom6B8io9UVFLtuQlHPEvqJFIt5muj7\nNJv6dbj3o+IdivpDwXY9+obnRT0tlfOzP44K+J3AX4JtH3AMcCxw7Gc/y3HHHcexxx7blAFGJtxG\nrniC8m8Bpd4IwjeDSZS2/Cstn0R7vwX4zlgv4pspiL7Ph6K/Hf1esoQ9TkdQu5GxLmnHOnSU57oo\n/xD69nAUsIyCmC8L9rPcSFuAO1ARvz1Jd6F++eOBJwFP+s53OOGEE+ruPzfhNkYl3hVUyrrflVLm\n64f5Q1Qn9KWOtcNbwADqhtjMUMt9M0MfALspzGVTLs1jZP55l1zzIQpi/mCQ1qFTLCxP0rIgvxyd\nqyV8yPSjsej/G6Q/J+UnJunkG2/k5JNPrqmYm3AbxgjpI931U8odlPZmsAe1BqsV/bSUp76Agwy1\n4MNtKPSb0YfbPIrFPBZ3v51Fddb8ADpY6QEKYn5/sl2LulGOAlag8eQrKCwKvZRCOO4jwB+BPyTp\nd0m7/wp4KvC0W27hpJNOGvZsi6NauE9Hf+CHRduJJcorqRPvt3b/tNEqeN9zmvjHln98PO4fcFQv\n/Gl1ptJav1//HcXC7rdborJdqHiHYp6V5lJ+HMQ+VNTXooK+lsKc549SWLjCD+PvQaNWZqDW/G+A\nW4Fb0PjzU4DTgOf85jecdNJJLT8qtCWE+weoz/RAtA3TgZTytPpp+QOoNTGedGGPy7IeElkPjlLH\n0643gfxYYcbIiPsCKnkLCF1Avt4+tLOvWt//9JQ6zZgPpo+hAp+VtqPtno8K+fyU/Hx0kFGaJX8A\nFfNwXhYfUz6VQsfmccAJaNjjb4CfAj9O2vA84KybbuL0009nwoTWm0GnJYS7ES8ZDhXvNPHPeiCU\nylfz2XjbR3WiP9w3kFIPGr9vIzzzQT9D3wJiH3/aW0Dag2AsI+8IPhx9kNTjLaAfjTP3Qv5osn0k\n2H802d+NWugLUCFfkKSFydavguTXEtpAoVPzNuBP6GClE1HXyTNQa/1nwPWo2L8OeNvataxYsaIO\ndzs8Ro1wtxIDlBf6ah4IWQ+TUm8kfgulHxRZZdW4rcq9qZgrq3E49O9eqehndQTvRt8CplJZ+GdW\n3oeETmZ4b6FPoGK+KSVtpDDP+hgKg4l81EoH6k6Zi3Zm/hr4OeoHPxk4C3gKKuBXAGcC5z/yCPPn\nzx9GS2uLCfcop4/KXE/xw6Qat1UlD5rYlTUSV9VI3GETMVdWpfRT2biAcgPEdqG/w1ITxVU6cGw6\nQ5fwC9cvfYhCZ+dadL6VDWikygmo1f0UVPS/DvwQeA3wTuBq4CrgM8ArGtwZGWPCbbQE3hIcrltq\nf5nPVnPeQ6h4V/KWkFVW6kFSbRpH+z9I+qhskri0KSLihwQUhHwGKubT0Tj1I1C/uA97nJ/s70Kt\n7ltRP/dE4J+AlwOfBb4MXJl85mXAB4F/bKJ4m3AbRkSWK6sSV1Wph0+58rjOwSQ/QLGQTyBb5ONj\nE1LK47JwOyFjPyuNp/VcXAcouH18ejxJj1GYI953km5EO0Q7UIv7xeiApYtQMf9eUvflqBW+CHga\n8J1bb+WpT30qMErDAU24DSObfoYKu+9oPxiVxaKfdjzcZpWFx/qCY31RWR/a2TmegpCXyvs0LmUb\np7Fl0pgoSZRI2YZI0HYfbrkNtbBPAz4NXAZcjg7m+RHwDjRS5Xrgk8Dvk3OZcBuGkRsc6lrqS9LB\nlHxY1hfV9/n+YL8/KDuEvnH0R2kg2Lok74JEyjZmIEn7Uav8J8CNwJOBk1DBPhUNG/wcamWfDHwU\neBbasXm/v0aLCXdDosVeQvqTNM6XOz4mqlOv45XUK/eZtONp+dB6MIxWQyhY0a3GAXTI/3aK53HZ\nifq+XwO8G+1wfAD4EDog5zQ0tHAAvb/jUYF+GvAc4FfJ9mXAJxp2N9XREOF+PelP0oGoLD7un7Lx\nsYNBmavgHAMMfaqnXb/SemltSisLzxW305/XUfw6mPbAqPTBEZdJlZ+J21EuHz+AKrl+qfPUOkkF\nZWn7pcqytmmv8mlllaRm4KI0EG3j33CalewtaL/11nZssae5aw4EaT8aiui3e4LkOyr9qNUjgJkU\nT7I1A3VxvAON1/4YGs/9k2S7H/Vtnxjcu//el6PCDTpsvlVpiHC/tBEXySn+n6PUw6fUA6bcwyQ+\nXsnDJD5/2vXjz6bdh69zkPR/+jRByBKJ+JqVlLlhljmGtsml1Bso8VmXUVYqxWT5c8koyyLNtZCW\nD6+Z9eDLeoMN0/hkG/u5vT98YpKPO1UPo7CoyOQkPyVIYcjgjKR+2v0/DDwbeCMqvnuAvwe+igr7\nt1FXiV/ffktSDhqdsjvJd6acu1WwgXVNxnegtPbMCUYjSRN0Fx2Py8qR1pkX59vBZbcZWA28CXWT\nALwVeCXqzwa4AR1843kIHbgD+jAaSPJHohb9Y02O6U7DhNswWox2EdFy+IEz29Foj3DBiXAY/MNJ\nuaBLqv0j6gaJGQD+FngVBdG+HfhvdHIqULfN94Hzg/0N6GhLUHeO9+dPR8MOnXMtsyKPx4TbMIyK\nGUD9zvvQedb3oq4IP+e690PHi27HMdc7kvwkCgNmQj/1ItSdMR+dl2QOKtx3oG6QNOH+Auq//lBQ\ndj7wLtS9Ajr16+LknKAx3nMojMY8QGH5uPHJNQ8dOsT48a3VPWvCbRgtig/FSwvBi8PxSqU43rvU\n1AS+czAt7UvqHkax73kKGic9LWW7lMJQde+bnoG6IGZQfbTKU0l/G9kPnIfGYnu34+Po4JpLg3q3\noBNNedZRcJOQ3GO4gMR44ODBg+0n3CJyBvAp9Pu6wjl3fpmPGG1I2ElXLqV1hJbrFC11vFQEUKl8\npelQRj6MnkiLqMhKfRn7aTHQYxk6wCUc1Rhuw06/sCxtdOUU0ueAmZSynRxsD0P9wM1CKPigQ76F\nLml2fFD2I+Dp6EPC82fUYvdsQmcZ9OynWLjHAv39/SNocX0YkXCLyFj0gfZs1BX1exH5rnPu7rDe\n5ZTvTS93bGAEZdX08FcTHVBp3VKfT4taiPPxtlyEw3CPxZEZpfbjMigO/xOGhgT6smrCCdPi3sem\n5OOytNDKrP1SaQKFUX5po/3GB/lxUVk4SnB8dGx8lA8jL8YF120tz2rzCTsPQ74NvDoquwV4ZlR2\nH/D2YH8b6irxhK4S0O+/0YNvKmGkFvcpwFrn3EMAIvI1dCqAIuH+LYWndKnYVcmoB0NjYqspy4qp\nTbtWNfFHGxU0AAAaC0lEQVS4ldbN+vzYEnVLxQ/HscZZMcXDKRtbot7YMmWGUQ27UaHYhrpyJlAI\n+5uNLqwQzwQ4hoLBEP7mfgd8JKp7N/CCqGwT6j/37KQwlzcUZrEM2zhjxoyWE++RCvdCtFPWsxGd\nNbGIK1M+mPU1ZFkYVt/qW/32qf9a1L2xL6P+YjS65DhUlD1erOPw2XEU+6pJPn9GyrmXBO3xPnvP\nh5Pth6iMrGiTLKGvtn4WIxXu1noMGYaRC+5FreHrM46/CZ1P+xsVnu8QQzs6D1XwOUfruKN6e3vp\n7e2tqO6IJpkSkb8CznPOnZHsnwsM+A5KP8nUb8l2aRAdS3NzjDSlXavS4cjl6hmGUT13AjejS4dt\nRi1f32k6FXWVLEAnhIpf4cdQmLXQMxN9GIT+6qcCF6JzkHhmA3dRGCl5Hup28Zb2B9AHwL8n+/5/\nvN0mmfoD0CkiR6Huo1ehc7sU8TbSOxEJyqvpsKw2pV1rJJ2V4Z+wUv92mn86y19drqyUL7zcPBpp\nc3fE5eX247LY951VljZkOqujstrJutLKKumYzOqoTOuIjDspza8/fI5OUrX4/8fYVXIMGjFyelC2\nHB14Ewr3XIqHuE9FhcszHn0otDojEm7n3CEReTuF8Mkr44gSKMxp204MJ8okK18qwiQr6iQr+qSS\neTT6y9SL5wApVebD6+JzZM1BEofx9QXH4vNUM1dLVlkc2ldNaGAc+heH//kHTBxJEkeUZKUwiiTe\nj+e3jue+TlvwIG0hhKwFFcLl3vI03ULcKel5DhqzHQr3Sejq7m8IypahMwEek+zPQhcX9kxEBxG1\nOiOO43bO/QD4QQ3akivMVWKEgh7OL50Wyx3WyZqrOo7jTksH0VjjXRQWOvADbeJFELIG4MQDb8ZS\nEHEfv+1juH2anJLCwTfhwBuf/PJiw10oOI0s4f47NMTtI0k7QEX84ugzx6Krv7842Z+HDqv3TEK/\nX89kYPPu3bQaNnLSMIZJODnYpGY2ZAQ49GERj5yMR1D6aVb3UjzkfQtDp12Nh70/gYqpXxvSj5yc\nmSQ/5N2vEzkXFdRppAt+Wtly4HnoIggfS8qOTs59M/DcpOwp6NB4z2KKw+ImUxzp0g+MG9d6Mmkr\n4BiGUVf8QsHh2pA7KCyC4JOfZGpzkgbQGfoWoELu0Njsiag/O+ZRdOTkt4C/Tsq+BHwR6EUF/zHg\nKPSBc1jSrnnoQ0iAa9FJqa5Lrj8e6OvvZ8yYxvZotMQKOIZhjF7GU7Cuq2EP6sbYhAqtoDHYJ2XU\nnw9cgU7Z+lOgB40X/wTwHXQlrpmou+RnqIXu5/h+BH1ATETfEEDFfBI0XLQrwYTbqJi4c7OSOUgq\nmZek1DmrTfG5XEadeD+trJLpAyotC7+/uC5BGSn14zJS8uWQlG1ayG0lkUZpUwqU65RNm2cl7DA9\njGLf+mTUT95JdQsavBC15E9DY8Cfhs7JcTY6udRMdDX3G1HhBnWzPIgKt1D4XreirpZWpCWEOxaE\nctu0yIDhHM+qN9JzjuT8tThezWfSRDbtM6GwxPOCpIXvlSoThgqAlKlXSdmYlOvG9dPKpES9+Bql\nyrJCNNNEMg7zJCgjpX5cRko+izTxj8uyHjrxAyjt99yH+r/DyJt4YeB4RsOwEzWcnfAJCv50KHSE\n+s5PvwKO95WHPvI5qNvjVOB1wN8At6LzlbwMHdRzAyrcHwM+j/4dF6ITLUHxXCXrUAu/FWmIcM8m\nW8x8guzJg7LideNY3ErrScZnS01UVGmbfH58ldeXCq9faQxzpWWx+Ka1aUxUzzAaQR9D5/zeiUbU\n7KLgK78fHeQXLhY8HbgI+CU6EOcB1IXyJXRJs7noupNPRiNj9ibX3ISKP8nxW1B/c7vNVVIRd1Ne\nJE0QDGN04C34eDrcrLfTAVQjfLjhkRVew6Hi+8/oaMkpqMUPKtxbk/xm1LgEHXNyZpL/0XBurkFY\nVIlh5JhDFMdpp23DeO64rFSKF2nIWsghTvGc4vHc437wUtZ0uVlvhrGbqRz+GhNQl0oHOqT9CuAq\ndPKq3Wj8931J+1aiw/B3oZb6tuRc7Tbk3TBGDQOkD2BJG9iSdazUSjRZg2Ti8nAfihdISFs0Ie4I\nTFtkISybFJTFCzPEIzizUtghOTba1muecYcKrl/DcgtqTW9EXSV/QDs6n4W6QG5HV3//D9SafwHw\nZtRn/hrgHODcOrSzFpjFbbQ0XiwPkC525cpL1UnrGCslvn0UxDBOh1V5bEJ0fEKJz8efC+u2g+Xl\n0A5Jvz5luFblrii/M8qH61juRDsyZyVpLhoiuAAd6r4CfXj0olEl+1E/+Aq0M3MZ8GXUtbIR+GFf\n3+CSZa1mcZtwG0NwZFuGaSIY58uJazWfDcUyXGIrTQzjpbhqdcyX+8VjDf2NHGCosJbK747Kw9GV\n49FoER81kpafTnFUSZj8Wpbj0IfAQ+i0sPehfWx3oH7uFWiUiV9g4XLgf9D5t09DOy5nADfu3cvk\nyZNDAa3Bt1Y5Jtw5op90ARuJcJb6bNaxgxR8g6FwhcKWJXZZ9X2alHKOUp81sawthyiI5k6qE944\nP5ZssfWhe9NTymIhnkrliwY/gY6QfDhJG5O0AQ3fW4da30tRke5O0rFotMifgB8D30/234CGD34W\nDRX8MPCWJoyUjDHhLkM5v2W1ollOQEuVDVCdEKZZnuF+tSIZ5ltvrNjoxVGYWCoreTHdSbHIxnUO\nki203pINRTVNfP12Qg3vbzcqyI+ioxjDvE+bkvuYj8ZeL0CXIVuEzjmyFI27PhIV79tQof4jOkvg\nDnTx4Gcn2zuAr6DD598M/OvWrcye7eNLmktLCnf4Kl6pL7KUJTmc42FHTyl/ZClxq+R4JSLpt7Y4\nbHvRT7aIplmypeqMZ6irIEt0Q4s2tohrOVNfKXxH4eYSyYvzZtR6n4+K7rxkOz8oOxIV6jkUGxXb\n0SiQu1FXyF3oNK270GXPTkBjtU9J6v8EnYvk16hr5DVf/SoveclLOOywcAGz5tMSwr2IYjH1C3JW\n4m+s5PV8uP5P3ylkYml4HPobrUZYsyzfA2jscCigcT62aNP8u4dTuSuhnvRTHK0Rbn1+c1A2HhXh\nrOSFeR76PWWxCx1ksyZI9yWpD52TZCWwKtkei2rOPWj0yC/RDskx6CyBz7vhBk4//XSmTZtGq9IS\nwr2OoWJpr+JGLRlAR9dVY8XuJF2gId26jQU2rcMsLJtCa//OHYWpWUMBjpM/tgPtuPPTroZbn58X\nlE2usB0DqOX9QJLuD7ZrkzauQOOw/dwlXUmahxpem9Bwv9+hoyh/j1rnTwP+5vLLeeYzn8ny5csz\nF+ttNVpCuFvZx200lzTrNvbTlvLZ+ryfyS1LaEOfrd9Ps2wPRw2MvHIQHRHoRTfMpyVQ8ZvDUCGe\nG5XNZvjhhzvQiZweQkX5wSA9hH7vPmRvOSrSy5P9Iym8Fbuk/p/RUZF/StJBdNbAk4Gn3HQTT3nK\nU5gzJ1yBMl+YcBs1p5+CdRtar7HAViK+kG7JlnItpLkZptHa1u1w6Uf9uLEIZ+X3oCLshThMfhKm\nUKRLuSiq4XFUULPSACrMy9D5sJcH+8sy2rEH9VffjnYk3pbkp6C+6ycBT/7mNznhhBNYunRpXaxp\nCwc0mko8yKGScK8sy3cfpSMOqolWyLN1OxwGUJELrV6/gECaIHv3hBfbLEH2onwEte+zcahve12Q\nHoq2/aggH4UK8dIoP7NEu/rQFdr/gq7+fkeSHkH918cBx110EccffzzHHXdcQy1pE25jWHhXQjx4\nIc2KLTf4YQKlBzZUKsRTaU/rdjg4VIhD4c0S4S2o9TyVodbvnIyy2dR/Md9+1Ee8nmJx9sK8Hn3A\nLg3SUdG2lDB7+tCOxbtQgfbpATSM7xi0Y/Ho66/n2GOPpaOjo+nLhplwjyL8Uk3DHVEW5h1Dow1C\nX22adZsmyK0SmdDqOPRNIxbfNEHeilqikxhqEWcJ8mxqF/9cKfsoiPJ6hgr0pqRdS6O0JMgfXsX1\nDlAYsXhXkB5Aoz2O9unaa1m1ahU9PT0tF47nMeFuccKY22qENS3fR7bAVmvdtubPOT94IS4lvnE6\njGLRTRPiMN9Md88Aeh9ekP0IwlCgd6Mi7IU4zC9FxXQ497ATFec4bUR92CuTdPR113H00UfT1dXV\nsgKdhQl3HUjrJCslsllW8G7USphKugWbNYw3a0hvowY5jEZCH7GPKQ4FOSzzrglvEaeJ7+yUslby\nu+9DxTgU5Dg/jYIYh+K8FB1ROJfhu7YcOrTcD3K5J8jvQoeTe4FeeeONrFy5ko6OjsHJmfKOCXdC\nKLZpQlqpS2E3+qP2SxrFboEsSzbLCjaxbQ796MrbsdWbZQ2HPuIwzc0om01rCXGIj18OreX1UdqN\niq8f0r2YYlFeQuXx0qU4gMZM3xOle9Hv2wt0z6c+RU9PDytXrmTRokVNn9OjXWkJ4e6kILb70XCe\nUlZtlnUbH2/1AQ6jkT5UXEuJb5geR//Gfs3AckI8i8b7iIfLToaKcbi/CY0CWUJBhH2+FtZyjEMH\n09ybpFCcH0Y7GleiIt3zxS/S09NDT08PM2bMqFELjEppCeG+GxPbvHKA0sK7jWKR3oNGGHixjV0R\ncZpFPueUfoLCrHShGG8Itv0UuzAWp2zr8TawH7We701J4yjMmNdz/vn09PTQ3d3N8uXL28a90Q60\nhHDnoXNyNOBQYc0S37R0kKHWcJzC4zPJ/4O5H7VMQxGOtzvQEX2xEIf5esRUewZQK/leNILjnmR7\nLxr/vIyCQHdfeeWgQM+aNatOLTJqiQl3GzOACkia+GaJ8TjShThLnA+nvXz/DvWpb4hSKMyPoKIb\nC3EozvOpf3w1aJ+Ot5bvC/Jr0L9NNzpnR/eFF9Ld3U13dzfLli1revyzMTJMuHNEH8WCmya+Ydlj\nFHfUVSLItejIamX2MFSUQ2HegMazx0IcivRww+OGy0E0xtkL831Bfg/aR+St566vfIXu7m66uro4\n/PBqoquNPGHC3ST8zGtpQpwlyHtQn28lAjw7SaPJK3mAgl859C+HwvwEQ8U4FuZmTOY5gHZGhqLs\n8xvQh0UXiThfdtmgOC9cuDA3M9q1IxYOmHNCt0SWCMeCLBSLbjlBnkH+/cPDpQ8VtliQw/Q4Opl+\nKWGeRXNdOzsoCHIo0mvQB4Z3bXRdcAHd3d10dnayYsUKJkzIS6zM6MKEu8V4gsrEd1uSHkP/8UIR\nLifEtZp5Le+EnX1x2ohay1vR6JMsQV6E+pVb4cG2j8JCtPdRmNj/XtTt0U3BvdF57bV0dXWZayOn\nmHDXEb9MUiWWsM8fQEV2FtniG4rwLEaXW6JSHPp9ZonyBgrxyqEPORbmBbRWWGAfOld0uNqKF+mt\n6JDuQXG+/PJB63nevHnm2mgjTLiroI+ho+m85ZvWUbcdHZQRCm2WT9jvT6e9oiXqQVYERujOeBjt\nIPVW8ZIg70V5Ia05enEAvZdQnH1+Hfow8X7nzksuGbSclyxZwtixjYgpMZrNqBbuXWQLb5o17Adx\npFnDWWX5mrqm+fjpSGOfctj5t5FCBEaa68JvWzlSxaGDg7wghyJ9P9qv0EWyHNYFF9DV1TXod544\nsRUfN0YjaUvhFpFPAC9AXXv3A290zu0MjjvQf+wsF0RcNht9rW4FX2ae2Ul2BIYvG0OxEIcWs7eU\nW3c51WIep1iYw+04AnEGOr/2tUGBnjp1apNabBjZ1Fu4nwP81Dk3ICIfTy703uB4y3ZO5pndDLWO\nY4EeYKjLYlG0nd7oho+QPRQ6BcPVvtegfRbhIrKd11xDV1cXHR0dNlrQyB0Nc5WIyEuBlzvnXhuU\nmXBXyR6KBThNnPsoLchelPPow9+PvrqtidJ9qFW9AhXmDqDriisGLWfrFDTaiUYK903AV51z1wVl\nJtwBeyltJW9AfU6xEC+iWJTrOQdGIzhIIWIjFufN6Cx1nUnq+uxn6ezspLOz06YRNUYNIxZuEbkZ\nDaGNeZ9z7qakzvuBE51zL48+O2qEu5Qo+46+/QwV4VikK1m7Lw8cQiMz0sR5I4WRgh1A1yWXDIrz\n0qVLbZ4NY9RTTrjL/oc4555T5gJvAJ4PPCurznlBfnWS8sReNOwttpA3Bvv70c68UJCfhPbcelFu\n9qi+WjOA3rv3M6+lOJxuHoHlfPHFnJ6I87Jly2ykoGFE9Pb20tvbW1HdkXZOngFcBJzqnNuWcrzl\nLW5vKZeylkNRTvMnL6L9RNnj59iILec16MRIs0isZqDzggsGLefly5czadKkJrXaMGpDu4YDrkHH\nvjyWFN3qnHtrcLypwr2XoZZxvD1AtvvCi3W7irLHoUtopYnz/ej0oR0E1vMNN9DR0UFHRwdTptjA\nfqN9aUvhruDidRPuMPoiy6cci3KcHw2i7AkHooRuDZ+fREGYO9BYZ289T5uWl2huw6gtJtxVsIuC\n+IYp9C3H0Rex1byQ0SPKHoeOPo2F2afxFMS5E+i87rpBcZ4+PW9R3YZRf0y4E/yIvtg6DvcPURDf\ntNC4dgiJGwnbKRbkUKShMEqwE+j48pcHxXnmzJnNaK5h5JZRLdynUxDlAYZaxrFfeTSLsucxhroz\nfH6Ags+5A+i8+upBcZ41a5YNRDGMGjGqhfv7FMQ5ryP66sHjDHVneIHuoxCt0QF0fulLg+I8e/Zs\nE2fDaGNaQrhbORyw3uwkXZjXUphfoxMdyt151VWD4jx37lwTZ8MYpZhwN4CdDO0M9PsHKLg1VgCd\nV145KM42v4ZhGGmYcNeIXQy1mL1I70PFeTDW+cor6ejooLOzk/nz55s4G4ZRFSbcVbCL9A7Btehg\nnlCcO664YtByNnE2DKOWmHBHhJZzLNJ7KB4huOLyywfndD7yyCNNnA3DaAijUrhjt0Yo1Gni7C1n\nE2fDMEJGdThgPa5QqTh3AJ2JW8MsZ8MwqsGEexjE0Rprg7SXJEqDgjj7DkETZ8MwaoEJdwZhnHPs\nd97P0GiNFStWmDgbhtEQRrVw7yA9lG4tKs6DVjPQYaF0hmG0CKNauKcSiXMyQrCjo8MGoRiG0bKM\nauEeGBgwcTYMw6iQcsLdkCWzTbQNwzBqR0OE2zAMw6gdJtyGYRg5w4TbMAwjZ5hwG4Zh5AwTbsMw\njAxEpCWDK0y4DcMwcoYJt2EYRs4w4TYMw8gZJtyGYRg5w4TbMAwjZ4xrdgMMwzBalUZPLlUpZnEb\nhmHkDBNuwzCMnGHCbRiGkTNMuA3DMHKGCbdhGEbOMOE2DMPIwOYqMQzDMGqCCbdhGEbOGLFwi8g5\nIjIgIjNr0SDDMAyjNCMSbhFZDDwHWFeb5hiGYRjlGKnFfTHw7lo0xDAMw6iMYQu3iLwY2Oicu72G\n7TEMw2gZnHMtOV9JyUmmRORmYH7KofcD5wKnh9WzznPeeecN5levXs3q1auraaNhGEbb09vbS29v\nb0V1ZThPExE5BvgpsC8pWgQ8DJzinNsS1HPQujNsGYZhtCI+dtw5l2oQD0u4Uy7yIPBk59xjUbkJ\nt2EYRpWUE+5axXGbMhuGYTSImljcmSc3i9swDKNqGmVxG4ZhtB02V4lhGIZRE0y4DcMwcoYJt2EY\nRs4w4TYMw8gZJtyGYRg5o+SQd8MwjNFMq4Yym8VtGIaRM0y4DcMwcoYJt2EYRs4w4TYMw8gZJtyG\nYRg5w4TbMAwjA5urxDAMw6gJJtx1oNLlh9oJu+fRwWi851bEhLsOjMYft93z6GA03nMrYsJtGIaR\nMxqyAo5hGIZRPbYCjmEYRptQV4vbMAzDqD1mcRuGYeQME27DMIycUTfhFpGzROROEekXkROjY+eK\nyBoRuUdETq9XG5qJiJwiIr8TkT+LyO9F5ORmt6neiMg7RORuEfmLiJzf7PY0ChE5R0QGRGRms9tS\nb0TkE8nf+DYR+aaITG92m+qFiJyRaNQaEXlPs9sTUk+L+w7gpcAvwkIRWQW8ClgFnAFcJiLtaPlf\nAHzAOXcC8O/JftsiIs8EXgQc55w7BriwyU1qCCKyGHgOsK7ZbWkQPwaOds4dD9wHnNvk9tQFERkL\nXIpq1CrgNSKysrmtKlA3wXTO3eOcuy/l0IuBrzrn+pxzDwFrgVPq1Y4m8gjgrZEZwMNNbEsjeAvw\nMedcH4BzbmuT29MoLgbe3exGNArn3M3OuYFk97fAoma2p46cAqx1zj2U/Ka/hmpXS9AMS3cBsDHY\n3wgsbEI76s17gYtEZD3wCdrUMgnoBJ4hIr8RkV4ROanZDao3IvJiYKNz7vZmt6VJ/D3w/WY3ok4s\nBDYE+y2lUyNac1JEbgbmpxx6n3PupipOlcuYxBL3/37gncA7nXPfEpGzgKvQV+rcUuZ+xwFHOOf+\nKvHnXw8sb2T76kGZez4XCPtoWm8auWFQyf+1iLwfOOicu66hjWscLa1JIxJu59xwhOhhYHGwv4ic\nuhFK3b+IfMU59+xk9wbgisa0qn6Uud+3AN9M6v0+6ayb5Zzb3rAG1oGsexaRY4BlwG3JtJ+LgD+K\nyCnOuS0NbGLNKfd/LSJvAJ4PPKshDWoOsU4tpthT0FQa5SoJLZHvAq8WkQkisgx9xf5dg9rRSNaK\nyKlJ/jS0I6ed+TZ6n4hIFzAh76JdCufcX5xz85xzy5xzy9B/6hPzLtrlEJEzgHcBL3bOHWh2e+rI\nH4BOETlKRCagARXfbXKbBhmRxV0KEXkpcAkwG/hvEfmzc+55zrm7ROR64C7gEPBW157DN/8J+IyI\nTAT2J/vtzFXAVSJyB3AQOLvJ7Wk07fgbTuPTwATg5uRN41bn3Fub26Ta45w7JCJvB34EjAWudM7d\n3eRmDWJD3g3DMHJGO8ZPG4ZhtDUm3IZhGDnDhNswDCNnmHAbhmHkDBNuwzCMnGHCbRiGkTNMuA3D\nMHKGCbdhGEbO+P99Dd/wc0BBIgAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x10ab5b9d0>"
]
}
],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"CNLF.plot_imex_stability_region(N=200);"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYXFW1/vHvIhDmJAQSpgSDCDLKIEblCrQKyGCCAyBc\nZvwpV0RQmQQEghN6AQcE9QqCoCCijFGZBBtQEBFJIBCGoExhCEMmCJCh1++PtU/qdKWqu7q7qrpO\n9/t5nnq6u+r0qV2Vznt2rbP3PubuiIhIcSzX3w0QEZGeUXCLiBSMgltEpGAU3CIiBaPgFhEpGAW3\niEjBKLilJmbWbmafrfLYBmY238ysfFszO9DMbm5mW3vDzE42swsbtO9vmdnLZvZ8I/bf38xsczO7\nr0H7vtfMNm/EvotMwd0AZvaUmX0kfX+YmXWY2ffLttk73X9J+nlc+nl+2W3f9PjVZvbzsn1ca2bn\nVWnDCDO72MxeMLN5ZvaYmZ2Ue7zDzN7Zg5fl6bbsA+7PuPvqXpoUsHRbd7/c3T/Wh+ddyszacu/R\nPDN73Mw+35t9lXP3s9z9c/XYV56ZbQB8FdjU3der437NzP5tZg9XeKzdzN5M79MrZna9mY1Jj/2f\nmf0kt+0KZvZGlfvGm9km6fdnmdmrZnaTmW1S9pTfBM7O/f7I9Lf5evq/cEAXr2NLM7s5Hdg6Kmxy\nDvCN2t+ZwUHB3Rhe9v2TwL5mNiR3/6HA4ywbhsNTCGa336X7jwI+ZWZtAGb2GWAb4GtV2vADYBUi\nMIYBE4EZZdtYz15W3fTleWem92UYcCzwEzPbok7taoQNgFfd/dWe/qKZLd/FwzsBKwKjzGz7sscc\n+KK7rw5sBKwEZB2HO9LvZrYHngZ2LLvPgfuB4cB1wCbA2sA/gOtzbVwXaEvbZC4A3gJGAwcCP+2i\n17wQuBKo+GkOmAx82MzWrvL4oKTgbjwDXgQeAj4G0SMBPgjcQI0h5u4vAccBF6Ze3I+AI919QZVf\n2R74jbvPTb//mLtfnZ7/zrTN1KxXn3rof0g9q9fMbLKZrV+2z3elj65zzew6M1sj7S/7tLDM31P6\nxHFXheedZ2b7mdlDZvbx3PYrpF7i1jW8JzcCrwKbpd81M/uamc1I+/ht1sb0+CFm9nR67Otln4wm\nmdmvcttONLOHzWy2mf3FzDbNPfaUmR1nZlPNbI6ZXWlmK1Z47bsAtwDrpff54hr3faKZPQjMr/Se\nJocCVxMhemgX79HctE12cLsL2Cz9DQJ8iAjOVc1szXTfjsDd7r7E3e9z90vcfY67LwZ+CLw7977u\nCtzv7gtT+1cFPgWc5u4L3P1v6fkPrtK+x939EuCRKo+/RRxAPlbp8cFKwd08vwIOSd/vT/wxv11h\nu6pB7u6XEr33fwE3uvstXTzf34Fvp+DcuGw/WY/rPble/XLAL4ge4gbAm8D5Ze06BDgcWBdYDFQs\n03TR/vzzDnP3q4DLgINym+1J9KqndrUvM1vOzCYSPcIH0t3HEJ8sdkptnE30/kg9vguAA9Jjw4F8\n6WLpJ59UCrgi7W8t4E/A5FwP2IF9iTDZEHgPcFiF1/tnYA/g+fQ+H1HDviH+PvYARrj7MuUDM1sF\n+DTwW+AqYH8zW6F8s7TtmkSQ3pva9Cyde9g7EWF+d9l9d1LZTsAL7j47/bwV8Fju8U2Axe6e/3Q3\nldKBozemA90eyAcTBXfzXAu0mdkwovdxaZXtXkk9sez27rLH/wqMBH7dzfN9CbgcOBp42MyeMLPd\nq23s7q+5+7Xu/pa7vw58B9g5vwlwmbs/knr5pwH7mVlfyy2XA3uZ2Wrp54OJg1w165nZbGAB8Z4e\n7O5PpseOBL7u7s+7+yLgTGCfVKLaB7jB3e9Oj51O5zJV/nV8BviDu9/m7kuIOuvKwA65bc5z9xdT\ngE0mylaVlL8/3e3b075nunulAztEEM9Lvdnb0317lT3neWY2B3gZWA34Yu7xO4Cd07/deOAeIrx3\nSvftkLbp/EKiTn4+UbPPDAdez/28GjCv7FfnA6tXeS21mA+M6MPvDzgK7iZJH/n+SATeSHe/h8q9\n6zXdfY3cbWlvJvWcjyN6jt/vqgaaAvgsd98eWJPomf3OzCr+BzCzVSxOXD1lZnOJ/7jDy4L52dz3\nzwArEL3GXnP354G/EQE7AtidCPNqnnf3NYBhRLnolFwbxwHXZgc94uP3YqI2uy7wXO553yTKLJWs\nR7y+bFsnXnu+dPRi7vs3icCqxbo17PvZ8l8qcyhwTfr9JUR9OV8uceBL7j6C+DTwDuKTTOZOoue8\nFfDv9Lf5t9x9K5N66BkzG0WUfS5w99/mHppN51B+nfi3yRtOhG9vDUvPI4mCu7kuI3or3fWWl5HC\n6SLipOMxwBvASV3+UuLu84GzgFWJj/aVHEd8zB3v7sOJ3rbR+eCyQdn3i4BXan8VVV1KlEv2JWqr\nL3T3C6mmehIRClkJ6hlg97ID3yrp4PACMCb7fTNbmTigVTKTCLtsWwPGpvsrNqe79uY8X8O+q+4v\n9Xo/AhxqMWLoBWA/YM9c3RrSv5u7TyM6C9/N1cvvIkoPe6XvAR5O7dgL+EdWs07PuQYR2te5+1ll\nTXqQ+LvJPA4sb2bvyt23NTCt2muqwWZEuUUSBXcTufsdwC7Aj7vYrFrp4QtEieQ7qZf2WeDECqWU\n2InZaWa2vZkNNbOViBEYsynVI18iRhxkViN6jnNTAJxRoV0Hmdlmqcb6DeB3uSGAtSp/XoiSx3bE\nAemyWneUSh7nAiemu34GfMfi5C1mNirVwQF+D0wwsw+a2VBgEtXf698R5ZuPpNrxccQoiburbN+T\nctFVPdx3uYOBR4mw3DrdNiE+Tfx3ld+5lBhhtC9Aqj/PIv4m7kz3OdHLXnofQCrt3Qz81d1PqbDv\nPwPbpfcUd3+D+DTwjfQp7kPABHLlL4sT2Tvlfl4JGJq+XzF/ojc9th1waw3vzaCh4G68TuOf3f0v\n7j6n0mPJHOs8jvvLKYi+DXw2ndnH3acTofVzKusALiFqnDOBjwJ75UahTAIuTWWFfYjRAisTPei7\ngRvL2uZEqP6S6L0OJYI2/3i3r7/C82ZlpGuIUsc1VfZT7XkuBkangP4RMVLnFjObR9Rux6fneISo\n+19J9HrnE+H1dm6/2djzx4hPAD8m3r+9gAnZe1/Da6zaZnd/vIf7LncI8BN3n5W7vUQctA7JbZd/\nzkXEe3Ni7vE7iDLX33L33QWMovOJyU8SI5QOz/1Nzks9/2y00+3AJ3K/cxTxtzSL+HT5P+nvFTMb\nS7z3D6WfxxHnK6alNr9JnIzMTAD+4u750tSgZz3vMInUn5mdBmzs7od0u3F9nm814hPIu9z96WY8\n50BlZpsBl7r7+Bq2PRDY3N1PrXHffweOSAdeSRTc0u9SaeZ+YoTIXxv4PBOA24jSxrnA+9z9vY16\nPpFGUalE+pWZfY44qXhjI0M7mUiUjWYSdfb9G/x8Ig2hHreISMF0tRZCn5mZjgoiIr3k7hVHLKlU\nIiJSME0JbncfVLczzjij39ug16zXrNfc91t/ZVh31OMWESkYBbeISMEouBugra2tv5vQdHrNg8Ng\nfM2tqKHDAbNRJY18DhGRRskWnmx2huWeV6NKREQGgoaO4xYRKbJWrRaoxy0iUjAKbhGRglFwi4gU\njIJbRKRgFNwiIgWj4BYRqcLMlo6pbiUKbhGRglFwi4gUjIJbRKRgFNwiIgVTl+A2syFm9oCZTa7H\n/kREpLp69biPBR4BWnNiv4hIL9R6RZpm63Nwm9kYYE/gIqD1xs2IiAww9ehx/wA4Aeiow75ERKQb\nfVrW1cw+Dsxy9wfMrK3adpMmTVr6fVtbm66iISJSpr29nfb29pq27dMVcMzsO8DBwGJgJWAYcLW7\nH5Ie1xVwRER6qLsr4NTt0mVmtjNwvLtPyN2n4BYR6aFmX7pMCS0iA0arrlWiiwWLiFShiwWLiEhd\nKLhFRApGwS0iUjAKbhGRgunTBBwRkYGsVQdWqMctIlIwCm4RkYJRcIuIFIyCW0SkYBTcIiIFo+AW\nEamiVdcqUXCLiBSMgltEpGAU3CIiBaPgFhEpGAW3iEjBaK0SEZEqtFaJiIjUhYJbRKRgFNwiIgWj\n4BYRKRgFt4hIwSi4RUSq0FolIiJSFwpuEZGCUXCLiBSMgltEpGAU3CIiBaO1SkREqtBaJSIiUhcK\nbhGRglFwi4gUjIJbRKRgFNwiIgWj4BYRqWLArlViZmPN7C9m9rCZTTOzY+rRMBERqcz6Ok7RzNYB\n1nH3KWa2GnA/8Al3n25mDq07FlJEpCtZb7vZGZZ73ord/T73uN39RXefkr5/HZgOrNfX/YqISGV1\nrXGb2ThgW+Deeu5XRERK6hbcqUzye+DY1PMWEZEGqMtaJWa2AnA18Gt3v6788UmTJi39vq2tjba2\ntno8rYhIQzWztt3e3k57e3tN29bj5KQBlwKvuvtXyh7TyUkRkR7q7uRkPYL7Q8CdwINAtrOT3f0m\nBbeISM81PLi7eXIFt4hIDzV8OKCIiDSXgltEpGAU3CIiVQzYtUpERKS5FNwiIgWj4BYRKRgFt4hI\nwSi4RUQKpi5rlYiIDEStOnlQPW4RkYJRcIuIFIyCW0SkYBTcIiIFo+AWESkYBbeISBVaq0REROpC\nwS0iUjBNmYBzeAt+1Bjs9C8iraa3U128m++97NaRbkuAN4FX01cH1gQ2BzY65xw23HDDXrao8Zpy\n6bKLG/YM0hutORdMpPcdiuz33gbmpdvc3Nc5wGxgFjA/bdeRfm8MsAawIvAE8H5gY+A/wA1pv82e\nQdka15xs2DOIyEC2mOgRzwJeSl/Lv8/fOoC1gdHAykT2zCfCexawOrBJenwuMA04HPh2er7TgKHp\nK5QOCK0W3FqrRESaxoHXWTZwy4M4+3k2MBIYRSmQR6fv35m+rgK8DDwHPAI8BDwMDCfCe366jQRe\nI0J5l9SeHwH/7qK9w4Bn5szp68uuOwW3iPTJIuAVug/il9PX5SiFbz6Q3wl8sOy+tYAh6XmcCOcp\n6XZT+voisBXwbmAh0WPO2rUS8Gdgy/S8B6Y2FJ2CW0Q6caI23F2vOLvNJXqz5UE8mqgVl9+3ag1t\nWAhMB26mFNRTiFDeJt0+TZQ5pgJ/IOrRGwEzgXHAycDBwHt69S60NgW3yCCwkOjxlgdxFsbljw0l\nAjfrGWdfNwF2TD9ngTySUq+4N2YT4ZsP6MeI8N0a2BY4kQjrtYC7geuAM4iTjDumNjhxoDiP6Lkv\nAQ7rQ7tamYJbpICcGCmR7wlXC+GXiLryWpR6xVkQr00Mf8vfN4qoGzeizf+hFM5ZWL9GlDq2Af4L\nOIoobWRteJ3oeZ8I/AkYC3wC+BpwTXpsHLAn8AVghwa0vdUouEVaxFuU6sCVShLlt5XpHMTrpK9b\nAG10Ll2sQXNn271FjNjIAjq7DSN60dsQZYxziPJGedueBy4jyh9/BT4A7A2cSZx4/B7wDLAuMfJj\nAvG+/YjogQ90Cm6RBukgepPd1YqzsH6TZXvEWSBvVfbYaGLccSuYRan3nH39N1G2yEL6E+n7tars\nw4EHgclEWM8A9gAOBX5DlG4uB/YiTjjuCDwJ7EScfFwN+D1wZZ1f2zxgxIgRLXclHAW3SA8soPaT\ndq8Q44YrjaDYpuy+UcAIWntG6xIiUMtDegGlgP4ocBxRfunuwLIQuIsI6uuJXvdE4LtEMK9AlIN+\nStSttyXGW08myiNHAscQoT3YKLhlUFtCaYJHV0Gc3beYykE8BngvnYN4FKWhaUXzOjEeOitxTEk/\nj6YU0p8jwnQDaj/gzCGG8d0A3Eic7JxIhPGWuf08A/wQ+CXRy76ROBAeTvTk3wbOT227prcvssAU\n3DKgOPAGtQfxbKKnW35ybm1g+wr3rU5r94p7yol6cnkv+jlgM0ohfUD6fngvnuMpIpivB/5BlDcm\nAucSNeq8KcDZRLgfltozAjie6GXvSBxALiPC+7u9aM9AoOCWlreYzhM8ugriWel38oGbfT+OWIci\nC+G1iUWFBst/gkXEMLvyoXdQCuiJxMzCTen9+9IB/IsI6uuIf5uPA0cDu7LsOG4n6tT/S4zd/hJR\nHhkG3Eb0sj8KnE6MzZ5K1P3/3Mv2DQSD5W9WWki2fkStQTyXGBWR7/lm32/EshM8BmPNs9xc4mRf\nPqCnE0PpspD+avo+G5nRF28BtxNhPZkI3b2B/yMOlpXGeS8iTiaeQ5Ssjgf+mygvzSdKMTcBFwGP\nAt8gyiLr9LGtA4GCW+piEaXxw+WLAVUa4rY81Wfa7VB235r0bYLHQOZEPTgf0FOJ93grIpjfB/y/\n9HM9D2qvUJqxeBsxQ3FvoJ2oXVczD7iQGLq3ETG072OUDh43A58n1hN5hAj37xHnECYAvyN64M2g\ntUqkUJzSimr5CR6VJndkS2WuybIz7UYTH7vziwM1aoLHQPc2EWT5gJ5KjOfOpoEfQNR930VjDnaP\nEUF9A9Gj34UY6vdzqg/1y8wkwvoXRMnkWiKMM3OJESm3EEtB75K+fhP4FbAf8EXiRGWzgrtVKbgH\nkbfpfoJHfvbdiiwbxKOIk1Y7lT02El1OqZ5eZdkZhk8QCzFlIf3x9HVUA9uxhJhino2vnk/0ek8B\nPkyMqe7ONOJE5PXEpJt/AuWXKLgR+B9i7PY0oqd7PtHTvoyoc59OfLJ7sy8vaIBQcBdYBzG8qpYg\nnkWMt81P4Mj3gjencxCPprb/lNI3HcREkvwJw6lEOWHrdNuZGK+8Jc35N3md6PVOJkoh6xMnLX8N\nbEdtB2gH7iBGiNxPtH8GcYDPm03Utm8neuLZcqsXEGF/OxH2XwCOJerhUofgNrPdiSGXQ4CL3P17\nfW7VIPYmXU9zzgfxy0TNslIQb8my6xe3+gSPgW4BnaeBTyHKDWsR9eGtgc+mr+No7ieYmURITwbu\nJKaYTyAWchrXg/0sBq4mAnY+Ufq4msoHnD8SveyJlKbDA/yA6G3flm7LEycnZ/egHQNdn4LbzIYQ\n7/EuxL/9fWZ2g7tPr0fjBoIllKY917IGxUKWPWE3mjjzvw2dJ3iMprgTPAYyJ9aIzk9emQI8TdT7\ns1Ed+xGBvUY/tfEhSvXqGcDuRO/2cno+XvsNoh79A6KH/nUi+CsdfF4jRrTcRUywyderzwZ+BvyF\nCPJvEZ84fgm8A40YyvS1xz0emOHuTwGY2ZXEieUBHdzZBI9aShSvEX+A5XXi0XSeaZcF9DDUKy6S\nxcDjdO5FTyEO2FmpYw9i/PGm9O+BdhHRm87COlucKT/FvKdeInpuP0v7uJxYUrWayUQv+5PE+5UP\n4rOAS4kSyxiirv1xohy4GTHksNnBPVDXKlkfeDb383PEsM1OKgVRtbehWmi12vYb0Hk88drpvh9V\n2PY14mRTI9uj7ft3+32IkD6amAa+ftq2Vdr/FKUrxOSdl2692f99wG7A/sA9xEiWrra/kAjnK4he\ndLX9jyVGr9xInJDdJe0b4IT09dyyNlZr/5IKz3EFcSWcSr6Vvp5e5fFy2bUhy1UL+p5uX01fg7u1\nDkNNdCSVa8vH9WejpN98jih/jO7vhlTxIhFi9fQcsB4x5vpd3WwLUSoaT9TPuzOb6BCtRoyGGgwn\nytvb22lvb69p274G90ziAJkZS/x7dtKTdO/pkaCZ2/+V+OPLxjPPoPPwupeI8cnlgT6aOHtbft9a\nTW6/tq/v9i9QurTWWcRH/5WIAM/KJNOJSUW1jqluVPvfT4yFvo4YMbIpcVJwIjGiqDf735P4PzAx\n7e94YpRMtV730cQBbjxwCTFCJdv/K8QknI8QU98XEWWo64myyTPEUMiziYPQOcSU9+/S9dT3IWk/\nmQNTu6v5OpWv8l5NT3vKXW3f1tZGW1sbAGeeeWaX++lrcP8T2NjMxhFr1XyGmAMwIH0o3arJrmCd\nhXk+1P8D/J3OswuzBY7Ke+2Vgn80A2+Bo6Jbl86XxiqfxXgVcCrxb70lpUDfhvrPYuyOEf8xDyDK\nJe1EvXkPIqj2JurdH6L2UFiRKF0cS6yZfUK6nUiceC2v569D1NZ/lZ7380RQrkh0Ym4hwvsYonxz\nIbAvsTrg34iLQ0joU3C7+2IzO5roeAwBfjGYR5QYEa6rE72D7mRLimYBnz/Z+U+WnTq+hGXHYVc6\n6Zl935uTTdJ7Rox8eAcRhJm5xAiOB4h/118QMyDHUJpMkwX6ejT+4DyUqE3vRgTkVKJnezzRwdiD\nmA25G6Uhet3t71DgEGJtkXOJy4p9mQjn/D4sbbcrUW58L3FC8r3EzNvbiNEtxxDniz5F/B+Zk9p1\nX69e8cDT53Hc7n4jcR5BemgIpQDeoobt84v453vzzxGrsVVbxL9ajz5/YnU46s03ynCW/bS2mDgB\nl/XOf0gEqFMK8ezrpjTuIGyUDh5nEH9LNxALOx1BjBDZmyiHjKlhX3uk27+IssaGaT9fJk7YZtYl\nDhZXEKWLI4BJxHt1U7rvSKJHvj3xHm1N84Nba5VIn61CTIYYV8O2HUQpJj+FPft+CssOX3yL6rMq\nK5VtWuWyWUW1PHGw3oLSCIds/Hc2tPAPxCiHZ4nwLu+d92Zt7O6MIS7WexQxgeZmImBPJ05gTSR6\n49vQ9YF+O6J88jQxtnsrohRzAlE2Iv3+gcSokSOJ0Ti/JGrgNxM9/jPT7++d7pOg4B6gliM+eq5J\n1yefMm/ROeCzry8SH/PzJ2BfofpJ2EpB3+wL1RaVET3RdYlea+YNSlejmUIE4kPEgTYL8ez2Dur3\nyWl1YpjjPsQnhHuIEN+XGOkxgQjyD1P9QP4O4tPEGcQa27umNp9AnIg04m/mWmLVvwlEKeWbwK1E\nYP85bb8/1YfxDTYKbgFiNMQG6dYdp7RGSvmiVdNYNvzfIE4+1Rr0K9fpNQ0UqxJD6PLD6JbQeY2T\nC9PXNyiNaMnCfAv6/glpeWKCzY5ECeRRoqTyTSJQdyNCfE+WXY8E4uB9CjFc9tfECJNVibr6Pmn/\n+xEHgaOIHvtFxIFiT+KE6vNEgH+nj69lIFBwS48Z8R9xDeIjfHfeJnrplYJ+eoX7slUJawn6kQzO\ntbqHEGteb0L0gDMvUwrz24kyw5PEOOvy2nl3y7BWY8RMxs2Ak4h/uz8Sa5J8kSh5ZEMNy8d3r0is\nx3J4+p2ziUA/hlgzfBQxGuca4sTkwen7TxIllYOIg8Rgp+CWhluRODG1fncb0nkd8PITsY8T61vk\n75tHaR3wWoK+/LJZA80oIuB2yd33FqV1vKcSPeUpRCkkP+Z8W+LCBj0ta61NnFw8gjiBflt6jrOJ\nA2sW4u/P7Xs5oiwygRgmey5Rzz+S6I1/mujdH0uc1P1Buv2UGDI42Cm4paUYMbZ9BF1fRSWTv/JO\nedDPYNmr8CxH12Ply6+8MxD+g6xElB62y93nxDT4LMx/Q4y/fpU4kZgP9K2o/YC3CqVA7iAuDjyZ\nGBY4K/fYrpQupvEBor49gwjnLYhPEV8hRp38njgoHESsY3IGcdKyGQbqWiUi/WoFYuzzejVsm02Q\nqrQ42JPEybd8fT6bIFUe9NWCfzWKM6TSiKF6GxJliMwcSmF+L3HNyEeJcx/lo1rWoevXuxyl2vy3\niff4BmJ89sFEPXsisZDU2kRZ5QJiWOD5xHom7yfq2v8ihhS+Soy2uZpYWXGwUnDLoJGfIFXL2hpL\nKF1dPt9rf4mYqFJes++g6zHz5csdtOIEqRHEDMW23H2LiPDOxpyfQwT7ckSAb0spzDeheqhsRPSi\nv0IsvPYnojd+PHGuJBsvvhnRoz6JmJxzBPF+fY2YRXkC0WN/lSilDEYKbpEqhlBaAXKrGrbPlvvN\nB/zLxDT4f+Z+nkWEzjCqB3v5z/253O8KxOvfiugpQ3x6mUnpROg1RE95JjH8ND+q5T0sOwNzJFH6\nOIgYMXIHMYJkd2Im5kQiyD9HlFmuIWrgc4jgv4kI/43q+1ILQ8EtUierUio/dKeDCJ5KI23+ReeQ\nf5k4wVitLl8p6Bu97rcRk3XGEL3gzOuUpvdPIdYlmZbat03ZbWzaz1CiB70r8GNKJ1C/Qkzg2ZMI\n8tvTPs8mTrbuQJy4PIgYZtjRoNfaihTcIv1gOeLjf61D8t5i2Z589v2DdA7+8glSlWry+fvWoH69\n+dWIafL5iyksIU48ZmH+0/T1bTr3zLcmeuvZz6cTU/AnU5qC/19EiH+FuOIOxMWMf55eQ7Ze90Cn\n4BYpgJWIHurY7jYkyhiz6bxAWRb00yrct4A4gNQ62qana2MPIS7i8G5isk7mJUqllluIK7r/h6iT\n53vmnyEuFjwvbXcdsYbJO4nx4C8TJz7biBOeTv0ORFqrRESawoga8khqnyBVPts1+/oInYN+FhHc\nXZVt8kE/kurjwtemtEphZgHwMKWRLVcTnyhGUAryTxHB/TzRG78rPT6XCPLViKA/pIbXXlQKbpFB\nbkVK9eru5CdIldfnHyOuaZm/5up8ShOkaunRrwK8L90yHcC/KfXOL0nfzyXKK3ul332BKBM9TVxg\n4avU/6o/rULBLSI168sEqfKVKp+gc23+JSKQuhpO+WGidDI67f9BIsSz+vmTxMJWqxOB/wqx2uAO\nxEHhPjqfTC0qBbeINExPJ0jNp9STf5FS6D9OXAUnH/RziVJMFupbEZN2jBiaOYpST/1JSitlLiau\nAPQacTnCWtrWahTcItISjDgZOIy4Tmd3FtO5VJPv2b+RtlmFCPdXiLLJPGKGbAdRqx9JnOS8khjO\nOY2Ycv8wcaJ0Yd9fVkNYI+fgm5mDjg4i0hhdjR6x3ONWdj/ESdNsGyfCfFG6f1y6ZRdvaPZaJWaW\nPW/Fl9iUTF2wsFWPWyJSVF2FqbsvfTz/taOjA3dnyZIlLFmyhIULF7JgwQKWLFmCmbHWWmsxYsSI\npcGZfW01TQnuFVZoxVUZRESKSVeUEhEpGAW3iEjBKLhFRApGAz5ERKpotSvfZNTjFhEpGAW3iEjB\nKLhFRApGwS0iUjAKbhGRglFwi4hUYWYtOe1dwS0iUjAKbhGRglFwi4gUjIJbRKRgFNwiIgXTp+A2\ns7PNbLqUEZwMAAAGOUlEQVSZTTWza8xseL0aJiLS3/IXZGglfe1x3wJs4e5bE9fzPLnvTRIRka70\nKbjd/VZ370g/3guM6XuTRESkK/WscR8B/KmO+xMRkQq6XY/bzG4F1qnw0CnuPjltcyqw0N2vqHP7\nRESkTLfB7e67dvW4mR0G7Al8tNo2kyZNWvp9W1sbbW1ttbZPRGRQaG9vp729vaZtrS9nTM1sd+Bc\nYGd3f6XC4w6texUJEZGuZOuUNDvDcs9bcaGUvgb3E8BQ4LV01z3uflTucQW3iBTWgAzuGp5cwS0i\nhdWqwa2ZkyIiBaPgFhEpGAW3iEjBdDscUERksGrV83PqcYuIFIyCW0SkYBTcIiIFo+AWESkYBbeI\nSMEouEVEqjCzpbMYW4mCW0SkYBTcIiIFo+AWESkYBbeISMEouEVECkZrlYiIVKG1SkREpC4U3CIi\nBaPgFhEpGAW3iEjBKLhFRApGwS0iUoXWKhERkbpQcIuIFIyCW0SkYBTcIiIFo+AWESkYrVUiIlKF\n1ioREZG6UHCLiBSMgltEpGAU3CIiBaPgFhEpGAW3iEgVWqtERETqQsEtIlIwfQ5uMzvOzDrMbGQ9\nGiQiIl3rU3Cb2VhgV+Dp+jRHRES609ce9/eBE+vREBERqU2vg9vM9gaec/cH69geEZGW4e4tuV5J\nl4tMmdmtwDoVHjoVOBnYLb95tf1MmjRp6fdtbW20tbX1pI0iIgNee3s77e3tNW1rvTmamNmWwG3A\ngnTXGGAmMN7dZ+W2c2jdFbZERFpRNnbc3St2iHsV3BWe5D/Ae939tbL7FdwiIj3UXXDXaxy3kllE\npEnq0uOuunP1uEVEeqxZPW4RkQFHa5WIiEhdKLhFRApGwS0iUjAKbhGRglFwi4gUTJdT3kVEBrNW\nHcqsHreISMEouEVECkbBLSJSMApuEZGCUXCLiBSMgltEpAqtVSIiInWh4G6AWi8/NJDoNQ8Og/E1\ntyIFdwMMxj9uvebBYTC+5lak4BYRKZimXAFHRER6TlfAEREZIBra4xYRkfpTj1tEpGAU3CIiBdOw\n4Dazfc3sYTNbYmbblT12spk9YWaPmtlujWpDfzKz8Wb2DzN7wMzuM7P39XebGs3MvmRm081smpl9\nr7/b0yxmdpyZdZjZyP5uS6OZ2dnp33iqmV1jZsP7u02NYma7p4x6wsxO6u/25DWyx/0Q8Engzvyd\nZrY58Blgc2B34CdmNhB7/v8LnObu2wKnp58HLDP7MDAReI+7bwmc089NagozGwvsCjzd321pkluA\nLdx9a+Bx4OR+bk9DmNkQ4HwiozYHDjCzzfq3VSUNC0x3f9TdH6/w0N7Ab9x9kbs/BcwAxjeqHf3o\nBSDrjYwAZvZjW5rhC8BZ7r4IwN1f7uf2NMv3gRP7uxHN4u63untH+vFeYEx/tqeBxgMz3P2p9Dd9\nJZFdLaE/errrAc/lfn4OWL8f2tFoXwPONbNngLMZoD2TnI2Bnczs72bWbmbb93eDGs3M9gaec/cH\n+7st/eQI4E/93YgGWR94NvdzS+VUn645aWa3AutUeOgUd5/cg10VckxiF6//VOAY4Bh3v9bM9gUu\nJj5SF1Y3r3d5YA13/0Cq518FvLOZ7WuEbl7zyUD+HE3rLSPXC7X8vzazU4GF7n5FUxvXPC2dSX0K\nbnfvTRDNBMbmfh5DQcsIXb1+M/u1u++Sfvw9cFFzWtU43bzeLwDXpO3uSyfr1nT3V5vWwAao9prN\nbEtgQ2BqWvZzDHC/mY1391lNbGLddff/2swOA/YEPtqUBvWP8pwaS+dKQb9qVqkk3xO5AdjfzIaa\n2YbER+x/NKkdzTTDzHZO33+EOJEzkF1HvE7MbBNgaNFDuyvuPs3d13b3Dd19Q+I/9XZFD+3umNnu\nwAnA3u7+Vn+3p4H+CWxsZuPMbCgxoOKGfm7TUn3qcXfFzD4JnAesBfzRzB5w9z3c/REzuwp4BFgM\nHOUDc/rm54ELzGxF4M3080B2MXCxmT0ELAQO6ef2NNtA/Buu5MfAUODW9EnjHnc/qn+bVH/uvtjM\njgZuBoYAv3D36f3crKU05V1EpGAG4vhpEZEBTcEtIlIwCm4RkYJRcIuIFIyCW0SkYBTcIiIFo+AW\nESkYBbeISMH8f6G21xYyd8NSAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x109e22190>"
]
}
],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"SBDF.plot_imex_stability_region(N=200);"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvXe4XNV57/9ZiF4kkBBCQgIVJCFUsISQUD9qqBxJFINx\n4m7fJDfNN4nzc42vjROXGydOYoPtYMem2AZj0xGm61CFJJoqIAkhUO8V1XPO+v3x3ct7zT5z+pQ9\nM+v7POvZe2b2mVkzZ+az3/2utxhrLUFBQUFBpaMTij2BoKCgoKDWKYA7KCgoqMQUwB0UFBRUYgrg\nDgoKCioxBXAHBQUFlZgCuIOCgoJKTAHcQS2SMabGGPO5Rh670BhzwBhjkscaYz5mjHm8kHNti4wx\nXzHG/CxPz/0vxpgdxpjN+Xj+YssYc6kxZkmennuRMebSfDx3KSuAOw8yxqw3xkyJ9j9tjKk3xvwg\ncczV0f2/jG73jm4fSIwbosfvNcbcmniO+40xP2xkDmcbY35hjNlijNlvjHnbGPMl7/F6Y0zfVrwt\nG42GD1j7vrX2LBsnBfzxWGvtr621M9rxun+UMabK+4z2G2NWG2P+vC3PlZS19rvW2j/LxXP5MsZc\nCPwDcIm1tkcOn9cYY9YZY1ZmeazGGHM4+px2GmMeNMb0jB77b2PMj71jTzLGfNDIfaOMMQOiv99u\njNlljHnMGDMg8ZL/DHzf+/u/Mca8Yow54r7fzbyXv4++p/uMMf9jjDnZe/jfgG+1/JOpDAVw50c2\nsf8OcIMxpoN3/6eA1TSEYacIgm78Lrr/r4DrjDFVAMaYG4EPAV9uZA7/AZyOgNERmAesTRxjWve2\ncqb2vO6m6HPpCPwf4MfGmME5mlc+dCGwy1q7q7V/aIw5sYmHJwKnAF2NMSMTj1ngr621ZwH9gFMB\nZzg8G/2t00jgPWBC4j4LvAp0Ah4ABgDdgMXAg94cuwNV0TFOmxDMf9GC9zgD+BIwBbgI6Avc5B3y\nMDDZGNOtueeqJAVw518G2AosB2YAGGM6A2OAh2ghxKy124AvAD+LrLj/Av7CWnuokT8ZCdxlrd0X\n/f3b1tp7o9d/LjpmqbPqIwv9kciy2m2MedgYc0HiOS+OLl33GWMeMMacEz2fu1po8H2Krjiez/K6\n+40xHzHGLDfGzPGOPymyEi9rwWfyB2AXMCj6W2OM+bIxZm30HL91c4we/6Qx5r3osX9KXBl90xhz\np3fsPGPMSmPMHmPMAmPMJd5j640xXzDGLDXG7DXG3G2MOSXLe58GPAH0iD7nX7Twub9ojFkGHMj2\nmUb6FHAvguinmviM9kXHuJPb88Cg6DsIMB64GzjDGNMlum8C8JK1ts5au8Ra+0tr7V5rbS3wn8BA\n73OdDrxqrT3mveb91toH0f+mOX0K+Lm19k1r7V5kXX/ae64j6AQyI/ufV6YCuAunO4FPRvsfRT+m\no1mOaxTk1trbkfX+GvAHa+0TTbzey8C3I3D2TzyPs7iGeVb9CcD/IAvxQuAwcHNiXp8EPgN0B2qB\nrG6aJubvv25Ha+09wB3Ax73DZiOremlTz2WMOcEYMw9ZhK9Hd38eXVlMjOa4B7glOv7SaP9Posc6\nAb7r4o9XPpEr4DfR850LPAo87FnAFrgBwaQPMAwPNt77fQqYBWyOPufPtuC5Qd+PWcDZ1tr6LO/9\ndODDwG+Be4CPGmNOSh4WHdsFuA5YFM1pA5kW9kQE85cS9z1Hdk0Etlhr90S3hwJvN3JsS4ySSwH/\nf70M6OafcIE3gWZP5JWkAO7C6X6gyhjTEfgEcHsjx+2MLDE3BiYefwHoDPyqmdf7W+DXwN8AK40x\na4wxMxs72Fq7O7KUjlhrDwLfASb5hwB3WGtXRVb+14GPGGPa6275NVBtjDkzuv0JdJJrTD2MMXuA\nQ+gz/YS19p3osb8A/slau9laexxdcl8fuaiuBx6y1r4UPfZ/yXRT+e/jRuARa+3T1to65Gc9DRjr\nHfNDa+3WCGAPI7dVNiU/n+ae20bPvclam+3EDgLxfmvti8Az0X3Vidf8oTFmL7ADOBP4a+/xZ4FJ\n0f9uFLAQwXtidN/Y6JjMNyI/+c3IZ+/UCTjYyDxbUgjpTGCfd3t/tD3Lu+8AcHYLnqtiFMBdIEWX\nfPMR8DpbaxeS3SLpYq09xxt/tGYiy/kLyHL8QVM+0AjA37XWjgS6IMvsd8aYrD8AY8zpRgtX640x\n+9APt1MCzBu8/feBk5DV2GZZazcDLyLAng3MRDBvTJuttecAHZG76KveHHsD97uTHrAKXRl0Q1b2\nRu91D9P4pXwP9P7csRa9d991tNXbP4wA1BJ1b8Fzb0j+UUKfAu6L/r4O+Zd9d4kF/tZaeza6GrgI\nXck4PYcs56HAuui7+aJ332lEFrqTMaYrcvvcYq39rffQHjIhm/FnzbwPEPQ7erc7RdsD3n0do9cJ\nihTAXVjdgayV5qzlBorg9HO06Ph54AO0qNOsrLUHgO8CZ6BL+2z6AlqAGmWt7YSsbUPmj+/CxP5x\nYGfL30Wjuh25S25AvtUtzf1B5FP9EvqhOxfU+8DMxInv9OjksAXo6f7eGHMaOqFl0yYEO3esAXpF\n92edTnPz9bS5Bc/d6PNFVu8U4FNGkRhbgI8Asz2/NUT/N2vtCmQsfM/zlz+PXA/V0T7Aymge1cBi\n32cduS2eAB6w1n43MaVl6HuTTS35XFaSebVyGbDNc8WA1jCadJ1VmgK4Cyhr7bPANOBHTRzWmJXy\nl8hF8p3ISvsc8MUsrhQ9iTFfN8aMNMacbIw5FUVg7CH2R25DEQdOZyLLcV8EgG9kmdfHjTGDIh/r\nt4DfeSGALVXydUEujxHohHRHS58ocnn8O/DF6K6fAt8xWrzFGNM18oMD/B6Ya4wZYxRu9k0a/6x/\nh9w3UyLf8ReAI8gPnE2tcRfd08rnTuoTwFsIlpdFYwC6mvjTRv7mdhRhdAOAtXYtsB19J56L7rPI\nyv7jfQCRa+9x4AVr7VezPPdTwAjjhfAZYzpE37kTgQ7GmFOMF1FltJDt1jvuAD4Xfa/OQSeZX3rH\nnoq+G082/9FUjgK486+M+Gdr7YJo9bzBY5H2msw47r+LQPRt4HPRyj7W2jcRtG4lu+rRD2AHsuam\nAtVeFMo3gdsjt8L1KFrgNGRBvwT8ITE3i35ktyHr9WQEWv/xZt9/ltd1bqT7kKvjvkaep7HX+QVw\nXgTo/0KROk8YY/Yj3+2o6DVWIb//3cjqPYDgddR7Xhd7/ja6AvgR+vyqgbnus2/Be2x0ztba1a18\n7qQ+CfzYWrvdG9vQSeuT3nH+ax5Hn80XvcefRW6uF737nge6krkweS2KUPqM953cH1n+LtrpGeAa\n72++jtYgvhS918PA1wCMMb3QZ788+vvHgX8FFgDr0eK7bzTMBRZYa33XVMXLtN5gCgrKvYwxXwf6\nW2s/2ezBuXm9M9EVyMXW2vcK8ZrlKmPMIOB2a+2oFhz7MeBSa+3XWvjcLwOfjU68QZECuIOKrsg1\n8yqKEHkhj68zF3gauTb+HbjCWnt5vl4vKChfCq6SoKLKGPNnaFHxD/mEdqR5yG20CfnZP5rn1wsK\nyouCxR0UFBRUYmqqFkK7ZYwJZ4WgoKCgNspamzViKbhKgoKCgkpMBQG3tbaixje+8Y2izyG85/Ce\nw3tu/ygWw5pTsLiDgoKCSkwFAfe3jOHVV1+lvr5BobOgoKCgoFaqIODeC3xs5Egu6NCBzxrDvffe\ny759+5r9u1JVVVVVsadQcIX3XBmqxPecRuU1HNBFlbhXeAeVx3sU5dmORPm+1atWcckll9D+CqFB\nQUFBuZNjUj452czrZoViQcHt6wNU4MCB/ERUd7L60UepqqritNNOy9u8goKCglqiAO4mZIEVxBB/\nA7XiqL7lFqqrq7noooua+vOgoKCgslJJgDupPaj473xUoq4bkTVeU8PYsWM56aRkl6agoKCg8lFJ\ngttXHfAKsTX+DupOWn3bbcyaNYvzzjuv3fMMCgoKSpNKHtxJbUFW+KOogvsAtMA5e/FiLr/8ck44\nIYSmBwUFlbbKDty+jqHOuY8ii3wPao09+557uOqqq+jUqVNTfx4UFBSUSpU1uJNaR+xSeYE43HD2\nypUMGjQohBsGBQWVhCoK3L4+QL2Q5kejA9EC5/z5TJ48OYQbBgUFNasQDlhEWdRK2kH8DWAiUP3j\nHzN79uwQbhgUFJRVAdwpUlPhhuPGjePEE/NapjwoKKhEFMCdUvnhhvOBd4nCDW+/nZkzZ4Zww6Cg\nClYAd4nIhRvOR11lLyGyxl95heHDh4dww6CgClJZg9sY0wEZrhuttXO9+y3AF4AqYDxwdrtfrXA6\nBjxPbI3vR+GG1b//PdOnT6djx47FnF5QUFCeVe7g/gfgcuAsa+08734LcBPwLLAYGIggPgnVIykl\nkK8ljhlfCFxBVN3wrbcYMGBACDcMCgrKifIObmNMT+A24NvAP2SzuN0rHAWWADXRWIQyHychmE8A\nzmnXbAqng8iV4kB+KpFL5bHHmDRpEqeeemoxpxcUFFTCKgS4fwd8B+gI/GNT4E7qGJkgfxnoTwzy\niZQGyC2wjNilshzNv/qnP6W6upqePXsWcXZBQUGlpryC2xgzB5hlrf1rY0wV8IVs4P6G9zdV0cim\nY8hRXoNcKy8BF3t/MwHo3ObZFk67gMeBR6JtT2AOUP3ii4wePZoOHToUc3pBQUEpVE1NDTU1NQDc\ndNNNQP7A/R3gE0At8hZ0BO611n4yerxdUSUO5M8imC8E+hKDfCLpB3ktcgk5a3wTMAOo/vWvmTlz\nJp07p/0dBAUFFVoFCwc0xkyila6S1uo4sUVeQybIJyGQd8nRa+VLG5Bf/FH0HoYSLXAuXcrQoUPD\nAmdQUFDBwf2FbFEl+QqkOQ68Sgzyl4A+ZFrkaQb5EXQ1MR+5VWqJFjgfeoipU6dy+umnF3N6QUEV\nr7IOB2zixQuagONA7lwrLwEXkQnycws0l9bKAm8TQ/wVFPc++4c/ZM6cOfTp06eY0wsKqkgFcBdB\nx4HXiBc7X6R0QL4P1VNxbpUuRC6VBQsYN25caN8WFFQABXCnQLXEFvkCYovchR9OIp0gr0cWuIsZ\nX0to3xYUVAgFcKdQtcQWeQ2ZFvkk0gvyZD2VgUTWeKinEhSUUwVwl4BqgdeRNf4s6qJzIZmula5F\nmltjStZTOYDqqcy57z6mTZvGWWedVczpBQWVtAK4S1DZLHIf5Gm0yNeQWU/lSqLknzVruPjii4s5\ntaCgoBYqgDuH8i3yGtIP8gPAU8R9OM8icqk89RQTJkzg5JNPLub0goKCGlEqwD0TGIWq6V2BOs6U\ng5rykVeRrqiVetSyzYUbvg1MBeb84hfMnj2bbt3K5b8SFFT6SgW470fFpNw4i0yQXw50ytssCicH\nche14lvkk0kXyLcTL3A+iWrCVANzlixhxIgRYYEzKKiISgW4/VewKJzNB/kbqBDTFd74EFDqfdh9\ni9wtdvoW+STSkdl5HM3tEQTyfYSGEUFBxVTqwJ1NtcAqYpAvBt5CbcMcyEcBlwKl3MbXgdz5yF8C\nepO+FP13iF0qC4HRRAucq1fTv3//Yk4tKKgiVBLgzqbDwFIE8Vei7UZkifuW+cVAqZZl8jM7XUJQ\n2qofHiRe4JxPWOAMqiyFcMAcaB/KfPQt84PASDIt8x45er1CK1vRrH5kZnYWs7GERVE1zqXyNjAN\nqA4LnEFlqgDuPGkbArjvMz+FTJCPpDQ66STlyti6xc6FCORVpKND0Da0wPkIsspDBmdQuSmAu0Cy\nwHoyYf46CkF0IL8CGA6UWtHU48St3p5FIPc7BE2keM2Xj5G5wHkAlaidc//9TJs2jTPPPLNIMwsK\narsCuIuoOrTY6WC+GC2GDiDTMh8MlFLNPdchaEE0XPPlKhR+OIHihVmuIfaLvwyMRQucc9atCyVq\ng0pGAdwp01G0+OlAvgR4H7iMhoufpXLBfwy9lxoE8sXIfTGZuGdnMQL79iNXyiMog7MzUZRKTQ1j\nx44NJWqDUqsA7hLQfhTl4Vvm+5GP3LlYSmnx8yh6D84iX4JCKqsQzMejKJFCqh4twD6MrPF3iXpw\n3nkns2bNokuXNAREBgUVVwHc7dQ2Mhc+lwAnU5qLn0fIBPkrwBBiH/l4oNCe6M3ICn8EeAYYRmSN\nL1vGkCFDQg/OoIpUAHeO5RY/fRfLa0B3Gi5+pj3z8zDyP7uolVcROKuiMQ44o4DzcT04H4lGPZFf\n/NFHmTx5MqeeemoBZxMUVDwFcBdAdWRmfi6Jbg8khnkpZH4eRpEqLrPzdZTwVBWNsRQuEseiz9BF\nqbwRzWHOf/831dXVXHDBBQWaSVBQ4ZUKcO+i+BmAhdYRBBvfMt+ILHHfzdKX9GZ+HkKFsmoQzJcB\nI4h95GOAQtnAu4DHEcgfR6UCqoE5ixYxcuTIEDMeVFZKBbjPQpfcg71xabQtBd9wrrQX+ZV9y/wQ\nmQufaS57e5BMkK9A/v3J0RiNkp/yrVqUVepcKruJYsbvvZfp06eHrj9BJa9UgLse2ACs9MaqaJxF\nJtAd1IuVSFJobUEAXxRtX0GfiQ/ykRQ++qMlOoCSbtxi51tozi78cBRayM231hG7VPyuP3PeeYe+\nffsWYAZB5aoQDphFFsVOr0LWmwP6myhxxFnlPtDLoW53U/LL3jqYL0XlYEd5YxiFgWJrtA/1v3Qg\nX4MgOhmYguqu5zti23X9eZjMmPG5zz3HmDFjOPHENK8yBKVNAdytkLPQHczdeBO5VobQEOhptEhz\npePos3AwX4xKrw4lE+b9SVey0B7gOWKQv4siVZxrZQTQIY+v78eMP4yMhJnA3LvuYsaMGZxzTiU5\n6oLaogDuHKgeheKtIhPob6E61oOJoX5pNMq1QsZBFIboFj8XI1D6yUKjSVey0E4U7leDQL4RZXM6\nkF9Gfk88G4ldKs+iE8ccYM6bbzJw4MAQMx7UQAHceVQdArrvP3dAP4+GPvRBFDY+uVDaQbzo6dws\nJ5O5+DmS9KwfbCcuYftMdHsicqtMRv+rfIH8UPSaLoPzNCKIhzrjQZ4CuIugOrRw5cN8Jaoj3YPM\n6BYH9LQnzbRGyUqJi1Bsdk9i94prE1eIaJDmtIUY4jUoCmcSsY/8EvITOmlR6ObDyCJfA1wFzA1p\n+BWvAO4UqRb5iH2Yr0CLgj0RzH0/+kAKF6+cb9Wi9+u7WFaj9zma2DIfSPH95RuI/eMLUO2VKmLX\nSr66H21FC5sPE6fhzwXmrFzJoEGDgkslKO8K4G6FjiN4+/7zFchqv4hMd8sQVEK1HC6oD5FZXGsR\nSnhx/nI3ip2r+C6ZIDfEIJ+CknJyrSPRazlr/CSiKJUnn2TixInBpRKUFwVw50DHkFW6MjHeA/qQ\nCfPByBIs9UKlO4ktcmedO3/5KGSdj6Q4ZWJB36k1xAudC5Cby0F8Mrk/0ViUPeog/jYwHZhz++3M\nnj2bc889N8evGFSpCuDOo46iH28ybHEjgrcP88Go7Vg+w9/yKecvd+GIrrPQhaQjvtyicFHnH18A\nnEsM8iq0UJ1LbUMLm48ATxNXNpwbXCpB7VQAdxF0CEW0+IlFK5HvdAAxzN22N8X3J7dFzl/uw3wt\n2ePLC42wemQdP4Mg/jxav/BBnsso7iPohOFixk8mgvhTTzFx4sTQLCKoVQrgTpEOIqvQ95+vRLU2\nLqEh0HuR3gJUjcnFlzuYL0bZjG7Rc3S0zbX125xq0RXCM9Fw/TqnRGMCuUvisijb9RHgIeTSmQHM\n/dWvmDVrFp07V1rJtaDWKoC7BLSfhguiKxEEk/7zwaj2dykBfSuZWZ+LUSy5D/IRFDa2/lg0Jwfy\nJehKwVnk48hdaOgW5FJ5CFnlw4F5wLzVq+nfv3+OXiUoHyrbcEBjTC/gDmREWeBWa+0Po8csKCa2\njzd6o3KmXSgtABVau2kI8xUoPj0J9CFA1+JMs9WqR1aoD/PlyI3kYD4axdUXak3gCKo4uAD5q5eh\nxVdnkeeqYNZhdKJ4CLlUOqJQw3nPPcfYsWPp0KFUV0HKU+UM7vOB8621bxhjzkTlIa6x1r7pwD0f\nLWytQyFdbtQSQ7xPllGu6ert1XYawnwlimRxEHdjMOnJlGxKR1ESzGJimG9Fhal8y/wCCnOyP4j8\n4s5Hvho1knAgH077TyrJWiobUXnaeb//PVdddVUoT5sClS24s7zgA8CPrLVPN+cq2UsmyP2xHl06\n90Fg7+3t90HRDCGCNpZF/RuTMF+FrDrfMh9CadRx2U0cirgoGsUKSdyNCmY518omlJ4/FYF8MO0/\nobyP/OIPIh/8WGDeLbcwd+5cevXq1c5nD2qLKgLcxpjeqH7PYGvtwfb4uC2yuJIwd1b7ZuB8GrfY\nu1OakRq5Vj0Cgg/z5SiMsRsNgX4J6c0STYYkLkJWeh8yXSxDyH+LuK3E8eNPIwvd+cen0v7ORvuB\nJxDE/4AMlXnAvFdfZfjw4SHUsEAqe3BHbpIa4F+stQ9E9+VtcfI4Sol2QH+XTFfMPpTt6Fvpfbzb\npeA+yKfqUNr/CmSVLyfOEr2QTFfLEBTSl8aAtuPIH+3D3LWI82Ge7wid95Al/nS0PYkY4lNoX5XG\nWtR56GEE8iNEEH/sMaqqqjjllDRUmilPlTW4jTEnoau8P1hr/9O73wJ8wzu2Khr51gdkB7obHcgE\nub9/Eem1OvOtY8gadxa6s9I3IngPSYzepO/KZh+ZvvJFaI6jvXEF+avhbtFn6EBeg1buHciraHsP\nVotyBB6Kxkq0+D/vzjuprq4ONcZLWDU1NdTU1ABw0003AfldnDTA7cAua+3fJx5LZTigRT5LB/V1\nZLphNqCIl2xQ74Osp0pb+3dJRc4ydy6XPSj6Iwn0HqQnYsgii3iRN5aik85o1KVnFLq6yIeLpQ5d\nFTwdjRfRSdCBfAJtD4XcTuwXX4AWc68Grl63jj59+rRv4kFFUyGiSsajdZtlxIz+irX2sbSCuznV\nIR/6O8QWum+170buhGxQ74My8tICrXxrHzHEfSv9OA1hPpS2W5q5lu9icWMTAp+Deb4aURxDVwIO\n5K+hOPap0Whr6OEh1LbtQQTzbgji17zyCiNGjAh+8RJSSMDJgw4jC87B3Af8uwja2aBeSW6Y7WRa\n5m7/DBoCfTDpaGyxB0WxvEwM89PIBPnlwOk5ft0PUNPlpxDI1wLjiUE+jNa7o+rQ/B8EHoheYx5w\nzeOPU1VVFaoaplwB3AWWRQDIBvV1yA3TlUyY+5Av52gYi3zlDuQriCNcutPQOi922VyL/n+LiGG+\nEtUqdzC/ktz3+txFHK3yNPo+ObfKNPR9aa3eQgB/MNqfAVxz113Mnj2bjh2LVeMxqDEFcKdMdcTR\nML4Lxvna9xOHOGaDezmmZLjGFj7MV6Crmn4I4kO8bW+Kd3I7gmqe+DDfh9wbDuSjyK1L6H1iiD+N\nrtimRWMKrc+Y3YoWNh9ESUZjgWt+8hPmzZtHjx5p6lJauQrgLjF9QCbUfct9PbpM92Huj57kP365\nkDqCinL5MHcLoi5M0Yd6t+JMk63ErpWFKBuyO4L4mGibq9hyi8I3n0aulefQicyBvLULnQeAxxDE\nH0VXOdcC17z1FgMHDszBjEtbZR0O2MSLB3DnUBbVgPatdAf1dciv3JPsUO9LbsuYFlN7yfSdL4/G\niTSE+RAKnyFaF83LWeQLkYtoBJkwPz8Hr1WL/PJPIpi/irJJHchH0vITxjGUPXc/AnlHBPFrFy9m\n5MiRFbm4GcAdlHcdJV40db71dd5wses+zPtF2wtJZ4JNS2VRFb4kzN9EgPRh7vznhXy/e1EkyUJi\nN0tHMkE+nPb79F2NlacQzDeguHEH8gG0LOKpHngFQfx+dCV4DXDdM88wYcIETjyxnK7tGldFg3sk\n+vE0NdIQVVDOsmjRKwlzNzajy/sk0N3oTGmGONahKI0k0DcgiA1NjJ4U5n1aVLhqIYL4y6hi4mXE\nIB8Tzac92oYs8SejcQJK2JmGFjtb6h9fRQzx91BFw2sfeojp06dz6qnlGydV0eB+GX2BtiB/4Jbo\ntr9/IgJHc4A/j/Ly46ZFx9EP0rlenPvFgR4aul5K2Vo/hKzxZcQwX46uWpIwHwJ0KsCcDiK3h4P5\nQrQQ6YN8BNDWBHeXdems8WfR/3B6NMbTslDV91CEyn0okWkGcN3ddzN79uyyq2hY0eBu7hUsWpnf\nSgx4H/T+/i5k/fkwzwb87uhStBStxLTJzzR1QPfBvhklqjiQ+9Z6P0qrLswOGsJ8Fepf6UA+jMK4\nW1w44suoVvjLKHRyGIL4GBQR0tamyMeR2+ZJVNBqRfR8V0VjCM3/frYjf/h90RyrgOtuu4158+aV\nRfp9AHeOVIs6kDuQ+1Z78r5aGoK9e5b7ulF6FmOadIxM37q/XYf8tkkr3W17kv7yAfXofSxHUHdg\n34jgPYwY5sPQdypfBsNB5HteiEC5ELkZfZB/iLZ9n/ei+PEnonEYWeJXRdvm2s3tRRmbv0d1WsYA\n1996K9dccw1du5ZKm490qOzA3Rp9QEOr3W19t81OdCncEsh3IljxrZFFn6/vfvG3u1A2qQ90N/qQ\n+yzFXOoQim7xYb4MfT+SMB9M7lqh+bLIN76QGObrUIanA/kY2tYd6R3gcWSRL0D/nxkI5ONoeiH1\nIAovvBeFG14O3PDjH3PdddfRrVuxAjdLRxUN7paqntiK98G+Jct9vhWfDezu/uCLb5kOo/j0d7KM\n9cgt5kDug/1i0tn6zkW3JGG+GsVbO5BfFm0vJPfvYR9ygTiL/GV0VTk2GuNQYbDWJDEdR1Exj0fj\nLdRIYgYwE/0/GtPh6G9+h7phfQi4Afjwli2cf34ugiLLTwHcOdYHxJa6b7Un73O+eB/uSdi7kWar\nspiqR4WfskH9nejxfo2MtLlgjiHYOZgvjfYPIZhfRgzzIeQ2yqoO+elfRDB/Ca1ZXIkgPg5le7bm\ne7gTLXI6kJ+GAD4TNZRoLHb+CHLD3IMgfhnwkVtu4frrr+e885pzxlSOAriLpFq00NUc4Lcgf2R3\nmoZ7dyot5SqgAAAgAElEQVSr6mBLtJvY5bI2sd2FLFzfQvddMGlpPeAWQ5cRw/wt1PjBWeZu5LIZ\nxDYE8BejsQy5c8Z5o3sLn8uik9HjqFvPElTLZSYwC7XJyzbvI9HfOIiPBG689Vauu+46unTp0qb3\nVS4K4E65XETNFhSd4YM9CfojZAd6cnQlXdZmMXSIOAomCfUN6HNyQL+YTLgXO6fgOHKtLE2MozSE\n+WByU23yMFr0fBFVKnwJGQrjvXEJLTtxHECLk48hkFsE8Nmotko2a/ww8onfE/3deOCjd9zBNddc\nU3Yhhi1RAHcZ6TCZUN9MdsDvRfBOAr1HYr8blemHr0VRMEmor0Ww70QM9OQoRDx3Y9pOQ5ivQSeb\ny1Dm5WXIh9zeGI56ZPm/gGD+PCqANg5BdQKKKW8u09OiePk/IDAvRi6a6mj0z/I3B1GI4d2oFssM\n4E/vv59Zs2YVvE1bCAcMKpiO09By35xlfyda4MsG9x7e/vlUTrhkPfp81ibGGgT408gO9P4Up0nE\nUeS/fgOB3G1PRwD3Rz/aV1VxMwL4C9FYi9wbE9BC5Riav1o5gHzjjyL3yJnAHATxCTQ8EexGi5q/\nQXHmHwY+/uyzjB8/nhNOyH+NyADuoNTJ98NvJgZ6EvI7UBJNNqj38Ea5x8O7Il9JmK+Jxklkgtxt\nCw1116rtjcTYjVwtw4lhPoS2+/v3IZfK88gyfgO5biZGYzxNFzazqETuIwjiq1G8+FzkVkl6ud8H\n7gLuRFb5x4FPrV5N//7Z7PbcKIA7qGRVRwz4TTQEvIP+DgSoHs2M8yg/H7yLV1/jDR/wJxJD3Af6\nAArnftlNbJW/Ho13ojl8CLk+HNTb4lU+jFwhz6N0+peRlT8pGhNRBmpj2oos8YdQ3PiHiPpnRs/j\nZKP3cQeyxPsCn7n1Vm688cacN4WoaHA/jSy2c6JtJ8q3y0slq4445n0TmVD3x24E7wvIDvYLolEu\nUTQWndR8oK/2bp9OJsjd9mLyv1B6GLkg3kC9L1+LbvdEIHcwH0HrrxqOozKzz0bjRRS3XhWNSTQO\n8sNogfNBBPKuqDrhtdF83PfiOIpM+SXizNXAnz3/POPGjctJGdqKBncVWjDbE20PoDO6D/NzWrGf\nllCuoLbpGLEP3lnwviW/KRpHiEGe3F7g3c5HRmKh5BJ21hDD3G3XIbANyDJ6kz+3VC2qifIqArmz\nzrsggF/ubVuzCFobPd+zQA3yk1+EIk0mI5Bnq2tTj6x3V9iqDrgO+buvJDYCdyAr/OfR7f8NfHL3\n7nbVTKlocCdfoQ75x3yYu23yvuT+HnSZ7cO8M5lg90fy/tMpDyuuEvQBsZXuAz4J+zPIhLk/ekbb\ncym9/7trc7cGgdRtV6P3fRGC+EBvOxCtNeT6vdajKwUHc7ftiBYo3bichr7pxlQbPc8zyDWyEGV0\nuibJ42h4UnYx4/ehRcv9COA3Ioib6JgXgJ8i18uHgc8vXcqwYcNa96aLqFSCuz2y6DJqT2Ik4e7f\nt9vbr0Uw70wm0JPwz/b4aZTej7/c5XzLSbBv9PY3oZNAd2KY98yy353SWVw9Suw/f5sY6G+jK5oB\nKO56oLftT25ivp3qkY/8VRQD7mDeBUH8imhcTst85keRZe16ay5DMHaFrrJ1u1+FAP5bFLv/UbRo\nOSR6fBvwM+An6DP4x/nzmTVrVuq7+ZQduNurozQEvIN7tn3/PksmyDuTCfnGtudQOkAoVx0mE+Qb\nE9sN6FK7CzHM/dGLGPBpd9XtRgB/ixjqbyPXSw8Ecn8MIndXJPXoBPIKyqBcghYSL0Jp9aOj7TCa\n/03sRy4VV3Z2LwL4zGjru2ksAv1daMGyM/Ap4GNoPeUYSu75fnT8V+++m+uvv54OHdK5TB7AnUMl\nLf3djWz3oJRr39o/nexgb24EK79wqkUW2kZvbCAG+0Zk1Z9NDPLk9kIE9/a2IMuHahG83/LGm9H2\nBARwB/JLo9GL9gcSHEcLnou98S5KFhqNrOoraT6l/120EPkYcq1cguK/55C5YFmPgH87WtycAvwF\nstwNcp/8Mwop/PYDDzBv3rzUWeAB3CmQRQuyPsx3J/Z3efft8rYgK/CcaJsN7l28x9x+KS/YpVn1\nCO4O5tnGVvQ/6OWNC6Ph9ruRnsgqi7IyHcRXRfur0FpUEuaXohC89tiqB5BV/rI3TkQAd1UML6fx\nE+AxFHY4H3gYGVXzUOTJZGJrfj/KwPwJAvXngc+gdZH5wJeR5f6j5csZMmQIaVEAd4nrEDHEfcj7\n9/lbt29oCHN/3799rnc7nReOpaU6BG8H8ve9rdvfhyzzC7OMi6KRhpPvPmKYr0L1x1cht9JAlHDj\nxhA077ackCyyqF1N8ZeQ/34EcZr9OBrvpvQWsq7vj/5uLlqwnIYgbqPn/A8U1fI3wN8hgN8KfCO6\n72vHj6eiEXIAdwXKIuAnYe5vd0b7O71j9qEoAQf2cxNbf79rtN+Z4L9viw4j18t7COTvEYPd3deR\nGOL+6I0qHOY21aR1OoCs8pXRWBFt9yKIu16drv54U4k3jWk/qiv+IsrMXIJi2yehEOOJZI8t34ga\nONyNXEM3Ap9DrhmQD/47qH7K14G/RCGZn0XupN9u3/7Hjj0hHDAo9aojtuQd3JP7/nZndPyZ6Id5\nLjHQz23kvq6ELkItkXPJvJcY673tycQQ753Y70NxqhzuQQBfnhink1nVcBiy2Ftj2x4jTuipQRZ0\nf2RVT0NWefIq5R0U2/1LtDD7dyg88CS0mPn56HnvQJ/ZPyHoP7t5M927dw/gDipP1aMfqwP6jmh/\nR+K2u28HSqxxED8vsc12XwB9Q7kwyPXIxbA+y34n4l6f/uiHQh8L5WO36ErCVTR0tcc3Ievcz84c\nSstDFo+hhU7XtX45gvecaFzoHVuHaqL8AC0w/wvwkWhuPwK+izI0R0WP3Q+8ePgwp52mU0EAd1DF\n6wgxyLcTA317I/cdRQDvFo3zGtnvhk4IwU+vE+oW4kYTrjb5umh/H7IwkwWxLiY3kSQt0QEEcJeZ\n+RpyY/QnjgEfiWDekiidvQjg86NxIXKTfJRMiD8N/COKAroDLfw/glwlL6LP4Eb0+fxr9DcB3EFB\nrdQRBPDtaNFvO3IjbEvsb0M/3nNovC+ov6288vyxPiCzsqG/vxtZ5gMT4xKarvaXCx1BFrkfB/4u\ngrff/LhnM89Ti6JO7kZd50chN8lV6OrtOPAlBPglaL3gh2iB82lklQ9BV5MQwB0UlFfVIms+2Sou\nWZ98C7Iq/Trkfoq82+9OOmOy8ykH9WQCz9vIjXEJCgn0o0ny2S3yIPJtu4iThciXPY64YNUAGnen\nHUbJN99DJ+2fE1cb/Av0PfgJ+u70QXHeQ1EG5q+j4wK4g4JSIBdb7yDup8z7WZXbkJXpp8gnozzO\nJz0x2fmUK4jlwgNdJMkKtMg4JBrDozGY/Jz0LEr3d+VjFyDoTkHdcmaRPYqlFvgv4D+RlX0++v9e\nEm1PRoWpBgH/B0H7PuDeAkMbAriDgtqlOuSO8VPjXcieG3sQ0HsTh+v5/SzLve2tA7qLJnH+6neR\ni2U4SqYZg6JJch0lbZHv/mkU4vcMOoF8GPhTBGhf/xudhL8e3b4ALXJegKD+LgL8CuAG4M0Ugrsg\nkeZ9aVjbwx/nkhlJcAYhiiAoHepA7E65opFjDhPHYq+PxsPERaA6kLkAOABZo5eQ/ronLZEhrqU+\n3bv/EILf68jCvQV9TiMRxF3Ls/aGLRp0ouwH/DlazH4GFZ8ahCzwHxADfDSy0p3qvf2z0JUYaJF2\nUzvnli8VxOJeS8NUbr+mx07ihaYd0d+e18zojqycswmQD0qvXBMFB3FXb3sFsuz6ESepDEWLaPn0\nFxdbe1BSzUvIzfEaAukMVHfk0hy/3l5kWa9B1rgBvoZcW/+M3CdnIFifjPzfLwG/QP+7DkBtXV1B\n+lv6yrurxBgzE11hdAB+bq39f95jbXKVHERfdh/m2xL3bUFn73oa1oTw93tTHlZNUPnpCMo+XI4i\nKZahS/ZuyBqdgLIEexdpfoXQfmT9Po7iqM8HPo1cHLnq03kcGXgbouccC3wLJe28hxY5N0bH3oJc\nPj+Obp8G7PrgA04//fQczaZlyiu4jTEd0GLzNHRVsQT4E2vtm9Hjefdx70M+R1fJzR+uLkQf5PNy\nK+BD0GVr8SsSBAVlqg5Z489Howb5zT+J4pHbkjpeKqpDfurbUBnXL6JFwvYaXkeRC+Qocmt1RVf5\np6HP+EvIygb4N2QU/nt0+yxg0759Oe9l2Zzy7eMeBay11q6PXuxu1PbtzXY+b4vVCV1iDm3k8aPo\nzOIWTu6MtjtRdtWfotjOSgv5CkqnOhCnhf8NupR/CiWKfB34JkrTLkf3YAf0W7wKuTb+DsVVP0P7\n4L2DuH/pSrTG4FLjt5G5eHmEzMzNg0CnTp0KHg7YnNprcV8PzLDW/ll0++PAaGvt30a3G33yxh5o\n7AuZ7+N3kn31v1jzCceH45PHr0OREh9DmX/Fnk++j7eoVOsY4Kt5ms8vka/9tuaOb4STjdXxbu/x\n+ba403Uaaof2Uv5hW0GlrT7IP7uxuQPLRAaFEm7P42vUk54Y/JqaGmpqalp0bHvBvQmtAzr1Isv3\nqjV0b+2ZoKXHu0I3jyAf4iLkT5sEfII4k6pQ8wnHh+Nberzr6PID9IO7r8jzKcTxFrlJfkJm6F5b\nnv81dJXyJnK7fMt7/DTk93a6CS1m/nN0uzmXVGs9Fk0dX1VVRVVVleZx001NPk97wf0K0N8Y0xsl\nnt0I/Ek7n7PNqo0m4S9Qvk2c4dWReJHyWhT6k6uV66CgXOogcePc36C1nE8Bf01uG/6mTfWo9vb3\nUMTHfcSNf9uqAchoA8XOr0BAd81GdnrHnow++7SrXeC21tYaY/4GRfN0AP7HRZTkUvWo8M12b2wm\nM5Lk/ej+88gMDxyJvvCDCZAOSqfqkP96KarD8TxKKf8QSlB5iLgJQDmqHr33B1CfyI6ohsifk5sm\nHWdEr/EBShI6E6XtD0K5IL6L4DS0QJl2tTsizlr7BxTb3mIdI7Ph7k4yoZwcu1BYTjfi7MruCMyX\nE0O6B6EbS1B6VY8g4RJwlqEIp1XoOz0UJaP8GwrXKmfLeiNaFHwchf6dhaK87kcnrFxGzbyEokqc\nL3sKcpkMQqGW7xFb4EnXyYnAB0eP5nA2uVFBQpmnEvdL3I0+GNf89hziovrdUNLMSGJIuxrLAchB\npSDnrltLZrbkWmRVn4NyCIYgSP+vaL+YbcjyraPEVxOuut9hlGB0FfIr98nxa1q0jnUrsuRvIw4B\nvAKl4YOs71MQl7ogINZGjx1BMD/55PQFCxcE3F8mhnRndHZNy0puUFBrdIS4LonfSsyNLcjg6Etc\nn+Tj0X4/BIpy1n7gDeLmCK+jK4wBKKxvDvBt9FnkOhb9MFrEfQTVijkdZWG+jf4nTocSr+3qk3RB\n3gCH6W2kt/xAQcA9vflDgoKKKouycP2yrtn299OwtGsVcWXAXlRGMtdh4mbBbvF/JXJtDkVhfOPQ\nYuoQ8tOx3rUuq0GRJ4uj150DPIZcIcmTw3uo8t+d3nPsIW4QsZE4IedtdOJNo0LWd1BZy6I1Er+B\nQramClvQj9xvpHABWhScHe33QpZbpVwtWgTit4mbKrjtJnQV4Zop/K9o25f8tY7biSJtXorGq6hs\n7CTg79FCbmMupw9QaOG/ogzUSdH996ITTafo9vPA/+ftX5nTd5A7BXAHlZws8kluSwy/jdnWaGxD\n7gm/ZZlb2B5FZgeccvYzN6U9ZPri3fZtdDIbRNy6bGK034/8rjsdRi4X175sEfpfjkIul68gqHZq\n7AnQ9+RV4FfRqELWuatA+DoC/gPR7dfQe56O/Ny/Ae5evDg3byjHCo0UgoouiyyiphoGOyi7+88g\ns0lwtubB3aNtOUdntET1aMF0nTfeIV5APU5cKzxZN7xrlufLtT5AETauAcMryC9+CVpIHBWNS2ne\nmq9DVvlDyJoGJd98hrjKYi2yvr8Vba9HJ6+JwD9Ex96KOuA8F/1NaF0WVPaqRyUEdnpjh7fN1tXd\nEDfT6EpmYw23dWA+j1CqN6m9aKH0XW84SK9HZU37IEvZbR2gu1K4olXbUITJUuJFzPXIqh8OjEAh\nvpfR8hPuJpSo9BgKLeyJ/NzXRs/n3tsx1FzhW9ExP0ZXD5tQLfApqCrgchQJ98yyZQwbNgwI4A4q\nMVm0IOcaXuzy9n0g+2M3soi7olBOF+7ZhYZg9rseBWWXRWD2o1feJe628y6yNPsgq9Jt+yGfcx8K\n//n6tcaXE8P6OGoacRmK1x6OLOnWLOjuRJbwMwjY24HJqBnDTDJrcICigH4B/AxdRXyZOGDiN8jK\n/ntU3nU56pjzH8BHrPUB2ooZtl8B3EF/VC1x16HdZEI4G5jdcacg6J6b2Drwuvt8QIe4+5arHvnj\n/TDDZLghZEayJCHdmeKUeq1FbhcXXbI82q5HFv1gYlAPQ5Zua+e5CXgBJew8hz6n8WiBcSo6ASRd\nKHuQ7/pXyLL/KMrEdBmoL6F63wdR6YuRqBP8XwM3AzdGXAzgDsqZ6lHomkto2pVlm9zfhb6kZ6Mf\neZcs49xG9oNbon06QNzU431vuNsb0f/lQgRlt3WjN8Vv0efS8lcmxmq06DsERWe47UDaFhZZFz3v\ni9F4AX1vxyFQT0SgzhZVsQPFb98b/d0UVG9/LnK7WFTb/HvIt38TKjC3GzVsWAL89tVXGTFixB+f\nM4A7qIHq0CVwEsC7ybSMk2MvuvRNAvgcb79zYtsF/fgrJZStUDpM9u5L/jhOZjs9Ny7y7k/LAupR\nFFXyJkrFd9u1aI3hUuIQwMHIN90eN8xOFDGyMBpL0KLyOG8MJPtJy6IoEJdw8wZqxXU98lm7KKED\nyPK+JXqeL6AFy3q0OPltVM/oW1lalAVwl6ksysTygbsncXt3I/cdRF8uP6vUB3BT94U4zvzrELpM\n34Dg7MYGb3uAOMY7W9/TXsTdV9KkvQjKb0XD7b+PLPxB0bg0GgNpv5/8CILrYgTrRcg/7UL8roxG\nU3XxD6GQvj8Aj6IFx2pkVU8lPgFa4qa/90WP/RXyhdeiBJxvoZPPv61axaBBg9r57nKrAO4W6ggx\ndH3YNnWfA/GJ6MfpoNqZTNCek+W+zigGNVjAhZfLknRZkRvJzJB0YD5EDGW37UkM5AtId0JOHfKP\nJ5Nn3kIheJdEY5C3fzG5yfysQyeDJd5YhRYHR6E6LaOj12wqxM8iv/kT0ViIFjRno0XEYWSeFN9B\nYXx3ot/lp1G/zu7o//lL4PvofX7z+ecZP358e99qXlQx4LbosnUvMVTdfnI494R/Xx0xXM9J7Dd2\nn7udlsvcIP0ftxJDOFvq+ib0Y3fw7elt/f1zSZ+lnE27EZBXR1s33kEnFgflgcSJNBeQu/dWj1wp\nr3jjdeT7voI4Fns4qh/SnDagaJGnonEmigKZgSzmZNLNRuD3wF0owuZGVB9mFHqPm1Do38+Q6+VL\nCxdy5ZVpzYmUSgrc9ejS08F1TxP72aBsaAjbc4gX5PzbyWNOpzR+pJUqi74bDrybaQjnTejS+1wy\n09b94cBcalmShxGIVxMD2m2PEUN5gLffn9yHATpIv4ayEt04B0VmXI5AfTlx/Y/mtBO5P55GIX67\n0cLiVOSz7pvlbzaiRcjfIUv+ahQ5MhVZ2hZFodyC4P9x4G/ffpsBAwa04t0WT6kA92/IDuAkhPcj\ngDqwJrdJGCf381HIJij/OopqhTggb6YhnDdHx/oQ7pHldndKNxSxFoXRrUFQdtvV6CqiD3FGow/o\nbuTH6KhDJ4bXiQH9BvqtXU6cLDMSnSxbqt0orG9BNN5D4X0O1sPI7n5ai/zV96HPZi7wEQR3597Z\nCdyBrOsTgL8EPrlvHx07ltapOhXgvpEYsG5kA3AnwqJbOakWZcptbmJsQpZ0N2IQ+1sfzqX108uu\neuQKWJMYqxHAuhPD+WJiK/oi8vvbOIrC8Fwp1tdQGvr5CNDDiWHdGkiDoqOeQxbws+jKYQxye0yO\nnjfbe6tHJ4sHUUz2LmRZXxf9nTtB1wFPooXIJ1Bn+D977jnGjx/faFf1tCsV4C6FxcmglqsOuSS2\nkB3G7v6d6Efu4OuKOfkw7h4dk9YFvraoHp2QknBeg3ywXYjrgvijH4VZL9lPZm2QN5Bl3RcB2kH6\nMmRYtVZbUWU9B+r3UHf6iQi4I2n8qugIcpc8hEL8zgKuQcAeTeb35E3U6uxO9F36NPCne/Zw9tlt\nmXV2hXDAoNTLAdmHb7btDnSF5FwTPRL7DsjdKN8rKNeGLNnlZg1xpxu/YJNfwKklC3S5kEX/rzfI\nbG6wBSXKDPfGUNrmarQIzM8hWD+HvkMuYWYSstKb+h5sAeajeOwFyFVyNXKFDEwcux34LYL1RuS7\n/uTy5QwZ0t6WwtkVwB1UNB1HLgtXd3pzln1nIXdBl8eNQdmVRy1VP3JrVIvimtcmxhrki+5Mw6p6\nbhS6Nkgtcre8kRgWZRq6uiDDo7m29YRaj1wqz3ujFlnTE1E7siE0HeJXjyJP5kdjHWphNgeF+CXj\nuA8iC/zXKJtyDvCJxx5j6tSpnHhifk2DAO6gnOsIDRsD+FB2YzcKC/Phm9x3QC5XC7kxHUXui7Vk\nljp9B0G7G5lA7kfs1iiU5ZzUPuKiTa7K3ioUMXMZgrMr4tSD9i1cHkGQfSEaLyGwTvBGS9qQ7UF+\n6Pmoil8XlDhTjazzpCFwDDUSvgsl2owFPvarX3H11Vdz5pmFawAXwB3UIrmwt2wdWpLjA+LmAEkI\n+/d1JX9dSUpBB4mh/A6ZgN6KMh37kQnl/iiDsJgx+vXopLI0MXYg18YwBGdXwCkXONuJ4PwCsm7f\nQJmT4xCkxxG39mpu7m+gDMc/IJ/6BATqWWRvDlyLXCW/Rd3eLwX+BLhh+3a6di1EZfCGCuCucNWj\nH0Vj7bP824aGHVuyjS6E2HPQ92sHMZT9sRadCPsS16D2txeSDrfPAeLyp8ui7Qq0OHhZYvQjNydi\ni9wrrqDTS+iKbTQKzxsX7bf0hLALWdWPIWu5I4L0LORGyeZDr0PulntQXPZFKArtI++/T69eyQKt\nhVcAd5nqCDF0G9tuQYsqnWgI3/MT+z0o/07gbVEtWozKBud3kIvHWc0Xe/v9aL+7IJeqRz7dZYmx\nBVmYznp225YmsbREh1F4nYP0S8jd44o5jUWWfEvdZbUolf1xBOtVaDFyFspy7NfI39VFc/gdyng8\nH8Vjf2TtWvr1a+yvKksB3G2QRT65pmDstoeI22SdT0Mr2d3uRmV0/26PDhG31UqO91HDBQdj34Lu\nixYK06a9ZMLZ1aruTAznodG2P7l3Z20hE9LL0clhbDTGIb94a7QBgfpxlOnYi7iBwTgaLwFch9wv\nv0OW9XkI1jeUUDZjIRXA7ekocXSFD+Dk/lZ0WZe0hrPdDu6Klst1DU/C2d3ei3yfDspu69ptpbUm\nzHHkclhODOil6OTv6lMP80buoowz57AMAXphtD2AEl0cpEfS+miXD1CI3+MouWU7qhsyM9r2aOJv\na1Ec9++Rz7o7KrkaYN28yh7cFkVNNAVhd99BZPn6EPYt4u7efkifb5uOoVA5B2N/uw7B14eyD+ce\npDsRx8VFu3ZczpJejSxXH9BD0XvL1/vZgZriOki/ivz1zpoeQ+N1rJuSW1R8EoF6EUrGuQpZ1iNo\n+j0dRQk096KMx94I1h9es4aLL764lbOpXJUsuA/RELxJGG9FFvQZNLSKz0+M7ugSNc1gKAW5E2US\nyO72VgSxvjSEc18aVnZLqw4in62zoN3WEAPabQeT39DAWuRiWeiN7Wjh0EF6NG235DciUD+JCjKd\ng0A9HWU6ntXM3x9CFvm9KNzvUuDDwIfXr+eiiy5q46wqW6kCdy2yFLIBeAuZboxjZIdv0p98PqG1\nVq51DPmU19EQzOsQvBoD84WUVix4LYo8WZ4Ym1H50ySk81XQyddOBGdnUb+CrkYcpMegGtpt9Ykf\nQNX4HKx3oOJO0xGwL2zBc+xFkL4/eo4rgGtvvplrr72WHj2acqAEtUSpAPcwBOPdyOrNBuLkthPB\nd5wvWRS6tQ7FCSfhvAX9D5I+ZgfnNC4ENieLYOzA7BrbvoXe61AyIX0xhTkB1SJr/mViUDtr+kpi\na7o9n7mL/nCgfgOBdno0mnN/OG1DGYz3oUXPScCHf/lL5s6dS5cuTfWtKV1VdDjga8TFhErJGitl\nHUU1JLKBeR36ofpA9hcF0xLb3FbtQ2nZvgW9AlmofkPbYeiyvpDhl1uJremX0W/jQuK2Xe21piHu\nxegaEdQgX/M0BOoJtNy1sw5Z1fejz3AmcO3ddzN79mzOOqs5J0rpq6LBnZaoknKSRRbQuzSE87ro\nsZ7EYPbh3IfStJqTOoosZt+CXoGuJi4lE9JD0dVcoef3OjGkX0aV+XxIX0Fuoky2ktk1xiBIT41G\nS9+7RRExDyBYb0UFn66ZP5+pU6dyyimV5ZgM4A5qtQ4RQzkJ53dR5EvSanZw7kX5XN241O+km+Nd\n9F4dnN22D4VfhLYomuZlFInxcjTHAWS6PAaQGxfgARRq9xQC9iagithX3b8Vr1OLXB8PRMMA1wLX\nPv88Y8aMoUOHyi2YEMAd1EB16AeXBLPb34cucR2M+5JpNZdDYwFfLtxuBZkW9JsoXt6H8xC0eFgs\n++8A8hs7S3oRcm+MJoZ0W+KmG9Ox6HUcqJehnorTEKwvp3XulUPI3/0AKqfaC9W9vmbpUoYOHVqy\nDQhyrQDuCtUeGgLZ7W9AQHIg9gHdF60LlGv44h5iP/QKb5xAQwt6MMU9SdWh0MBFxNb0OlTgyUH6\nSuSaymUD3qUI0k8ji3ggsetjPK3PNdiBIP0AKuZ0BZEb5L33uPDClsSSVJ7KEtzGmO+j8rjH0LrX\nZ4PbAnEAABQLSURBVKy1+7zHyx7cR9AioA9mf2vJtJL9bW/Smw2YKx1C0FuRGPsQmAdHWwfpQvuh\ns2krsRW9CIXjnY/g7CA9jNwu4Fr0A3IW9QJ0UncNc6to27rEWpQI8yA6EUwHrr7jDqqrq+ncuRxW\nOspT+Qb3dOBpa229MeZ70Qt92Xu85MFdj8LIktayu70TXWYmwez8zudQGWGNx1AG4QpkSTtAb0L+\nVh/Og1EVuDRcTbjCS4uJYX0AAdpBehT5WczdQmxRP4N8zVO90do6IqDv6xJiWO9GnWSumT+fKVOm\ncOqp5W4qlIcK5ioxxlwLfNha+3HvvtSD2xWU8mH8rjfeQ/BNLgC60ZPKqnVdR7xQ6AP6HRTW5sO5\nkPHQLVE9Orks8sZbKPzOt6YvJj8n2z0oNM+BeiuypKehDudtSVEHnXz8Po1nIxfI1QsXMmrUKE44\nIQ2nyKDWqJDgfhi4y1r7G+++VID7MFrxbwzO9TT0M7v93hSv00kxZZEP3ofzCgS6rsQuDgfpQaTP\n7bMdWdIO0ovRSXi0N4aTv7o0H6CKeM8gWL+Nij1NicZw2n7S34kyFx+MnvsyIliHmiBloXaD2xjz\nJNmbXnzVWvtwdMzXgBHW2g8n/rYg4K5D9RaSQHZjF7o09y1lH86dqQx3Rja5eHAf0CujcTqZfuih\nKD46jWkXh1HMtG9N70ULcA7SV5BfH7qL/HgmGq+hrMQpyPUxmvaV9l2NrOoH0RXPVGDebbdRXV3N\nueee245nDkqbmgN3s1ex1trpzbzAp4HZ6HuUVd/09qui0Rq5DifOWl5PpuW8EVmBPpCnE4O5B5Xl\nzmhMe8iEs9u3xHAejjpnD6Fh09a0KOnyWIxCBi9BcJyFvnMDyK8fvQ6dLJzrY2H0mlOBr6HIj/aE\nA9ZFz/lQNA4gf/VXH32UyZMnB391mammpoaampoWHdvexcmZwL8Dk6y1O7M83mKLez8NgezGehSv\nm81a7oOs6crK52paB8iM5HCQPogA7YZzc5xPuq84tiI4O7fHEnSV5CzpUeTX5eFk0WfpXB/PARcQ\nW9QTaX/HmoOonOpDyBVyATAPmLt4MSNHjgzx1QVWuYYDrkFXf7ujuxZaa//Ke/yP4PZrZ/i+5vXR\n9ghxskk2P3OplAMtpA4hS9OB2bk4diDrMxlu14t0AxoErteIIb0YnYh8l8codIWVb1kUTudcHwuQ\nm2gqKnc6mZY1zm1OG9Ci4sMoXvtK4Oqbb2bOnDmhLGqRVZbgbsGLW5DVsIO4dkZy9EU/xLRDpVg6\niha2fOt5JXGonW89D0afaSm4hurQ+/Ct6bXoJOMAPYrWpW+3V++TCWqLLOrJ0TYXGK1HJycH6/eR\ne2fePfcwY8YMOnYst5zY0lVFg3s9gndawsLSquPAGjLhvAJ9fn3ItKAHo7C1Uqni56JUnKtjEYJX\nDzIhfRmFdXttQ4B2sN5PXPNjMrk7aRxG7pWHUfbimUQukGefZezYsZx4Yvh1pFEVDe5ihwOmTXUo\n7tl3b6xE1mZPMuE8BC14lZoPfy8CtO/ysGRC+gpy28W8JdqFYqkXRGMz8k27EL3B5G5BcyuC9EPR\na45Ai4tzQ8/FklEAdwWqHlnLSUC/jXyjg1F4nYP0IEqz1+VRlE7tuzw2owVDB+rRFMfHvg8tIjqr\n+l3UScaB+kPkzq3kSqI6F8ga1Kdx7q9+xaxZs0KKeQkqgLuMZZGfMrlI+BaKfvAjOYYgQBeyeH8u\n5ULxFntjBboquIIY0pdSHNeYn/SyAEXXjCZ2fYwkt+6lI9HrOBfIyciqnvPUU0yYMIGTT25P5HZQ\npSqAO4eyaEEwCeg3EYh9F4cbpb7MtIVMd8cryL3h3B0u+zBX5Utbq8Mo1tlZ1EtRiVO3mDia3LuZ\ntqFQvYei17yMyAWyahWXXHJJCNkLarcCuNsgi/yTSRfHKpTWPTjLKLSvNh/aj8DsLOklCIzOknZ+\n6fOKNUHklllE7KN+BUXUuMiPseS+RIFF9a+dC+Rt1FR37h13MGvWrJC1GJRzBXA3IYvqWazMMjqQ\nHdDl8hM9hmDkuzzeQ9azi5m+AoVqFtN+PI7g7ED9MopRr0KgnkB+UvCTLpCTiFwgTz7JxIkTgwsk\nKK8K4I60g+yAduneyVFMqzLXqkcLZb4lvRxVO3SQHoXed7HDC+tQmKAD9UsoFNK5PiaQmx6N2bQF\nuUAejl47uECCiqWKA/cusgP6GLEP+lJvvxvll/izmUxL2vmlfUt6BOlYIK1Dlr8D9fMoJNJlJk4i\nfzVTLKo14qzqtSgKZM6ddzJr1iy6dElrtZagclfZgns3mb5nt3+YOMzOXyzsTvkBGhTu5vzSLrHl\nKJl+6UKliLdE9ej/5ED9LLq6mUzs/sjn1c4HKBHmEWRdn0HkAnnmGcaPH89JJxX7miMoTQrhgG3U\nXhoCegX6ATo4++MCyhPQIL+ri5deEo0NKBbZh3Qf0vMZWBR144O6E7FFXYWyJ/OpDQjUjyCLfiTq\ntzcnJMIENaMA7ma0F/3AXSSHg/R+BGiXqOJgXQoFk9qjOhS94Ls8VqF4aT/CYzDpKiVgUZz3ApQt\nWIOSinxQ98rzHOrQSc3BeiOqBTLn7ruZMWMGZ5+dLy95ULkpgDvSfjJdG27sRYkpvvV8KWqHVe6N\nl1wCj0sRX4L6IJ5HDGhXujRt3XhcBb0aYlifSCaoexdgHvtROdRHgEfRZzcXmPPCC1x55ZV06FAK\nZbeC0qaKBvc/EAN6DwrnSro4KgHQTruIIe1AbYgh7UYal8YsSuN3ro8a5Lf2fdSFCiFcS2xVL0Jt\nweb86EfMmTOH3r17F2AGQeWuigb394gBfRGVA2iQL/51YkAvRv0CLyfTmu5Jel0/7xFDegGK0Kki\ntqrz1Vw3qeOoXvUjKBJkP/JVV99/P9OmTePMM9MQJxNUTqpocKchjrsQOo589L41/Q46YfmQHki6\nT14byPRRf0AM6ira3o28LdoJ/AFFgDyOYs/nAtVLljBixIjQwTyoLBXAnSe5pBYX3bEYxSP3JtPl\nMYz0l2TdRKZF7WpSVyFYD6JwoLYoOWg+sqxXRHOY87OfUV1dTffu3Qs0k6Cg4imAO0faRKYl/QrK\n4PMhfTnp7ICe1GYyFxP3oESXKgTJwRTWbXM4mouLrT6ByKp+7DGqqqo45ZS0n/qCgnKrAO42aDcC\nsx/lcZwY0GkottQabSV2eyxA7oeJxFb1UArvutlIbFU/i2LR5xDSy4OCIIC7WR1Ci4e+y2Mrsp79\nCI/epHfxMKltxKCuiW5PRFb1ZOS+KTSo69Bn62C9EZgJVP/mN8yYMSM0GQgK8hTA7ek4Ckn0XR5r\nkGvAh/QgSqPZrtN2ZLXWIIt6CyrG5BYTh1Gc97MPLSg+ghYYzweqgernnmPMmDGhz2JQUCOqWHDX\nozhf392xFIUj+n7pQjenzYV2IFA7H/UmYDwxqHPZjqs1sijbc340XonmVX3zzcyZM4eLLrqoCLMK\nCmq7QjhgHuU607hFw8Uo87ATmX7pEZRmR5odxH0Ta5CbYRzxYuJwineFcDSam3OBHCGKrX7oIaZO\nncrpp6ct1zMoqOUK4M6hdhP7pN2oJdPdUUqLh0ntRDCsicZ7CNTOoh5OceuTbEFp5Y+g1l2DiWD9\nxhsMGzYsLCwGlY0CuNuoD1BhfR/S25H17Ls8LqJ0Fg+T2kUM6gXEoHaLiSMoLqjr0ZWMc4G8g1p3\nzbnjDmbOnEnXrmkpGhsUlFsFcLdAx1DyhR/hsQ5VBfQt6YGU1uJhUj6oa4B3Ua9EZ1FfTvEr/rmi\nTfPRwmJnooXFBQsYN25cqFsdVBEK4E6oHi1k+YuHK1At6WTmYal399tNpo/aB/UkBOpiY9CiCBuX\nBLOYuGjT7Nmz6du3bzGnFxRUFFU0uOvR5b/v7ngNNd71Lem0tNNqrxyoa6Kxjtj1UUU6QA2ZC4vz\nUUx7NVD9wANMnTo1FG0KCiqSUgHuc9Glv591OJJ0li1ti7KBegxx1EdaQA3xwuJ8tLA4iAjWr73G\nhz70obCwGBSUAqUC3Bsor5ZiDtQu6eUdYlBXoZNSWkDd2MJi9e23M2vWrLCwGBSUQqUC3GnJnGyr\nfFAvIN0WNShj0V9Y7EIUrldTw9ixY8PCYlBQyhXA3QY15vpwi4lpsqghM2PxEeKMxTk338zs2bPp\n06dPMacXFBTUSgVwt0DZQD2WeDExbaAGZSjWIFg/ikIpq1HG4pQpUzjjjDOKOLugoKD2KIA7i5qy\nqKtIn+vDaSPxwmINKsdaDcwOGYtBQXlRRYcDFhvcjSW8VJFuUNehJrhuYXEDMAOY8+tfM2PGDLp0\nKZe4nKCgdCqAu4DKBuo0xlFn025UCtX1WOwBzAaqn3+eK6+8MpRCDQoqoMoW3MaYLwDfB8611u5O\nPFYQcLuiTC48r1QsatBns4LYql6KTjBzfvpTZs+eTa9evYo5vaCgilZawd0u880Y0wuYjhIjC6ad\nCNIO1H5RpltRBmZaQQ3KUHyGGNYdkK/6q48+SlVVFaeddloxpxcUFJRytfe6+wfAF4EHczCXRuXq\nUddE433ietQ/p/jV81qid4lB/QK6CqgGHlu5kkGDBoWFxaCgoBarzbwzxlwNbLTWLss1dFyHF5fw\nsgHFJVcB/0NpgPo48CIxrHehHoufvece7po+nbPPPruY0wsKCmqBCu0iaama5J8x5knUKjCprwFf\nQdnTfzy8sef5prdfFQ1ffs/EGuJWXJOAX1L8xgEt1TaUqfgo8CTQD1nVty1axMiRIznhhEK36A0K\nCioV1dTUUFNT06Jj27Q4aYwZAjyN3LUAPRFvR1lrt3vHZV2c3EYM6meJQV2FYF0qoK5HVQ6dVb0a\nmAZU/+IXzJw5k+7duxdzekFBQSWqgoQDGmPeBS5vLKpkC5mLia4LuQvPK1Zz27ZoP7KmXR2QTkQZ\ni08/zfjx4zn55FKvHh4UFFRs5TWqxFOT9B+EQF0FfI7SArVfB+RR1GBgLIL119aupV+/fkWcXVBQ\nUCWqIAk4tZQOqCGuA+LSy10dkNkPPsiUKVNCg4GgoKC8qlAWd5MqBWhvILMOyDCUsXj/0qUMHTo0\nhOsFBVWgipWA05zKMuW9JaoFFhLDejMK15sd6oAEBQVFSmvmZEWBeyfwGAL1E8CFyKqe89JLjBo1\nig4dSuHaICgoqFAK4C6CLPA6sVW9CpVurb71VmbPns0FF1xQpJkFBQWVggK4C6QDwFPEUSBnElXX\ne+IJJk6cyCmnnFLA2QQFBZWyArjzJIsSX5xVvQg1RZgNVK9eTf/+/fP46kFBQeWsAO4c6ghK5nGw\nPoJAPfv++5k2bVoI1wsKCipppSIcMBfaSOz+qAGGIFj//vXXueyyy0K4XlBQUMUotRZ3LfAycR2Q\nEK4XFBRUKSopV0kyXK8XUR2QF19k9OjRIVwvKCioIpRqcNcDbxC7QFYBU4DZIVwvKCiogpU6cLvq\neo9G4ywiq/rJJ5kwYUII1wsKCqp4pQLcbxJb1UtQuF41MCuE6wUFBaVYFR0O2JMorvrBB5k6dSpn\nnHFG3l4zKCgoKFdKK7gL0kvr/fp6/tta5s2bVxHQbmn7oXJSeM+VoUp8z2lUQcBdaTHWlfjlDu+5\nMlSJ7zmNCt1rg4KCgkpMBfFxBwUFBQW1XkX1cQcFBQUF5U55tbiDgoKCgnKvYHEHBQUFlZgCuIOC\ngoJKTHkDtzHmBmPMSmNMnTFmROKxrxhj1hhj3jLGXJWvORRTxphRxpjFxpjXjTFLjDFXFHtO+ZYx\n5m+NMW8aY1YYY/5fsedTKBljvmCMqTfGdC72XPItY8z3o//xUmPMfcaYTsWeU75kjJkZMWqNMeZL\nxZ6Pr3xa3MuBa4Hn/DuNMZcCNwKXokqtPzbGlKPl/6/A1621w4H/G90uWxljJgPzgGHW2iHAvxV5\nSgWRMaYXMB14r9hzKZCeAAZbay9Dzae+UuT55EXGmA7AzYhRlwJ/YowZVNxZxcobMK21b1lrV2d5\n6GrgLmvtcWvtemAtMCpf8yiitgDOGjkb2FTEuRRCfwl811p7HMBau6PI8ymUfgB8sdiTKJSstU9a\na+ujm4tQRYty1ChgrbV2ffSdvhuxKxUqhqXbAzW0cdoIlGP91i8D/26MeR/4PmVqmXjqD0w0xrxs\njKkxxows9oTyLWPM1cBGa+2yYs+lSPosqh1XjroA2ODdThWn2tW6zBjzJHB+loe+aq19uBVPVZIx\niU28/68Bnwc+b6293xhzA/ALdEldsmrm/Z4InGOtvTLy598D9C3k/PKhZt7zVwB/jaYsaju05Hdt\njPkacMxa+5uCTq5wSjWT2gVua21bQLQJNbdx6kmJuhGaev/GmF9Za6dFN38P/Lwws8qfmnm/fwnc\nFx23JFqs62Kt3VWwCeZBjb1nY8wQoA+wNKrF0xN41Rgzylq7vYBTzLma+10bYz6NCn5OLciEiqMk\np3qR6SkoqgrlKvEtkYeAjxpjTjbG9EGX2IsLNI9Caq0xZlK0P4X/v707RokgiIIw/BfCeAHB1A30\nAJ5CI/EChgaLqRdRDN0DGIlgIB7ARDBQ9CiKLJRBLxuJJq7Da+rLJquBmYKm+9FtI6dn17T3RNIO\nMFQv7Z/YfrG9aXtie0L7qXerl/ZvJO0Bp8CB7Y+x86zQI7AtaUvSQDtQcTNypqWV3fIu6RA4AzaA\nW0lPtvdtv0q6ot1UNgem7nN88xi4kLQOvC+eezYDZpKegU/gaOQ8/63Hb/g758AA3C9WGg+2p+NG\n+nu255JOgDtgDbi0/TZyrKWMvEdEFNPj+emIiK6luCMiiklxR0QUk+KOiCgmxR0RUUyKOyKimBR3\nREQxKe6IiGK+AJ6488gc9t9OAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x10afa54d0>"
]
}
],
"prompt_number": 9
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment