Skip to content

Instantly share code, notes, and snippets.

@dmitry-kabanov
Last active January 18, 2017 11:50
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save dmitry-kabanov/1009b117306a72770cc8c5617f32d531 to your computer and use it in GitHub Desktop.
Save dmitry-kabanov/1009b117306a72770cc8c5617f32d531 to your computer and use it in GitHub Desktop.
Examples of time series for Slava
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Examples of time series for Slava"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example 1: one non-oscillatory mode"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"t = np.linspace(0, 50, num=5001)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f39de2283c8>]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAi0AAAFkCAYAAADsVgtLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3XmUXVWd9//3lyQEAiTaIAljAwuZVAIppjDJjIAPoIJY\nSOtPUB8lNDzh6YX62C0N3WiDbQLa2A6IgEI0YDOFIQidQIAEQlVIACMghJnEhCFhCJCk9u+Pfau5\nKepW6lZu1bm36v1a66yqe+4+e+/slQWf7LPPPpFSQpIkqd6tU3QHJEmSusPQIkmSGoKhRZIkNQRD\niyRJagiGFkmS1BAMLZIkqSEYWiRJUkMwtEiSpIZgaJEkSQ3B0CJJkhpCj0JLRIyLiAURsTwiZkXE\nnmsof1BEtETEOxHxRER8uZMyJ0bE/FKdcyPiqGrajYjBEXFhRMyLiDcj4sWIuDIiNutQx/SIaCs7\nVkXET3syDpIkqe9UHVoi4iTgR8C5wO7AXGBqRGxSofw2wBTgLmA0cAlwWUQcXlZmX+Aa4JfAbsCN\nwA0RsUsV7Q4rXXte6fvPADuW6iqXgF8AI4FRwGbAOdWOgyRJ6ltR7QsTI2IW8EBK6azS5wCeB36c\nUrqok/IXAkellHYtOzcJGJFSOrr0+XfAsJTSsWVlZgJzUkqn96TdUpk9gAeAv00pvVA6N61U79lV\n/cElSVKhqpppiYghQBN51gSAlFPPncDYCpftU/q+3NQO5cd2VaaH7QJ8iDyz8nqH81+MiMUR8UhE\nfD8i1u+iDkmSVAcGV1l+E2AQsKjD+UXkWzGdGVWh/PCIGJpSereLMqN62m5EDAX+DbgmpfRm2VdX\nA88CLwG7AhcBOwAnVKhnY+BI4BngnU7/hJIkqTPrAdsAU1NKr6xtZdWGloYQEYOBa8mzLKeXf5dS\nuqzs42MR8TJwV0Rsm1Ja0El1R5KDjiRJ6pkvkteurpVqQ8sSYBV5EWu5kcDCCtcsrFB+WWmWpasy\n7XV2u92ywLIVcEiHWZbOPAgEsD3QWWh5BuC3v/0tO++88xqqUq2MHz+eiRMnFt2NAcUx73uOed9z\nzPvW/PnzOeWUU6D0/9K1VVVoSSmtiIgW4FDgJvifBbGHAj+ucNlMoOPjy0eUzpeX6VjH4e1luttu\nWWDZDjg4pfRaN/5Yu5NnZF6u8P07ADvvvDNjxozpRnWqhREjRjjefcwx73uOed9zzAtTk+UVPbk9\nNAG4ohQiHgTGkx83vgIgIn4AbJ5Sat+L5WfAuNJTRJeTg8YJwNFldV4CTI+Is4FbgGbywtuvVdHu\nYOAP5MeePw0MiYj2mZlXS8FnO+Bk4FbgFfIj2BOAu1NKj/ZgLCRJUh+pOrSklCaX9kY5n3x75mHg\nyJTS4lKRUeRbM+3ln4mIY4CJwJnAC8BpKaU7y8rMjIiTgQtKx5PAcSmlP1XR7hbksELpO8i3fRJw\nMHAP8B5wGHAWsAH5kelrS21KkqQ61qOFuCmlnwKd7iKbUvpKJ+fuIc+cdFXnH8gzJT1t91nyE0Zd\nXf8CcFBXZSRJUn3y3UOqS83NzUV3YcBxzPueY973HPPGVvWOuANNRIwBWlpaWly8JUlSFVpbW2lq\nagJoSim1rm19zrRIkqSGYGiRJEkNwdAiSZIagqFFkiQ1BEOLJElqCIYWSZLUEAwtkiSpIRhaJElS\nQzC0SJKkhmBokSRJDcHQIkmSGoKhRZIkNQRDiyRJagiGFkmS1BAMLZIkqSEYWiRJUkMwtEiSpIZg\naJEkSQ3B0CJJkhqCoUWSJDUEQ4skSWoIhhZJktQQDC2SJKkhGFokSVJDMLRIkqSGYGiRJEkNwdAi\nSZIagqFFkiQ1BEOLJEnqFfPn17Y+Q4skSeoVt95a2/oMLZIkqVf85S+1rc/QIkmSesVTT9W2PkOL\nJEmquVdeyUctGVokSVLNPfZY7es0tEiSpJp75BEYPLi2dRpaJElSzc2bB9tuW9s6DS2SJKnm5s2D\nj360tnUaWiRJUk21teXbQ4YWSZJU1xYsgLfeMrRIkqQ6N29e/mlokSRJdW3ePPjIR2DjjWtbr6FF\nkiTV1Lx5sOuuEFHbeg0tkiSpph55JIeWWjO0SJKkmnnrrfyixE98ovZ1G1okSVLNPPYYpORMiyRJ\nqnPz5sE668Auu9S+bkOLJEmqmXnzYIcdYP31a1+3oUWSJNVM+5NDvcHQIkmSaiKlOgwtETEuIhZE\nxPKImBURe66h/EER0RIR70TEExHx5U7KnBgR80t1zo2Io6ppNyIGR8SFETEvIt6MiBcj4sqI2KxD\nHUMj4tKIWBIRb0TEdRGxaU/GQZIkve/FF+G11+ootETEScCPgHOB3YG5wNSI2KRC+W2AKcBdwGjg\nEuCyiDi8rMy+wDXAL4HdgBuBGyJil7Iya2p3WOna80rffwbYsVRXuYuBY4DPAQcCmwN/qHYcJEnS\n6tq376+b0AKMB36eUroqpfRn4BvA28CpFcp/E3g6pXROSunxlNKlwHWletqdCdyWUppQKvM9oBU4\no7vtppSWpZSOTCn9IaX0ZErpwdL1TRGxJUBEDC+VH59SujulNAf4CrBfROzVg7GQJEkl8+bB8OGw\n9da9U39VoSUihgBN5FkTAFJKCbgTGFvhsn1K35eb2qH82K7K9LBdgA8BCXi99LkJGNyhnseB59ZQ\njyRJWoOHH4bRo2u/fX+7amdaNgEGAYs6nF8EjKpwzagK5YdHxNA1lGmvs+p2S3X/G3BNSunNsnbe\nSyktq6L/kiSpG1pbYffde6/+wb1XdXEiYjBwLXmW5fRa1Dl+/HhGjBix2rnm5maam5trUb0kSQ3t\n8ssn8eSTk9hgAzj22Hxu6dKlNW2j2tCyBFgFjOxwfiSwsMI1CyuUX5ZSencNZdrr7Ha7ZYFlK+CQ\nslmW9nbWjYjhHWZbuuo/ABMnTmTMmDFdFZEkacDaYYdmoJkrr3x/IW5raytNTU01a6Oq20MppRVA\nC3Bo+7mIiNLn+ytcNrO8fMkRpfNdlTm8vUx32y0LLNsBh6aUXutQZwuwskM9OwJbd+iPJEmqwpw5\nMHQo7Lxz77XRk9tDE4ArIqIFeJD8VM8w4AqAiPgBsHlKqX0vlp8B4yLiQuBycmA4ATi6rM5LgOkR\ncTZwC9BMXjT7tSraHUx+dHk34NPAkIhon5l5NaW0IqW0LCJ+BUyIiNeAN4AfA/eVnjaSJEk9MGcO\nfPzjMGRI77VRdWhJKU0u7Y1yPvm2ysPAkSmlxaUio8i3ZtrLPxMRxwATyY82vwCcllK6s6zMzIg4\nGbigdDwJHJdS+lMV7W5BDiuUvgMI8rqWg4F7SufGk281XQcMBW4HxlU7DpIk6X2trbBnl1vNrr3I\nTw6rkogYA7S0tLS4pkWSpE68+y5suCFccgmcXvb4S9malqaUUuvatuO7hyRJ0lp57DFYubJ3H3cG\nQ4skSVpLc+bkDeV6a/v+doYWSZK0VubMgR13hA026N12DC2SJGmtzJkDfbHs09AiSZJ6bNWq/M6h\n3l7PAoYWSZK0Fp58Et5+29AiSZLq3Jw5+aehRZIk1bU5c2DrreFv/qb32zK0SJKkHpszp29mWcDQ\nIkmSeigleOgh2GOPvmnP0CJJknrkqafg9dcNLZIkqc7Nnp1/GlokSVJdmz0btt0WNtmkb9oztEiS\npB6ZPRv23LPv2jO0SJKkqq1aBa2tfXdrCAwtkiSpB+bPzzvhOtMiSZLq2uzZEAFNTX3XpqFFkiRV\nbfZs2Gkn2GijvmvT0CJJkqrWl5vKtTO0SJKkqrz3Hsyd27frWcDQIkmSqjRvXg4uhhZJklTXZs+G\nwYNh9Oi+bdfQIkmSqvLQQ/Dxj8P66/dtu4YWSZJUlb7eCbedoUWSJHXbW2/BY48ZWiRJUp2bPRva\n2mCfffq+bUOLJEnqtlmz8oZyu+zS920bWiRJUrfNnAl77QWDBvV924YWSZLULSnl0DJ2bDHtG1ok\nSVK3LFgAixcXs54FDC2SJKmbZs7MPw0tkiSprs2aBR/9KGy8cTHtG1okSVK3FLmeBQwtkiSpG95+\nO7/Z2dAiSZLqWksLrFxZ3HoWMLRIkqRumDkTNtggvyixKIYWSZK0RrNm5fcNDR5cXB8MLZIkqUtF\nbyrXztAiSZK69OyzsHBhsetZwNAiSZLW4P77809DiyRJqmv33gs77gibblpsPwwtkiSpSzNmwAEH\nFN0LQ4skSerCa6/Bo4/C/vsX3RNDiyRJ6sJ99+WfzrRIkqS6NmMGbL45bLtt0T0xtEiSpC7ce2++\nNRRRdE8MLZIkqYLly2H27Pq4NQSGFkmSVMGDD8KKFYYWSZJU5+69F4YPL/YlieUMLZIkqVMzZsB+\n+8GgQUX3JDO0SJKkD1i1Km/fXy+3hqCHoSUixkXEgohYHhGzImLPNZQ/KCJaIuKdiHgiIr7cSZkT\nI2J+qc65EXFUte1GxGciYmpELImItojYtZM6ppe+az9WRcRPezIOkiT1V/PmwRtv1Memcu2qDi0R\ncRLwI+BcYHdgLjA1IjapUH4bYApwFzAauAS4LCIOLyuzL3AN8EtgN+BG4IaI2KXKdjcAZgDnAKnC\nHyEBvwBGAqOAzUrlJUlSyYwZsO66sGeX0xJ9qyczLeOBn6eUrkop/Rn4BvA2cGqF8t8Enk4pnZNS\nejyldClwXamedmcCt6WUJpTKfA9oBc6opt2U0m9TSv9KDkhdPVH+dkppcUrpr6XjzWoGQJKk/u7u\nu2HvvWG99YruyfuqCi0RMQRoIocCAFJKCbgTGFvhsn1K35eb2qH82K7K9LDdrnwxIhZHxCMR8f2I\nWL8HdUiS1C+1tcH06XDwwUX3ZHWDqyy/CTAIWNTh/CJgxwrXjKpQfnhEDE0pvdtFmVFr0W4lVwPP\nAi8BuwIXATsAJ1RZjyRJ/dK8efDqq40fWhpeSumyso+PRcTLwF0RsW1KaUGl68aPH8+IESNWO9fc\n3Exzc3Mv9VSSpGJMmwZDh8I++3T/mkmTJjFp0qTVzi1durSm/ao2tCwBVpEXsZYbCSyscM3CCuWX\nlWZZuirTXmdP2u2uB8nrX7YHKoaWiRMnMmbMmLVsSpKk+jdtGuy7b3XrWTr7h3xraytNTU0161dV\na1pSSiuAFuDQ9nMREaXP91e4bGZ5+ZIjSue7KnN4e5ketlvp6aGOdi+Vfbmb5SVJ6rdWrYJ77qm/\nW0PQs9tDE4ArIqKFPEsxHhgGXAEQET8ANk8pte/F8jNgXERcCFxODhonAEeX1XkJMD0izgZuAZrJ\nC2+/1t12S21/GNga2II8e7JTKdwsTCktiojtgJOBW4FXyI9gTwDuTik92oOxkCSpX5kzB5Yu7Seh\nJaU0ubQ3yvnk2zMPA0emlBaXiowCtior/0xEHANMJD/a/AJwWkrpzrIyMyPiZOCC0vEkcFxK6U9V\ntAtwLPBr8sxJAtpvrp1Xuu494DDgLPKeLs8D15balCRpwJs2DYYNg732KronHxT5yWFVEhFjgJaW\nlhbXtEiS+r2jj4aVK+GOO9a+rrI1LU0ppda1rc93D0mSJABWrMg74R5ySNE96ZyhRZIkAfDQQ/Dm\nm/W5ngUMLZIkqWTaNNhoI6jhU8o1ZWiRJElADi0HHACD63TrWUOLJEli+XK49144tOOuaXXE0CJJ\nkrj3XnjnHTjiiKJ7UpmhRZIkcccdsNlm8LGPFd2TygwtkiSJO+7IsywRRfekMkOLJEkD3MKFMG9e\nfd8aAkOLJEkD3p2lF+scdlix/VgTQ4skSQPcHXfAbrvBppsW3ZOuGVokSRrAUoI//rH+bw2BoUWS\npAHt0UfzmhZDiyRJqmt33AHrrQf77Vd0T9bM0CJJ0gB2xx3wyU/m4FLvDC2SJA1Qy5fDPffA4YcX\n3ZPuMbRIkjRATZ+et+4/6qiie9I9hhZJkgaoW26BbbaBnXcuuifdY2iRJGkASimHlmOOqe+t+8sZ\nWiRJGoDmz4dnnsmhpVEYWiRJGoBuuQXWXx8OOqjonnSfoUWSpAHollvgkENycGkUhhZJkgaY11+H\ne+9trFtDYGiRJGnAueMOWLXK0CJJkurcLbfAxz8OW29ddE+qY2iRJGkAaWuD225rvFkWMLRIkjSg\nzJ4NixfD0UcX3ZPqGVokSRpAbroJ/uZvYN99i+5J9QwtkiQNINdfD//rf8HgwUX3pHqGFkmSBojH\nH8874R5/fNE96RlDiyRJA8SNN+bN5I44ouie9IyhRZKkAeKGG+DII2HYsKJ70jOGFkmSBoCXX4aZ\nMxv31hAYWiRJGhBuugkGDYJPf7ronvScoUWSpAHghhvgwANh442L7knPGVokSernli2Du+5q7FtD\nYGiRJKnfu/VWWLHC0CJJkurcf/0XjBnTeC9I7MjQIklSP/bWWzBlCnz+80X3ZO0ZWiRJ6sduuQWW\nL4cTTyy6J2vP0CJJUj82eTLssQdst13RPVl7hhZJkvqpN97IMy394dYQGFokSeq3pkyBd97pH7eG\nwNAiSVK/NXky7LUXbLNN0T2pDUOLJEn90LJlcNttcNJJRfekdgwtkiT1QzffDO++CyecUHRPasfQ\nIklSP/T738PYsY2/oVw5Q4skSf3MK6/A7bf3r1tDYGiRJKnfufZaaGuDL3yh6J7UlqFFkqR+5je/\ngSOOgJEji+5JbQ0uugOSJKl2nn4a7r8frrmm6J7UXo9mWiJiXEQsiIjlETErIvZcQ/mDIqIlIt6J\niCci4sudlDkxIuaX6pwbEUdV225EfCYipkbEkohoi4hdO6ljaERcWirzRkRcFxGb9mQcJEmqN1df\nDRtuCMcdV3RPaq/q0BIRJwE/As4FdgfmAlMjYpMK5bcBpgB3AaOBS4DLIuLwsjL7AtcAvwR2A24E\nboiIXapsdwNgBnAOkCr8ES4GjgE+BxwIbA78obt/fkmS6lVK+dbQZz8Lw4YV3Zva68lMy3jg5yml\nq1JKfwa+AbwNnFqh/DeBp1NK56SUHk8pXQpcV6qn3ZnAbSmlCaUy3wNagTOqaTel9NuU0r+SA1J0\n7EhEDC+VH59SujulNAf4CrBfROzVg7GQJKluzJ4NTz4Jf/d3Rfekd1QVWiJiCNBEDgUApJQScCcw\ntsJl+5S+Lze1Q/mxXZXpYbudaSKv4ymv53HguSrrkSSp7vzmN7DZZnDwwUX3pHdUO9OyCTAIWNTh\n/CJgVIVrRlUoPzwihq6hTHudPWm3Ul/eSyktW8t6JEmqKytWwO9+ByefDIMGFd2b3uHTQ900fvx4\nRowYsdq55uZmmpubC+qRJEnvmzIFliyBL32pmPYnTZrEpEmTVju3dOnSmrZRbWhZAqwCOj75PRJY\nWOGahRXKL0spvbuGMu119qTdSn1ZNyKGd5htWWM9EydOZMyYMVU0JUlS3/nVr2CPPWDXDzw32zc6\n+4d8a2srTU1NNWujqttDKaUVQAtwaPu5iIjS5/srXDazvHzJEaXzXZU5vL1MD9vt7OmhFmBlh3p2\nBLbu0B9JkhrGiy/mNzqfdlrRPeldPbk9NAG4IiJagAfJT/UMA64AiIgfAJunlNr3YvkZMC4iLgQu\nJweGE4Cjy+q8BJgeEWcDtwDN5EWzX+tuu6W2P0wOIFuQnx7aqRRuFqaUFqWUlkXEr4AJEfEa8Abw\nY+C+lNKDPRgLSZIKd+WVMHQo9PcVC1WHlpTS5NLeKOeTb6s8DByZUlpcKjIK2Kqs/DMRcQwwkfxo\n8wvAaSmlO8vKzIyIk4ELSseTwHEppT9V0S7AscCvybMsCWi/uXZe6TrIYWcV+bHrocDtwLhqx0GS\npHrQ1gaXXw4nnggdll72O5GfHFYlETEGaGlpaXFNiySp7kyfnh9xvvtuOPDAonuzurI1LU0ppda1\nrc8XJkqS1MB+9Sv46EfhgAOK7knvM7RIktSgXn8drrsOTj0V4gP7wPc/hhZJkhrU1VfnTeWK2pul\nrxlaJElqQCnBpZfC8cfD5psX3Zu+YWiRJKkB3X03zJ8Pp59edE/6jqFFkqQG9NOfwk479d+XI3bG\n0CJJUoN56SW4/vo8yzIQFuC2M7RIktRgfvnLvAPuQFmA287QIklSA1mxAn7xCzjllP6/A25HhhZJ\nkhrITTfl20Pf/GbRPel7hhZJkhrIT34C++0Ho0cX3ZO+15O3PEuSpAK0tORHna+9tuieFMOZFkmS\nGsTEibDttvCZzxTdk2IYWiRJagAvvAC//z2cdRYMGlR0b4phaJEkqQH85CcwbFh+OeJAZWiRJKnO\nvfkm/Pzn8PWvw0YbFd2b4hhaJEmqc7/+dQ4uZ55ZdE+KZWiRJKmOrVwJF18Mn/88bLVV0b0plo88\nS5JUx669Fp5+GiZPLronxXOmRZKkOtXWBt//PnzqU9DUVHRviudMiyRJdermm+HRR+E//7PontQH\nZ1okSapDKcEFF8CBB8L++xfdm/rgTIskSXXorrtg9myYOrXontQPZ1okSapDF1yQ17EcfnjRPakf\nzrRIklRn7rsPpk+H//oviCi6N/XDmRZJkurMuefCJz4Bxx1XdE/qizMtkiTVkWnT8nqWG26AdZxa\nWI3DIUlSnUgJ/vEfYY894Nhji+5N/XGmRZKkOnH77XD//XDbba5l6YwzLZIk1YH2WZb99oMjjyy6\nN/XJmRZJkurADTdAa2te0+IsS+ecaZEkqWArV+ZZlkMOgYMOKro39cuZFkmSCnb55fCnP8FVVxXd\nk/rmTIskSQV64w343vfglFN8k/OaGFokSSrQD38Ir7+et+1X1wwtkiQV5MUX4d//HcaPh623Lro3\n9c/QIklSQf7pn2CDDeDb3y66J43BhbiSJBVgzhy44gr4yU9gxIiie9MYnGmRJKmPtbXBuHHwsY/B\n//7fRfemcTjTIklSH7vqKpg5E6ZPh8H+n7jbnGmRJKkPvf46nHMONDfDJz9ZdG8ai6FFkqQ+9L3v\nwfLl+akhVcdJKUmS+sjcuXDppXDhhbD55kX3pvE40yJJUh9YtQq+8Q3YcUc466yie9OYnGmRJKkP\nXHopzJoFM2bAkCFF96YxOdMiSVIve/ZZ+H//D04/Hfbfv+jeNC5DiyRJvSilfFvowx+GH/yg6N40\nNm8PSZLUi66+Gm6/HW6+GYYPL7o3jc2ZFkmSesmiRfB//g+cdBJ8+tNF96bxGVokSeoFKcFXvwrr\nrAM//nHRvekfvD0kSVIvuOwymDIFbrwRNt206N70Dz2aaYmIcRGxICKWR8SsiNhzDeUPioiWiHgn\nIp6IiC93UubEiJhfqnNuRBzVk3Yj4vyIeCki3o6IP0bE9h2+nx4RbWXHqoj4aU/GQZKkzjz1FIwf\nD6edBsceW3Rv+o+qQ0tEnAT8CDgX2B2YC0yNiE0qlN8GmALcBYwGLgEui4jDy8rsC1wD/BLYDbgR\nuCEidqmm3Yj4FnAG8HVgL+CtUpl1y7qUgF8AI4FRwGbAOdWOgyRJnVm1Cr785Ty7MnFi0b3pX3oy\n0zIe+HlK6aqU0p+BbwBvA6dWKP9N4OmU0jkppcdTSpcC15XqaXcmcFtKaUKpzPeAVnIAqabds4B/\nSSlNSSk9CnwJ2Bw4vkOf3k4pLU4p/bV0vNmDcZAk6QMuvBDuvz+/yXmjjYruTf9SVWiJiCFAE3nW\nBICUUgLuBMZWuGyf0vflpnYoP7arMt1pNyK2Jc+clJdZBjzQSd++GBGLI+KRiPh+RKxfoe+SJHXb\njBn5hYjf+Y6byPWGahfibgIMAhZ1OL8I2LHCNaMqlB8eEUNTSu92UWZUFe2OIt/66aoegKuBZ4GX\ngF2Bi4AdgBMq9F+SpDVasgSam2HffeG884ruTf804J4eSildVvbxsYh4GbgrIrZNKS2odN348eMZ\nMWLEaueam5tpbm7upZ5KkhpFWxt86Uvw7rswaRIMHnD/d4VJkyYxadKk1c4tXbq0pm1UO6xLgFXk\nRazlRgILK1yzsEL5ZaVZlq7KtNfZnXYXAlE6t6hDmTkV+gbwYOm67YGKoWXixImMGTOmi2okSQPV\nv/873HZbPrbYoujeFKOzf8i3trbS1NRUszaqWtOSUloBtACHtp+LiCh9vr/CZTPLy5ccUTrfVZnD\n28usod32MgvIwaW8zHBg7y76BvlJpAS83EUZSZI6NX16fhnid74Dn/pU0b3p33oygTUBuCIiWsiz\nFOOBYcAVABHxA2DzlFL7Xiw/A8ZFxIXA5eRQcQJwdFmdlwDTI+Js4Bagmbzw9mvdaPfXZWUuBv4x\nIv4CPAP8C/AC+RFqImI74GTgVuAV8iPYE4C7S08bSZLUbc8+CyeeCAcdBOefX3Rv+r+qQ0tKaXJp\nb5TzybdeHgaOTCktLhUZBWxVVv6ZiDgGmEh+tPkF4LSU0p1lZWZGxMnABaXjSeC4lNKfqmiXlNJF\nETEM+DnwIWAGcFRK6b1SkfeAw8iPRm8APA9cW2pTkqRue/ttOP542HBD+P3vB+Y6lr4W+clhVRIR\nY4CWlpYW17RIkoD8XqHm5vzm5pkzYdddi+5RfSpb09KUUmpd2/rMhZIkVenCC/PsyuTJBpa+5Fue\nJUmqwuTJedHtd7+b17Oo7xhaJEnqpnvvzfuxnHyyC2+LYGiRJKkbHn8cjjsO9tkHLr8c1vH/oH3O\nIZckaQ3++lc46igYORKuvx6GDi26RwOTC3ElSerC66/DkUfmR5z/+7/hwx8uukcDl6FFkqQK3nwT\njj46byJ3992wzTZF92hgM7RIktSJd97Jm8c98gjcdRd84hNF90iGFkmSOlixAk46Ce67D26/Hfba\nq+geCQwtkiSt5t134fOfz29svvFG+OQni+6R2hlaJEkqWb4cPve5vOD2xhvzE0OqH4YWSZLITwcd\nd1y+JTRlChx2WNE9UkeGFknSgLd0aQ4sDz2Ubwt5S6g+GVokSQPayy/n20DPPgtTp8J++xXdI1Vi\naJEkDViPP543jlu5EmbMgI9/vOgeqStu4y9JGpBmzcqzKhtsADNnGlgagaFFkjTgTJoEBx8MO++c\n39y81VZF90jdYWiRJA0YbW3w3e/CySfDiSfCH//ou4QaiWtaJEkDwhtvwCmnwM03w0UXwT/8A0QU\n3StVw9BJ9D9EAAAM8klEQVQiSer35s/PMyvPP5/3YDn66KJ7pJ7w9pAkqV/77W9hzz3zraFZswws\njczQIknql5Yvh69/Hf7u7+Czn4XZs/PCWzUubw9JkvqdRx7J61eeeAIuuwxOPdX1K/2BMy2SpH5j\n1Sr44Q9hjz3y7w88AKedZmDpLwwtkqR+YcGCvPfKt74FZ56Z3yO0665F90q15O0hSVJDW7UK/vM/\n4TvfgY03hmnTfOFhf+VMiySpYc2bl7fi//u/hy9+MX82sPRfhhZJUsN5+2349rdhzJi8ady998LP\nfgbDhxfdM/UmQ4skqWGkBL/7Hey0E1x8MZx3HsyZk2db1P8ZWiRJDeHBB2H//aG5Oc+wPPpofo/Q\nuusW3TP1FUOLJKmuPfMMfOlLsPfe+VbQnXfCDTfA9tsX3TP1NUOLJKkuvfQSjBsHO+wAU6fmNStz\n5sChhxbdMxXFR54lSXVl8WK48EK49FJYf304/3w44wzYcMOie6aiGVokSXXhuedgwgT45S9h0CA4\n5xw4+2wYMaLonqleGFokSYV69FG46CKYNAk22gj+7//NO9pusknRPVO9MbRIkvpcW1tep/If/wG3\n3gpbbZXfGfTVr3obSJUZWiRJfebVV+HXv87b7j/1FOy+O1x5ZX6MeciQonunemdokST1qpTgvvvg\n8svzLaC2NjjxRPjNb2CffXwDs7rP0CJJ6hXPPANXXZWPp56Cv/1b+Kd/yreANt206N6pERlaJEk1\ns2RJ3vjt6qth+nTYYIM8q3LZZXDggbCOu4NpLRhaJElr5a9/heuvh+uug2nT8u2gT34yr1X57Gdd\nWKvaMbRIkqr25JNwyy1w001w99353MEH56eBPvMZGDmy2P6pfzK0SJLW6J13cji59dZ8/OUv+UWF\nBx+ct9c//nj4yEeK7qX6O0OLJOkDVq6E1tZ8u2faNJgxA95+O++ncvTReefaQw7Ja1akvmJokSSx\nYgXMnZtnU6ZNg3vuyW9UHjYM9t8fzj03h5WPfcxHlFUcQ4skDUAvvQSzZuVj5kx46KF8C2i99WC/\n/eBb38q3fvbc003fVD8MLZLUj6WUA8rDD+eZlDlz4IEH4Pnn8/dbbZU3eLvggvyzqQmGDi22z1Il\nhhZJ6ieWL89P9Tz6aA4nDz+cjyVL8vcjRsDo0fCFL+SAsvfesMUWxfZZqoahRZIazJIl8Oc/52P+\n/Pd/X7Agz6wAbLMN7LYbnHFG/jl6dN6R1vUoamSGFkmqMytWwHPPwdNP5yBS/vPpp/NLByHvLrvt\ntrDzznkTt513hp12yj8//OFi/wxSbzC0qC5NmjSJ5ubmorsxoDjmfWPlSli4EF58EX73u0lsu20z\nL76YP7/wQn5fz/PP55cKQg4mW2+dw8no0Xk/lB12yMFk++3zwll1n3/PG1uPQktEjAP+ARgFzAX+\nPqU0u4vyBwE/Aj4GPAdckFK6skOZE4HzgW2AJ4Bvp5Ruq7bdiDgf+CrwIeA+4Jsppb+UfT8UmACc\nBAwFpgKnp5T+WtUgqFf5H5a+55j33MqV8MoreTv7xYvzz46/tweThQvfDyQwiaFDm9liC/7nGDsW\nttsuh5TttssLZX16p3b8e97Yqg4tEXESOYB8HXgQGA9MjYgdUkpLOim/DTAF+ClwMnAYcFlEvJRS\n+mOpzL7ANcC3gFuALwI3RMTuKaU/dbfdiPgWcAbwJeAZ4F9LZXZOKb1X6tLFwFHA54BlwKXAH4AD\nqh0LSf3Lu+/Ca6+9f7z+euXPr776fih59dX315K0Gzo0v8l4003zTrGf+AQcddT74WTLLfNjxbfd\n5joTqbt6MtMyHvh5SukqgIj4BnAMcCpwUSflvwk8nVI6p/T58YjYv1TPH0vnzgRuSylNKH3+XkQc\nTg4gp1fR7lnAv6SUppTKfAlYBBwPTI6I4aXyX0gp3V0q8xVgfkTslVJ6sAfjIamPpZQDxvLleW+R\n5cvz8eab+Xjjjep/vvZarqszQ4bkNSLlx5Zbwu67vx9M2sNJ++8bbrjmMLLuugYWqRpVhZaIGAI0\nAd9vP5dSShFxJzC2wmX7AHd2ODcVmFj2eSx5FqVjmeO6225EbEu+bXRXWZllEfFAqcxkYA/yn7m8\nzOMR8VypjKFFIoeCtjZYtSofK1fmY8UKeO+9fHTn9+6WKw8e5b9X+vzOOx+c2ejMoEGw0Ub52HDD\n1X9uueXq5z70odVDSfnn9dc3XEj1oNqZlk2AQeTZi3KLgB0rXDOqQvnhETE0pfRuF2VGVdHuKCCt\noZ6RwHsppWVdlOloPYBrr53PAw+s/h/KSr+X63h+ba7vbl19cX1XdVV7fWfXPv74Uv75n1t7fH13\nrik/2tre/779947flX/urFxX31Wqo9J3a6pz5cr8uTxYtP/e2bnufP/uu0sZMqT1f+qttSFDYPDg\n/LP9GDw4H0OHfvDYaKM8c7Huuu+fW2+99z+X/97+ef3187bz7cfazGS89x4sWpSP3rJ06VJaW1t7\nrwF9gGPet+bPn9/+a02WjPv00JptA/Bv/3ZKwd0YeM47r6noLgw4bW29N+YrVuRj+fJea6IhNTX5\n97yvOeaF2Aa4f20rqTa0LAFWkWcsyo0EFla4ZmGF8stKsyxdlWmvszvtLgSidG5RhzJzysqsGxHD\nO8y2dNX/qeSFwc8AFe54S5KkTqxHDixTa1FZVaElpbQiIlqAQ4GbACIiSp9/XOGymeSndcodUTpf\nXqZjHYe3l1lDuz8plVkQEQtL5+aVygwH9iY/IQTQAqwslbm+VGZHYOsO/Sn/M79CfrJJkiRVb61n\nWNr15PbQBOCKUohof/R4GHAFQET8ANg8pfTlUvmfAeMi4kLgcnJgOAE4uqzOS4DpEXE2+ZHnZvLC\n2691o91fl5W5GPjHiPgLeWbkX4AXgBvhfxbm/gqYEBGvAW+Qg9J9PjkkSVJ9qzq0pJQmR8Qm5I3g\nRgIPA0emlBaXiowCtior/0xEHEN+WuhMcog4LaV0Z1mZmRFxMnBB6XgSOK59j5ZutktK6aKIGAb8\nnLy53AzgqLI9WiCHnVXAdeTN5W4HxlU7DpIkqW9F6s5zg5IkSQVbp+gOSJIkdYehRZIkNQRDyxpE\nxLiIWBARyyNiVkTsWXSf+ouIOCAiboqIFyOiLSKO7aTM+RHxUkS8HRF/jIjti+hrfxAR34mIByNi\nWUQsiojrI2KHTso55jUSEd+IiLkRsbR03B8Rn+pQxvHuRRHx7dJ/XyZ0OO+410hEnFsa4/LjTx3K\n1GS8DS1dKHtJ47nA7uQ3S08tLQjW2tuAvKD6dPJuxqspewHm14G9gLfI479uX3ayHzmAvEXA3uQX\nlw4B7oiI9dsLOOY19zz5RbBjyE9E/jdwY0TsDI53byv9I/Pr5P92l5933GvvUfJDMqNKx/7tX9R0\nvFNKHhUOYBZwSdnnID/9dE7RfetvB9AGHNvh3EvA+LLPw4HlwOeL7m9/OMivx2gD9nfM+3TcXwG+\n4nj3+jhvCDwOHAJMAyaUfee413aszwVau/i+ZuPtTEsFZS9pLH+5YiK//LHSyyFVI5VegAm0vwBT\na+9D5BmuV8Ex720RsU5EfIG8v9T9jnevuxS4OaX03+UnHfde89HSrf6nIuK3EbEV1H68ffdQZT15\nOaRqpzsvwFQPlXaUvhi4N72/H5Jj3gsi4uPkHbfXI29o+ZmU3y4/Fse7V5TC4W7AHp187d/z2psF\n/H/kma3NgH8G7in93a/peBtapIHpp8AuwH5Fd2QA+DMwGhhB3g38qog4sNgu9V8RsSU5kB+WUlpR\ndH8GgpRS+XuFHo2IB4Fngc+T//7XjLeHKuvJyyFVO+UvwCzn+K+liPgP8ms0DkopvVz2lWPeC1JK\nK1NKT6eU5qSUvkteFHoWjndvaQI+ArRGxIqIWAF8EjgrIt4j/wvfce9FKaWlwBPA9tT477mhpYJS\nQm9/SSOw2ksaa/byJ3UupbSA/Be6fPzbX4Dp+PdQKbAcBxycUnqu/DvHvM+sAwx1vHvNncAnyLeH\nRpeOh4DfAqNTSk/juPeqiNiQHFheqvXfc28Pda3Ll0Nq7UTEBuS/2FE6tV1EjAZeTSk9zxpegKnq\nRMRPyS8jPRZ4KyLa/+WzNKX0Tul3x7yGIuL7wG3Ac8BGwBfJ/+o/olTE8a6xlNJbQMc9Qt4CXkkp\nzS+dctxrKCJ+CNxMviW0BXAesAL4XalIzcbb0NKF1I2XNGqt7EF+FDGVjh+Vzl8JnJq69wJMdd83\nyOM8vcP5rwBXQbdfOqru25T893kzYCkwDzii/YkWx7vPrLYPlONec1sC1wAbA4uBe4F9UkqvQG3H\n2xcmSpKkhuCaFkmS1BAMLZIkqSEYWiRJUkMwtEiSpIZgaJEkSQ3B0CJJkhqCoUWSJDUEQ4skSWoI\nhhZJktQQDC2SJKkhGFokSVJD+P8BcbIO0FEGeFgAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f39e0673390>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"amp = 1e-8\n",
"alpha = 0.2\n",
"y = amp * np.exp(alpha * t)\n",
"\n",
"plt.figure()\n",
"plt.plot(t, y)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example 2: one oscillatory mode"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f39de1e6860>]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgcAAAFsCAYAAABYeho/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3Xd81dX9P/DXIWHLULYkslcYgQz2CqgoKIi4otZCh19b\nrZYua+3vq9aq/bZValvbWrVWReNWkC1LSgAZGeyEGZClBgh7Jef3xzufkhtDcm/yWedzX8/HI48r\nyb3nHD+Ee9+fc97nfZTWGkRERESWWl4PgIiIiPyFwQERERGFYHBAREREIRgcEBERUQgGB0RERBSC\nwQERERGFYHBAREREIRgcEBERUQgGB0RERBSCwQERERGF8G1woJQappSaqZTap5QqUUqNd7i/Wkqp\nJ5VSO5VSp5RS25VSv3ayTyIiIj+K9XoAlWgIIAfAKwA+dKG/XwL4HwD3ANgMIAXAv5VSR7XWf3Wh\nfyIiIl/wbXCgtZ4HYB4AKKVU+Z8rpeoAeBrAHQCaAtgA4Jda68+q2eUgADNK+wWAPUqpOwH0r2Z7\nRERERvLtskIYXgAwAMBtAHoDeA/AXKVUp2q2twLAaKVUFwBQSiUCGAJgjg1jJSIiMoZvZw4qo5SK\nBzAZQLzW+mDpt59TSl0PYAqA6uQK/A5AYwBblVLFkMDpUa312zYMmYiIyBhGBgeQmYIYAPnllhzq\nAPgaAJRS3QBsAaABfGNZovT7/6e1/lXpn28HcCdkmWIzgL4AnldK7ddav+HI/wUREZEPmRocXAbg\nAoAkACXlfnai9HEHgO5VtFNY5r9/D+AZrfV7pX/epJRqD+ARAAwOiIgoapgaHGRDZg5aaa0zK3qC\n1voCgPwI2mwAoLjc90pgdl4GERFRxBz94FNK3aeUylVKFZV+rVBKXRfmaxsqpRKVUn1Lv9Wx9M/x\nWuttAN4C8LpSaqJSqr1Sqr9S6peleQfV8QmAXyulxiql2imlJgKYCne2URIREfmG0lo717hS4yB3\n49sg6/6TAfwcQF+t9ZYqXjsCwBJIbkBZr2mtv6OUioEkHt4DoC0k12AVgMe01puqMdaGAJ4EMBFA\nSwD7IQHIk6WzEERERFHB0eCgwg6VKgTwM631q652TERERGFxLedAKVULUpOgAYCVbvVLREREkXE8\nOFBK9YIEA/UAHAcwUWu99RLPbQZgDIDdAM44PTYiIqIAqQegPYD5WuvCKp5bKceXFZRSsQCuAtAE\nwC0Avg9geEUBQmm54jcdHRAREVGw3aW1fqsmDTg+c1CazLez9I/ZSqn+AB4C8IMKnr4bAKZPn44e\nPXo4PTQqNXXqVEybNs3rYUQVXnP38Zq7j9fcXVu2bMHdd98NlH6W1oQXdQ5qAah7iZ+dAYAePXog\nKSnJvRFFuSZNmvB6u4zX3H285u7jNfdMjZflHQ0OlFJPA5gLYA+ARgDuAjACwLVO9ktERETV5/TM\nQUsArwFoA6AIwHoA12qtFzvcLxEREVWTo8GB1vp7TrZPRERE9uO5AYT09HSvhxB1eM3dx2vuPl5z\nc7leIbEySqkkAOvWrVvHJBYiIqIIZGVlITk5GQCStdZZNWmLMwdEREQUgsEBERERhWBwQERERCEY\nHBARUaAdPw585ztA48ZAv37ASh79VyUGB0REFFgXLgATJgDvvw9MnQrUqweMHg3k5no9Mn9jcEBE\nRIE1bRqwbBnwySfAE08AixcD3boBd98NFBd7PTr/YnBARESBdPgw8NRTwH33ASNGyPfq1wf++U9g\n40bg3//2dHi+xuCAiIgC6e9/B86dA/73f0O/n5oK3Hwz8OyzgI9K/fgKgwMiIgqc4mKZIbjzTqBl\ny2/+/MEHgS1bZJmBvonBARERBc68ecCePbKkUJHhw4EePYDXXnN3XKZgcEBERIHz7rvy4Z+SUvHP\nlQLuuAOYMQM4c8bdsZmAwQEREQXK+fPAzJnApEmVP+/WW4Fjx4D5890Zl0kYHBARUaAsWQIcPSpJ\nh5Xp0QNISJBAgkIxOCAiokCZNQu46iqgb9+qnztmDLBgAXctlMfggIiIAmXhQuDaayWvoCrXXgt8\n8YXsXKCLGBwQEVFg7N8vH/RXXx3e84cPB+rWldkDuojBARERBYZVt2DUqPCe36ABMGQI6x2Ux+CA\niIgCY9EiIDERaNEi/NcMHQqsWMG8g7IYHBARUWAsWwakpUX2miFDgMJCIC/PmTGZiMEBEREFwpdf\nAjt3AoMGRfa6gQOBWrWAzExnxmUiBgdERBQIq1bJ48CBkb2ucWOgd28GB2UxOCAiokBYtQpo0waI\nj4/8tUOGACtX2j8mUzE4ICKiQFi1SmYNwqlvUF5KiuQcHD9u/7hMxOCAiIiMV1wMrF4d+ZKCJSlJ\ndivk5to7LlMxOCAiIuNt2gScPFn94CAhQYohZWXZOy5TMTggIiLjrVsnywlJSdV7fe3aQJ8+DA4s\nDA6IiHxCa2D6dOC224Af/QjYtcvrEZkjJwfo3Bm47LLqt5GUJEEGMTggIvKFkhLgW9+Sr337gA8+\nkEp/zKAPT25ueKcwViYpCdi8GTh1yp4xmYzBARGRDzz1FPDWW0BGhuy3z8uTD7sJE6S4D12a1jJz\nYEdwUFICbNhgz7hMxuCAiMhjeXnAk08Cv/oVcMcd8r1GjYD33pMPvp/8xNvx+d2ePUBRkcy01ETP\nnpK3sHGjPeMyGYMDIiKPPfIIEBcHPPpo6Pdbtbo4o7B+vTdjM0FOjjzWdOagfn2gUyfZ+RDtHA0O\nlFKPKKVWK6WOKaUOKaU+Ukp1dbJPIiKT5OUBH38sswb163/z51OmAB06AM884/7YTJGTAzRvDlx5\nZc3b6tWLMweA8zMHwwD8BcAAAFcDqA1ggVKqgn8CRETR57nngNatJRGxIrVrAw88IAmKBw+6OzZT\n5ObKkkJ1KiOWx+BAOBocaK3Haq3f0Fpv0VpvADAZwFUAkp3sl4jIBKdOSQLivfdKAZ5LmTwZiI0F\nXnnFtaEZxY5kREvPnsCBA8Dhw/a0Zyq3cw6aAtAAovyyExHJcsLx48A991T+vMsvByZNktwDCnXi\nhNSD6N3bnvZ69ZLHaM87cC04UEopAH8CsFxrvdmtfomI/Or114GhQ4GOHat+7u23yx58TnmHysuT\nx4QEe9rr2lVmaaL9Ors5c/A3AAkA7nCxTyIiXzpyBFi4ELjzzvCef801QJMmwDvvODsu02wuvdXs\n3t2e9urUkQAh2mcOYt3oRCn1VwBjAQzTWh+o6vlTp05FkyZNQr6Xnp6O9PR0h0ZIROSuefPkJMEb\nbwzv+XXrAuPHAzNmSE0EElu2yDbQRo3sa7NnT/8HBxkZGcjIyAj5XlFRkW3tOx4clAYGEwCM0Frv\nCec106ZNQ1J1T88gIjLAzJlSkS8uLvzXjB0LvPEG8MUXkb0uyDZvtm9JwdKtm1Sp9LOKbpizsrKQ\nnGxPvr/TdQ7+BuAuAHcCOKmUalX6Vc/JfomI/Oz8eWDuXJkJiMS11wK1aslrSWzZAvToYW+b3boB\n+/dLsmO0cjrn4D4AjQEsBbC/zNdtDvdLRORbmZlS7jfcJQXLFVcAAwcyOLCcPQvs2GF/cNC1tFRf\nfr697ZrE6ToHtbTWMRV8ve5kv0REfrZokVT0q87e/Ouuk9cXF9s/LtNs2ybXwe5lBQYHPFuBiMh1\nS5YAI0fKEkGkRo4Ejh3jWQuALCkA9s8cNG0KtGx5cZtkNGJwQETkopMngc8/B9LSqvf61FTZubBs\nmb3jMtHmzTID07y5/W137cqZAyIicsny5cCFC9UPDurVAwYMYHAAyMyB3UsKlm7dOHNARBSRw4eB\nqVPljXngQODll4GSEq9HZYYlS+SgpZoU7Rk+XIIDre0bl4m2bLGv+FF51sxBtF5jBgdEFJGCApna\nfvVVWf9u0wb4/vflcCAmyVXts8+AESNqdoLg8OHA119fXHOPRloD27dfTB60W9eucu7FoUPOtO93\nrlRIJKJgOH0auOEGeWPOyQHat5fvZ2QAd98NdOgAPPGEp0P0tbNngays8EsmX8qgQZLMuHKlc9Pq\nfrd/v5xq2aWLM+136yaPeXky0xNtOHNARGH71a/kbm3mzIuBAQCkpwOPPQb89rfAmjWeDc/3cnKA\nc+fkw70mLrtMgoLVq+0Zl4m2b5dHp4KDjh0lAIvWpEQGB0QUlrw84K9/BR5//OKxtmU9+qjUpJ86\nNXrXaauyapUkFPbpU/O2+veP7uBg2zZZmunQwZn269YF4uOBnTudad/vGBwQUVgefRRo2xZ46KGK\nfx4TA/zxj1L9b/Zsd8dmipUrgeRkOfmvpvr3BzZskKWeaLRtG3DVVRJsOaVjRwYHRESXtH078OGH\nwK9/Xfmb8bXXyu6FZ591b2wmWbVKro8dUlMlATQ72572TLN9u3NLCpZOnaQ8czRicEBEVfrzn4Fm\nzSTpsCo//SmwdKkk3tFFBw7ITg+7goPevWXqO1qXFrZtcz444MwBEdElHDsG/OtfwA9+EN4U7k03\nyXHCL73k/NhM8vnn8mhXcFC7thz5HI0JoCUlMnPQubOz/XTsCBw5Il/RhsEBEVXqgw9ky9i994b3\n/NhY4FvfAt5+GzhzxtmxmWTNGtkSFxdnX5upqdE5c3DggORauDFzAAC7djnbjx8xOCCiSk2fLqV+\nI/lQ+/a3gaNHgU8+cW5cpsnOljt9O/XrJ3fQx4/b267fbdsmj27kHADRmXfA4ICILumLL6Tcbzi5\nBmV16ybZ9G+95cy4TJSdLR/mdkpMlMcNG+xt1++2bZMaBE5tY7RcfjnQpEl05h0wOCCiS3r7bUl6\nu/nmyF87aRIwf74sSUS7AweAgwftnzlISJBlnNxce9v1u+3bZRtj3brO9qNU9CYlMjggokv6+GNg\nzBi5e4rUTTfJuvCnn9o/LtNY2w3tnjmoW1cOHoq24MCNnQqWjh25rEBE9F9ffQWsWAGMH1+913ft\nCvToIQFGtMvOBpo2DS05bZe+faMzOHB6p4KlUyfOHBAFzrlzwOuvA9/9LvDDH8o0N0v7hmfOHHkc\nN676bdx0kyQlRvtpjdnZ8iFek5MYLyUxUXIOouXI7JISuZN3c+Zgzx7g/Hl3+vMLBgcUWHv2yJ7y\nb39bDrxZtAi47jrgllu4Dh6OmTPl+rVqVf02xo4FCguBdevsG5eJnEhGtCQmAidPRs/U96FDslxl\n7SRwWseOEtzu3etOf37B4IAC6cAB2X535IjsL1+3Dti6FXj/fZk9mDgx+u4EInH2rFynG2+sWTsD\nBgCNGgELF9ozLhMdPSrT0nYnI1qsHQvRsrSwe7c8tmvnTn9WrYNoCb4sDA4ocIqLgbvukruLpUuB\nlBT5vlKSQT9zJrB4sRwxTBVbsULuRq+/vmbt1K4NjBwZ3UmJOTny6NTMQcuWUlwpWoKDggJ5dCs4\nuOoq2TYZbYWQGBxQ4Dz/vAQFb75Z8RvIqFHAb34D/O530VldLhyLFgHNm9tztPA118hJjSdP1rwt\nE+Xmyq6Cbt2c66NPH2D9eufa95OCAtk907SpO/3Vrg1ceeXFoCRaMDigQPnyS+CJJ+QcgLS0Sz/v\nF7+Qg2t+/GMmKFZk0SK5frVseIe45hpZwlm2rOZtmWjjxov1CJzSsyewZYtz7fvJ7t3uzRpY2rVj\ncEBktMcfB2JiZGagMjExwJ/+BKxcKUcR00XHjkmexujR9rTXrRvQtq0s5USjjRuBXr2c7SMhQdbE\no+Esi4ICZ7aEVobBAZHBDhwAXnkF+NnP5HjhqqSlyXr4M89w9qCszz6TvA27ggOlgGHDgOXL7WnP\nJFq7FxyUlAD5+c724wcFBZw5cAODAwqM55+Xtd0f/jD81zzyiOxkiOZs+vIWLpQkLDu3ig0bBqxd\nG31bSPfsAU6ccD446NFDHjdvdrYfr2nt3bLCvn3RtcOJwQEFwvHjwN//LoFBJIlK11wjW8FeeMG5\nsZlm8WKZNbCzYM/QocCFC9GXALpxozw6HRxcfjnQpk3wg4PCQgkwvVhWKCmRACFaMDigQHjnHblD\nu//+yF6nFHDvvcCsWcD+/c6MzSRHj8oH2vDh9rbbq5dkmEfb0sLGjVLnIT7e+b4SEoIfHLi9jdFi\n9RdNSwsMDigQXn5Zqh9W5034rruAOnWAV1+1f1ymWblSHocMsbfdWrWkzf/8x952/W7TJtlJ4ETZ\n5PKiIThwuwCShcEBeerAAWDGDKlF/9VXXo/GHBs2AJ9/Dnzve9V7fZMmwB13SDJjtCcmZmYCLVo4\nc6jNsGFSXOnCBfvb9is3khEtCQlyING5c+7054WCAqBBA6nB4aYGDeTfBYMDctWRI8B3vgPExclB\nNePHy9avH/1Ipsqpcv/6l1SJu+GG6rdx991SAW3NGvvGZaLMTLnDd+JOd+hQ+X221uGDrrhY7uTd\nCg569JDAa/t2d/rzgrVTwY2ZmPKibccCgwOP7d0LDBoke+2nTZOElz17gCeflA+9kSM5i1CZkhLg\n3Xflzr927eq3M3y4BBjvvGPf2Exz/rzMwNi9pGBJSpL6EtGSlLhjh5xR4ebMARDsYkhe7FSwMDgg\n1xQVSbb8mTOyzevBB6VMZ3w88PDDche3dy8wYYK8ydA3rVghiYS3316zdmJj5bTG996LnqNvy8vJ\nkfMonAoOGjSQD8poCQ7c2qlgadFCptuDnHfgRQEkC4MDGymlhimlZiql9imlSpRS453szyRaA9/6\nlhw/+umnFa/x9u0r+QdZWVLml77p3XdlCWbgwJq3dfvtEoytWlXztkyUmSl1Ipw6PRAA+vePnqWb\njRvlw7plS/f67NEj+MGB1zMH0XLz4PTMQUMAOQB+CCDKU71CvfGGfPC//jrQpculn9e/P/Dcc8A/\n/hG95WcvpaREjmC+9VZ7zgAYOlTeyGfMqHlbJsrMBFJTJUBwSmqqfGhGwyFMGze6t1PB0q2bJCUG\nUVGRbLX1Mjg4d05u6KKBo8GB1nqe1vp/tdYzAHiQQuJPX30FTJ0K3HkncOONVT//vvuAESPkMZoq\ndFUlM1N2eNx6qz3t1aoFjBsnNQ+i0cqVzi0pWPr3l6AuK8vZfvxg69aLeQBu6dpVSigHcdeNNaXv\n5bJC2XEEHXMOPPDkk5LJPG1aeM+vVQv4y18kC/mll5wdm0lmzJCqcHYsKVjGjZNp2Wg7u33/fkmG\n7d/f2X569gTq1w/+0kJJidzBd+3qbr9du0q10CDe3XpVAMnC4IAcVVAgSwQPPxzZWmTv3sA998ip\ng9zeKGbPBsaOtWdJwXLNNbLrYfZs+9o0wdq18pia6mw/sbGS0xD0pMS9eyXRuFs3d/u1gpEgHsC0\ne7cUK2vd2pv+mzaVapfREhw4eMJ49U2dOhVNmjQJ+V56ejrS09M9GpF9nnhC6qA/+GD1Xjt9umxx\nrM7rg2TnTpm2ffppe9tt3FiWcGbNAh54wN62/WzNGqBVK6m14bT+/YGPP3a+Hy/l5cmj2zMHHTtK\nsJyXZ38JbK8VFMiBYHbeDERCKX/tWMjIyEBGRkbI94qKimxr35fBwbRp05DkZMq0R/btk0TE//s/\noGHDyF/frp3s53/uOTlgKNaXf3vumDNH7vCvvtr+tseNk5mdkyer9/dkojVrZNbAjeS5lBRZUiss\nDO9obRPl58tdrtvr43XrSp9BnDnwcqeCpV07qUPjBxXdMGdlZSE5OdmW9rms4KIXXpD11uqW+QWA\nn/9c/pG895594zLR7NlyZ9Sokf1tjxkjWcnRckiQ1hIcpKS401+/fvKYm+tOf17Iy5PtyTEx7vdt\nJSUGjZcFkCxxcbJkFA2crnPQUCmVqJTqW/qtjqV/duGMMn85eVJyDb73PZm6rq7ERCAtTY4njlYn\nTwJLlki+gRO6d5dEx0WLnGnfb3btAg4fdj7fwNK1qwTJ2dnu9OeF/Hz3lxQs3boFMzjwsgCSJT4e\n+OILb8fgFqdnDlIAZANYB6lz8CyALABPONyv70yfLvt07cgVuPdeOd1u69aat2WiJUukYuS4cc60\nrxQwenT0BAfWzgG3goOYGKBPH6nIGFR5ee4nI1q6dpXSzUE64OrUKdkC7vXMQXy8LIedOuXtONzg\ndJ2Dz7TWtbTWMeW+vuNkv3700kvyYWZH5DtxoqzVRuu2xgUL5Do6eWc2erTc2RYWOteHX6xZI2+6\nLVq412e/fsGdOTh9WtalvZo56NpV6qH4JXHODl5vY7RYCbvRMHvAnAMX5OYC69YB3/2uPe3VrQtM\nngy89lp0nrmwaJF8eDuZPDd6tKzFL1niXB9+sXate7MGlr59Zebr9Gl3+3XD9u3yu+PlzAEQrKUF\nrwsgWeJLF8QZHJAtXn1VtonZuUY+ZYrc1c6bZ1+bJjh4UIoUjR7tbD/x8fImG/SlheJiCVzdSka0\n9OsnfQfx+GbrQ9mrmYO4OKBevYvbKYOgoECWo9q29XYc1sxBNCQlMjhw2Nmzkm9wzz01O1K4vJ49\npTBSuW2ugWfdyaelOd9XNOQd5OVJUS23Zw569ZL96kHMO8jLk1omzZt703+tWnJeS5BmDnbvlsDA\n6+3b9evLki6DA6qxWbPkDn/KFPvbTk8HZs6MroqJixZJYORGlbQRI6QE7sGDzvflFasyok1bo8PW\noIHsCgli3oG1U8HNA5fKC9p2Rj/sVLBEy44FBgcOe/ttKRfbo4f9bd9xh6zZzpxpf9t+tXgxMGqU\nO30NHSqPQa53kJ0t+/HLFSR1Rd++wZ058CrfwNK1a/CWFbxORrTEx3PmgGroxAkp1nP77c6036ED\nMGhQ9Cwt7NolX24FB23byjUOcnCQkyMf0l7o10+SdYuLvenfKV7WOLB07iwVWc+c8XYcdvFDASRL\ntBRCYnDgoFmz5M7eriOFK3LbbbK17/hx5/rwi8WLZT115Ej3+hw2TGpKBJHW3gYHffvKfvFt27zp\n3wlffy0FpbyeOejcWf5+g3C66LlzcjQ7lxXcxeDAQe+8I4fMdOjgXB8TJsg/nmjYtbB4sSzRNG3q\nXp/DhskHaBCDr4IC4OjRi+WM3WYFJRs2eNO/E7zeqWDp1Eket2/3dhx22LtXAh2/zBzExwNHjkil\n1iBjcOCQY8eAuXPlzt5JHTpItbkZM5ztx2tWzQG3lhQsQ4cCJSXAypXu9usGa73fq5mD5s1li2+Q\nggNrnb9LF2/H0aaNbGfcscPbcdhh92559EtwEC2FkBgcOGTmTNnG6OSSgmXCBMltOH/e+b68snu3\nTC1aSYJu6dZNKgcGMe8gJ0f+39q08W4MvXsHq9ZBfr4cK1y/vrfjqFVLZg+CEBxYBZCuusrbcVis\nQkhBzztgcOCQjz4CBgxw5xd6wgSZHg7q2jgAZGbK46BB7varlAQkQby22dkya+DllrtevYI3c+B1\nvoElSMFBmzZSGdYPrEJMDA4oYufOAZ9+Ctxwgzv9JSXJVFeQlxYyM+VN14vCMkOHAp9/HryZmZwc\n7/INLL17ywdYUNZv/bBTwdKpUzByDvy0UwGQ5ZoWLbisQNWwYoUksF1/vTv9KQWMHy/Bgdbu9Om2\nFSuAIUO86XvwYNl1sn69N/074fBhORzIq3wDS69e8ju7ZYu347BDcbF8GPtp5mD3bvO3ivqpAJIl\nGmodMDhwwNy5QMuW7t6VjRsn/4iCVPjEUlQkU89eBQd9+0rp61WrvOnfCbm58uh1cNCzpzwGYWlh\nzx7JM/LLzEHnzjLbZfqHmJ8KIFkYHFC1zJkjswa1XLy6I0bImtz8+e716ZZVq+Tu0qvgoF49CfQ+\n/9yb/p2QnS1Jc15/kDVsCHTsGIykRCsw99PMAWB23sGFC/Ih7LfgIC6OywoUob175Y3OrSUFS8OG\nwPDhwax3kJkph514+UE2YECwZg5ycmS9PybG65HIOIIwc5CfLwG6lc3utXbt5O/X5LyD/ftlWYTL\nCu5jcGCzuXNlxuDaa93ve8wYYOlSWR8PkhUrZN3fy6z6gQOlkl9hoXdjsJMfkhEtQdnOmJcn9Q38\nEHABshTWrp3ZMwfWNka/zRzEx8tyZxCLo1kYHNhs7lzZbnf55e73fd11Ukt92TL3+3bKhQtyx+7V\nkoJlwAB5XL3a23HY4cwZSQD0Ot/A0quX1LAwPfDKy/N+maY807cz+q0AkiUaCiExOLDRuXPAwoXA\n2LHe9J+QIL+0Qco7WL9etrl5HRx07CjbKIOwtLB5swRdfgkOeveWR9OXFvLz/ZNvYDF9O2NBgfy7\na9jQ65GEioZCSAwObLR8uZzE6Ha+gUUpmT0IUt5BZiZQpw6QkuLtOJSS2YMgJCVmZ8v/j/Wh7LUu\nXWQK3OSlhZMn5YPCrzMHpm5x9uNOBSA6CiExOLDRnDlSycvLO7IxY2TKeM8e78Zgp8xMIDlZdgx4\nbeBACQ5KSrweSc3k5Mgdrl/uxmrXBnr0MHvmwLo799vMQefOErh8+aXXI6kevxVAstSpI+eCMDig\nsMydK3fuXibOXX21JEQuWODdGOxkJSP6wYABUqba9COGvTym+VJ69TJ75sDaxujHmQPA3LwDPxZA\nssTHA/v2eT0K5zA4sElBgazlepVvYGnaVO60Fy/2dhx22LtXvrzON7CkpsqjyXkHJSVSAMlvwYG1\nY8HU6e/8fNlu26yZ1yMJ1bGjPJqYd1BSIjOgfpw5ACS/izMHVKW5c2UL09VXez0SYPRoCQ5MfaO1\nWIct+WXmoGlTmf42Oe9g1y7ZfpWY6PVIQiUkyDHn+/d7PZLq8eNOBUCWjlq3NnPm4NAhqTjp5+CA\nuxWoSnPmyB1u06ZejwQYNUr+YZlerz4zU9ZMW7XyeiQXmZ6U6JeyyeVZZZQ3b/Z2HNXlx50Kls6d\nzQwOrBoHfl1WYHBAVTp7Vu7UvdqlUN6QIZIws2iR1yOpmcxM/8waWFJSJHHu7FmvR1I9ublyopyf\nAi5APgDq1TMzONDavzMHgLm1DvxaAMkS9EJIDA5s8J//SEaw1/kGlgYNpBCTyXkHJ07IB5lf8g0s\nqalymI2pJzTm5sqSgpdJsxWJiQG6dzczOPjqK/mQ8OvMgam1DnbvBho39sdsbEWsQkhBTUpkcGCD\nOXNk36uR/0qSAAAgAElEQVRf9o0DsrSwdKm5x7VaWwb9Fhz06QPExgJr13o9kupZv95/+QaWhAQz\ngwO/7lSwdOoEfP21eXe4ft6pAAS/SiKDAxv4YQtjeaNGyba77GyvR1I9mZkXEwD9pF49CQJNDA6O\nHZOERD8HB5s2mZdIm58v//Y7d/Z6JBUzdTujXwsgWa68Uh4ZHFCFdu0Ctm71z5KCpX9/yVQ2dWnB\nyjdw89jrcKWkmBkcWEshffp4O45LSUgAjhyRZFqT5OXJh5gfCnVVxNrOuHOnt+OIlF8LIFnq1ZP8\nHQYHVKG5c2Wa2Q9bGMuqUwcYNszM4KC4GFi50n/JiJaUFLnDPXXK65FEZv36i9UI/SghQR5NW1rw\n804FQM4maNTIrJkDrf2/rAAEu9YBg4MamjMHGDpUEmf8ZtQoSZY8d87rkURm0yZZH/VbvoElNVUC\nmJwcr0cSmdxcCQzq1PF6JBXr1EmCF9OCAz/vVABkycO0HQuHD0uSt59nDoBgb2dkcFADZ874awtj\neaNHy92tafvyMzNlNqZ/f69HUrGePYG6dc1bWrB2KvhVbKzcgZsUHFy4IB+6fp45AGRpwaTgwK9H\nNZfH4KAGlFL3K6V2KaVOK6VWKaVSne7TLZ99Bpw+7b98A0tiInD55ebVO8jMBPr1ky2ZflSnjlxb\nk4KD4mKpz+DXfAOLaTsWdu+Wra1+njkAZObApJwDvxdAssTHMzioFqXU7QCeBfAYgH4AcgHMV0o1\nd7Jft8ydK5GjVd3Nb2JigBEjgCVLvB5JZDIz/bukYDEtKXHnTplF8vPMAWBecJCfL49+nzno1Ek+\ncM+f93ok4SkoAOrXl3wJP4uLkyUQ0/KPwuH0zMFUAC9qrV/XWm8FcB+AUwC+43C/rpgzR2YN/LSF\nsbxRo+SgIFN+effvl7sxvyYjWlJSZJeKKXvHrbLJJgQHX30lXybIy5MPMWvPu1917CizR6Yc5W7t\nVPDzeysQ7EJIjgUHSqnaAJIB/HdSW2utASwEMMipft2yfbsc3evXfANLWpokJK5Y4fVIwmON0+8z\nB6mpklGdleX1SMKTmysH8LRs6fVIKmftWDDlXJD8fKBLF39uuS3LqnVgytKCCTsVgGAXQnLyV7o5\ngBgA5XctHwLQ2sF+XTF3rmRWjx7t9Ugq17On7MU1ZWkhM1PeFKwCI37VvbvkRJiytLB+vf/zDQD5\noI2JMWdpwe87FSxXXSXX1ZSkRL8XQLK0bSuPQQwOYr0eQEWmTp2KJk2ahHwvPT0d6enpHo3om+bO\nlToCjRp5PZLKKSWzByYFB36fNQAks75fP3OCg9xc4LbbvB5F1erUkQDBlOAgPx+YPNnrUVQtNlY+\nbE0JDnbvNuP3tUED4IorvKl1kJGRgYyMjJDvFRUV2da+k8HB1wCKAZQ//60VgIOVvfDBB6dh0qQk\np8ZVY6dPy4ftb3/r9UjCM2oUcP/9sj7u52Dm1Ckp92zCmy0geQezZ3s9iqodPSp3Yn7PN7CYkpR4\n4oSsNZswcwCYs2Ph2DH5nTVhWQHwbjtjRTfMWVlZSE5OtqV9x5YVtNbnAawD8N+Jd6WUKv1zpSvg\nfr/LXbpUahz4Pd/AkpYmyUjLl3s9ksqtWSP7xk2YOQAk72D7din562cbNsgjgwN7bdsmj37fqWAx\npRCS349qLi+o2xmdTqN5DsD3lVL3KKW6A/gHgAYA/l3Zi/weHMydK2t4fi1DW16XLrKG7/dSypmZ\nUmmyVy+vRxKelBR59HtSYm6uTNeb8iGWkAAcOOD/oMvvpzGWZxVC8vvBVqYUQLIEtRCSo8GB1vpd\nAD8D8BsA2QD6ABijta50o9Lmzf6+2CZsYSxLKVla8HvQlZkJDBwoiVMm6NJFlmnWrPF6JJXLzZUP\n3Nq1vR5JeEzZsZCfL8m+l1/u9UjC06mTLIV8/bXXI6lcQYH8rrZp4/VIwsPgoJq01n/TWrfXWtfX\nWg/SWleZwhUTA8yY4fTIqmfbNom+TVlSsKSlyR2uX+/GSkrksCVTlhQA2b6WnOz/pES/l00ur2tX\nubZ+X1rIyzNnNgYw5+jmggKZmfX79lBLXJzU5ThzxuuR2MuXlz8lBfj4Y69HUbE5c2SKdtQor0cS\nmVGjZDpx2TKvR1KxLVskcDEpOADkd9XPMwfFxcDGjWZsY7TUry9T4H4PDvLzzVlSAC4e3ez34MDv\nRzWXZ9U62L/f23HYzZfBQVqaJP358S537lxg+HDgssu8Hklk2reXL7/mHSxfLjNGAwZ4PZLIpKZK\n1Tm/VvTbvl1215g0cwBIfQ4/Bwdamzdz0KiRLIP4fceCKQWQLFZwELSjm30ZHIwcKVnrftsmdvKk\nBC1+PWipKn7OO1i+HOjb17ygy0pK9OvsgSllk8vz+46Fgwdla7BJMweAGTsWTCmAZAlqlURfBgct\nWsgd5EcfeT2SUIsWAWfPAuPGeT2S6klLk21tfrzLzcwEhg71ehSR69ABaNbMv8HB+vWyU8XvB9iU\nl5Agd2LHjnk9kopZBy6ZFhz4/ejmU6eAL780Kzi47DKgaVMGB66ZOBGYN0+mRP1i9mygc2fz3hAs\naWnyuHSpp8P4hv37gV27zAwOlPJ33kFurln5BhZrx8LWrd6O41Ly8iRhzkryM4XfCyFZB0OZtKwA\nBHPHgm+Dg5tukijy00+9HonQWpIRTZ01AKQOeNeu/ltayMyUR9OSES2pqRIc+HH/uGk7FSzdu0vg\n5delhbw8mTWqW9frkUSmUycJxv1001WWaQWQLAwOXNStmxQZ8suuhfXr5S/f5OAAkNkDvyUlLl8u\n052m7GsuLzVVpkL99uZw5IhMzZsYHDRoIHePmzZ5PZKKmZaMaLF2LPh19mD3bpmRsQ40MgWDA5fd\ndBMwc6YkJ3pt9mxZWxo+3OuR1MyoUfLG5qdtN6YctnQpfk1KtJIRTVxWAGRpwc/BQffuXo8icn4/\nurmgQD5oTSnYZWFw4LKJE4HCQn+cCTB7NnDNNeZNI5Y3cqQ8+mVp4cQJICfHzHwDy5VXypffgoPs\nbKkZYOIdLiDbGf0YHJw9Kx+uJl7XNm2AevX8m5Ro2k4FS1wccOgQcO6c1yOxj6+Dg5QUmV7yemmh\nsBBYtcr8JQUAaNlSzi7wS3Dw+edSqMfkmQPgYt6Bn2Rny6xBrC8PZq9az56SoHb8uNcjCbVjh1T0\nNDE4UMrfOxZMK4BkiYuTnCM/zcjWlK+DA6VkaeGjj7xN9po3T94MTK1vUF5amn+Cg+XLpTa9KYdY\nXUpqqpRRLinxeiQXZWUB/fp5PYrq69lTHv2WlGgduGRicAD4e8eCaQWQLPHx8hikpQVfBweALC3s\n2SNTz16ZNQtISjI3Ya68UaPkzcHKDPZSZiYweLA5ddQvJTUVKCqSioR+cPq0bANMSvJ6JNXXo4fc\nIPhtaSEvT04PbdXK65FUj18LIZ07J3feps4cAAwOXDV8uNxZelUQ6exZ2cI4YYI3/TthxAh50/V6\n9uDCBTlsyeR8A4vfkhI3bJDlGpNnDho0kO2CfgwOunUz51TW8jp2lLoixcVejyTU3r0yQ9yhg9cj\niVzjxlKemsGBi2rXBm64AfjwQ2/6X7xYqrTdfLM3/Tvh8sulVLHXWxqzsyUhcdgwb8dhhyuukDdd\nv5zQmJUlZ1X06uX1SGrGj0mJpm5jtHTqdPEu3U9275ZHE5cVgODtWPB9cAAAt9wibxBevEl8+KFU\nRbTWP4PCOmfBy1yOpUvl7jA11bsx2MlPSYnZ2bIVsF49r0dSMwwO7OfXo5t375bZGGv93jQMDjww\nZozUrs7IcLff4mJgxgyZNTB1CvFS0tLkF9nLNfIlS2SXQp063o3BTqmpcsfuh7ocWVlm5xtYevaU\n39OiIq9HIr7+Gjh82OzgoH17eT/zY3DQtq257wcMDjxQty4waZIEB27e6WZmyiFFQVpSsAwbJtPO\nXuUdXLgA/Oc/F897CILUVEkE9Dq7/vx5yTkwOd/A4rcdC9ZOBRMLIFnq1pUPMr/tWNi929wlBYDB\ngWfS0+WX2c1p2w8/lOI2QZn2LqtxY0mi8yrvICtL8g2sokxBkJQkuy68XlrYulUSaYMQHHTvLtfU\nL0sLeXly1925s9cjqRk/7lgwPTiIjwcOHPDHzKEdjAkORo4EWrd2b2lBa9khMXGi+dvsLmXUKFn3\n9yLvYMkSoGHDi1n+QXDZZbL9zuvgICtLHvv29XYcdqhfXxI9/RQctGsn4zKZHwshmR4cxMVJnZMD\nB7weiT2M+diLiQFuuw145x13tuCsXCn1FSZNcr4vr6SlScnPLVvc73vpUtnCaFoN9ar44fjm7Gy5\ns23c2Ntx2MVPSYmmJyNa/FYI6dw5YN8+84MDIDhLC8YEB4AsLRw4ACxb5nxfb70lyTEjRjjfl1eG\nDJEPZ7eXFs6fl8qIQVpSsKSmynr/mTPejSE7OxhLChY/BQdbtwYnODh8GDh61OuRCKvGAYMD/zAq\nOBgwQKbDXn/d2X7On5cZivT04C4pALKNcNAgYOFCd/sNYr6BJTVVfn/Wr/em/5ISCQ6CsFPB0rOn\n7Mn3+oPs/HmZig9CcGAd3eyXpQXTaxwAsqOuQQMGB55QCpgyBXj3XSlM5JRPP5UtS3fd5VwffnHd\ndcCiRZLA5pbFi2V9PjnZvT7dkpgoszFeLS1s2yYHFQUtOAC8nz3YtUuSzYIQHPjt6GbTaxwAMv4g\n7VgwKjgAgMmTZcr27bed6+OttySxLDHRuT78YuxYuYt381js+fMl3yFo+QaAbBNLTJTTJr2werU8\nBmmHTbdu/tixYPqBS2VdcYXc6fpp5sDkGgcWBgceiosDrr8eeOUVZ9o/dkx2Kdx1V/AKH1WkTx/Z\nrjlnjjv9HT8u9SPGjHGnPy8MGiQJrV5Yswbo0kVKZAdFvXqSYOl1rYMtW6R+ftu23o7DLn7asWD6\nTgVLXJzkTwSBccEBAHzve3KH5MS6bkaGzEx8+9v2t+1HSsnsgVvBwZIlMjUb5OBg8GCpPPnll+73\nvXo10L+/+/06zQ9JiZs3S0nqoNw0+GnHQlCCg/h4zhx4atw4OS71pZfsbVdr4MUXpX0r8zQajB0r\nWdhuvFHMny93LKYXkanM4MHy6PbswblzkozI4MAZVnAQFH4qhBSU4CAuTpJn/XbiZXUYGRzUrg18\n//vAv/9tbwbzunXy5nrvvfa1aYLRo+WaujF7MH9+sGcNALl7aNsWWLHC3X7Xr5cAIUj5BpaePWUb\n85Ej3vRfUhK84KBjR5kCP3fO23EEocaBJS5OAoNDh7weSc0ZGRwAwP33yy/Vyy/b1+aLL17MaYgm\njRvLWQtOBwc7dshX0IMDpWT2wO2ZgzVrgNjYYFRGLM/rHQt79wInTwbrdNZOnSTosbYReiUINQ4s\nQap1YGxw0Lq11CH485/tqWV96BDwxhvAffdJNcZoM26c5AOcOuVcHwsWyIdXkA5bupRBg+TD2s27\nstWrJcHU9NK+FenaVf5dehUcWMmQQZo58Mt2xiDUOLAwOPCJqVMl6nzvvZq39de/ygfXD35Q87ZM\nNH68JGLOm+dcHzNnygxFUMr6VmbwYLmeOTnu9RnUZERAtoh26eJtcNCwodn78MuLi5PlRK/zDoJQ\n48DSrJn8rjI48FhiohTxefLJmiWAnDgBvPCC5DFccYV94zNJ585y1/n++860f+yYFFu66SZn2veb\nfv3kTcKtvIPjx2WrXRDzDSy9eklpai9s2iS1T4JUMTUmRu7W/RAcBKHGARCsQkjG/6o/+aS8Kb75\nZvXbeOEFeXP98Y/tG5eJbrkFmDXLmXMB5s2T8rMTJtjfth/VqSMf1G4FB2vXyrptUGcOALkZyM31\n5hTRzZuDlW9g8cN2xqDsVLDExwej1oFjwYFS6ldKqUyl1Eml1GGn+klJkWOVH3useuu7hw8DzzwD\n/M//yFGs0WzSJAmSPv3U/rY//lgS5aLpGg8e7F5wkJkJNGkSrDXx8hITZbeC23dlWgdvp4LFD4WQ\nghYccOagarUBvAvg7w72AUBmD/buBf74x8hf+9RTktD4//6f/eMyTUKCTJ1+8IG97Z47B8yeHT1L\nCpbBg2WL1p49zve1fLn0F6Rp7/Kscua5ue72u2+fBM1BDA6smQMvZmMsu3YxOPAjx95KtNZPaK2f\nB+D4KmHPnsBPfiJBwvbt4b8uOxt4/nngkUekqBLJ7MGMGfZm2S9eLDkH0bKkYBkyRB4/+8zZfoqL\nZduk1V9QxcfLeQBuBwdWEmRQg4NTp7zbl3/6tARf1s6JIIiLk/+nkhKvR1IzgbnPeOwx2d44ZYqs\nbVfl7Fkpw5yQAPz8586PzxS33SaFpebOta/Nt96Sw2qi4SCrspo3B3r3BpYudbafjRsl+Bo61Nl+\nvKbUxbwDN23eLNtDg3R3a/H66GZrG6M1jiCIi5PPoK++8nokNROY4KBhQ2D6dLmDeuSRqp//0EPy\npvrqq8HIkrVL795y3O9rr9nT3qlT0XWQVXlpaVI/wkmZmbIlLcg7FSxeBQdB26lg8To4sPoN2swB\nYP7SQkS/7kqpZ5RSJZV8FSulujo12KoMGQI8+6x8/e53FT9Ha1l+ePFF2aWQnOzuGE0weTLwySf2\nRL6ffCJbRdPTa96WiUaOlDXVggLn+li+XAK6Bg2c68MvEhOBbdukWqFbNmwI5k4FQG6qWrf2bsfC\njh1y6mabNt7074SgBAexET7/jwBereI5Nf41mzp1Kpo0aRLyvfT0dKSH8Qnz0EOS0fzII3L++h/+\nINO7AFBYCPz0p3JX/NRTsqxA35SeLtcpIwN48MGatfXmm8CAAcE+aKkyI0bIjMlnnwH33ONMH5mZ\nwK23OtO23yQmSoC/caP8XjmtuFiCgyBfXy93LOzYAXToEKxZmRYtZDba6e2MGRkZyMjICPleUVGR\nbe1HFBxorQsBFNrW+yVMmzYNSUlJ1X7944/L+uADDwDvvCOlbAFZcoiJkeDAqTfqIGjeHLjxRuBf\n/wJ+9KPqLwfs2yfnNfz5z/aOzyRXXCHFpZYsceZ3bu9e2Q0R9GRES8+e8m84N9ed4GDnTlkaC3K+\nTKdOkSVy22nnzmAtKQAS6LRt6/zMQUU3zFlZWUi2aTrcyToH8UqpRADtAMQopRJLvxo61WdZkydL\nssvjj0uGc9OmwBNPyBQvA4Oq3XuvvAEvX179Nl5+WaYM777bvnGZaORI55ISrXyGoCcjWurVk+RW\nt/IOrH769HGnPy94WQhpx47gBQdAMLYzRrqsEInfACj7MZxV+pgGYJmD/f5X8+bAL37hRk/Bc+21\nkoQ1bZqchxCpCxeAl16SwCAazlKoTFqabJl1otjLokVSXKpFC3vb9TM3kxLXr5c1+ZYt3enPCx07\nylbGEyeAyy5zr9+SErlZC9JOBUsQggMn6xxM0VrHVPDlSmBANaOUlJP++OPq3VW8/74sK9x3n/1j\nM83w4TLVuGiRve1qLW2OHm1vu36XmCgf2m7sI8/NDfasAeDd6Yz798uWcs4c+FOA0kDIbnffLWvm\nf/hDZK8rKQGefhoYM0buaqPd5ZfL+ridtSMASbjdty86g4Pjx+Wu02nr1wc73wDwLjgI4jZGixUc\neFl5sqYYHNAlNWgAPPyw5A5EkrA0a5ZkeD/6qHNjM83YsXJmRTgFusK1aJEcM16dZR+TWbnKWVmV\nP6+miopkKSjoMwctW8qWRrd3LOzYITOUQSwuFRcnsyKFjqfvO4fBAVXqgQfkzSPcsyfOn5eAYuTI\n6PvQqsz110sVQzsPYlq0CBg40N11Yj9o2VJKKa9d62w/1vHQQQ8OlPJmO+POnZLVX6+eu/26IT5e\nHk1eWmBwQJWqXx/47W+Bt98GFiyo+vl/+QuQnw/86U/Oj80k/frJh5pdSwvnz8uZFVdfbU97pklN\nBdascbaP9eul8mT37s724wde7FgI6k4F4GIhJJOPbmZwQFWaPBkYNUq2Nx45cunnbd4M/PrXwA9+\nEPx12kjVqiWzB3Pm2NPe8uUy7X3DDfa0Z5qUFGDdOmeTEnNzZcdONJRX79RJKk+6aceOYO5UAORG\nIDaWMwcUcEoBr7wiSWC33VbxiY2FhVJFrkMH4Pe/d3+MJrj+epmqtuNu4pNPgCuvvLj+Hm1SUmSZ\nxsniPevWRU959e7dJb/izBn3+gzyzEFMjPz7ZHBAgde+vWxPXLZM7lbLHvG6fbvMLHz1FfDhh9FR\n4786rrtO7kI//LBm7WgtwcG4cdF5mBVw8UPbqaWFs2dlWSElxZn2/aZ7d5mFcatSYlERcPhwcIMD\nwPztjAwOKGxpacC8eXJH1bkzMGmSlFlOSJCDcBYvlup1VLEmTaS41Hvv1aydvDx5E7/xRnvGZaIr\nrpAPFqeSEjdskLyOaAoOAGDrVnf6s5YwgnzmCoMDiippaZJw+LOfSf7BuXNyAmZODtCrl9ej879b\nb5WDkvbtq34bH34oiaLRVt+gvJQU54KDtWtlzTjoOxUszZtLwOVWcJCXJ49dPTvD13mmBwdOlk+m\ngGrWDHjsMa9HYabx4yUD/oMPqn/iZUYGMGECl29SUmR5pbhY1njttGYN0Lt3MLfZXUr37hc/tJ2W\nny9lqYNcWj0+/mIhJBOX/zhzQOSipk1laaHcSath27hRvsI4vTzwUlPlxMQtW+xve+3a6FlSsHTv\n7t7MQX5+8Jcg4+Lk9/PoUa9HUj0MDohcNnkysGoVsGlT5K/NyJAAY8wY24dlnKQk2SL6+ef2tnvq\nlPzdRMtOBYsVHLhR8jcvL9hLCoD5tQ4YHBC5bPx4WeN95ZXIXnfhAvDGG8AttwB16zozNpM0aiQ5\nAZmZ9rabmytLFdE4c3DihByI5CStZeYgWoIDU/MOGBwQuaxOHeDb3wZefz2yfeWzZ8tdCE+6vGjI\nEPuDg9Wr5e8o2hJs3dqxsH+/7G4K+rJC69Yys8XggIjCdu+9ss/7jTfCf83f/w707x99092VGTpU\n7kK//NK+NjMzJZ8h2mZnOnSQZFmng4P8fHkM+sxBbCzQpg2DAyKKQNeuwMSJUk2yuLjq5+fkAPPn\ny0FYdNGQIfJo1+yB1tKW1W40iY2VugNO71jIz5fdJUEtnVyWydsZGRwQeeSRR6SY0dtvV/3c3/xG\niv5wl0Ko+HjgqqvsCw4KCmTaOxqDA8CdHQt5eRIY1K7tbD9+wOCAiCKWkiL1Cn75S1mDvZTMTOCj\nj+RQq1hWJvmGIUPkICo7WO0MHmxPe6ZxIziIhmREi1XrwEQMDog89NxzcibF//5vxT8/e1YSEPv3\nB771LXfHZoqhQ6Wk96lTNW8rM1M+IJs3r3lbJureXZJeT5xwro9oqHFgiYuT6+nG9lC7MTgg8lDH\njsBTT0mQ8P77oT/TGrj/fnkz/ec/7a8CGBQjRsg2TztmD6I138DSo4c8OjV7cP48sHNn9MwcxMVJ\noHXsmNcjiRyDAyKP/eQnwB13SD7BX/8qswWFhVIs6ZVXgBdfBBITvR6lfyUkSFb4woU1a6ewUKpP\nDhtmz7hMlJAgjxs2ONN+fr4k4FrbJoMuPl4e9+zxdhzVwRVMIo8pJVsamzUDfvQj4Kc/lTvhBg2A\n114D7rnH6xH6m1LA1VcDn35as3YWLZLZmmg+0KphQ5nN2rjRmfatqqA9ezrTvt+0by+Pu3bJWR0m\nYXBA5AOxsTJr8MMfyodUgwZyJHPLll6PzAzXXCMB1ldfAS1aVK+NhQvljtaqbBetevd2buZg0yb5\nnY6WnI7WraVexu7dXo8kcgwOiHwkIeHi1C6Fz7rbX7RIlmgipbXMPNx4o73jMlHv3sDLLzvT9qZN\n0TNrAEiFxPbtZebADatX29cWcw6IyHhXXikfOgsWVO/1O3fK3d3VV9s6LCP17g0cPAh8/bX9bUdb\ncAC4Gxy89559bTE4IKJAGDsWmDUrvIqT5c2bJ0s7I0faPizjWGvjducdnD0LbNsWfcFBhw7uLSsc\nOGBfWwwOiCgQJk6UnIMVKyJ/7ccfA2lpQOPG9o/LNJ07y8FTducdWDsVojE42LXLnVoHBw/a1xaD\nAyIKhAEDZEvjRx9F9rojR4ClS4GbbnJkWMapXVvqHdgdHETbTgVL+/ZS5+DIEWf7OX3a3j4YHBBR\nINSqJR/wH30U2V3arFmydXTCBOfGZpreve1fVti0SbL3r7jC3nb9rkMHeXR6acHuWgoMDogoMCZN\nkjfhVavCf83778usQ9u2jg3LOFZwYOdUeDQmIwIXgwOnkxILCuxtj8EBEQVGWpqc0vivf4X3/C+/\nBObMAe66y9lxmaZ3b+D4cXs/0HJzgT597GvPFM2aSXEpBgdERB6pVQuYMkWOwa7spEvL9Onymjvv\ndH5sJklKksd16+xpr6hItov262dPeyZRyp0dC7t2Aa1a2dcegwMiCpTJkyUweOutyp9XUgK89BIw\nfrzc3dFFrVpJpUi7goOcHHmMxuAAcKfWwa5dUu/DLgwOiChQ2reXxMQ//KHymgezZsnpgz/+sWtD\nM0pysn3BQXY2UK9e9By4VJ61ndFJu3bZW/rbkeBAKdVOKfWyUmqnUuqUUmqbUupxpVRtJ/ojIirr\nV7+Sgjvvvlvxz7WWo7KHDInuI5orYwUHdiQl5uRIHkNslBbst5YVnKx1YMrMQXcACsD3ASQAmArg\nPgBPOdQfEdF/paTI1sSf/1wS68p7802pQ//kk+6PzRTJybJv3o473uzs6F1SAGQ26/RpSYB1wokT\n0radO24cCQ601vO11t/VWi/SWu/WWs8C8EcANzvRHxFRec8/Lx9uDzwQesdWUAA89BBw662yu4Eq\nlpwsjzVdWjh7Fti8ObqDA6e3M1rJjibMHFSkKYDDLvZHRFGsXTvgn/8EXn9djsIuLJTZgtGjpUzy\nP/7h9Qj9zUpKXLu2Zu1s3ChFpvr2tWdcJnI6OLDa9f3MQXlKqc4AHgDAf45E5Jq77pIA4dVXgebN\npZ/Vg6oAAA78SURBVNhR/frAkiXRV6mvOpKTax4crFolJZmjOTho0kR2xGzf7kz7u3YBdevK77hd\nIkoPUUo9A+DhSp6iAfTQWueXeU1bAHMBvKO1Dqs0ydSpU9GkSZOQ76WnpyM9PT2S4RIR4fvfB268\nEVi0SAKCa66J3sS4SA0aBPz2t3LnX91rtnKlLCnUq2fv2EzTpYskydolIyMDGRkZAKT6ZGws8NOf\nFtnWvtIRpE8qpZoBqGpH8E6t9YXS518JYAmAFVrrKWG0nwRg3bp165BkVeEgIiJPLF8ODBsmeQfV\nfUvu3BkYN05yQKLZPfdIcLBypf1t33QTcOYM8PTTWUiWZJFkrXVWTdqMKBbUWhcCKAznuaUzBosB\nrAHwnciHRkREXkpJkeObly+vXnDw5ZfAjh0yAxHtunQB5s51pu2dO4GhQ+1t06k6B1cCWAqgAMAv\nALRUSrVSStlY3JGIiJxUrx6QmgpkZlbv9Z9/Lo8DB9o3JlN17gx8/TVw9Ki97WotOQdW0qNdnEpI\nvAZARwCjAewFsB/AgdJHIiIyxNChMnNQnQI+K1fKMc3t2tk/LtN06SKPduYdALIL58QJQ4IDrfVr\nWuuYcl+1tNYxTvRHRETOGDoU2L+/egcHLVsmFSiVsn1YxnEqONixQx47drS3XZ6tQERElzRkiJxc\nuWRJZK87cUKWFUaPdmZcpmnSBGjRwv7gwGrPCj7swuCAiIgu6fLLJe9g/vzIXvef/8gWyFGjnBmX\nibp0sb/WQX6+LN00amRvuwwOiIioUmPGAJ9+Wvkpl+UtXiwV+7p2dW5cpunc2ZmZAyeuMYMDIiKq\n1Jgxck5FJNUSFy6UWQPmG1xkdyEkQGYOGBwQEZHr+veXNfNw9+nv2SPHNI8b5+y4TNOlC3D4sHzZ\nQWsJNuzONwAYHBARURViY4GxY4EPPwzv+TNmyHkK11/v7LhMY32I5+dX/rxwHTokR5Jz5oCIiDxx\n++3Ahg3Ali1VP3fGDFlSaNzY+XGZpFs3edy61Z72nNqpADA4ICKiMFx3nXzYv/NO5c87cABYuhSY\nONGVYRmlYUOgfXtg82Z72svPl5yOTp3saa8sBgdERFSlunXlA3/6dKCk5NLPe+MNWVK4/Xb3xmaS\nhAT7goNt24CrrnLmxEsGB0REFJb77pOKfPPmVfzzkhLg1VeBm28GmjZ1d2ymsDM4cGqnAsDggIiI\nwjRgAJCcDPzpTxX/fNYsWU//wQ/cHZdJEhKkFPXJkzVva8sWoHv3mrdTEQYHREQUFqWAX/5SCiIt\nXhz6s+Ji4IkngBEj7D8+OEh69pQtiHl5NWvn7FlZVujZ055xlcfggIiIwjZpkhzBfP/9cn6C5S9/\nAbKzgaef9m5sJujRQx5rurSQny8BGYMDIiLynFKSV/DFF5KgmJcHvPwy8LOfAQ89BAwe7PUI/a1R\nIyA+vubBwaZN8uhUcBDrTLNERBRU3bsDM2cCt9xycc17yhTgD3/wdlymsCMpcdMmoE0bORjLCQwO\niIgoYmlpwK5dwPLlcifcu7fXIzJHQgLwySc1a2PTJudmDQAuKxARUTU1bixllRkYRKZPHzm6+fjx\n6rfB4ICIiChA+vWTx9zc6r3+zBkJLhgcEBERBURCAlCnDpCVVb3Xb94sBad69bJ3XGUxOCAiInJR\n7dqyFJOdXb3XZ2UBtWoBiYn2jqssBgdEREQu69evZsFBjx5Agwb2jqksBgdEREQuS0qSpMKzZyN/\nbVaWvN5JDA6IiIhc1q8fcOECsHFjZK+7cEESGRkcEBERBUyfPkBMDLB2bWSv27JFdiswOCAiIgqY\nBg0koXDFishet26dPFrbIZ3C4ICIiMgDgwcDmZmRvWblStkK2aiRM2OyMDggIiLywJAhwI4dwKFD\n4b8mM9OdI7EZHBAREXlgyBB5DHf24PBh2eFgvc5JDA6IiIg8EB8vX+EGBytXyiNnDoiIiAJs+HBg\nyZLwnrtsGdC6NdChg7NjAhgcEBEReea666RS4sGDVT93wQLg6qsBpZwfF4MDIiIij4wZIx/2CxZU\n/ryDB4GcHAkm3MDggIiIyCMtWgDJycDcuZU/b8ECCSKuvdadcTkWHCilZiilCpRSp5VS+5VSryul\n2jjVHxERkYnGjgXmzav8nIUZM4CUFAkm3ODkzMFiALcC6ArgZgCdALznYH9ERETGuf124OjRS88e\nHDsGzJ4tz3OLY8GB1vp5rfVqrfVerfUqAL8DMFApFeNUn0RERKZJSJByyNOnV/zzjz6SWYXbbnNv\nTK7kHCilrgBwF4BMrXWxG30SERGZ4p57gJkzgS++CP2+1sALLwCjR0tNBLc4GhwopX6nlDoB4GsA\n8QBucrI/IiIiE333u0DDhsCzz4Z+f9kyYM0aYOpUd8ejtNbhP1mpZwA8XMlTNIAeWuv80udfAeAK\nAO0APAbgmNb6hkraTwKwbvjw4WjSpEnIz9LT05Genh72WImIiEzyxBPA00/LyYu9egHnzwMDBgC1\nagGrV8ujJSMjAxkZGSGvLyoqwrJlywAgWWudVZOxRBocNAPQrIqn7dRaX6jgtW0B7AUwSGv9+SXa\nTwKwbt26dUhy+rBqIiIiHzlzBkhKkqDghReAV18FPvgAWL4c6N+/6tdnZWUhOTkZsCE4iI3kyVrr\nQgCF1ezLSkSsW83XExERBVa9esCsWcANN0hxpAYNJEkxnMDAbhEFB+FSSvUHkApgOYAjADoD+A2A\nbQBWOtEnERGR6Tp2BNavBzZvBtq1A8qtsLvGqYTEU5DaBgsBbAXwEoAcACO11ucd6pOIiMh4sbFA\nnz7eBQaAQzMHWuuNAEY70TYRERE5i2crEBERUQgGB0RERBSCwQERERGFYHBAREREIRgcEBERUQgG\nB0RERBSCwQERERGFYHBAREREIRgcEBERUQgGB0RERBSCwQERERGFYHBAREREIRgcEBERUQgGB0RE\nRBSCwQERERGFYHBAREREIRgcEBERUQgGB0RERBSCwQERERGFYHBAREREIRgcEBERUQgGB0RERBSC\nwQERERGFYHBAREREIRgcEBERUQgGB0RERBSCwQERERGFYHBAREREIRgcEBERUQgGB0RERBSCwQER\nERGFYHBAREREIRgcEDIyMrweQtThNXcfr7n7eM3N5XhwoJSqo5TKUUqVKKX6ON0fRY7/gN3Ha+4+\nXnP38Zqby42Zg98D+AKAdqEvIiIiqiFHgwOl1PUArgHwMwDKyb6IiIjIHrFONayUagXgnwDGAzjt\nVD9ERERkL8eCAwCvAvib1jpbKdUuzNfUA4AtW7Y4Nyr6hqKiImRlZXk9jKjCa+4+XnP38Zq7q8xn\nZ72atqW0Dj8VQCn1DICHK3mKBtADwHUAbgEwUmtdopRqD2AngL5a6/WVtH8ngDfDHhARERGVd5fW\n+q2aNBBpcNAMQLMqnrYLwLsAbij3/RgAFwC8qbWeUkn7YwDsBnAm7IERERFRPQDtAczXWhfWpKGI\ngoOwG1UqDkDjMt+6EsB8AJMArNZa77e9UyIiIrKFIzkHWusvyv5ZKXUSslthJwMDIiIif3OzQiLr\nHBARERnAkWUFIiIiMhfPViAiIqIQDA6IiIgohG+CA6XU/UqpXUqp00qpVUqpVK/HFCRKqWFKqZlK\nqX2lh2CNr+A5v1FK7VdKnVJKfaqU6uzFWINAKfWIUmq1UuqYUuqQUuojpVTXCp7Ha24TpdR9Sqlc\npVRR6dcKpdR15Z7D6+0gpdQvS99fniv3fV53myilHiu9xmW/Npd7To2vty+CA6XU7QCeBfAYgH4A\ncgHMV0o193RgwdIQQA6AH6KC5FCl1MMAHgBwL4D+AE5C/g7quDnIABkG4C8ABgC4GkBtAAuUUvWt\nJ/Ca224vpEhbEoBkAIsBzFBK9QB4vZ1WekN3L+T9u+z3ed3ttxFAKwCtS7+GWj+w7XprrT3/ArAK\nwPNl/qwgJzn+wuuxBfELQAmA8eW+tx/A1DJ/bgw5E+M2r8cbhC8AzUuv+1Bec1eveyGAKbzejl/n\nywDkARgFYAmA58r8jNfd3mv9GICsSn5uy/X2fOZAKVUbEuUvsr6n5f9oIYBBXo0rmiilOkCiz7J/\nB8cAfA7+HdilKWTG5jDAa+40pVQtpdQdABoAWMHr7bgXAHyitV5c9pu87o7pUrpEvEMpNV0pFQ/Y\ne72dPHgpXM0hpZUPlfv+IQDd3B9OVGoN+eCq6O+gtfvDCRallALwJwDLtdbW2iCvuQOUUr0ArISU\nkT0OYKLWOk8pNQi83o4oDcL6Akip4Mf8PbffKgCTITM1bQA8DmBZ6e++bdfbD8EBUdD9DUACgCFe\nDyQKbAWQCKAJ5PC315VSw70dUnCVlsr/E4CrtdbnvR5PNNBazy/zx41KqdUACgDcBvn9t4XnywoA\nvgZQDEmuKKsVgIPuDycqHYTkefDvwGZKqb8CGAs5ofRAmR/xmjtAa31Ba71Ta52ttX4Ukhz3EHi9\nnZIMoAWALKXUeaXUeQAjADyklDoHuWPldXeQ1roIQD6AzrDx99zz4KA02lwHYLT1vdJp2NEAVng1\nrmiitd4F+cUp+3fQGJJpz7+DaioNDCYASNNa7yn7M15z19QCUJfX2zELAfSGLCskln6tBTAdQKLW\neid43R2llLoMEhjst/P33C/LCs8B+LdSah2A1QCmQhKJ/u3loIJEKdUQ8gukSr/VUSmVCOCw1nov\nZGrw10qp7ZAjs5+E7BiZ4cFwjaeU+huAdADjAZxUSlmRfJHW2jqOnNfcRkqppwHMBbAHQCMAd0Hu\nYq8tfQqvt8201icBlN9jfxJAodZ6S+m3eN1tpJT6A4BPIEsJbQE8AeA8gLdLn2LL9fZFcKC1fre0\npsFvINMfOQDGaK2/8nZkgZIC2WKkS7+eLf3+awC+o7X+vVKqAYAXIZn1/wFwvdb6nBeDDYD7INd5\nabnvTwHwOgDwmtuuJeT3uQ2AIgDrAVxrZdDzersmpI4Kr7vt4gC8BaAZgK8ALAcwUGtdCNh3vXnw\nEhEREYXwPOeAiIiI/IXBAREREYVgcEBEREQhGBwQERFRCAYHREREFILBAREREYVgcEBEREQhGBwQ\nERFRCAYHREREFILBAREREYVgcEBEREQh/j9j+FYpmavPcgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f39de278f98>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"amp = 1e-8\n",
"alpha = 0.023\n",
"freq = 0.68\n",
"phi = np.pi/2\n",
"y = amp * np.exp(alpha * t) * np.cos(freq * t + phi)\n",
"\n",
"plt.figure()\n",
"plt.plot(t, y)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example 3: three oscillatory modes"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f39de1f4a90>]"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAFkCAYAAAAt0UHnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl41OXV//H3Ya2o4E4UUUQtYBU1kc0FqQioqKBANWpx\nQa0Fl8b66K+L4tJqxSqgdauo4BaLgLigUoIbyqYJUFzADSrKIiCCCsh2fn/ckzrETNaZ+SaTz+u6\n5uLJ5Hzv+2QeKif3au6OiIiISCapF3UCIiIiIsmmAkdEREQyjgocERERyTgqcERERCTjqMARERGR\njKMCR0RERDKOChwRERHJOCpwREREJOOowBEREZGMowJHREREMk7KCxwzG2Jmi8xsg5nNNLMO5cR3\nM7NCM9toZh+Z2fmlxAwwsw9jbc4zs5Mr26+ZDY218Z2ZfW1mU8ysY4mYxmZ2r5mtMrNvzWycme1V\n1c9CRERE0iOlBY6ZnQXcCQwFjgTmAZPNbI8E8a2AF4GpwOHASGCUmfWIizkaeAp4CDgCeA6YaGaH\nVLLfhcAQ4FDgGGAx8G8z2z0uZgTQG+gHdAX2AcZX+oMQERGRtLJUXrZpZjOBWe5+VexrA5YAd7v7\nsFLibwdOdvf2ce/lA83c/ZTY108DTdz99LiYGcAcdx9clX5jMTsDa4Hu7v6amTUFVgJnu/uzsZg2\nwIdAZ3efXZ3PRkRERFInZSM4ZtYQyCGMxgDgoZoqALokeKxz7PvxJpeI71JWTFX6jT3zG+AbwmgP\nsTYalGhnIfB5GfmLiIhIDdAghW3vAdQHVpR4fwXQJsEzWQnim5pZY3f/oYyYrMr2a2a9gaeBJsBS\noIe7fx2XyyZ3X1dGXz8Rm+LqRZjy2pgoTkRERH7iZ0ArYLK7r65OQ6kscGqDVwlrffYALgGeMbOO\n7r6qGm32Ap5MRnIiIiJ11LmE9bZVlsoCZxWwFWhe4v3mwPIEzyxPEL8uNnpTVkxxmxXu1903AJ/F\nXrPN7CNgEHB7LLaRmTUtMYpTVv4QRm544oknaNeuXRlhkkx5eXkMHz486jTqFH3m6afPPP30mafX\nhx9+yHnnnQexf0urI2UFjrtvNrNCoDvwPPxvsW934O4Ej80ASm757hl7Pz6mZBs9imOq2G+xekDj\n2P9dCGyJPRe/yHi/EvmUtBGgXbt2ZGdnl9OdJEuzZs30eaeZPvP002eefvrMI1PtJR6pnqK6Cxgd\nKzhmA3mE9S6jAczsNmAfdy8+6+YBYEhsN9UjhOKiP3BKXJsjgdfN7GpgEpBLWBB8SSX6bQL8iVAA\nLSNMUV1O2Ab+DIC7rzOzh4G7zGwN8C2hQHpbO6hERERqtpQWOO4+Nnb2zM2EqZ25QC93XxkLyQJa\nxsUvji38HQ5cCXwBDHL3griYGWZ2DvDX2OtjoI+7f1CJfrcCbYGBhOJmNfAOcKy7fxj3I+TFYscR\nRnZeIZydIyIiIjVYyhcZu/t9wH0JvndhKe+9SRiRKavN8ZRz4F45/f5AOLyvTLG4K2IvERERqSV0\nF5VkhNzc3KhTqHP0maefPvP002dee6X0JOO6yMyygcLCwkItTBMREamEoqIicnJyAHLcvag6bWkE\nR0RERDKOChwRERHJOCpwREREJOOowBEREZGMowJHREREMo4KHBEREck4KnBEREQk46jAERERkYyj\nAkdEREQyjgocERERyTgqcERERCTjqMARERGRjKMCR0RERDKOChwRERHJOCpwREREJOOowBEREZGM\nowJHREREMo4KHBEREfmfZcvgiivg4IOhVSu44AL45JOos6o8FTgiIiICwOuvwyGHwNNPw0knwYAB\n8OqrcOSR8OKLUWdXOSpwREREhOnTQ1Fz1FGwcCHccw/ccQd88AGceCL06wdvvhl1lhWnAkdERKSO\n+/JLOPNM6NgxjNTsttuP39tpJ/jXv+DooyE3F1avji7PylCBIyIiUoe5w2WXQb16MG4cNG7805hG\njeCJJ2DjRrjmmvTnWBUqcEREROqwsWPDqM3998NeeyWOa9EC/vIXGDMGiorSl19VqcARERGpozZs\ngKuvDtNTffqUH3/JJdCuHfzxj6nPrbpU4IiIiNRR994LX30Ft99esfgGDeDPf4bJk2HOnNTmVl0q\ncEREROqgtWvhtttg0CA46KCKPzdgALRuXfGiKCopL3DMbIiZLTKzDWY208w6lBPfzcwKzWyjmX1k\nZueXEjPAzD6MtTnPzE6uTL9m1sDMbjez/5jZd2b2pZmNMbO9S7Txuplti3ttNbP7qvN5iIiI1AT3\n3Qfffw/XX1+55xo0gN/9DsaPhxUrUpNbMqS0wDGzs4A7gaHAkcA8YLKZ7ZEgvhXwIjAVOBwYCYwy\nsx5xMUcDTwEPAUcAzwETzeyQSvTbJPbsTbHvnwG0ibUVz4F/As2BLGBv4NpKfxAiIiI1yKZN4Zyb\ngQPD4uHKOu+8UOiMGZP83JIl1SM4ecCD7v6Yuy8ALgPWAxcliP8t8Jm7X+vuC939XmBcrJ1iVwIv\nu/tdsZgbgCLg8or26+7r3L2Xu49394/dfXbs+Rwz27dETuvdfaW7fxV7fVedD0RERCRq+fnhSoar\nr67a87vuCv37w6hRYZt5TZSyAsfMGgI5hNEYANzdgQKgS4LHOse+H29yifguZcVUsV+AXQgjNt+U\neP9cM1tpZvPN7FYz26GMNkRERGo0d7jzTujdG9q2rXo7F18MH38Mb72VvNySKZUjOHsA9YGSM3Qr\nCNM9pclKEN/UzBqXE1PcZqX7jbX9N+CpEiM0TwLnAd2AW4FfA48nyF1ERKTGe/NNmD+/6qM3xY47\nDvbdN5xyXBM1iDqBqJlZA+AZwujN4PjvufuouC/fN7NlwFQzO8DdF5XVbl5eHs2aNdvuvdzcXHJz\nc5OTuIiISBWMHg0HHgi//GX12qlXL+yoeuopGDkS6tev3PP5+fnk5+dv997atWurl1ScVBY4q4Ct\nhAW68ZoDyxM8szxB/Dp3/6GcmOI2K9xvXHHTEjihAutrZgMGHASUWeAMHz6c7OzscpoTERFJn+++\ng2eegeuuA7Pqt3fWWTB8OEybBt26Ve7Z0n7pLyoqIicnp/qJkcIpKnffDBQC3YvfMzOLfT09wWMz\n4uNjesbeLyumR3FMRfuNK25aA93dfU0FfqwjCSM9yyoQKyIiUqNMmBC2hv/618lpr2NH2G+/mjlN\nlepdVHcBl5jZQDNrCzxA2KI9GsDMbjOz+E1mDwCtY2fUtDGzwUD/WDvFRgInmdnVsZgbCYuK/1GJ\nfhsA44FswhqbhmbWPPZqGItpbWZ/NrNsM9vfzE4HxgBvuPt7yfuIRERE0mP06DA11apVctozC9c8\nvPBCzdtNldICx93HAtcANwNzgPZAL3dfGQvJIkwPFccvBnoDJwJzCdu9B7l7QVzMDOAc4NJYzJlA\nH3f/oBL9tgBOBfaNtbGUMCqzlB93Wm2K5TEZ+BC4gzDic3r1PhUREZH0++9/4bXX4PyfHJ9bPb17\nw5dfwrx5yW23ulK+yNjd7wNKPf3X3S8s5b03CSMyZbU5njACU9V+/0vYaVXW818Qdk+JiIjUeo8/\nDjvuCP36Jbfdrl1hp53CjeRHHJHctqtDd1GJiIhkOPdw6nD//qEYSaZGjaBnT5g0KbntVpcKHBER\nkQw3fTp88knyp6eK9e4Ns2bBypXlx6aLChwREZEMN2YM7L8/HH98ato/5ZQwSvTKK6lpvypU4IiI\niGSwDRvCNu6BA8PhfKmQlQWHHw5Tp5Yfmy4qcERERDLYxImwbl0ocFLphBNCgVNTtourwBEREclg\nY8bAMcfAQQeltp/u3eGLL8IFnDWBChwREZEM9eWXMGUKXHBB6vvq2hUaNIBXX019XxWhAkdERCRD\nPfFE2MY9YEDq+9p553B1Q01Zh6MCR0REJAO5h6sZzjgDmjVLT5/du4fTkrdtS09/ZVGBIyIikoHe\neQcWLEjP9FSxE06A1atrxrUNKnBEREQy0Jgx0KJFGFVJly5doHFjeOON9PWZiAocERGRDPPDD5Cf\nD7/+NdQv8+bF5GrcOKzDeeut9PWZiAocERGRDPPCC7BmTequZijLcceFAifq83BU4IiIiFSCOyxZ\nAosXw5YtUWdTujFjwkhK27bp7/vYY2HFCvj00/T3HU8FjoiISAX88AP89a9hXct++8EBB8Cee8Ll\nl8PSpVFn96Ply+Hll6MZvYGwDscMpk2Lpv9iKnBERETK8eWX0KkT3HQT9OkTpoBefhkGDw5rXQ49\nFCZMiDrL4PHHw4F7ubnR9L/LLnDYYdGvw2kQbfciIiI127Jl4ZTezZth9mw44ogfv3fSSfD738Ol\nl0K/fjByJFx5ZXS5usOjj4azb3bdNbo8jj0WCgqi6x80giMiIpLQpk3Qvz9s3BimXOKLm2K77QbP\nPBMKnauuggceSH+exWbNgg8/hIsuii4HCAuNP/oIvvoquhw0giMiIpLA0KHhwLw334T9908cZwZ3\n3BHW6QwZEtbpnHZa+vIs9uij0LJlOHAvSsceG/58++0wmhQFjeCIiIiUYu7cULTccAN07lx+vBmM\nGAGnnw7nnguffJL6HOOtXx/WA11wQXrPvinNvvuGgjDKhcYqcEREREpwD7uj2rWDa6+t+HP164ct\n2llZ4YLLDRtSl2NJEybAt9+m92qGshxzDMyYEV3/KnBERERKmDQpTK/ceWe4jbsymjYNa3IWLIBr\nrklNfqV54AH45S+hdev09VmWTp2gqChM20VBBY6IiEicbdvgT3+C44+HHj2q1sbhh4fprfvugylT\nkptfaebMCQXZ5Zenvq+K6tw5LNKO6uJNFTgiIiJxnn8e/vOfcKifWdXbGTw4LPa96CJYuzZ5+ZXm\n3nvDupfTT09tP5VxxBHhbqqZM6PpXwWOiIhInOHDw/qRY46pXjv16sEjj4Ti5ne/S05upfn6a3jy\nSbjssnDAX03RqBEceWTYuh4FFTgiIiIxhYVhS3heXnLa23//sLNq9Ohw+nEqjBoVptUuuSQ17VdH\n584awREREYnciBHQqhX07Zu8Ni+8EHr3DgXIqlXJaxfCAYTDh8N558FeeyW37WTo1Ak++wxWrkx/\n3ykvcMxsiJktMrMNZjbTzDqUE9/NzArNbKOZfWRmP7kuzMwGmNmHsTbnmdnJlenXzBqY2e1m9h8z\n+87MvjSzMWa2d4k2GpvZvWa2ysy+NbNxZlYD/wqJiEh1rVkTdj/99rfJPUfGDB56KCy4HTIkee1C\n2JK+YkXltrKnU/H5QVFMU6W0wDGzs4A7gaHAkcA8YLKZ7ZEgvhXwIjAVOBwYCYwysx5xMUcDTwEP\nAUcAzwETzeyQSvTbJPbsTbHvnwG0ibUVbwTQG+gHdAX2AcZX+oMQEZEaLz8ftmyBgQOT3/bee4cd\nVWPHwr/+lZw2t2yBYcPCVRJt2iSnzWTbf39o3jyaaSpz99Q1bjYTmOXuV8W+NmAJcLe7Dysl/nbg\nZHdvH/dePtDM3U+Jff000MTdT4+LmQHMcffBVek3FnMUMAvY392/MLOmwErgbHd/NhbTBvgQ6Ozu\nsxO0kw0UFhYWkp2dXeHPSkREotWhQyhEnn8+Ne27w1lnwdSp8N57oa/qGDUqTHsVFYXFvDVVnz7h\nlOWKbJcvKioiJycHIMfdi6rTb8pGcMysIZBDGI0BwEM1VQB0SfBY59j3400uEd+lrJgq9guwC+DA\nN7Gvcwh3dcW3sxD4vJx2RESklpk/H959N7WXVJqFUZwGDcLt49UZX/j++3CFRG5uzS5uIExTzZ4d\nFkKnUyqnqPYA6gMrSry/AshK8ExWgvimZta4nJjiNivdb6ztvwFPuft3cf1scvd1lchfRERqoTFj\nYM89w2LgVNpjj7Ae58UX4eGHq97O3/8Oq1eHs3pqus6dYd26cLJzOtX5XVRm1gB4hjB6MzjidERE\nJM3cw9qYAQOgYcPU93f66XDxxeHU4cLCyj+/YAHceitcfTUccEDy80u2o44Ko1fpXoeTyiOBVgFb\ngeYl3m8OLE/wzPIE8evc/YdyYorbrHC/ccVNS+CEuNGb4n4amVnTEqM4ZeX/P3l5eTRr1my793Jz\nc8nNzS3vURERSaNZs2DJklDgpMs994QrDM48M0yN7blnxZ7bsgUGDQpb2W+4IaUpJs3OO8Ohh4YC\nJ34KMD8/n/z8/O1i1ybxyOeUFTjuvtnMCoHuwPPwv8W+3YG7Ezw2Ayi55btn7P34mJJt9CiOqWi/\nccVNa+CX7r6mRL+FwJbYc/GLjPcrkU+phg8frkXGIiK1wDPPhJ0+xx2Xvj5/9jMYPx5ycuDUU8MC\n3KZNy3/uj38MBdnrr8MOO6Q8zaTp1OmnW8VL+6U/bpFxtaV6iuou4BIzG2hmbYEHCFu0RwOY2W1m\nNiYu/gGgdeyMmjZmNhjoH2un2EjgJDO7OhZzI2FB8D8q0W8DwnbvbOA8oKGZNY+9GgLERm0eBu6K\nnc2TAzwCvJ1oB5WIiNQu7qHA6dcvuWffVETLlvDKK7BwYShyvvmm7PgHHwwXeA4bBscem54ck6VT\np7Bz7Lvvyo9NlpQWOO4+FrgGuBmYA7QHerl78ZmGWYTpoeL4xYRzZ04E5gJ5wCB3L4iLmQGcA1wa\nizkT6OPuH1Si3xbAqcC+sTaWAstif8bvkMojnMszDng99v1+Vf9ERESkJolieipedja8/HLYxdW5\nc7jks6Rt28Ji4ssugyuuSN41EunUqVP4Oaqy5qiqUnoOTl2kc3BERGqPa66BJ56AL79M/whOvI8/\nhjPOCAuIzzsvjCjttRd88EG4KbyoCK6/Hm68sXo3nEdl61bYZZfwM5R16nIyz8GpQfeOioiIpNcL\nL4RdTVEWNwAHHxyKmH/8I7zGxBZvmMHxx8Mbb6R3jVCy1a8fdlPNTuMCjzq/TVxEROqmjz4Kr9NO\nizqToFGjsPX700/hv/+FuXPD5ZyvvVa7i5tiHTum904qFTgiIlInvfhi2M3UvXvUmWzPDPbbDw4/\nHHbbLepskqdTJ/jiC1i6ND39qcAREZE66YUX4IQToEmTqDOpGzp1Cn+maxRHBY6IiNQ533wD06bV\nnOmpuqBFi/BSgSMiIpIikyeHnT2pvntKtlfagX+pogJHRETqnBdegCOOCIftSfp06hSupti6NfV9\nqcAREZE6ZevWcLjeqadGnUnd06lTOM34ww9T35cKHBERqVOKiuDrr+Gkk6LOpO7JyYF69dIzTaUC\nR0RE6pSCgnDDdceOUWdS9+y0E/ziFypwREREkm7KFOjWDRo2jDqTuildC41V4IiISJ2xfj28/Tac\neGLUmdRd6bpZXAWOiIjUGW+9BZs2qcCJUrpuFleBIyIidUZBAey9N7RrF3Umddchh4S1OKm+eFMF\njoiI1BkFBWH0xizqTOqu4pvFU70ORwWOiIjUCatWwZw5mp6qCdJxs7gKHBERqRNefTX8WdNuD6+L\n0nGzuAocERGpEwoKwtqbFi2izkTScbO4ChwREcl47uH8G01P1QzpuFlcBY6IiGS8zz6DxYuhR4+o\nM5FiqT7wTwWOiIhkvIKCsHvn+OOjzkSKpfpmcRU4IiKS8QoKwj+oTZtGnYkUS/XN4ipwREQko23d\nGnZQaf1NzZLqm8VV4IiISEabOxe+/loFTk2T6pvFVeCIiEhGKyiAHXf8cWuy1BypXGisAkdERDJa\nQUFYXNyoUdSZSEmpvFlcBY6IiGSsDRtg2jRNT9VUqbxZXAWOiIhkrLffhh9+UIFTU6XyZvGUFzhm\nNsTMFpnZBjObaWYdyonvZmaFZrbRzD4ys/NLiRlgZh/G2pxnZidXtl8zO8PMJpvZKjPbZmbtS2nj\n9dj3il9bzey+qnwOIiKSfgUFsNdecOihUWcipUnlzeIpLXDM7CzgTmAocCQwD5hsZnskiG8FvAhM\nBQ4HRgKjzKxHXMzRwFPAQ8ARwHPARDM7pJL97ghMA64FPMGP4MA/geZAFrB3LF5ERGqBgoIwemMW\ndSaSSKoWGqd6BCcPeNDdH3P3BcBlwHrgogTxvwU+c/dr3X2hu98LjIu1U+xK4GV3vysWcwNQBFxe\nmX7d/Ql3/wuhmCrrr/56d1/p7l/FXilYCiUiIsm2ejUUFel6hpquY8dws/iXXya33ZQVOGbWEMgh\nFBAAuLsDBUCXBI91jn0/3uQS8V3Kiqliv2U518xWmtl8M7vVzHaoQhsiIpJmr70WLtns3j3qTKQs\nRx8d/nz77eS22yC5zW1nD6A+sKLE+yuANgmeyUoQ39TMGrv7D2XEZFWj30SeBP4LLAXaA8OAnwP9\nK9mOiIikWUEBtGkDLVtGnYmUJSsLDj447HY76KDktZvKAqfWc/dRcV++b2bLgKlmdoC7L4oqLxER\nKV9BAZx0UtRZSEUcdxy8+SZceGHy2kxlgbMK2EpYoBuvObA8wTPLE8Svi43elBVT3GZV+q2o2YT1\nOgcBZRY4eXl5NGvWbLv3cnNzyc3NrWYKIiJSnkWL4NNPNT1Vk+Xn55Ofnw/AkiXwn//AFVesTVr7\nKStw3H2zmRUC3YHnAczMYl/fneCxGUDJLd89Y+/Hx5Rso0dxTBX7TbSLqqQjY7HLygscPnw42dnZ\nFWxWRESSacqUcJHjL38ZdSaSSPwv/Z99BgceCAMGFDF9ek5S2k/1FNVdwOhYwTGbsLupCTAawMxu\nA/Zx9+Kzbh4AhpjZ7cAjhKKkP3BKXJsjgdfN7GpgEpBLWFR8SUX7jfW9K7Af0IIwKtM2Vggtd/cV\nZtYaOAd4CVhN2LZ+F/CGu79X7U9GRERSpqAg7M7ZZZeoM5GKOOAAaNEC5sxJXpspLXDcfWzs7Jmb\nCVNEc4Fe7r4yFpIFtIyLX2xmvYHhhO3gXwCD3L0gLmaGmZ0D/DX2+hjo4+4fVKJfgNOBRwkjMg7k\nx96/KfbcJuBE4CrCmTlLgGdifYqISA21dStMnQpDhkSdiVSUWViHk8wCx8IOakkWM8sGCgsLCzVF\nJSISgXffhQ4d4I03oGvXqLORirr/frjiiiK2bs0ByHH3ouq0p7uoREQko0yZAjvuCJ07R52JVEbX\nrmH0LVlU4IiISEYpKIBu3aBRo6gzkcpo1w6aNk1eeypwREQkY6xfD2+9pesZaqN69WDYsCS2l7ym\nREREojVtGmzapAKnturQIXltqcAREZGMMWUK7LNPmO6Quk0FjoiIZIwpU+DEE8O2Y6nbVOCIiEhG\nWLEiHPev6SkBFTgiIpIhpk4Nf554YrR5SM2gAkdERDLCyy9D+/aQlRV1JlITqMAREZFab+tWeOUV\n6N076kykplCBIyIitd6778KqVSpw5EcqcEREpNZ76SXYdVfo1CnqTKSmUIEjIiK13qRJ0KsXNGgQ\ndSZSU6jAERGRWm35cigshFNOiToTqUlU4IiISK32yivhYL9evaLORGoSFTgiIlKrvfRSuMNor72i\nzkRqEhU4IiJSa/3wA0yerN1T8lMqcEREpNZ69VVYtw7OOCPqTKSmUYEjIiK11oQJcNBBcOihUWci\nNY0KHBERqZW2boXnngujN7o9XEpSgSMiIrXS9OmwciWceWbUmUhNpAJHRERqpQkTYO+9oWPHqDOR\nmkgFjoiI1DrbtoUCp29fqKd/yaQU+mshIiK1zvTp8PnncPbZUWciNZUKHBERqXWefBJatoRjj406\nE6mpVOCIiEitsmkTjB0L55yj6SlJTH81RESkVpk8Gb7+OhQ4IomowBERkVrlySfDwX7t20edidRk\nKS9wzGyImS0ysw1mNtPMOpQT383MCs1so5l9ZGbnlxIzwMw+jLU5z8xOrmy/ZnaGmU02s1Vmts3M\nfvI/FTNrbGb3xmK+NbNxZqbr3EREIrJqFTz7LFxwQdSZSE2X0gLHzM4C7gSGAkcC84DJZrZHgvhW\nwIvAVOBwYCQwysx6xMUcDTwFPAQcATwHTDSzQyrZ747ANOBawBP8CCOA3kA/oCuwDzC+oj+/iIgk\n1+jR4c/zf/Krr8j2zD3Rv+1JaNxsJjDL3a+KfW3AEuBudx9WSvztwMnu3j7uvXygmbufEvv6aaCJ\nu58eFzMDmOPugyvbr5ntDywCjnD3/8S93xRYCZzt7s/G3msDfAh0dvfZCX7mbKCwsLCQ7OzsSnxa\nIiJSFndo0wY6dAjTVJJ5ioqKyMnJAchx96LqtJWyERwzawjkEEZjAPBQTRUAXRI81jn2/XiTS8R3\nKSumiv2WJgdoUKKdhcDnlWxHRESS4LXX4OOP4bLLos5EaoNUTlHtAdQHVpR4fwWQleCZrATxTc2s\ncTkxxW1Wpd9EuWxy93XVbEdERJLgzjvhsMN09o1UTIOoE8hUeXl5NGvWbLv3cnNzyc3NjSgjEZHa\na/58eOklePxx3RyeKfLz88nPz9/uvbVr1yat/VQWOKuArUDzEu83B5YneGZ5gvh17v5DOTHFbVal\n30S5NDKzpiVGcSrUzvDhw7UGR0QkSYYNg/32g7POijoTSZbSfumPW4NTbSmbonL3zUAh0L34vdhi\n3+7A9ASPzYiPj+kZe7+smB7FMVXst7SV1oXAlhLttAH2K5GPiIik0EcfQX4+XH01NGwYdTZSW6R6\niuouYLSZFQKzgTygCTAawMxuA/Zx9+INfw8AQ2K7qR4hFBf9gVPi2hwJvG5mVwOTgFzCguBLKtpv\nrO9dCcVKC8CAtrFCaLm7r3D3dWb2MHCXma0BvgXuBt5OtINKRESS749/hBYt4De/iToTqU1SWuC4\n+9jY2TM3E6Z25gK93H1lLCQLaBkXv9jMegPDgSuBL4BB7l4QFzPDzM4B/hp7fQz0cfcPKtEvwOnA\no4TRGweKJwJvij0HoTDaCowDGgOvAEOq9aGIiEiFzZwJ48fDmDHws59FnY3UJik9B6cu0jk4IiLJ\nsXlzOPOmXj145x2oXz/qjCTVknkOjnZRiYhIjTRsGLz3HsyereJGKk8FjoiIVMnmzTBvHixbFgqQ\nli2hXTtokIR/WaZNgxtvhGuvBQ2GS1WowBERkUr59FO4/XZ46in4/vvtv7fTTtCjB5x3HvTuDY0b\nl95GWRaakr4TAAAgAElEQVQtggED4Jhj4KabkpOz1D0pv01cREQygzsMHw6HHgovvAD/938wYwYs\nXQpLlsAbb4QdT59/Dv36wT77QF5emGaqqE8+gW7dYOed4V//0rZwqToVOCIiUq4tW+CSS8JZNJdd\nFkZxhg6Fzp1h771h332ha1f4wx/g3Xfhgw9g0KBwKeZhh4W4UaPg229Lb989jAgddVQY9XntNWhe\n8rhWkUpQgSMiImVyh9/+NmzVfuyxMIrTpEnZz7RrFxYJf/FF2Oa9225w6aWhGBo0KLQ1dSpMmgS3\n3QZHHAHnngs9e4YdU/vum56fTTKX1uCIiEiZbrstjL6MGQO//nXlnm3UCM48M7yWLIHRo8N9Uo88\n8mNMkyZw8skwcmSYnhJJBhU4IiKS0LRpcP314TVwYPXaatnyx7a+/x5WrAhrbFq0CGfdiCSTChwR\nESnVunVh2uiYY+CGG5Lb9o47QuvWyW1TJJ4KHBERKdXQofD112EUJxln24ikk/7KiojIT8yfD/fc\nA7feCvvvH3U2IpWnWU8REfmJvDw4+GD43e+izkSkajSCIyIi23n99bCF+9lnwy4okdpIIzgiIvI/\n7mGXU04O9OkTdTYiVacRHBER+Z+CAnjrrXAAn1nU2YhUnUZwRETkf26/PVyXcPLJUWciUj0awRER\nESDsnJo6NdwJpdEbqe00giMiIgCMGBHugOrfP+pMRKpPBY6IiPDVV+Hm78svD9cniNR2KnBERISH\nHgr3QV1ySdSZiCSHChwRkTpu2zZ4+GH41a9gt92izkYkObTIWESkjnv9dVi0CB57LOpMRJJHIzgi\nInXcqFHQpk24NVwkU6jAERGpw9asgQkTYNAgbQ2XzKICR0SkDnvySdi6FQYOjDoTkeRSgSMiUoc9\n8UQ4tbh586gzEUkuLTIWEamjPv0UZs2C/PyoMxFJPo3giIjUUU8/DTvuCKedFnUmIsmX8gLHzIaY\n2SIz22BmM82sQznx3cys0Mw2mtlHZnZ+KTEDzOzDWJvzzOwn18JVpF8zu9nMlprZejObYmYHlfj+\n62a2Le611czuq8rnICJS0+TnQ58+ocgRyTQpLXDM7CzgTmAocCQwD5hsZnskiG8FvAhMBQ4HRgKj\nzKxHXMzRwFPAQ8ARwHPARDM7pDL9mtl1wOXApUBH4PtYTKO4lBz4J9AcyAL2Bq6t0ochIlKDzJ8P\n778PublRZyKSGqkewckDHnT3x9x9AXAZsB64KEH8b4HP3P1ad1/o7vcC42LtFLsSeNnd74rF3AAU\nEYqVyvR7FXCLu7/o7u8BA4F9gL4lclrv7ivd/avY67sqfA4iIjVKfj7suiv07Bl1JiKpkbICx8wa\nAjmE0RgA3N2BAqBLgsc6x74fb3KJ+C5lxVSkXzM7gDAiEx+zDphVSm7nmtlKM5tvZrea2Q4JchcR\nqRXcQ4HTvz80alR+vEhtlMpdVHsA9YEVJd5fAbRJ8ExWgvimZtbY3X8oIyarEv1mEaafymoH4Eng\nv8BSoD0wDPg50D9B/iIiNd4778DixZqeksymbeJlcPdRcV++b2bLgKlmdoC7Lyrr2by8PJo1a7bd\ne7m5ueTqvygiErEJE2DPPaFr16gzkbosPz+f/BJnFKxduzZp7aeywFkFbCUs0I3XHFie4JnlCeLX\nxUZvyoopbrMi/S4HLPbeihIxcxLkBjA79txBQJkFzvDhw8nOzi4rREQk7dxh/Piwe6p+/aizkbqs\ntF/6i4qKyMnJSUr7KVuD4+6bgUKge/F7Zmaxr6cneGxGfHxMz9j7ZcX0KI4pp9/imEWEIic+pinQ\nqYzcIOzIcmBZGTEiIjXW++/DJ5/AGWdEnYlIaqV6iuouYLSZFRJGP/KAJsBoADO7DdjH3YvPunkA\nGGJmtwOPEAqQ/sApcW2OBF43s6uBSUAuYVHxJRXo99G4mBHAn83sE2AxcAvwBWHbOWbWGjgHeAlY\nTdi2fhfwRmzXlYhIrfPss7DzztC95K+JIhkmpQWOu4+NnT1zM2H6Zy7Qy91XxkKygJZx8YvNrDcw\nnLAd/AtgkLsXxMXMMLNzgL/GXh8Dfdz9g0r0i7sPM7MmwIPALsA04GR33xQL2QScSNhOviOwBHgm\n1qeISK00YQL07g2NG0ediUhqWdhBLcliZtlAYWFhodbgiEiNsmgRtG4NY8fCgAFRZyPyU3FrcHLc\nvag6bekuKhGROuLZZ8PIzck/udxGJPOowBERqSMmTAgnF++0U9SZiKSeChwRkTpg+XKYPh3OPDPq\nTETSQwWOiEgd8NxzUK8enHZa1JmIpIcKHBGROmDCBDj+eNh996gzEUkPFTgiIhnum2/g1Vc1PSV1\niwocEZEMN2kSbNkSrmcQqStU4IiIZLjnnoOjjoJ99406E5H0UYEjIpLBNm6El1+Gvn2jzkQkvVTg\niIhksFdfhe++U4EjdY8KHBGRDPbcc3DggXDIIVFnIpJeKnBERDLUtm2hwOnbF8yizkYkvVTgiIhk\nqFmzYMUKTU9J3aQCR0QkQz33HOy5J3TpEnUmIumnAkdEJENNnBiuZqhfP+pMRNJPBY6ISAZasAAW\nLtThflJ3qcAREclAzz0HTZpAjx5RZyISDRU4IiIZaOJE6NULdtgh6kxEoqECR0QkwyxbBjNnanpK\n6jYVOCIiGeaFF6BePTj11KgzEYmOChwRkQwzcSJ07Qq77x51JiLRUYEjIpJB1q6FqVM1PSWiAkdE\nJIO88AJs2gT9+kWdiUi0VOCIiGSQZ56Bzp2hZcuoMxGJlgocEZEMsW4dTJ4MAwZEnYlI9FTgiIhk\niBdfhB9+gP79o85EJHoqcEREMsS4cdCxI+y3X9SZiERPBY6ISAb47jt4+WVNT4kUS3mBY2ZDzGyR\nmW0ws5lm1qGc+G5mVmhmG83sIzM7v5SYAWb2YazNeWZ2clX6NbObzWypma03sylmdlCJ7zc2s3vN\nbJWZfWtm48xsr6p8DiIiqTRpEmzcqN1TIsVSWuCY2VnAncBQ4EhgHjDZzPZIEN8KeBGYChwOjARG\nmVmPuJijgaeAh4AjgOeAiWZ2SGX6NbPrgMuBS4GOwPexmEZxKY0AegP9gK7APsD4Kn0YIiIpNHYs\n5OTAAQdEnYlIzZDqEZw84EF3f8zdFwCXAeuBixLE/xb4zN2vdfeF7n4vMC7WTrErgZfd/a5YzA1A\nEaFYqUy/VwG3uPuL7v4eMJBQwPQFMLOmsfg8d3/D3ecAFwLHmFnHqn8kIiLJ9c03YYHxOedEnYlI\nzZGyAsfMGgI5hNEYANzdgQKgS4LHOse+H29yifguZcVUpF8zOwDIKhGzDpgV19dRQIMSMQuBz8vI\nX0Qk7caNgy1b4Oyzo85EpOZI5QjOHkB9YEWJ91cQiovSZCWIb2pmjcuJKW6zIv1mAV5OTHNgU6zw\nqWj+tdLGjbB6ddheKiK1zxNPQPfusM8+UWciUnM0iDqBTNW3bx6HH94Msx/fy83NJTc3N7qkAHeY\nPx/Gj4e334bCwjC8Xax58zCPf+qp4SyNPfeMLlcRKd/nn8Mbb8Bjj0WdiUjl5Ofnk5+fv917a9eu\nTVr7qSxwVgFbCSMh8ZoDyxM8szxB/Dp3/6GcmOI2K9LvcsBi760oETMnLqaRmTUtMYpTVv7/s2TJ\ncDp1yuaxx2CHHcqLTr2NG8NveSNHwnvvwa67htuGr7kmHOm+007w7bfwyScwYwZccQXk5cEFF8D1\n10OLFlH/BCJSmiefhCZN4Iwzos5EpHJK+6W/qKiInJycpLSfsikqd98MFALdi98zM4t9PT3BYzPi\n42N6xt4vK6ZHcUw5/RbHLCIUKfExTYFOcbkVAltKxLQB9iuRT6n+/nd46SXo1g1WlJwIS6PNm+Hu\nu8POiksvhYMOCosRly+HiRPhT3+CgQPhzDPh/PPhllugoCB8/4YbYMIEaNsWhg+HrVuj+zlE5Kfc\nwy8uffuGX1JEJI67p+wF/Iqwe2kg0BZ4EFgN7Bn7/m3AmLj4VsC3wO1AG2AwsAk4MS6mC/ADcHUs\n5kZgI3BIRfuNxVwbe+804DBgIvAx0Cgu5j5gEdCNsHD5bWBaOT9zNuCFhYX+7rvue+/tvv/+7vPn\ne1pt2+Y+aZJ7mzbuZu4XXui+cGHl21mzxv3yy0MbPXq4r1iR/FxFpGreftsd3CdPjjoTkeQoLCx0\nwhrZbK9mDZLSbeLuPha4BriZMPXTHujl7itjIVlAy7j4xYRzZ04E5hK2ew9y94K4mBnAOYTza+YC\nZwJ93P2DSvSLuw8D7iEUP7OAHYCT3X1T3I+QRziXZxzwOrCUcCZOheTkwOzZsMsucPTR8MorFX2y\net5/H04+GXr3DosO58yBRx6Bn/+88m3tsgvccw9MmQLz5kF2dljDIyLRe/DBMDp74olRZyJS85iH\nUQdJEjPLBgoLCwvJzs4GwhHqubnhpNE//hFuvBEapGD106pVoe0HHoBWreDOO+H009luoXN1LF0a\niqZFi+CFF+C445LTrohU3po14ReYG26AP/wh6mxEkiNuDU6OuxdVpy3dRZUGO+0U1rv89a/wt7+F\nxb2ffZa89jdtghEj4OCD4fHHQx/vvw99+iSvuIHwH9M33ggjUz17wtSp5T8jIqnx+OPh7JsLL4w6\nE5GaSQVOmtSvH37LmjYNli2DX/wiFDzVOXtm61Z4+mk47DD4/e/DIV8ffxx2RjVuXP7zVdG06Y+L\np08/PWw1F5H0cg/TU336QFZGncolkjwqcNKsS5ewhuXKK8N0Utu2YUpp48aKt/Htt/DQQ6FIys2F\nAw+EuXPh/vthrzRcBdq4cThHp2NHOOUUKKrWIKKIVNbkyfDBB3D55eXHitRVKnAisNNOcPvtYdFu\np04weHCY/vnNb8LallWrto/ftg0WLw7bQc87D/beGy67LBRH77wTRlQOOyy9P0OTJvD889CmTTgU\ncMmS9PYvUpf9/e9w1FFw/PFRZyJSc6nAidAhh4QppgULQsEyeXKY9tlzzzAV1Lo17Lsv7Lhj2Cnx\n61+HtTXXXhsKnokTw3/korLzzqEga9QoLD5eV/JSCxFJurlzw/q3a65J7ho7kUyjXVRJVtouqopy\nD4XL7NlhRGT16lA87LJLGK3Jzg5XKdQ0H3wQtsF36RIKnlTsEBORIDc3nDb+ySf635pknmTuotL/\nPGoQszBSc8ABUWdSOYccEtbknHRSuN7hnnuizkgk+dxh5cpwynfjxrD//vCzn6U3h/nz4V//Cuvt\nVNyIlE1TVJIU3bvDP/4RXvfdF3U2IsmzdClcd11YzN+8ORx+eBhR3XnnMHL5j3+kb3r2+uvDL0AX\nXZSe/kRqMxU4kjS/+U3YHXblleE+K5HabOvWsJj3wAPDluyTToJx42DmTHj99XBx7R57hFHLVq3g\n1lurd+xDeWbMgOeeC7svGzZMXT8imUJrcJKsOmtwMsGWLWFX1axZ4VWV6yFEovbtt/CrX4WF/7/7\nHQwdCs2alR775ZfhcM0HHwwbA/75z3CYZzJt2QIdOoTztGbNCn+KZCKdZCw1VoMGYY3A3nuHQmfN\nmqgzEqmcr78O26+nTw8Fzl13JS5uAFq0COvO5syB3XcPz/6//webNycvp/vvD8dK3H+/ihuRilKB\nI0nXrFnYTbV6NQwYkNz/0Iuk0vffh8L888/DqeM9elT82V/8Ijxz++3hHrjjjgv3tlXXggWhYLrs\nsjCKIyIVowJHUuLAA8POqjfegKuuijobkfK5h4M058+Hl1+G9u0r30a9euGcqrfeghUr4IgjID+/\n6jlt2BC2he+3H9xxR9XbEamLVOBIynTrFobU778f7r036mxEynbHHeHwzCefrP5ISadO4UC+U06B\nc84Jh3RWdqfV1q2h4Fq4MBRJO+5YvZxE6hqdpCApdfHF4SDAq64Kt5337Bl1RttbsyacDr1yZVjb\nsOee4Td3/WNSt7z1VrgM9w9/CKeJJ0OzZvDUU+GU78GDw8W0o0dXbAHy5s0waFAouCZODCNBIlI5\nGsGRlLvjDujVK+xKWbAg6mzCCbA33QSHHgq77RbWSpx5ZriZ+eijwzUZ3brBo4/Cpk1RZyuptn49\nXHghdO4MN9+c3LbNwijM3Llh4f3xx4cpp88+S/zMokXhXKmnnw6jSaedltycROoKFTiScvXrhyH2\nFi3CAs7Vq6PJo6gI+vULI0l33hnu8XrssbDm4quvwpqJoqJwUGGjRuE36DZt4JlnoslX0uPPf4Yv\nvggFbapOB27dOixAfvTRcIbOQQeFkaJRo8LVLPPnh8trL744HCL43/+Gs6TOPjs1+YjUCe6uVxJf\nQDbghYWFLtv79FP33Xd379bNfdOm9PX75Zfu55zjDu4HHuj+0EPu339f/nPz57uffnp4bsAA92++\nSX2ukl4zZ7qbud9xR/r6/P5793/+071TJ/d69cLfr+JXq1but93mvm5d+vIRqUkKCwsdcCDbq/nv\nsQ76S7K6ftBfed58E048MUxXjR6d2vt0Nm0Kp83efDPssAPcdhucf37l+nSHsWPDKc377BO2vx94\nYOpylvTZti1MS23eDO++G835MmvXhimpDRvCCGfLlrohXOo2XbYptVbXrvDEE2FnyaZNYY1BKo6d\nf/PNcG7IwoVw+eVhzc0uu1S+HTM466xw/9Bpp4U1Oq++Gs48kdrtscfgnXfC35WoDs9r1kwLiEVS\nRWtwJO1+9atwp8/EidC3b3IvKvz667CO4fjjQ0EzZ04YxalKcROvbdtwF1BWFvzyl/Dee8nJV6Kx\nbl04PO/ss8MicxHJPCpwJBJ9+8KLL4btuV26hJ1N1bFtW/iNvG3bUDw98EBouyqHtSWyxx5h9Gaf\nfcKusCVLkte2pNett4YiZ9iwqDMRkVRRgSOR6dkzXBy4ZQsceWQ4EHDbtsq389pr0LFjWF9zwglh\nK/pvfhNOlU223XeHV14J02qnnBLWUEjtsmxZGNW75pqw5kVEMpMKHIlU27ZhHcS554bD0Dp1Csfk\nl7f2fcuWMMXVtWsoaho0CNtwn346TCOlUlYWvPRSGME5++xw4qzUHrffDj/7GVx9ddSZiEgqqcCR\nyDVtGqaUXn/9x5GR1q3Db9jjx4cC6P33w0mwDz0EF10UiowzzgiF0LPPhvUxxx6bvpwPOSTsrpo8\nGf7yl/T1K9XzxRfh79rvf1/9dVkiUrNpF5XUGMcfH4qYt94K62nGjg0H8sWrVy8UFxdfHBYrR7kT\nv2fPsDtr6NCw3bhXr+hykYq57bZwDceVV0adiYikmgocqVHMwq6W4p0ty5eH1/r1YUttq1Y1656o\nP/0pjB6dcw7Mmwf77ht1RpLI55+HEcCbbw6jhiKS2VI2RWVmu5rZk2a21szWmNkoMyv3nyYzu9nM\nlprZejObYmYHlfh+YzO718xWmdm3ZjbOzPaqbN9m1tLMJpnZ92a23MyGmVm9uO/vb2bbSry2mlnH\n6n42UnFZWeGckKOPDmfP1KTiBsKI0uOPhzUdF1xQtUXSkh5//WuYlrr88qgzEZF0SOUanKeAdkB3\noDfQFXiwrAfM7DrgcuBSoCPwPTDZzBrFhY2Itdcv1uY+wPjK9B0rZF4ijGB1Bs4HLgBKXrXnwAlA\nVuy1N1BYzs8tdczuu8OYMTB1atidIzXPZ5/BI4/AddfBTjtFnY2IpENKChwzawv0Aga5+7vuPh24\nAjjbzMra43IVcIu7v+ju7wEDCQVM31i7TYGLgDx3f8Pd5wAXAscUj6yYWbsK9N0LaAuc6+7z3X0y\ncD0wxMzip+0M+Nrdv4p7ac+M/MSJJ0JeXjg8bv78qLORkv7yl1CI/va3UWciIumSqhGcLsCaWAFS\nrIAwItKptAfM7ADCKMnU4vfcfR0wK9YewFGEUZf4mIXA53ExnSvQd2dgvruviouZDDQDSh7C/7yZ\nrTCzaWZ2Wlk/tNRtt94KP/952PK+cWPU2VTexx/D3/4WdqdlZ4frKXr3hhtvDOuLaquPPw6L1v/w\nB2jSJOpsRCRdUlXgZAFfxb8RG/n4Ova9RM84sKLE+yvinmkObIoVPoliKtJ3VoJ+iIv5DrgaGACc\nArwFTDSzUxPkL3Xcz34W7tZauDAsPq4t3nknbM3/+c/DSMd334WDE4u33Y8YEdZBnXBCWFBd29xy\nCzRvHg5/FJG6o1IFjpndVsrC25KLcH+eqmTTyd1Xu/sId3/H3Qvd/Q/A48D/RZ2b1Fzt24eRnLvu\nCtc61GTffRe2S3fqBIsXh8XSX30FU6aEs2LuvRcmTQrvPfMMrFkDxxwTntmwIersK2bBglB0/ulP\noQAVkbqjstvE/w48Wk7MZ8ByoOTOpvrAbrHvlWY5Yc1Lc7YfXWkOzImLaWRmTUuM4jSPa7cifS8H\nOpTov3nc9xKZDfQo4/v/k5eXR7NmzbZ7Lzc3l9zc3Io8LrVYXl64Z+uCC+A//6mZB8otXAh9+oSt\n03feGYqWRDdqN2oE/fuHqat77glTPTNmhJOkW7RIb96VddNNIcdBg6LORERKys/PJz8/f7v31ibz\n/ht3T/qLsIB3K3Bk3Hs9gS1AVhnPLSUsIC7+uimwARgQ9/UPwBlxMW2AbUDHivYNnARsBvaIi7kU\nWAM0LCO/h4B3y/nZswEvLCx0qbsWL3Zv2tT93HOjzuSnXnop5NaunfvChZV/vrDQfd993ffZx33B\nguTnlyzz57ubuT/4YNSZiEhFFRYWOmG5SrZXsxZJyRocd19AWLT7kJl1MLNjgHuAfHf/3wiJmS0w\nsz5xj44A/mxmp5nZYcBjwBfAc7F21wEPA3eZWTczywEeAd5299mV6PvfwAfA42bW3sx6AbcA/3D3\nzbHcBprZ2WbWJvb6I2Er+d1J/8Ak4+y/f5jiefJJ+Ne/os7mR88/H0Zujj8+XHT68ypMKGdnh3U7\nu+wS2nn//eTnmQw33hj+/3DhhVFnIiJRSOU5OOcACwg7mF4E3gRKLvM7mLBzCQB3H0YoRh4k7J7a\nATjZ3TfFPZMXa28c8Dph1KdfZfp2923AqYSRnumEQmo0MLREO9cD7wIzgdOAX7n7YxX66aXOO/dc\nGDAgbE3+8suoswnTZv37hwJn/HjYeeeqt5WVFe4Oa948bJFfvDhZWSbHnDnhZxw6NNxvJiJ1j3l5\n1zZLpZhZNlBYWFhIdpQXJUmN8PXXcOih4fXKK+Hk4yi89FJYQ3PqqeHG9WT9o79iRThlumHDcI/Y\n7rsnp93qOu00+OijMLrUQBfSiNQaRUVF5OTkAOS4e1F12tJt4iIptNtuMHp02Jl0773R5PDKK6G4\nOeUUyM9P7ohG8+ah/dWr4fTT4Ycfktd2Vc2cGUarbrxRxY1IXaYCRyTFevaEK66Aa6+FDz5Ib9+T\nJ0PfvnDSSWEtUKNG5T9TWQcfHAqKwkIYMgSiHhS+/vowYnbWWdHmISLRUoEjkgZ/+1u4Cf3cc9N3\nhsy//x3W2/ToEc6xSUVxU6xTp3B2zsMPhz+j8sYbUFAQtodHNR0oIjWD/hMgkgZNmoTpoQULwqLj\nVI9yFBSE4qZ7dxg3LrXFTbELLgg3dV95JUyblvr+SnIPZ/RkZ4cpORGp21TgiKTJEUfAP/8Zbh6/\n777U9TN1alhk+8tfhp1EjRunrq+S7rorLDru3z/9O8eeeSYcQDhsGJilt28RqXlU4Iik0a9/HUY4\nfvc7eO215Lf/73+HnVLHHw8TJqT/eoKGDX+cDuvXL32LjjduhOuuC4Vd9+7p6VNEajYVOCJp9ve/\nh9GVPn2gqFqbILf30ks//gM/cWJ0dy/ttVcorubODYur0+Huu2HJkjB6IyICKnBE0q5hw1AAtGsX\ndjctWFD9Np98MuyWOvnkaEZuSurQAe6/Hx56KEzLpdLSpeEW9Msug7ZtU9uXiNQeKnBEIrDTTmHE\npXlzOO64cPVBVbjDLbfAeefBOefA2LHpWVBcERdeGBZUX355OJsmVa68EnbYAW6+OXV9iEjtowJH\nJCK77x62NR98cJiyKnGpbrlWrgxTUjfcEIqcRx+tOcVNsREjwmhOv36wfHn58ZX1wgthIfXIkeFQ\nRRGRYipwRCK0227hlOO+fcMIzMCBsGxZ2c9s2waPPBIOs5s1Kxyy9+c/18ydQ40ahW3q7uFerk2b\nyn+motasgcGDw7ScDvUTkZJU4IhEbMcd4fHHwwjMpElhRGfwYHjrrbA7CGDrVvjww7BAuU0bGDQo\nXHL5n/9A797R5l+evfcORc6sWXDVVck5A8gdLrkEvvsuHCxYE4s7EYmWbmoRqQHMwkF5ffqEs2Qe\nfjgs0jUL63XWrw9FTvH266eeClM/tcXRR4ezfy65JNxEPnRo9dr75z/D1NS4cbDffsnJUUQyiwoc\nkRpk113DepqhQ8PozNy5YSpmp52gdWs45phwKnJtdPHF8NVX8Kc/ham5qm4hf+21sHB58OBQ7ImI\nlEYFjkgN1KBBuHIgOzvqTJLrD3+Ar78OO582boRrrqnc9NK8eXDmmWFR9ogRqctTRGo/FTgikjZm\ncMcdYVv3tdfCF1+EdUUNG5b/7KxZYUHxgQeG05Ir8oyI1F1aZCwiaWUWpuHuvTesyznmGHj//cTx\nW7fCP/4BXbuGwxGnToVmzdKXr4jUTipwRCQSgwfD9OnwzTfQvn04rHDKFFi3Lnx/2bJwMWmHDmG9\nzsUXw6uvqrgRkYrRFJWIRKZDB3jvvbAravjwcOUEQL164bwfCHdrTZ8OXbpEl6eI1D4qcEQkUo0a\nhV1RQ4bA/Plh99j338Oee0KnTtCiRdQZikhtpAJHRGoEszBV1b591JmISCbQGhwRERHJOCpwRERE\nJOOowBEREZGMowJHREREMo4KHBEREck4KnBEREQk46jAkYyQn58fdQp1jj7z9NNnnn76zGuvlBU4\nZrarmT1pZmvNbI2ZjTKzHSvw3M1mttTM1pvZFDM7qMT3G5vZvWa2ysy+NbNxZrZXZfs2s5Fm9q6Z\nbfdxjLcAAAdsSURBVDSzogS5tDezN81sg5n918z+ryqfhaSe/iOUfvrM00+fefrpM6+9UjmC8xTQ\nDugO9Aa6Ag+W9YCZXQdcDlwKdAS+ByabWaO4sBGx9vrF2twHGF+Fvh14GHg6QS47A5OBRUA28H/A\njWZ2cVk/g4iIiEQvJScZm1lboBeQ4+5zYu9dAUwys2vcfXmCR68CbnH3F2PPDARWAH2BsWbWFLgI\nONvd34jFXAh8aGYd3X22mbWrSN/u/rvY9/YCSjs79TygITDI3bfE+jgSuBoYVa0PSERERFIqVSM4\nXYA1xQVGTAFh1KRTaQ+Y2QFAFjC1+D13XwfMirUHcBShKIuPWQh8HhfTubJ9J9AZeDNW3BSbDPz/\n9u4txK7qjuP492dQY9QQ2tQxqMV7QCqjjlKDl1qVeAEjQrEX86A+iChUfFFLC7EqCm0Nlaov3oNE\nEUHUBy+IlpomaUjiBTGJ2ERj1XiZwKSESmL892GtY3YOc+KcOWufI3t+H9iQ2Xvtvdf82Zn5zzp7\nrf9cSa5nbGZm9j1WVy2qQ4HPqzsiYpekrflYp3OCNGJT9VnlnCFgR058OrWZzL079WfjOPdpHRvr\ncN50gHXr1nVxK+vV2NgYa9eO+yqV1cQx7z/HvP8c8/6q/O6c3uu1ukpwJN0F3LyXJkF692UqOxJg\n4cKFA+7G1DMyMjLoLkw5jnn/Oeb955gPxJHA8l4u0O0Izl+AR76jzUZgC9A+s2ka8IN8bDxbAJFG\naaqjOEPAG5U2+0ma2TaKM1S57mTu3ak/Q237hirHOnkJuAL4APiqi/uZmZlNddNJyc1LvV6oqwQn\nIkaB0e9qJ2kFMEvSyZV3Yc4jJTD/6nDtTZK25HZv5+vMJL03c19utgb4Ord5JreZC/wYWJHbdH3v\nDlYAd0iaFhG78r75wIaI6PTxVCtGS7u4j5mZme3W08hNSy0vGUfEelL29YCk0ySdAfwNeKI6g0rS\nekmXVk79K/AHSZdIOhFYAvwHeDZfdxtpavdiSedIGgEeBv4ZEau6vPcxkk4C5gAHSBrO2765yVJg\nB/CwpBMk/RL4LXB32WiZmZlZaXW9ZAzwG+Be0gymb4CnSdPAq44Dvp2RFBF/kjSDtGbNLOB14KKI\n2FE550ZgV77e/sCLwPWTuPeDpPVxWlpvkR0FbI6IbZLmk0aPVgNfArdGxEMT+ebNzMxscBQRg+6D\nmZmZWVGuRWVmZmaN4wTHzMzMGscJTkGSrpe0KRfnXCnptEH3qSkknSXpOUkfS/pG0oJx2uy1UKt1\nR9LvJK2StE3SZ5KekXT8OO0c90IkXSvprVwoeEzSckkXtrVxvGsk6Zb8M2Zx237HvRBJi3KMq9u7\nbW16jrcTnELyLKu7gUXAycBbpEKhswfaseY4EHgTuI60oOQeJlio1bpzFmkG4k+B80m12V6WdECr\ngeNe3EekxVRPAUaAV4Fnc409x7tm+Y/Sa0g/v6v7Hffy3iGtLXdo3s5sHSgW74jwVmADVgL3VL4W\naYr7TYPuW9M20sy4BW37PgFurHw9E/gfcPmg+9uUDZidY3+m497XuI8CVznetcf5IGADcC7wGrC4\ncsxxLxvrRcDavRwvEm+P4BSQ184ZYc8ioEGapj6v03lWxgQLtVrvZpFGz7aC4143SftI+hUwA1ju\neNfuPuD5iHi1utNxr81x+ZWDf0t6XNIRUDbeda6DM5XMBqYxfqHQuf3vzpQzkUKt1gNJIi3EuSwi\nWp+VO+41kPQT0krq04H/ApdFxAZJ83C8a5ETyZOAU8c57Oe8vJXAlaQRsznArcA/8rNfLN5OcMxs\nIu4HTgDOGHRHpoD1wDBpEdRfAEsknb33U2yyJB1OSt7Pj4idg+7PVBAR1TpT70haBXwIXE56/ovw\nR1RlfElaXXm84pzdFPi0yakWaq1y/AuQdC9wMXBORHxaOeS41yAivo6IjRHxRkT8nvTC6w043nUZ\nAX4ErJW0U9JO4GfADZJ2kEYOHPcaRarv+B5wLAWfcyc4BeSsfw2pqCfw7ZD+eRQqGmadRcQm0oNf\njX+rUKvj34Oc3FwK/DwiNlePOe59sw+wv+Ndm1eAE0kfUQ3nbTXwODAcERtx3Gsl6SBScvNJyefc\nH1GVsxh4VNIaYBWpZtYM4NFBdqopJB1I+g+gvOtoScPA1oj4iN2FWt8HPgBup1Ko1bon6X7g18AC\nYLuk1l9UYxHxVf63416QpDuBF4DNwMHAFaTRhPm5ieNdWERsB9rXYNkOjEbEurzLcS9I0p+B50kf\nSx0G/BHYCTyZmxSJtxOcQiLiqbzmzW2kobQ3gQsi4ovB9qwxTiVN3Yy8taq6PwZcHRMr1GrduZYU\n67+37b8KWAITLpBrE3cI6ZmeA4wBbwPzWzN7HO++2WOtLce9uMOBpcAPgS+AZcDpETEK5eLtYptm\nZmbWOH4Hx8zMzBrHCY6ZmZk1jhMcMzMzaxwnOGZmZtY4TnDMzMyscZzgmJmZWeM4wTEzM7PGcYJj\nZmZmjeMEx8zMzBrHCY6ZmZk1jhMcMzMza5z/AwFvxb2bH4KSAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f39e06734e0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"alpha = np.empty(3)\n",
"alpha[0] = 0.029\n",
"alpha[1] = 0.057\n",
"alpha[2] = -0.01\n",
"\n",
"a = np.empty(3)\n",
"a[0] = 1e-5\n",
"a[1] = 2e-5\n",
"a[2] = 1e-5\n",
"\n",
"omega = np.empty(3)\n",
"omega[0] = 0.9\n",
"omega[1] = 0.125\n",
"omega[2] = 0.01\n",
"\n",
"phi = np.empty(3)\n",
"phi[0] = np.pi/2\n",
"phi[1] = np.pi/3\n",
"phi[2] = np.pi\n",
"\n",
"y = 0 * t\n",
"\n",
"for i in range(3):\n",
" y = y + a[i] * np.exp(alpha[i]*t) * np.cos(omega[i]*t + phi[i])\n",
" \n",
"plt.figure()\n",
"plt.plot(t, y)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python [conda root]",
"language": "python",
"name": "conda-root-py"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment