{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Lab 1: Identification of Damping Using Log Decrement\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "clear all\n", "close all\n", "clc\n", "imatlab_export_fig('print-png')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Beam Properties\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "wn =\n", "\n", " 108.2591\n", "\n", "\n", "wn2 =\n", "\n", " 219.7387\n", "\n" ] } ], "source": [ "l=21.75*0.0254;% length in meters\n", "h=0.5*0.0254;% height in meters\n", "w=1*0.0254;% width in meters\n", "rho=2700;% density in kg/cubicmeter\n", "E=7.31e10;% youngs modulus in Pa\n", "I = (1/12)*w*h^3; % moment of inertia (m^4)\n", "k = (3*E*I)/l^3; % stiffness (N/m)\n", "V = l*w*h;% volume (m^3)\n", "m = rho*V;% mass (kg)\n", "wn = sqrt(k/m) % analytical natural frequency of massless beam with concentrated mass\n", "wn2=(1.875)^(2)*sqrt((E*I)/(m*(l)^3)) % natural freqency of a uniform section beam\n", "%c_cr = 2*sqrt(k*m); % critical damping coefficient" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Acceleration vs Time plot\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA\nB3RJTUUH4QoZBDcBlmRaTQAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJ\nbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAyNS1PY3QtMjAxNyAwMDo1NDo1OX82MkoAACAA\nSURBVHic7Z19WFRl/v9vB08IConyoKQBZmBqGiJXiqJJsrG5PSB12aJb+y3zu5B++ZkPtZnakm2W\notmmKUZX9DVrLXXtW5epQT6Eu60PGAI2kk67AqaDuZIyjMMwvz/OdvY0j/fM3HPOfWber8vL6557\n7vP5vM89w/nM/dzDZrMRAAAAQG10agsAAAAACEFAAgAAwAkISAAAALgAAQkAAAAXICABAADgAgQk\nAAAAXICABAAAgAsQkAAAAHABAhIAAAAuQEACAADABQhIAAAAuKCn2gIAX3R3d3d0dNxwww033HCD\nWgI6OzsJIZGRkaoI8Iqamprm5ubw8PAHH3wwFPz6j3aVAyWwAWCz2Wy29evX33777dIXIzExccGC\nBUoKWLt27e233x4WFiYK6N2794IFC0wmk5IavCU/P58QEhsbGzgX+/fvf+qpp5566il5VSjg13+0\nqxyoBVpIgHR3d99///2ffvqp+LJ3795dXV2tra1lZWWrV69WTMbnn39+8uTJ+Pj41NTUr7766tq1\na2VlZVeuXNm8ebNiGjikoaFh/fr1hJCVK1dKmXfeeSchJCoqSjVZFGhXOVALBCRAXnzxRTEajR49\neuvWrcOHDyeEnD59+vXXXxcLdHd3//WvfzUajSaTqWfPnsnJyZmZmXILp0+f/vrrr7u6uiIiIm66\n6aaMjAyd7j/Dk998883Jkye7urpuvvnmCRMmuJIxYcKEpUuXig+sb775ZtSoURaLpbKyUgxIRqPx\n66+/FkXGxcU5teBRJyHk3Llzx44dM5vN4eHhOTk50dHRNG9R3oL7wnV1dRcvXgwPD8/Ozj527Njp\n06cnTZo0cOBAN5pPnDhx+vRpMb1///5evXrFxMRkZGT88pe/zMjICA8Pl/ttbGxsaGjo6uoaMGDA\n5MmTpY9A7vfEiRN6vT4mJuYXv/iFU/EnTpxoa2sTC0u1Wl1dTQhJTk4eOnQo8fRxe6Vc7k7szRs9\nevSwYcPEt06dOjVgwIApU6ZQ1jDQPGo30YD6xMbGEkLCwsKam5udFnj00Uftvjbp6enff/+9+G5J\nSYndu4WFheJbP/zww7Rp0+RvjR49+p///CeNqqlTp4qqrFarzWbbuXOnaGHnzp2uLnGv88qVKwUF\nBfJ34+PjPb7l/hbsOqAoC8+bN0+6F/ea7awRQqZOnerot62t7Z577pEXu/nmmw8ePGjnt6ioSCpw\nzz33OK3DFStWiNUuaZBqft++fe4/bjmUysVisbGxhYWFUsm33357/vz50ssHHnhAMuvPNwrwDwJS\nqNPQ0CD+Yefk5Lgqs2jRovnz52/dunXbtm3Lli0TBIEQMmfOHJvNdubMGfHyBQsWbNu27d133330\n0UeLiorEC8VnR+/evVetWvXuu+/efPPNhJARI0Z4VNXa2ip6EZ9iNrqA5EanzWa77777RAsTJ05c\nt27dvHnzJCVu3nJ/C04fr+4Li9x11105OTm7du1yr3nlypUZGRniJQUFBY888siKFSsc/UrRaM6c\nOaWlpfHx8YSQmJgYMahIfuPj43/729/eeuut8gDjWPPiMN7atWvlspOTkz1+3HIolUvRJTk5+be/\n/W1MTAwhJDw8PCIiYtasWSNGjBDfPXz4ME0NA62DgBTq7N69W/ybf+SRR9wUs1qtJ0+e3L1796ef\nfjpu3DhCSGJios1m+9vf/iZe/uabb7a2tsovOXXqlPjWH//4RzFn+/btYk51dbUbX99//704vSI+\nPv67774TM48ePfroo48++uijR48e9UGnJEbeMjCbzZRvuboF+eOVsjAhZNu2bTSaRcRhGELIjz/+\nKGU69fvoo4+K737wwQdizsqVK+V+GxoabDbb4cOHxZd/+tOfnNahGJ4zMjJsNtuVK1fEGCmGEzcf\ntyMeldtkAenMmTM2m01qGIna9uzZI77csmULTQ0DrYOAFOqIAzOEkGnTprkq8+abbzqOQkdFRdls\nNqvVmpaWJmWGh4ffc889YmeR9LBw5N1333Xlq6GhITExkRAyaNCgU6dOeXUvbnRKYsRHmxyat1zd\ngvzxSlm4d+/elJpFPD7WJb9bt24V3zWZTGLOrFmzpMKSzR9//FF896mnnnJajdL0llOnTr311ltE\n1oPn5uN2hD4gxcTEiC9Fd+SnJtHBgwflFejbNwpoCExqCHVGjhwZFRX1448/7t2798KFCwkJCXYF\nGhsbxbGHu+66q7i4WBCEV155RfqlrNPp/vrXv1ZUVBw+fPibb745derUnj17vvzyyx9++EEa6L7r\nrrvErhUJu5cShw4deuCBBy5fvnz77bfv2bNn4MCB9DfiUaeYuHbtmt2FNG/R3AJl4YiICErNlEh+\nu7q67G6kZ8///IHbzYBwQ15e3qBBg5qbm//3f/9XbE7de++94hfDzcftz8I1aa6/RGpqqmMxH75R\nQGOoHRGB+jz33HPil2HSpEnSvIa2trZly5bZZD9Ld+3aZbPZTCZTcnIykbWQ5L9/pc6i6urq5uZm\nMT1v3jy5u7///e/iPAU7du7cKT40c3Jy2trazD8hvvvll19OnTp16tSpX375pdO7cK9TGhoZPXq0\ntCxGHAx385bHW5D/3veqMI1mkTfffFMso9frpUy5qdbWVrGANAq4bt06MWfTpk2Ofj22kGw22/PP\nP08IEZuqhJDdu3eL+W4+bkcjHpXbZJMaxJdSC6mtrc3m0ELy9hsFNAcCErCZTKZJkyaJf+phYWHJ\nycmxsbHiM9pms508eVJ8Ky0traioKD09XRxUEB+aZrNZEIT77rtvzpw5RUVF4thPeHj4lStXbDab\nNK3rnnvuKSoqevTRR8UCUpiRI06rc0R8Anqc1OBep002Pezmm29++OGHf/nLX0pT6dy85f4W7B6v\nXhWm0WyTdaCFh4fHxsbOnz/fjd9x48bl5+eLn92QIUPE+OpDQPruu++k+h80aJCU7/7jtoNGuVcB\nyWMNA62DgARsNpvNYrGUlpYOGDBAegwJgvD444+L765atUrqVCkpKZGPSVgslttuu00ePwYNGiSP\nGcuWLbvxxhuldyMiIh5++GGnv2f9DEjudYqUlpaK87hERo8eTfOWm1twjDFeFabRbLPZ5s+fL40z\nzZ4929GU1Wp97rnn5J2Bv/zlL6VJBz4EJJvNlpOTIxZ74YUXpEyPH7cdHpV7G5Dc1zDQOj1sNpvT\npwAITYxGY319/U033TR06FD5aseOjo4jR45kZGT06dPH8aru7u7Tp0+fP39++PDhjqNQhJBvv/32\n3Llzw4YN82pYyAfc6xQ5e/bsP/7xj4yMDPnSV49veXULXhWm0UyD+ClcunQpMzMzoFsRevy4FUCx\nbxRQEgQkAAAAXMDFLLvu7u7jx483NzdbrVa7BfNNTU1btmwxmUy5ubm5ubmUbwEAANAcXJyHtGzZ\nsv/+7//eunXr8uXL5fl6vb6goCAhISEjI6O0tLSyspLmLQAAAFqEiy47i8UiCMKBAweeeuqp+vp6\nKX/OnDlDhgx59tlnCSEHDhwoKSk5duyYOALs5i0AAABahIsWkjjV1ZGamprx48eL6ezsbLPZLO16\n4uYtAAAAWoSLgOQUk8nU1dWVlJQkvtTpdJGRke3t7e7fAgAAoFG4mNTgFLEvUdy3WEQQBKvV6v4t\nOfJNtwAAAOj1erUluIPfgCT24zU2No4dO1bM6ezsFJf+uXnLDs5r3xVpaWlQrjBQrjxQrjD8/0bn\nt8tOEITExERpny7xSE3xwEo3b3kLzcgT5egUQ1MKu4NySqCcrSmF3UG5BlB1n4h/Y7Var1+//vnn\nn48YMeL69evXr18X89etW3f//fd3dnbabLbS0tIZM2ZIl7h5SyI1NVUR+eyBcuWBcuWBcoXhXzYX\nXXa7d+9++umnxfTIkSMJIfX19YIgFBUV6fX6zMzMqKio6Ojo8vJy6RI3bwUBn332mdoSfATKlQfK\nlUe7yjmHi4A0bdo06eBIOYIgSGd80b8FAABAi/A7hqQM2u3nhXK2phR2B+WUQDlbU5zDxU4NAUKj\nM2EIIQaDISUlRW0VvgDlygPlyqNR5fw/EkO9hQQAAIATgjwgyZuxSCONNNKhnNYAak/zCyA0cxxr\namqYlGFr6uzZs0q6g3IblEO5N6Y0qpz/ad8YQ+IRjfZQEyhXAyhXHo0q5/+RGORddgAAALQCAhIA\nAAAuCPWApN21AlDO1pTC7qCcEihna4pzMIbEIxrtoSZQrgZQrjwaVc7/IzHUW0gAAAA4AQEJAAAA\nF4R6QNJuPy+UszWlsDsopwTK2ZrinCAPSB5XL2dlZTEpQ+nXK23u01lZWUzKEELEYv6XoUyjzgnq\nHHVOXcYOVmX4BJMaeESjQ6YEytUAypVHo8r5fyQGeQsJAACAVgj1gKTdfl4oZ2tKYXdQTgmUszXF\nOeiy4xGNdggQKFcDKFcejSrn/5EY6i0kAAAAnICABAAAgAtCPSBpt58XytmaUtgdlFMC5WxNcQ7G\nkHhEoz3UBMrVAMqVR6PK+X8kBnkLickCN6SRRhrpIEjzD1pIPKLR318EytUAypVHo8r5fyQGeQvJ\nI9rt54VytqYUdgfllEA5W1OcgxYSj2j09xeBcjWAcuXRqHL+H4mh3kICAADACT3VFuCBL774Ys+e\nPVardeTIkY888kh4eLiY39TUtGXLFpPJlJubm5ubq65IAAAA/sN1C2nTpk3PPffcyJEjs7Ozd+zY\n8cQTT4j5er2+oKAgISEhIyOjtLS0srLSZxfa7eeFcramFHYH5ZRAOVtTvGPjmClTpmzZskVMf/vt\nt6mpqdeuXbPZbE8++eTLL78s5u/fv3/06NFdXV2Ol5OnqxSTypazZ8+qLcFHoFx5oFx5NKo8NTVV\nbQke4LqFNHDgwI6ODjHd2dnZs2dPscuupqZm/PjxYn52drbZbA6SXwcAABDCcD2G9MILL/z+978/\nc+aMIAgnT55cuXJlWFiYyWTq6upKSkoSy+h0usjIyPb2dnWlAgAA8BOuA9L3339/5coVQkjv3r07\nOztbW1sJITabjRASHx8vFRMEwWq1OrWQlpYmJj777DOnBWpra9PT093LoCnD1lRzc7OS7qCcQDmU\ne2NKW8rz8vI8uuAFtfsMXWK1WtPT03fu3Cm+vHjx4rBhw+rr669fv56amnrkyBGp5OjRo/fu3eto\nAWNIygPlygPlyqNR5RhD8h2z2Xzt2rWBAweKL+Pi4sLDw8+dOycIQmJiothaIoQYjUaTyTR06FD1\nlAIAAGAAvwEpIiIiISFh79694ssDBw6YTKZbb72VEJKfn19RUWE2mwkhGzduTE9P1+KqaQAAAHL4\nDUiEkNdee626unrMmDE5OTnz5s1btmzZLbfcQggpKioaNGhQZmbmhAkTDh8+vGrVKp9daHetAJSz\nNaWwOyinBMrZmuKcYN7LrseCaltZjtoqfEGjO2URKFcDKFcejSrHXnYAAAAAFQhIAAAAuCDUA5J2\n+3mhnK0phd1BOSVQztYU5wR5QPJ4rG9WVhaTMpR+vdLmPp2VlcWkDCFELOZ/Gco06pygzlHn1GXs\nYFWGTzCpgUc0OmRKoFwNoFx5NKockxoAAAAAKoI8IPVYUO2+gHb7eaGcrSmF3UE5JVDO1hTnBHmX\nHSFEi712Gu0QIFCuBlCuPBpVji47AAAAgAoEJAAAAFwQ6gFJu/28UM7WlMLuoJwSKGdrinMwhsQj\nGu2hJlCuBlCuPBpVjjEkAAAAgIpQCUg9FlQzWX2NNNJII63RtAZQ98DagEKerpJOMXd1nHlNTY1H\nOzRl2JqiOSAZytmagnIopzelUeX8H2EeKmNI2tpGSKM91ATK1QDKlUejyjGGBAAAAFCBgAQAAIAL\nQjEgyTe40+5aAShna0phd1BOCZSzNcU5oTiGxP94kkZ7qAmUqwGUK49GlWMMCQAAAKACAQkAAAAX\nhFBAcrpY7PDhwx4XlNGUIYRs3ryZSRnKtFyeP2Xs/venDGUadU5Q56hz6jKEXZ3zD8aQeESjPdQE\nytUAypVHo8oxhgQAAABQgYAEAACAC0I9INF0sFJ2wjI0pbA7KKcEytmaUtgdlPMP72NIVqt127Zt\ntbW1N9xww5QpU+6++24xv6mpacuWLSaTKTc3Nzc31+m1GENSHihXHihXHo0qxxiSX1gslpkzZ27f\nvn3UqFFJSUkff/yxmK/X6wsKChISEjIyMkpLSysrK9XVCQAAwH96qi3AHZs3b75+/fpHH32k0/0s\ncJaVlRUWFhYXFxNCBgwYUFJSMmvWrLCwMG/t899UAgCA0IHrFtL27dt/85vfGI3GQ4cOXb58Wcqv\nqakZP368mM7OzjabzT73n2q3nxfK2ZpS2B2UUwLlbE3xjrrHMbmhq6srNTX1d7/7XXZ29mOPPTZ8\n+PDNmzfbbLaOjo7U1FSDwSCVHDNmzCeffOJowdUBfR5P7VMdmuO/+ATKlQfKlUejyvk/oI/fLrvu\n7m5CyPnz56uqqgRBOHr06MyZM6dMmTJw4EBCSHx8vFRSEASr1erKTlpaGiGE/OpNg8EgZUppeSY/\nNDc3qy3BR6BceaBcebSlPC8vT20JtPAbkMLCwnQ63UMPPSQIAiFk7Nix0dHRDQ0NN998MyGksbFx\n7NixYsnOzs6IiAhXdsRZJT0WVMtmxUgzZAxD3jDweaSsFufwiEC58kC58mhIuTSz7t+/zjmG3zEk\nnU43dOhQedNHbDMJgpCYmNja2ipmGo1Gk8k0dOhQSrPyw5Aoy/DZz6vdHmooZ2tKYXdQTol2lauJ\n2n2G7qioqLj33nuvXbtms9mqq6uHDRv2z3/+02azrVu37v777+/s7LTZbKWlpTNmzHB6udMxJLtM\n94NMaqHRHmoblKsBlCuPRpVjDMkvHn/88dOnT48bN65v375Xr15dvXr14MGDCSFFRUV6vT4zMzMq\nKio6Orq8vFxtpQAAAPyF64BECFm5cuXKlSvtMgVBWL9+vT9m3YwY2fXX8Ta2BAAAwQq/Y0jM8Th6\n5KoAn/282u2hhnK2phR2B+WUaFe5ivC+l50/2O1lJ6alqCOlnRZQd+M7je6URaBcDaBceTSqHHvZ\nAQAAAFQgILE5rhhppJFGmvO0BlB7ml8AcTrDW0zI004LyI3U1NTQuKMpRmmKZlIpQ3dQboNyKPfG\nlEaV8z/tG2NIGENiCZQrD5Qrj0aVYwwJAAAAoAIBCQAAABcgIHlGHBX0uIwJqxwogXK2phR2B+WU\naFe5imAMiXYMST6YFOiBJY32UBMoVwMoVx6NKscYEgAAAEAFApI7aM6qAAAAwAQEJM/w2c+r3R5q\nKGdrSmF3UE6JdpWrCAKS5xXOWVlZrsr0WFDttBVFs1Laz9XXWVlZTMqIN8ikDGXaTX16VcYO1Lmb\nNOqcoM61ACY1uJvUYJd2OtPBaQE/0eiQKYFyNYBy5dGockxqAAAAAKgI/oDkqleNHoW7lRma4rOH\nGsrZmlLYHZRTol3lKhL8Acl/pB5qzZmidMfQlMLuoJwSKGdrSmF3CitXEfZHmF++fPn8+fPr16+/\ndOnSL37xiylTpiQlJel0iHwAAADcwTIg/d///d8f//jHH374gRDSp08fnU732muvvfLKK4SQESNG\nbNq0KS4ujqE7AAAAwQSbhovRaMzIyHjppZeeeOKJL7/8Uq/XHzt27MiRI3V1dXV1dXv27ImKipo4\nceKcOXOYuFMYHzpnpVEru+Er9FCzdQfllEA5W1MKuwudMSQ2LSSj0fjWW2+lp6c7vhUeHp6cnFxZ\nWWm1WtetW8fEncJgDInelMLuoJwSKGdrSmF3oTOGxKaFNHz4cKfRSE5YWNjTTz/NxB1b/FkEF2j7\nSCONNNKs0vwT/Atj5fi5MJY4rIGlWTnrwx7hGl12R6BcDaBceTSqPMgXxl66dEmv17e2tsoza2tr\n/ZOkAu4XKincrczQFJ891FDO1pTC7qCcEu0qVxHfA9Lx48ezsrJmz549ZcqUhQsXSvmFhYUshHEE\nxpDoTSnsDsopgXK2phR2hzEkzyxcuHDdunWHDh3S6/Xff//9rFmzGMoCAAAQavgekP71r39NnjxZ\nTG/ZsqV///6BaxsdP378o48+MhqNUk5TU9Py5csXL168b9++QHgM9ElIOGkJAADs8D0g9enTp729\nXXq5bt26xMTEQMQko9G4aNGiJUuW/OMf/xBz9Hp9QUFBQkJCRkZGaWlpZWUlc6dyMIZEb0phd1BO\nCZSzNaWwO4wheeb+++9fu3atPGf16tWDBg3yW5I9S5YsmTdvnjynrKyssLCwuLh4xowZK1asWLt2\nrdVqZe5XIisri1WDBj3UbN1BOSVQztaUwu4whuSZhQsXvvTSS3aZr7766rFjx/yT9DM+/vhjQsi0\nadPkmTU1NePHjxfT2dnZZrNZu78O5KEO/XgAgFDGr50awsLCHDMjIyPlLy0WiyAIvtm/dOnSa6+9\n9v7778szTSZTV1dXUlKS+FKn00VGRso7DwEAAGgRxrt9d3d35+fnv/nmm4mJiYSQwsLCY8eOCYJQ\nUVFx5513emuttLT0iSeeSEhIsFgsUqa4kjc+Pl7KEQSBssvOYDA4pp1mytM7duwgpC/NVR4zd+zY\nkZ6ePuQNg/urhrxhODvXw73U1tZ63B2DpgxbU83NzUq6g3IC5VDuqUxeXp5HF5zAOCC1tLS0traK\n0eiLL744duxYeXn5iRMnioqKjh8/7pWpr7766ujRo9OnTz9w4EBXVxchpLa29sYbb0xOTiaENDY2\njh07VizZ2dkZERFBYzMlJYUQg13aaaY8PX36dFJTLcskhEjrtA3OMn9WQJ45ffp0769ycy809+sZ\nhqZoSkI5W1M0JaGcrSmakvwol3ZnSEtLo3GnIowD0o8//hgVFSWmX3311fT09MmTJ0+ePHnDhg2X\nL1+OiYmhNxUWFjZixIj33nuP/NQqqqqqioyMvPXWWxMTE6XtIYxGo8lkGjp0KNsbcYO3Iz30OwYB\nAEAow/6APmlzvAsXLpSVlYnp8PBwsZVDz9ixY6U2kMViGTly5MKFC8Wc/Pz8ioqKe+65Jzw8fOPG\njenp6VrcVwoAAIAcxge5xsXFtba21tXV7dy589q1a2L3msViMZvNffr0YeWlqKho0KBBmZmZEyZM\nOHz48KpVq1hZdoq3U/jcNKGwyoGtOyinBMrZmlLYXeisQ2LcQoqLi5s5c+bDDz9MCPmv//ovccbd\n/v37e/fuTTnM4xRBEOSb1AqCsH79ev/VUpKVlUW2Yx0SlEM5Y1MKu4Ny/mHcQiKELFu27MSJE3V1\ndc8++6yYM2bMmABt8BMKuGpvYdESACDIYBaQxo0bN2fOnHPnzhFCIiIiwsPDpbf69+/fv39/Vo5C\nGQQhAEAQwywgvf/++x0dHdOmTbvtttumTZt2/Phx+eIhTcPnXnYS/g9ZoW+dEihna0phd1DOP8wC\nUkpKypYtW+rq6j7//POxY8cWFxePHDkyJyfn/fff52cbBafPbo9H/9J3znrVgnHvVzTl7bHEdmUk\n5e7tiMX8L0OZzsrKYlLGDn/KUKZp6hN1TlkGdU5Zxg5WZfgkgEeYX7p06ZNPPqmoqLhw4UK/fv2m\nTZv2/PPPB8iXU9yHB/lx5k4z5ceZezz4XJ5mcvC5bwVUR6NHOxMoVwMoV5ggP8LcPf3793/ssccO\nHjx44sSJV1999cCBA42NjYFzF8Rg6AgAEAoEMCBJREREZGdn79u3b/jw4Qq4Yw7nY0j+u0PfOiVQ\nztaUwu6gnH8CEpDa29sNBkPTzwmEI2VQeGmCzzg2pER3HhtYWJ9BCZSzNaWwOyjnH/ZbB/3P//zP\nnj177DJ1Ot2pU6eY+wLewtWAEwAAyGEckL777rs9e/a89dZb48aN8/kYJE7gbeSGRo/TGRAAAKAJ\nGHfZdXR0JCYmZmdnaz0ayVGmWzkQ8Y/PHmrt9q1DOVtTCruDcv5hHJASEhL4WXXECq2MIfnsDn3r\nlEA5W1MKu4Ny/mEckPr37//rX/968eLF7e3t1p/D1hHwE946JAEAgP0su6ysrF27dmVmZg6XMXLk\nSOaOAkGgm73+rAB3g2N0scthsvIcaaSR1nSafxjv1NDe3p6ZmTl58uRFixbpdD+LdrfccgtDRzT4\n0Ahwun1DTUGvCds77Qr4tlNDTUGvrKwsb6/yuYDkzuleD1LmYdnmK/Lak8+JcFrGEZoV7DSmKN0x\nNAXlUE5vSqPK+d+pgfEsu+bm5tjY2PLycrZm1YXJeUhihOB/DMnV3Dz0rbN1B+WUQDlbU5zDuMsu\nMjIyLCyMrU0+CaYxGI/3Ekw3CwDgFsYBKTk5OSYm5tChQ2zNAgAACHoYBySj0Xj16tXZs2dPmDDh\nbhm5ublsHSmJwksTGOLenc+HZfiJdtdnQDlbUwq7g3L+YT/Lrk+fPsOGDYuNje0jo3fv3swdKUbQ\nr0Py35RXsU27fetQztaUwu6gnH8YT2qIi4vbtWsXW5uaIERGWTxuRzTkDYOtTHvnxAAAeECJ4yeC\nlRAJQq4I8dsHADCHTUAyGAyXLl3yWOzTTz9l4i5wOO2HZdg5u3nzZqdmA7SoTbqc0pcUYxzL0Cvx\n6Ovw4cNMyhC6+vSnzu3SNPVJWeeOtepbGco06pygzrUAm4WxjY2N+fn5KSkpZWVlqampdjurXrp0\naevWre+88050dPQXX3zhvztKWC2M9ViA5oxzymWt/lvw5xB0Jmera3S7cY0eSk2gXA00qjxUFsYO\nHz781KlTK1aseOihh7q7u6Ojo6Oionr06GGxWH744QeLxRIbG/vaa69lZ2czcaddmD+jtbiEyGn0\nAgAAZmNIOp1u2bJlp06d2rNnT0FBQXh4uCAI/fr1W7t27eHDh2tqahCNmMBhgAEAACawn9SQnJz8\n7LPP7t69+7PPPvvLX/6Sm5vbv39/5l4ADQovqnCFPIi6CqhudoP1U5XClQDlQPWlJwAAIABJREFU\n9KYUdgfl/MP+CHOGNDU17du3z2AwREZGPvDAA2PGjJG/tWXLFpPJlJubq+lVtwGFyS58/zHlB+57\n5xyjlMfNXulVYWUJW3dQTol2lasI19O+CwsLDQbDnXfeKQjCb37zm507d4r5er2+oKAgISEhIyOj\ntLS0srKSlcfg7hDTxN3RNKcAAEEJ1y2kqqqq6OhoMX3jjTdu2LAhPz+fEFJWVlZYWFhcXEwIGTBg\nQElJyaxZs0JkU1f+CfRUBUyFACBY4bqFJEUjQkh8fHxXV5eYrqmpGT9+vJjOzs42m83B0X/KHKla\n/G9qKFzDlO5o7gujAmzdQTkl2lWuIlwHJAmLxfLuu+9Onz6dEGIymbq6upKSksS3dDpdZGRke3u7\nqgI5RZkOcTdRwbdA2GNBtT/K7Tr9MCrA1h2UU6Jd5SoSkC47q9UqLj+SZyYmJvpscMGCBbGxsWIf\nnbiSNz4+XnpXEASr1eqzca8wGAzuM6W000wFCjC8qseC6rNzU/zxpaRaOY4FhrxhkO7FKc3NzW7e\n5RkoVx5tKc/Ly1NbAi3sA9KaNWs2bdpkl6nT6U6dOuWbwYULF168ePHtt98WR4nEbSAaGxvHjh0r\nFujs7IyIiPBDshcMecPJ4y8lJYUQg136p4XcBn8KSD/2XVvw3y8hxOA002MBj1ex8/uzws4y5WNL\nHq9ygqa3hdXirgEiUK4A0u4MaWlp6irxCOMuuwsXLmzatGnp0qV/+9vfamUcO3bMN4OLFy/+9ttv\ny8vLIyMjxRxBEBITE1tbW8WXRqPRZDINHTqUzQ0EF5ysQ1LAnZu+QfemnF7oypp2RwWgnK0phd1h\nDMlH2traEhISZs2aFRMTE/lzfLD2/PPPnzx5cvPmzRERERaLReoDzM/Pr6ioMJvNhJCNGzemp6dr\n6NeKkuAkJ6emfJ7iod1RAShna0phdxhD8pHevXszHM758MMPCSETJ04UXwqCUF9fTwgpKirS6/WZ\nmZlRUVHR0dHl5eWsPAI+UWtBknwzWbt8x91mAQB+wjggJScnDx48uKqq6u677/bfmquNaQVBWL9+\nvf/2Ac8oH4T8XOHk1eVYTQWAI4y77IxGo9FoLC4unjBhwt0ysLuPD3CyeEiUwfkYEo0pP+uzx4Jq\n35Y9SVd5K0A0xSQwa3c8A8rZmuIc9uuQ+vTpM2zYsNjY2D4yevfuzdwR5/Cw7Q3GkNybCsRnJLqj\nt+xmt1kfRr9cBT/tjmdAOVtTnMO4yy4uLm7Xrl1sbQYxPAQtwA++9eO5mrDu/twpn49SRGcjDdo9\nqVJdArVTg9FoPHPmjMFgwB4KrqA5ulgtvwx37nFlX7y2x4LqQPhlcgS1DzjacezwZKLBrv+Q3qZc\nj1NtTo+xdzU53gdfknJXfl3truvUl1yDlHZV53IN9H7lBaT0kDcM9H4dfbmq50Cn+Yd9QLp8+fL4\n8eMnTpx477335uXlZWZmPvTQQ8y9cAXlGT92yJvYrtJ+inH6XaTx66jB/6+1V379dOeqDn27d6/8\nSsq99eWYph/9oqlPb8tM2N5J71feweimt9AuLS0zt6tzp9/hCds7HX2Rn+pcvMTRl5Qv2XTzWduF\nWPc4tSPXJtane780nwuNXx/K8An7gJSdnd23b9/q6mpxSewHH3zQ0NAgbkMHFAZjSGxNUbrztifW\nVXk+xzPc4xgY/MSpKR9GyPxx5wY3nzVD5aEzhsQ4IJ05c6Znz567d+++6aabxPWw6enpX331VUND\nQ0dHB1tfQEkw3EUDq1piFdIC7ZdbguZGQg3GAamzs7Nfv352mdHR0YIgXLt2ja0v4Aof/hqD8g84\nKG/KNwIUsehn/bH1C4IVxgEpOjq6paWlu7tbntnQ0GCxWPr27cvWF/AITYd4gP7U/TTr5xiSx4E0\nVsYdYVjnHCqnRLGtCH1zR2/Kn68x1iH5AOOANHjw4EGDBo0aNeqjjz46d+7cuXPn1qxZM3369Pvv\nv1/cpRsoiZLdyj786bqZC0ujnNKjb5UQ0P3ulDelsDv+R7+kSZ6uTPn/Qw3rkHyA/fETVVVVjz/+\n+JIlS/7toGfP3/72t7///e+ZOwI+gw4QtQjokhQtfqxYowPkBOSAvrfffpsQcvny5bCwMPkx5CGL\nFp8UcmxlOVq/BTvwHFQSf748+KRCigAeYR4TE4NopC6KdSszf2SoO57hz+2oqNwfvFqC45Eg2Pkw\noKacror1zZRXZfiHTUBqbGzMycn57rvvjEZjjjOYbP4dZPizSp/ShdNu5UCs7vb2D8ZjOkAd4gqs\nZne6tte3OvetEnz2696dm2+d++13A13n8rWlNLtOuPdLX+ced5GQr5J25ddNGW/r09syfMKshRQb\nGyslHImLi2PlKGgI0E4NTPxO2N7pVSvB20FXJe9dK3693ZWVrV//CTW/bjSIjWzHHSVorg1omn/Y\njCENHz5827ZtYlpKAPcE2ahMgMAQAgChA/vzkObPn2+X2d7e/thjj7F1FGQEKDhx1Vr3YSdWJqFI\nK+MZjqGX4ZoYGvgcifFIKIx+YQzJd/7+979TZoJAI7bW1Wph+ODX6U6U3lqw84vVPO6RasyVOx8a\nqfzUubcTRPlR7m0xbXXNuSKAs+wkzGazTqeEIx7w/2esJrry1BLp7cMxQBPWg28evJ+oUiHozg0+\nmMWJRYsW5eTkzJgxo62tzW6K3aRJk26//XZWjoAEnok+EKCnmEezTPx6fO5zGyl9EBbc8YbbT0pd\nmAWkG2+8sW/fvuKGdX1l9OvX7+WXX/7ggw9YOdI6Sn4L3XQrc/unLgrTbod40I9nKGxKYXeUWxGe\nnevkiN4AuaMspt0/GTnMAtLzzz+/Y8eOTZs23X333TtkfPTRRw8++CArL8ARP3eEc2+BFZwPQjBE\nc2NIgXCnaeVKtswwhmQH46GduLi4DRs2sLUJ2IKOAgAAn4TKXANN4NvvMn+a6qxWdPuzs0MgdoXw\n1i9zDXYjBEruiOGqPG9+A9EN5YNfp2M58jI1Bb08/obz53792dEjELuuqEtAAtLLL788adKkCTKy\ns7MD4SjUcBqx3DTVPX4XxWvFYv40+b291uNKcjfK5eVpHhY0fv3Rb4erHWJ88Ov/ZkLyDihWde4t\n/tQ5pV+nux/5U+f+6PTKr3zfI0q/rMrwCfuANGnSpB07dgwePLitrS0pKalXr15tbW1Dhgxh7gh4\nRP5ddPPg5vAry7Bv3T3MJzthJIbSlN3UAH8+iKCvc4wh+UhLS8uFCxeOHDny+uuvx8TEbN26taqq\naunSpSaTia0jLcLtxDbfCLLbCREw2xjwDOOA9K9//eumm24ihOh0OovFImbOmjWroaGho6ODra+Q\nIjgeItq9i+BeExMIuK0xtYRxWyFcwTgg3XDDDWKiV69eV69e7e7uFl92dXVdu3aNoaOmpqbly5cv\nXrx43759DM0GGQxXOShM0CtX2JTC7lRUzs9xVgzdYR2SjyQkJLS0tBBCIiIi+vXr94c//KGlpWX5\n8uWEEIYnUOj1+oKCgoSEhIyMjNLS0srKSlaWtY5v27hx2Psc9Mo9YivLmbC9k4kpSqSdD9kegaEA\nGEOiN8U5jI8wj46Ofu+998T0tm3bHnzwwQ8++EAQhPLycoZeysrKCgsLi4uLCSEDBgwoKSmZNWtW\nWFgYQxchi4qb02i3Q88HxGpEH04QIP1FYHzOf9jPsktJ+ffMmcGDBx87dqyurq6+vn7y5MkMXdTU\n1IwfP15MZ2dnm83m4Giu+gNXfwluxGjuEYynDACKwbiF1NDQMHv27L/+9a9STnh4OFsXJpOpq6sr\nKSlJfKnT6SIjI9vb29l6CQQGgyEQBeSZdgV27NhBSF83V4n/y4tJheUvHf93/5arAvRXUUqyu0Gf\nNdPIo3xLMeViYsgbhoMPCpP+YvH/fimV231kTgs4mnJqobm5mf5+z85Nkc8Ud1TO5Avgs3Lf/O7Y\nsSM9PZ14ora21mMxN2Xy8vI8uuAExi2k3r17szXoiM1mI4TEx8dLOYIgWK3WQPv1H6ntyLaAPFNK\niw2R6dOnu79K/F9eTCosf+n4v/u3XBVwc5XUcvJKkt0N+qyZ5qYo32KrfMgbBrFm3PgdNGgQk/ul\nVC4m3BdwNKWAciZfAJ+V++Z3+vTpKRTQFHNTRv8ThHsYB6Tk5OT4+PhPP/2UrVk5giAQQhobG6Wc\nzs7OiIiIwHn0Aenx6qqHSpmeK/Q1Ad9Q5jQNtqBzNQhg3GVnNBqvXr369NNPv/jii/LWkk6nYzU/\nWxCExMTE1tZWyaPJZBo6dCgT44Bb8LjhBHXPZBIDIb4JwQr7SQ19+vQZNmxYQkJCHxlsu/Ly8/Mr\nKirMZjMhZOPGjenp6VLTOPhQ4Hcoh1NCKFeWsFLO5wk37mH+XMbaL7amJNx8UliHZAfjFlJcXNyu\nXbvY2nSkqKhIr9dnZmZGRUVFR0eznVMegmRlZZHt1U6f+1498hhOfhUlsSqmJFDuD741gBhKcmUq\nQC0zrEOyg3FAkmhqarp27dodd9xBCBE3DYqMjGRlXBCE9evXs7IGGMKqncEWtfp50L8EgFew77Kr\nqqpKS0v71a9+NWPGDKPRSAj5/PPPNTTvEGgULT76Azo1QIsVohYcztEITRgHJKvVWlxcvHTpUr1e\nHxsbK2bm5uZeuHBB2msVqIjTvzrK3mcl/2g1NJ7h5iw+N/Cg3I6gV+7tUXsKgDEkOxgHpLNnz8bH\nx8+aNUueGRERIQiCJtau+oN2f2Fx2PvsVd86Vz9vNbELn9Ma04Ryp7BSrtb+gUyKcfi5+ADjgNTd\n3S2uE7LDarXqdEF+XDrbY4Ppn7CBO5raUUMgjqb2p3yg/XqrwSm+NaH89+vGjlp+xbT7vsTA+fWt\nPA9+cYS5j/Tr16+lpcWud66qqkqn08XExLD1xRtqHRscoKO4vfLF6shn3vx6q8FbO6Hml9VnHWp+\nWaX5h3FAiouLGz9+fEZGhrhZQ3Nz84svvlhcXDx//ny2jjhEW79ERMQfquoqd/pjORDjGcoM8vMw\nEuNbH6YCyimFeftJ8VDnvoExJDvYd6O98847999//+LFi9va2h555JHt27c/88wzs2fPZu6IN7T1\nS0SOP8p9e8R7vEoT4xkYibFDrcE8j8pFYUFZ596a4pyArENasWLFihUrfvzxR0JIVFRUIFwoDybR\nAgC8gqvpNpoggBMNoqKigiYasUX+HfXt+6rkVQrAiTC1fnNw6xc/wuxAjQUaNi2kxsbGuXPnuinQ\no0ePqqoqJr64Rbt9uBwqx6iA8kC58hw+fJimq42mGKUpzmHWQop1S1xcHCtH3KLdb4N75cybLzQG\nQ3kkRi2gXHkwhmQHmxbS8OHDt23bxsQUCDK024Oh1gAABh48gp6xYCVQY0hNTU0nTpwQ0x0dHeL+\nqsGNkovRDh8+LD2zAu2Xh4V4Pvj1+FhXckGuVvy61yBWKQ8Lcr1djKyYXx/Awlg52Fw1IPjZfPa4\nZ3ZWVpb0PfN/EZz7ryzzhXg+LKqgsR+IhaI+3K/HuxPL+1Pnvk1i9vZzYbW4UoHF15RtSrbfcybY\nDfx4+z33pwyfYHNVX3D67Vf4g2fojkPlfq4sUbHLKxDjGcr0UPmgXFtdZ8y/5z50rtrVGMaQ7MDm\nqsB3AvTc19ZjDriByUeJwbzQAZur+gXP31eP2ngWH8TgMQeAK7C5KjNc9VAH6OnDcKxS4WFPTpT7\nEBgYHqjDts6ZNES0u5qHErlyTlrhlBqwl52PhPLmqhyOxChvyil2z30NKffBnXbXxEzY3slwiRiH\neDtuFzgl3oIxJN8J2c1VgWJw9bCwAz1yAPgMNlcFAKgGJ11nCoOfLK7A5qrMCM2RGHXd+WzKtycC\nzXpMVqaYw6rpFgjlyjygHZVrJTBgDMl3zpw5Y5dz+fJlx8zgwNWXQIHV3VKXsf8rtLOystz8ZTJf\nGe5DZzeN/UDUuZ1NtnVO79ertD94a5/Sr8fnfoDu1+PgpT9+awp6Ufp1xJ97VOtvQRkYBySj0Xjv\nvffaZXZ1dTlmBgdMVnf78DNNyVXlHu3bLVBl5VdeLWod26y8X/Gumful/I6FTj3775fV7hVKpvmH\nfUC66aab7DLj4uJ69ux56dIltr6AHaF5tBI/SmjgVi23wkBIEYSLVdUiBEdiPOL0MacJ5YF2B+WU\nuOkZE+F23E5h5drqmnMF+3VILS0tdrsEVVVVdXV19e/fn62vwOHbr0V5FxZTOR7ccWVKYXeOpgJa\n+aI7Ji5Q5z67U92UD/2ffhI665AYT/uOi4sbN25cZmbmM888k5ubSwh56623PvjggyeffNIHa01N\nTfv27TMYDJGRkQ888MCYMWPkb23ZssVkMuXm5oqO/ITn6ac8a9Mc2u2bUmsvKNEsh6d3a/ejBK5g\n32VXWVl53333rVq1aurUqVOnTt25c2dJScnChQt9MFVYWGgwGO68805BEH7zm9/s3LlTzNfr9QUF\nBQkJCRkZGaWlpZWVlUzvAAB+H3bcCmMOqzvlZFNHp15C59OkJCALY1evXr169er29vYePXr4sxSp\nqqoqOjpaTN94440bNmzIz88nhJSVlRUWFhYXFxNCBgwYUFJSMmvWrLCwMCbifUbTowIe2/s1Bb0m\nbO9k5Y6JHUKnnCEBHYkJ6LNJu2NIzJX7EJ+8ah3KT86U0nKb8jTlF5immMJ/CwEigOuQoqOjxWjk\n8zokKRoRQuLj47u6usR0TU3N+PHjxXR2drbZbOZhQC+YRgUUc+fnU5jnOre7Ncc1MQHtYbODh5EY\nP4dm/cfOlMcmi9O0xwKu3NFI8qdYEEQjwryFJK5D0uv18kxxHZJdpldYLJZ33313+vTphBCTydTV\n1ZWUlCS+pdPpIiMj/T9syWAwOE27L0Cf6VWB5uZmxfyyvcqpcqeFz85NUUbY2bkpQ94w0CtXTBir\nApJy3yw43i9ljfn5+Z6dm2KnXO5XfrlHYVJh8ahlRwvyAk7/ul3lu0L+PecfDR3YzT4guVmH5H6i\nXXd3t9VqFdN2hyotWLAgNjZW7KOz2WyEkPj4eOldQRCkC30mJSWFEINj2n0B15nERZqqwKBBg1wX\n9t2s0wI/JeRpzwVcXeVUOaVZW1lKjwXVdgXsMn0S5p1yygqRC3OqnIkw2e9uL5Tbyv5T0lmNeaxG\nD8Jk9j1fJS8sFyYWcK2cOLvqP2mnJZVEdte8IzUG0tLS1FXiEY7WIe3duzf9J+QnKi1cuPDixYtv\nvvmmOEokxqrGxkapQGdnZ0REBL0jtl0l8i5jJnYoYbVrC6UpHgYhbGU5drfDSrnTDhnfTFH6ojHl\ncSGLb+5cdTRRrpuxU+60xrxaguN/nTM8yUnh7zmfylWEo3VIeXl59T8htZAWL1787bfflpeXR0ZG\nijmCICQmJra2toovjUajyWQaOnQo2xvxgWA6Vci3DnGf3QXClHQLDEdr3AxCeOvC21EBV/bdh1JX\nleBxnMmNWco6d1XGW7Q7EqNd5SrCPiCJ65Defvvtc+fOnTt3bvny5cXFxb6tQ3r++edPnjy5efPm\niIgIi8UiNZvy8/MrKirMZjMhZOPGjenp6Zw0n3mYxOnVUC1bvx69eBwZ9sqdD1d5ZcrjmDZbFx4L\nsxIDALewn/ZdWVm5cOHCVatWvfLKK4SQ8PDwkpIScfjHWz788ENCyMSJE8WXgiDU19cTQoqKivR6\nfWZmZlRUVHR0dHl5OTv5miRoVs76uVbDqxm9DM0q7AKAYCUgY0irV68+derUkSNHjh49WldX51s0\nIoTof44YjQghgiCsX7++rq6upqZm9+7dgwcPZqfdd7S+PoOhO/9HBSiHT2hWlvBZURjPoATK2Zri\nnABOapDWIYUIPK+J8dmU/6MCjv1L2u1bh3K2phR2B+X8E5CdGl5++eXdu3fLp2LrdLpDhw4FwpdW\nCII+GYYjQAAA4Aj7FtKkSZN27NgxePDgtra2pKSkXr16tbW1DRkyhLmjkCXQMQAxBgCgCowDUktL\ny4ULF44cOfL666/HxMRs3bq1qqpq6dKlJpOJrSNl8Pholvfbbt682Wl+gI6Xll5SHofsZr2UmCMf\ntnG0I//fK7+Oabl4f8oQujqnKUOZpqlzys+Fpj5R56hzEVZ1zj89xI0PWNHQ0DBv3rzq6urLly9P\nnTr12LFjYv6IESOOHDkirSVSBlcTz5zOSbPbHlFMO820KyA5ktL0ma4KnJ2bkpKS4qqAFFqcFnBV\nWBkMBgMnU/C9BcqVB8oVJi0tzZ8t3BSAcQvphhtuEBO9evW6evVqd3e3+LKrq+vatWtsfTFB4ec1\nesMAAMAVjANSQkJCS0sLISQiIqJfv35/+MMfWlpali9fTgiJi4tj60st+AkqnBz0AgAATGAckKKj\no9977z0xvW3btk8++SQnJ2f79u3crl312MFK/0yn6atl2J/rxpS32+vxo9xbd1BOCZSzNaWwO4WV\nqwj7ad9jx44VE4MHDz527JjZbA4PD2fuhRVZWVlkO5s9DmjXCjB15z5e8rnKQbvrM6CcrSmF3UE5\n/wR8t2+eoxGHuD/VzWmO+3wAANAKHB0/EYIgigAAgESoByQ+u4wVdgfllEA5W1MKu4Ny/gn1gORD\nxyvNojNWC9Oka21lOXY2JeVuFtbRlCE/VYL/ZSjTWVlZTMrY4U8ZyjTq3A7UuZs0h3XOP4wXxnIF\nzcJY92tg7VawulnWSugWxhLXS1zlaWlhrOb69DS6YJBAuRpAucKE3MJY4BuaCzwAAMCcUA9IHHYZ\nUwYn7fZQQzlbUwq7g3JKtKtcRUI9IPG57EDE/zVGfK5ygHK2phR2B+WUaFe5ioR6QFIRdNMBAIAc\nBCQAAABcEOoBic8uY4XdQTklUM7WlMLuoJx/Qj0gBa7L2GmPHHqo2bqDckqgnK0phd1hDAl4BwaE\nAADATxCQvIZ+tbYYpZis+kYaaaSR9jOtAWzBC3m6yuk/+Vs1NTWOmVJabsdVppi22Ww1NTVS2lGJ\nmBDLeOTs2bMey9CYonTH0BSUQzm9KShXWHlqaiqNOxUJ9RaSzx2vjn106KGmBMrZmlLYHZRTol3l\nKhLkAUkKG0zGeHw2ghEmAADwSJAHJOVB7AEAAN/QRkA6fvz4Rx99ZDQapZympqbly5cvXrx43759\nPpu1O9PBT7DKgRIoZ2tKYXdQTol2lauIBgKS0WhctGjRkiVL/vGPf4g5er2+oKAgISEhIyOjtLS0\nsrLSZ+N8dhkr7A7KKYFytqYUdgfl/KOBgLRkyZJ58+bJc8rKygoLC4uLi2fMmLFixYq1a9darVa1\n5AEAAGAC7wHp448/JoRMmzZNnllTUzN+/HgxnZ2dbTabmTdXMRQEAAAK01NtAe64dOnSa6+99v77\n78szTSZTV1dXUlKS+FKn00VGRra3tzu1kJaWRn71JiHEYDBImVLaYDDU1tYS0tdNAfrM2tra9PR0\n93dEU4YQ0tzc7LEMQ3dQTqAcyr0xpS3leXl5Hl1wAkcBqbu7W+p5EwSBEFJaWvrEE08kJCRYLBap\nmM1mI4TEx8dLOYIguOqy0+v14lHiKSkphPw7eEjplJSUlJQUUuO0ACFEOqXY4CzzZwX+bcoT9Mce\neyzJ0B2UU5aEcramaEpCuf+mpGPL09LSaNypCEdddnv37k3/CYvF8tVXXx09enTQoEEHDhw4ePAg\nIaS2trapqUmMVY2NjdKFnZ2dERERqukGAADAAo5aSHl5efKmZVhY2IgRI9577z3yU6uoqqoqMjLy\n1ltvTUxMbG1tFYsZjUaTyTR06FBVNItgwAkAAPyHoxaSHWPHji3/iQ0bNhBCFi5cOHPmTEJIfn5+\nRUWF2WwmhGzcuDE9PZ2+oW2Ht7Mh3MQerHKgBMrZmlLYHZRTol3laqL2ZnpUXL9+PTU19ciRI9LL\n4uLi22+/PSsrKy8v75///KfTq+x2RPW4d6rt51ujOi3gaD8Q0GzdyCdQrjxQrjwaVc7/5qocddm5\nQRAEaVxOfLl+/XoV9QAAAGCONgKSurjqpsPQEQAAMITfMSRl0G4/L5SzNaWwOyinBMrZmuKcHjab\nTW0NgaLHgmpbWY64DklszUhpp5nSJfLLVVFuMBh8nqahLlCuPFCuPBpVnpaWJh/74JAgbyGxOsqX\nh+OHkUYaaaT9SfMPWkhoIbEEypUHypVHo8rRQuIOuxij3X5eKGdrSmF3UE4JlLM1xTkh10Jy2ipy\nLEDQQvIJKFceKFcejSpHCwkAAACgAgEJAAAAF4R6QNJuPy+UszWlsDsopwTK2ZriHIwhYQyJJVCu\nPFCuPBpVjjEkDYCdgQAAgAeCPCDRLBCjyUcaaaSR1npaA6i72XhAcXrShJQQ0zU1NXb5jnbEMh6h\nKUZpimZze4buoNwG5VDujSmNKuf/+AmMITm5RHGl9mi0h5pAuRpAufJoVDnGkAAAAAAqEJB+Bg/N\nIwAACE1CPSBpd60AlLM1pbA7KKcEytma4pwQHUMi3AwXOUWjPdQEytUAypVHo8oxhgQAAABQEboB\nidvmEQAAhCahG5BEtNvPC+VsTSnsDsopgXK2pjgnyAOSx9XLWVlZTMpQ+vVKm/t0VlYWkzKEELGY\n/2Uo06hzgjpHnVOXsYNVGT4JiUkNmpjIIEejQ6YEytUAypVHo8oxqUFNUj8pUlsCAAAAWoI5INGg\n3X5eKGdrSmF3UE4JlLM1xTnB3GUntk/RZackUK48UK48GlXOf5ddT7UFeMBqtW7btq22tvaGG26Y\nMmXK3XffLeY3NTVt2bLFZDLl5ubm5uaqKxIAAID/cN1lZ7FYZs6cuX379lGjRiUlJX388cdivl6v\nLygoSEhIyMjIKC0traysVFcnAAAA/+E6IG3evPn69evbtm2bNWvWk08+uW7dOjG/rKyssLCwuLh4\nxowZK1asWLt2rdVq9c2Fdvt5oZytKYXdQTklUM7WFO+oexyTe3JycnanMxtSAAANbUlEQVTs2PH9\n998fPHjwhx9+kPKHDx++f/9+MW21WocNG3bw4EHHy8XTqORn7jk9f49DaI7/4hMoVx4oVx6NKuf/\ngD5+x5CsVmtzc/PevXvXrl07ZMiQI0eOzJ8/f/bs2SaTqaurKykpSSym0+kiIyPb29udGklLSyO/\nejMtLe2zzz4TcwwGg0I34AfNzc1qS/ARKFceKFcebSnPy8tTWwItHAWk7u5uqedNEITu7m5CyPnz\n56uqqgRBOHr06MyZM6dMmTJw4EBCSHx8vHShIAiuuuzEWXbSxBJbmWYmxmhxDo8IlCsPlCuPhpRL\nD8C0tDR1lXiEozGkvXv3pv+ExWIJCwvT6XQPPfSQIAiEkLFjx0ZHRzc0NIgvGxsbpQs7OzsjIiJ8\nc6rdfl4oZ2tKYXdQTgmUszXFO2r3GbrjV7/61TvvvCO9HDNmzK5du2w221133SUmbDbbxYsXU1NT\nnXbpOo4haQWN9lDboFwNoFx5NKqc/zEkjlpIjuTn52/btq2jo4MQ8sUXX3R0dKSnp4v5FRUVZrOZ\nELJx48b09HQNNZ8BAAA4haMxJEcef/zx06dPjxs3rm/fvlevXl29evXgwYMJIUVFRXq9PjMzMyoq\nKjo6ury83I0RTWzNAAAAgOsWEiFk5cqVdXV1Bw8ePH78+LRp08RMQRDWr19fV1dXU1Oze/duMUr5\nhnb7eaGcrSmF3UE5JVDO1hTnBP9edmqr8AWN7pRFoFwNoFx5NKqc/0ci7y0kAAAAIQICEgAAAC4I\n8oDk8Xjgw4cPMylDCNm8eTOTMpRpuTx/ytj9708ZyjTqnKDOUefUZQi7OucfjCHxiEZ7qAmUqwGU\nK49GlfP/SAzyFhIAAACtgIAEAACAC0I9INF0sFJ2wjI0pbA7KKcEytmaUtgdlPMPxpB4RKM91ATK\n1QDKlUejyvl/JIZ6CwkAAAAnICABAADgglAPSNrt54VytqYUdgfllEA5W1OcE+QByeOis6ysLCZl\nKP16pc19Oisri0kZQohYzP8ylGnUOUGdo86py9jBqgyfYFIDj2h0yJRAuRpAufJoVDn/j8QgbyEB\nAADQCqEekLTbzwvlbE0p7A7KKYFytqY4B112PKLRDgEC5WoA5cqjUeX8PxJDvYUEAACAExCQAAAA\ncEGoByTt9vNCOVtTCruDckqgnK0pzsEYEo9otIeaQLkaQLnyaFQ5/4/EUG8hAQAA4IQgD0hMVlwj\njTTSSAdBWgPYgpfU1FSPZWpqapiUYWvq7NmzSrqDchuUQ7k3pjSqnOaRqC4YQ+IRjfZQEyhXAyhX\nHo0q5/+RGORddgAAALQCAhIAAAAu4D0gffHFF88+++yiRYsqKyvNZrOU39TUtHz58sWLF+/bt88f\n+9pdKwDlbE0p7A7KKYFytqY4h+sxpE2bNr3zzjtPPfVUdHR0RUVFVFTUli1bCCF6vf7hhx/+3e9+\n179//zfeeGP27NmPPfaY4+X8d5i6AsqVB8qVB8oVhn/ZPdUW4I4///nPc+fOnTlzJiFkxIgR9957\nb0dHR2RkZFlZWWFhYXFxMSFkwIABJSUls2bNCgsLU1svAAAA3+G6y27gwIEdHR1iurOzs2fPnuHh\n4YSQmpqa8ePHi/nZ2dlmszk4mqsAABDKcN1CeuGFF37/+9+fOXNGEISTJ0+uXLkyLCzMZDJ1dXUl\nJSWJZXQ6XWRkZHt7u1MLaWlpCuplCZQrD5QrD5QDORwFpO7ubqvVKqYFQSCEfP/991euXCGE9O7d\nu7Ozs7W1lRAiDnrFx8dLFwqCIF0oh/PeUgAAAHI4Ckh79+5duHChmK6trQ0LCyspKVm2bNmDDz5I\nCHniiScmTZo0ceLE1NRUQkhjY+PYsWPFwp2dnREREWrJBgAAwASOAlJeXl5eXp700mQyXbt2beDA\ngeLLuLi48PDwc+fOjRgxIjExUWwtEUKMRqPJZBo6dKgKigEAALCD30kNERERCQkJe/fuFV8eOHDA\nZDLdeuuthJD8/PyKigpxWdLGjRvT09O1uI0HAAAAOVyvQzp+/PiCBQuuXLnSt2/ftra2Z555RpwC\nbrFY/t//+3+HDh2KioqKjo4uLy8fPHiw2mIBAAD4BdcBCQAAQOjAb5cdAACAkAIBCQAAABdwNMuO\nIU1NTVu2bDGZTLm5ubm5uWrLcYlHnUePHv3uu++kl2PHjk1OTlZMHj3d3d3Hjx9vbm62Wq0FBQVq\ny3EJjU6t1HlTU9O+ffsMBkNkZOQDDzwwZswYtRU5h0anVur8xIkT1dXV58+fDwsLu+uuu+SzgrmC\nRiefdR72wgsvqK2BMXq9/qGHHrrrrrtuueWWVatWhYWF3XHHHWqLcgKNzg0bNuzatau7u7ulpaWl\npeWWW26R5sFzxdKlS9euXdvS0vLhhx+KewzyCY1OrdR5Xl5e375977zzzsuXL5eWliYmJt52221q\ni3ICjU6t1PnWrVsvX76ckZHR1dW1YcOG8+fPT548WW1RTqDRyWmdq3pebUB48sknX375ZTG9f//+\n0aNHd3V1qSvJKTQ6lyxZsmTJEsWlec3169dtNtv+/ftHjBihthZ30OjUSp1fuXJFSr/++utTp05V\nUYwbaHRqpc7lfPzxx8OHD1dbhWdc6eSzzoNwDEkrW69S6jSbzYcOHWpoaFBWnXeIWz3xD6VOTdR5\ndHS0lI6Pj+/q6lJRjBsodWqizuV0dHTINzDjFjc6OazzYBtD8mrrVRWh17lv376Wlpb6+vqEhITy\n8nIsAVYAbdW5xWJ59913p0+frrYQD7jXqZU6r6ur27Zt248//nju3Lk1a9aoLcclNDo5rPNgayHZ\nqLdeVRdKnSUlJSdOnNi6deuxY8dSU1Pnzp2rqMqQRHN1vmDBgtjYWJ6H7kTc6NRQncfExNxxxx1x\ncXEXLlyoq6tTW45LPOrktM7V7jNkzPXr11NTU48cOSLljB49eu/evSpKcooPOuvr61NTUzs6OgKv\nzkf4H0MSodfJf50vWLBgxowZ165dU1uIB+h18l/nIl9//XVqaurFixfVFuIBGp381HmwtZAEQdDE\n1qs+6Lx+/TohpGfPYOtl5RnO63zx4sXffvtteXl5ZGSk2lrc4ZVOzutcQtxX8+zZs2oL8QCNTn7q\nPNgCEtHO1quudP75z3+urKwUy0jTHC5fvvynP/3p9ttv53P6QHd3t8ViEcerLRaLxWJRW5FzXOnU\nYp0///zzJ0+e3Lx5c0REBM917kqnFutc0mm1WletWhUbG5uZmamuJKe40sl/nQfhXnZa2XrVlc7F\nixdfu3Zt/fr1hJAJEyZcuXIlIiLi6tWrY8aMWbNmTUJCgtrCnfDpp58+/fTT8pz6+noevt92uNKp\nxTq3O7FUEIT6+nq1xLjBlU4t1nlubu758+cjIiI6OjqSk5NffvnlUaNGqS3KCa508l/nQRiQRNrb\n269cucJnKJLjUafFYqmvrx85ciSHz/dgBXWuPFqpc4vFcvr06aFDh4aHh6utxR00Ojms86ANSAAA\nALRFEI4hAQAA0CIISAAAALgAAQkAAAAXICABAADgAgQkAAAAXICABAAAgAsQkAAICLW1tXq93mOZ\npqYmZfQAwD8ISAD4xV/+8pf33nvPLvPy5ctz5szR6Tz8fbW1tRUXF3d3dwdMHQBaAgEJAL84fPjw\ngQMH7DIrKiruuOMOcV9LN+Tm5oaFhf35z38OmDoAtIT627sCoF3OnDlz8eLFjo4OcavKgQMHpqSk\nWK3W999//6WXXpKXbGpqOn36NCGkT58+EydODAsLE/Pz8/O3bt3661//WnnxAPAGAhIAvrNjx47j\nx48TQp555hlCyEMPPVRSUrJ///6Ojo67775bKrZmzZrNmzePGTOmd+/e33zzzbhx41599VXxrezs\n7DVr1hgMBj73pAdASRCQAPCdRYsWGY3Gf/3rX+Xl5VLm119/PWjQIGnDSovFsmnTpjVr1kybNk3M\nMRqNUuFhw4bpdLr6+noEJAAwhgQAY1pbWx3PWmxra5PScXFxUlqn0/Xp0+ebb75RSBwAHIOABABj\nrl+/Lp9fJwhCSUnJH//4x/Hjx8+fP3/nzp120+rCwsKuXr2quEwAuAMBCQDG9OrVy2QyyXOKi4sP\nHjz43HPPCYKwdOnSuXPnyt/t7Ozk/+AuABQAAQkAv3BcbDRq1ChxpoOchISE++6779VXX33xxRer\nqqqsVquYb7FYTCbTLbfcooRWAPgGAQkAv7jtttvq6uoaGho6OjrMZjMh5M477zSZTGfOnBELnDlz\nZvv27WKbqbu7++TJkzExMdK07/379/fs2XPixIlq6QeAHxCQAPCL/Pz89PT0GTNmpKenv/LKK4SQ\nW2+99bbbbvv888/FAt3d3evWrbvjjjvuuOOOUaNGffHFF6+99pp0eXV19b333svPGdIAqAiOMAeA\nPZ999tmrr75aXV0t5ZjN5tra2ptuukk+XNTe3j5x4sQPP/wwLS1NDZkA8AVaSACwJy8vLzY2dvv2\n7VJOeHj4uHHj7CYvVFZW3nPPPYhGAIighQQAAIAL0EICAADABQhIAAAAuAABCQAAABcgIAEAAOAC\nBCQAAABc8P8BbS577pHd6iYAAAAASUVORK5CYII=\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "load('Case1-2.mat')\n", "% Adjust data to, on average, center around 0g's, then convert to m/s^2 \n", "Time_chan_2 = Time_chan_2*(9.81);\n", "mean=sum(Time_chan_2)/length(Time_chan_2);\n", "Accel = Time_chan_2-mean;\n", "%plot acceleration data\n", "figure\n", "plot(Time_domain, Accel)\n", "grid on\n", "grid minor\n", "title('Case 2: acceleration vs time')\n", "xlabel('t(s)')\n", "ylabel('acceleration (m/s^2)')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Data Analysis\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Enter values from plots for calculations. Use data cursor to\n", "% obtain x and t values.\n", "x1 = 81.45;\n", "t1 = 0.007031;\n", "x2 = 15.44;\n", "t2 = 3.181;\n", "n = 106;\n", "time = t2-t1;" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Log Decrement Method\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "d =\n", "\n", " 0.0157\n", "\n" ] } ], "source": [ "d = (1/n)*log(x1/x2) % delta" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "z =\n", "\n", " 0.0025\n", "\n" ] } ], "source": [ "z= d/sqrt(4*pi^2+d^2) % zeta" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Td =\n", "\n", " 0.0299\n", "\n" ] } ], "source": [ "Td = (time/n)%damped time period" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "wd =\n", "\n", " 209.8375\n", "\n" ] } ], "source": [ "wd = 2*pi/Td % damped natural frequency" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "wne =\n", "\n", " 209.8381\n", "\n" ] } ], "source": [ "wne = wd/sqrt(1-z^2) % natural frequency in rad/sec" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "k2 =\n", "\n", " 2.1187e+04\n", "\n" ] } ], "source": [ "k2=(wne^2)*m % one way to use experimental zeta and nat freq to calculate C" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": true }, "outputs": [], "source": [ "c_cr = 2*sqrt(k2*m); % critical " ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "c =\n", "\n", " 0.5042\n", "\n" ] } ], "source": [ "c=c_cr*z % damping constant" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Using Vibration toolbox for comparison\n" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The natural frequency is 210 rad/s.\n", "The damping ratio is 0.0025.\n", "The damped natural frequency is 210.\n", "A= 0.00189\n", "phi= 2.41\n", "Current plot held\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHCCAIAAABpA/R+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA\nB3RJTUUH4QoZBQIhDvnSBAAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJ\nbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAyNS1PY3QtMjAxNyAwMTowMjozMvVCl1wAACAA\nSURBVHic7Z17eBRF1v+LJENIIMglEIkgCWLCJQGTGJEA4mKQLChKiCsgors/0eWyhovgviKwL7Kr\nCwuiCAgRXkHkftVlI0SyrJgRBAxCAkwQBjWE1SBZAiQkk8z8/ijpbfs21T3Vl+o5n8fHp9Jzpupb\nNUOfqTp1upr4fD4EAAAAANYmxGwBAAAAAOAfcFcAAAAAA4C7AgAAABgA3BUAAADAAOCuAAAAAAYA\ndwUAAAAwALgrAAAAgAHAXQEAAAAMAO4KAAAAYABwVwAAAAADgLsCAAAAGADcFQAAAMAA4K4AAAAA\nBgB3BQAAADBAmNkC7EBtbe3169fbtm0bEgLuHyGEvF5vTU1N06ZNmzZtapaAmzdvIoQiIyNNESBJ\nUVFReXl5eHj4448/HoiNNWFXOcAMPiAwHnzwwTFjxjzzzDMZGRmFhYVmyzGZZcuWJScnc9+u2NjY\n6dOnGyngzTffTE5ODg0NxQKaN28+ffr02tpaIzXIMWLECIRQdHQ0/vPAgQOTJk2aNGkSX57Axpqw\nqxxgGphdBcqnn36Kb44FBQWvvfbar371K7MVmYPX6x0+fPiePXvwn82bN29oaKioqFi0aNHf/vY3\nw2R8+umnJ0+ebN++fUJCwuHDh2/cuLFo0aKrV6/m5eUZpkGOPn36IISioqLwn6WlpcuWLUMIvfHG\nG3I21oRd5QDTgLsKFO6HfNOmTYP53+prr72GfVXv3r03bNjQo0cPhFBZWdnbb7+NDbxe7xdffFFZ\nWVlbWxsWFhYXF5eens6voays7Ouvv25oaIiIiLjjjjvS0tL4i6tnzpw5efJkQ0PDnXfe2a9fPzkZ\n/fr1mz17Nr51njlzplevXh6PZ+3atdhdVVZWfv3111hku3btJGs4ceLEjz/+GB4ePmDAgGPHjpWV\nld1+++2CXyGnTp0qLS1taGi4/fbbBw4cKFgEluvIr3/967S0tPDwcITQ8ePHy8rKsP2BAweaNWvW\nunXrtLQ0vo3f5vhSjx8/7nK5Wrdu/fDDD0v26/jx45cvX8bG3CdSWFiIEIqLi+vatavfj4Crh0Q5\nvzm8Tti7d+9u3brhl06fPi0eVUT8KQNBitnTOzuwbNmy3//+9xkZGW6322wtphEdHY0QCg0NLS8v\nlzQYN26c4LuXkpLy73//G7+am5sreHXMmDH4pStXrgwbNoz/Uu/evb/77jsSVZmZmVhVY2Ojz+fb\nuXMnrmHnzp1yb+EWtX73u99xLWZmZtbV1fl8vsuXLw8ZMoQv5s477/zss8+4tyt0hL9cJugRbsIn\nWlJTbo4znjBhAmcwZMgQyX7Nnz8fDwU35txoFBQUKCvnQ6gcm0VHR48ZM4azXLNmzdSpU7k/H3vs\nMa7aQD5lIEgAd+VrbGw8cuTIzp07t23bJniprKxszpw5M2bM2LdvH3fx/PnzW7du3bp1K+ecbty4\ncebMmddee+2ZZ54xSrW1KC0txbeYQYMGydnMmDFj6tSpGzZs2LJly5w5cxwOB0Lo+eef9/l8586d\nw2+fPn36li1b1q1bN27cuAkTJuA34rtY8+bNFy5cuG7dujvvvBMh1LNnT7+qKioqcCv4fupT464Q\nQomJibm5uQ888AD+809/+pPP5+Ocx/PPPz9v3rz27dsjhFq3bo19gHJH+Df0N954Iy0tDRuPHDly\n1KhR8+fP94lu+srNcVLbt2//7LPP3n333Xz3Ix4NvBLw5ptv8vXExcX5Vc6HUDnne+Li4p599tnW\nrVsjhMLDwyMiIsaOHduzZ0/8qtPpDPBTBoIHcFe+WbNmpaamPvHEE4J/G2fOnElOTl62bNmmTZv6\n9+///vvv4+tut3vnzp07d+4U//Tr0aNHTU2NQbqtRH5+Pr77jBo1SsGssbHx5MmT+fn5e/bsuf/+\n+xFCsbGxPp/v0KFD+O0rVqyoqKjgv+X06dP4pb/85S/4yvbt2/EV5Y0t//73v/Gmj/bt21+4cAFf\nPHr06Lhx48aNG3f06FG5N3I+gPt8u3TpghDq2LEjJ2bcuHH4pU2bNuErb7zxhnJHfKIbOg7/IISu\nXbsmaeO3OU5qaWmpz+dzOp34z6VLl0p27dFHH0UIpaWl+Xy+q1evYl+OnY2ycgF+lft47urcuXM+\nn4+bVGFte/fuxX+uX7/eF8CnDAQV4K589fX1Pp/vwIEDAnc1fvz4119/HZcPHDjQu3fvhoYG8Xtv\n3ryJy1999VXv3r3xolOwgQNCCKFhw4bJ2axYsUIc24uKivL5fI2NjYmJidzF8PDwIUOG4CUv7rYl\nZt26dXJtlZaWxsbGcj5GVV/wbbd169bclSeeeAK3uG3bNlzYsGEDfqm2thZfGTt2rHJHfOrdFdd3\nueawMR5Dn8937do1/OqkSZMku8ZthDl9+vR7772HeGuDysoFkLsrbhhxc+jWdOqzzz7jf4jaPmUg\n2LD/Vguv13vu3DlunQQhdPHixZYtW3K3TvwbU0xRUdFTTz2FywMGDKirq3M6nVyYGuPxeAYNGtSt\nWzePx/Pdd9+9+eabwZl6lZSUFBUVde3atX379v3www8xMTECg1OnTuH4yoMPPjhx4kSHw/HXv/6V\n+0UfEhLyxRdfrF692ul0njlz5vTp03v37v3888+vXLnCjeeDDz6IF4g4BH9yHDx48LHHHquqqkpO\nTt67d2+HDh009KihoUF8kdtWw71648YNXAgLC1PuiIYUNK7vcs1hBPsyFMjKyurYsWN5efkHH3yA\np2JDhw7FHxZd5RzciHEkJCSIzTR8ykAwYra/NIIlS5YUFxfjcnl5OTdn4iOYXdXU1CQkJPC3TqSm\npv7973+XrL+mpiY41wD5vPLKK/gb9cADD3C7LS5fvjxnzhwf7+fz7t27fT5fbW1tXFwc4s2u+L/T\nuSWvwsLC8vJyXP7DH/7Ab+7LL7+UnMju3LkT374HDRp0+fLlulvgVz///PPMzMzMzMzPP/9criPc\nChueB1y5cgVHjJKTkysqKvBLXIjurbfewldWrlyp3BGfaP6xYsUK/KrL5RK0jm38Nieo0O/syufz\nvfrqqwghPPVECOXn5+PrysoF+FXu4221wH9ys6vLly/7RLMrtZ8yEJwEhbvy3fJYcr7KJ3JXN27c\nSEhIuHHjBnelT58++FYLSFJbW8vtSggNDY2Li4uOjsY/rn0+38mTJ/FLiYmJEyZMSElJwZNa7K7q\n6uocDsejjz76/PPPT5gwAcecwsPDr1696vP5uG1vQ4YMmTBhwrhx47AB54T44K2AYvC9mHyrRWho\naPPmzR977DHu1z32EJyY+++/f8SIEbiDXbp0wQmzyh0R3NC5pbnw8PDo6OipU6eKbZSb0+CuLly4\nwI1Jx44duevKygWQKFflrtR+ykBwEizuyufz/eUvf5k1a5bcqwJ3VV9fn5CQcOTIEe5K7969+fsD\nATEej2fevHm33347d0N0OBy/+93v8KsLFy7kloZyc3P5cRePx9O9e3e+d+nYsSPfo8yZM+e2227j\nXo2IiHjiiSckf3fTclfR0dFLly7FE7XQ0NBXXnkFv9rY2PjKK69ERERwNf/617/m9iYod0T83Iep\nU6dyi9LPPfec2Ea5OQ3uyufzDRo0CJvhvY4kysX4Va7WXfnUfMpAcNLE5/NJ/vO2GW63+8MPP4yJ\niUlPT7/nnnvEBv/6178mTZpUUlLCXfnVr341derU4cOHI4QqKyv79+//ySefxMfHGyeaWSorK0tK\nSu64446uXbvyg3k1NTVHjhxJS0tr0aKF+F1er7esrOzSpUs9evQQR78QQt98883333/frVs3beEo\nQrKzs3fu3BkdHV1ZWdnQ0PDFF1+kpKQIBGOpP/30U3p6uji047cjalFujiLUlWvAmE8ZYJGgcFfY\nV+FV++XLl2dkZPA9ltfrbWxs/Oyzz3Jzc4uLi9GtzRdvv/32/v37t2zZEh4e/tprr5WWlnIL+oCN\n4bsrs7UAAPBfWN0ZePbs2YKCArfbHRkZ+dhjj6WmpspZer3e7du3Y1+FEJo4ceLy5cu7dOnSsmVL\nfCU/P3/atGm4nJSUhBAqKSlxOBwTJkxwuVzp6elRUVEtW7ZctWqVzn0CAAAAZGF1dpWenv7ggw/2\n6dPnzJkzGzdunD9/Prehiy7V1dVXr17t1KmTHpUDFuTmzZt447jkiiUAAGbBqruqrq7mpkdLly79\n6KOPCgoKzJUEAAAA6AerOa2cr0IItW/fXjKpEwAAALANrLorDo/Hs27duuzsbLOFAAAAADrC6mIg\nx4svvnj16tU1a9aIH/eCEOI/Bg0AAACQw+VymS3BH6ZmfQXK9OnTn3zySf6zJwQkJCT4raSoqIik\nLRIzilX5DBdPaAbi9WjR9uIJzUC8KjPjxZsLw7OrmTNnlpWVrVu3jh/HEpCYmMjATwYZQLxZgHiz\nAPFmwYR4VvOuXn311ZMnT65bty4iIsLj8SD5B6sDAAAANoBVd7V161aEUP/+/fGfDoeD//wke/DJ\nJ5+YLUE7IN4sQLxZMC2eCVjdGej6JYH4Ku4M1sDNKFZFCN0WKfaRVnN0zSiKp9ui7cWTm9Gqimnx\nhGbGizcXhmNXJDCxICuH2+1m94m6IN4sQLxZMC2eiVslq7MrAAAAIKgAdwUAAAAwgP3dFX/dVrLs\ndDr92iCE8vLyqNjwmyNpV7nM/7+l+khSJhwHin0kHIevvvqKln6KfSQcB8Ky8d9nivopftZM/5ul\n2EcmgNiVdWF6KRzEmwWINwumxTNxq7T/7AoAAACwAcHirppMLzRbAgAAAKAdVtOEKeJ0OjMyMqiY\nUayKELotUuwjCUyLV64Nnq0MWIf/+7//o/i1N5FgiV01mV7oWzTIbDnqYHopPMjFMxEJAIIBwq8i\nE9/YYFkMBAAAAJgG3BUAAADAAOCu2H6EF90WKfaRVnN0zegmmjCXtgIEJ7b5otrfXSmk0eHtghkZ\nGSSpdiR1ErbLhT0DTznEVfm1N76PJGXCcaDYR8JxIEwT9qvfBly/fr2hocHIFm/evHnz5k1Vb9Eg\n0vh+WQSmv6tBvdXC4vsvgny3gonYb6vF7Nmzi4uLuT/79u07a9YsScvNmzc/8MADHTp0wH+2a9du\n48aNmZmZmpsWVOiXyZMnI4Teeecd/kWsPyQkpGnTpnfeeefo0aPT09O5VzWIDLxfOqF2uPxip60W\nsJEdAOzPoUOH0tLSHnjgAfxnTEyMnOXkyZM3btzI3S43btyYkpISSNOCCrVx6NCh+++/f+DAgQ0N\nDV9//fXQoUOnTZv2P//zP7REWgcqw2VX7L8Y6Bc8I/abRwyxK1VV0WqOrlkwx65SUlKG3iItLQ1f\n/Ne//vX000/n5OS88MIL33///apVq65du7Z48eLx48d/8cUXCKF//OMflZWVCKEXX3yxqKhoypQp\nOTk5f//73xsaGl577bXs7Oy//vWvXBP79u175plnsrOzp0yZcunSJYSQuEKE0JIlS37zm9889dRT\nRUVF+EpVVdVLL72Uk5OzePFiOf1JSUmZmZlZWVkvv/zy/v37X3nllbKyMvwSJ1LQHU72xIkTc3Jy\nPvjgA3G1Ys0cS5YsGT169OjRo99991055YQjI/fGL774Ar9x9erVcsMVOGx9URWA2RUiTKAjMaNY\nFSF0W6TYR1rN0TWjmyxJWJtZT1QRLHRfv379p59+wuXw8PAWLVqUlZX95je/ycvLi4yMrKqq+umn\nn+67776mTZv269cvJSXlzjvvRAh98MEHQ4cO7dat2/vvv//555/Pmzfv2rVro0aNevTRRx9++OG0\ntLSXXnopNjb26aefRgj98MMPTzzxRLNmzQ4fPjxw4MCysjJxhaNHj/Z4PL///e+rq6tzcnI++uij\n9PT0nJycO++88/nnny8pKVm3bt24ceOUu9arV6+UlJQvvvgiISGBExkSEiLoTqdOnd5///19+/Yt\nWLAgJCQkNzcXIYSlcog14+ujR4++du3a5MmTQ0JC8DqqpHKEkN+RUXjjZ599Nm/ePITQCy+80KlT\nJ/FwUcEeOcIIIeSzNQkJCbiApu0XFBTKFuH8+fNmS9BOkIvnvngWITMzMzk5OfMWb731ls/nKygo\n6N2797Vr1/iW0dHRBQUF4j+joqJ2796NL44YMeLll1/G5aVLl44bN46zb2xsvHbt2rVr17p3737o\n0CFBhadPn46IiKirq8N/rlixYsyYMcXFxREREY2NjfjikCFDJk2aJNa/adMm/pURI0ZwZrgJye5E\nRUWtX78elzds2JCSkiLupljz6dOnw8PD+VVJKueaUBgZwjdOnTp1woQJ4vEPHMKvotW+sZIE4+zK\n4jssAEAPZs2a9eSTT/KvZGZmpqenR0dH9+vXb9iwYS+++GJYmNINoXXr1rgQHh7es2dPXL7tttvq\n6+u5JtavX5+SkhISEvLTTz9xkzmOM2fONDQ0cNtY6uvrBwwYcOHChQceeCAk5OfARKdOnUi6c/Xq\nVS4Up9wdbpqSnp5+5swZQT2Sms+cOdOnT58WLVooKycZGcI39uzZ89ChQyQdD2YgduVnYZdbyYHY\nlaqqaDVH1yyYY1eS5OXlXb9+/eWXX968efP//u//BlLVsWPH1q9ff/bs2V27du3YseO2224T2zRt\n2vSOO+64eIvKysodO3Y0bdr0u+++42w456fADz/8cPDgwT59+pB0p7a2FhcuXLjQvHlzEs3NmjU7\nf/68X+V+dQbyRorY4IuKAXcFsSt1ZiBep9qM54cffqivrw8LC3v44YcfeeQR7DM6dOhQU1OjoTb8\nLjxJ2rdv39mzZ/F1foUPP/xwfX09/379zTffZGZmVlRU4J0FlZWV//jHP5QbOnz48OOPPz506NC+\nffv67Q5CaOPGjbiwefPmJ554gkRzZmZmQ0PD1q1b8Z9VVVWSykmGRe0bNY+/Aqx/UTns764CT8WF\ncrCVCdOE/ZYtxahRo5rcYvDgwQihM2fOtGnTZujQob/61a8++OADnIn18ssvjx8/vmXLltxdnpAB\nAwb07t27R48egwcPXr58+aBBP6+38ysMCwv76KOPXnnlldTU1EceeaRVq1aFhYVNmzZdt27dr3/9\n68cff3zIkCHiORPm6aefbtasWWRk5JQpU0aOHCmeo0h2ByHk9XrvvffePn36HDlyRDCDlNMcFha2\nY8eOmTNn9uvXb+DAgQ899JCkcpJhUftGzeOvAHPfVVnMDp7pi+RWC7/bLviYuAUjyHcrmIj9tloo\ncOXKFcH2hABru3Hjhl+za9euXblyhX/F4/H8+9//piKA352oqKjPPvusrq7uxx9/VHiLnObLly8L\nXhIrJ0TzGwPETlst7D+78gvhTwyS7cgQu1IF0+Kp12YWrVu35u8pCLy2yMhIv2YtWrTgdhlgwsLC\nFJKXVQkQd6dp06bt2rVTeIuc5rZt2wpeEisnRPMbA8ceX1QEeVdIamFX89ZBiF2pgmnx1GsD9KC6\nutpsCeZjmy8qzK4AAAAABgB3FShmPbMAAAAgqAB3Zd1gjPEtQuxKFbYJCVDE3IM51LYeJCeP2OaL\nCu7KusEY41uE2JUqGAoJzJ49+xEef/7zn3VqKD4+/sCBA3rUvHnzZsEjaAlbx30fPnx4Tk7OtGnT\njhw5omyvoQmLIDdEDH1RlQmKrRbU1+vgMU4AW5AfIBIg+p3loflkDYXDR+x08ggKgsNH7D+7opU6\namRbUDa3bMs0YfEBIufPn3/xxRe5o3tnz56Nj7eQO3RD+fgM/KRzVQeOSNYpPlYDiU7WUDj1QxK5\nw0cUTh5RGAcOORnik0fkeqrfEPFrYO67KovZiV/6gnPfBKnBgjThoqIiH8GT2kke5Y6r8guhGUmy\nKt0WScxAPEltVku6zMzMfO+99y7fgkukffbZZ5977jmfz7do0aL+/fvji1FRUYmJibt37/7444+7\ndOmybt06n883atSokSNHFhQUbN++/fbbb//yyy854549e27ZsmXPnj2+Xz7BPSUl5eOPP96wYUPz\n5s1HjRq1Zs2aPXv2dO/eHVcoV2dUVFTv3r137969e/fu22+/fe/evT6fr7i4OCoqav78+Xv27Ckv\nL1+3bt3HH39cUFAwf/78u+++m+um5OPMxU9zT0lJef/99zl7l8vVvn373bt3FxQUbNmypbi4WGEc\n+E1Iyhg1atSwYcPy8/P37t37xhtvKPdUvyHiOpuQkEDytbfaN1aSoFgMVIZ8YdfvGmC/7Td9BJVB\n+EdVVdYUT15bWU4sxUbJSdhWwf/zrbfe2rRpEy4/+uijL774IkJo2bJlqampM2fOXLt27fHjxznj\n2bNnDx8+HCF07dq1hQsXpqen7969+z//+U/Tpk0RQj/++OOSJUs+/PBDbDxv3rzs7GyxgD/96U+P\nPPIIQmjr1q2dO3f+7W9/ixA6f/78p59++vTTT585c0auznnz5uHWDxw4sGvXrocffviee+4JDw/v\n06cPPq7+6aef9nq9NTU1999//4cffnj48GG5RzdJEhcXd+TIkWeeeQb/+d1333Xo0GHQoEHi/GLB\nOAjOyhLLuO2223bu3Hn58mVc1cMPP4wQUuipfkPEB2JXQQGVoBcEuoIcgdswC/EBIgihyMjIt99+\ne8iQIStXruTHPASHbiifgiH3wAjlA0cU6vR7rIbfk0qUERw+onCQitrDRy5duiQ4eYSwp9SHyJaA\nu9IL8FKA9fF6vXPnzh0yZMg777zz7LPP4l/xSHToBj4Fw+12U2xac53cqR9YLT5QmBx8+Mj8+fP5\nF/Py8lasWFFYWDh79uyqqqrXXnsNXyc5fIQvQ3zyCAqgp3oMO9NYcauF1+s9evTorl27tm/fLmdz\n9OjRbTwuXLiguTnj442QuqSqKmuKp16bKcyePbtt27affPJJcnLySy+9xF0XHLqh+fgMBTSfrCF3\n6gcJkoePyJ08gtQfPiI+eURDTzloHT5igy8qxoruas6cOS+88MKGDRvmzp0rZ7Nr166VK1d+eYsr\nV65oa6vJ9ELjF3Yh/KOqKmuKp16b3ogPEPn000/XrVu3du1ahNDKlSs/+eSTjz76CBsLDt3QfHyG\nAppP1igvL5c89UMB5cNH5E4eQeoPHxGfPKKhp4EPkeDwEba+qEqYvddDgvr6ep/Pd+DAgZ49e8rZ\nzJo1a9asWX6r8rszUHnjnypjhUr86pQkyM/gMJGgOkBEjMKhG3qcgqGtTsKTSlRVKDhIJZDDR8Qn\nj/hMOnzETgeIWDF25XA4SMzq6uoOHjzYpk0bLkQJAAAtJA/doHjUSIB1Uj+MQ65Cv4ePSF5v27at\n+KLm0dNj2FnEiu6KkIKCgosXL5aUlMTExKxatYrbPyMgMTERPbICIcSPWPLLO3bswJntcgaSZc3G\nCKHi4uKUlJQu77jPT5bWjCkvL1d4lV+VkWaEVdlePHltzAGHbmBsMw7cLU5MVlaWwWICgVV3lZub\ni/f2eDyeKVOmTJ48ec+ePZKWLpcL70e/5c/cgvKtlBE3z+G55Yy5spSxbCXcRbxdUO66WLycDyY3\noG5GWBWJJdPiVdUGACYimRWHcblcuJCYmGiUHO2w6q646bnD4Zg4cWJ2dnZtbW1ERIS5qjjgVBGA\niX//AMAQrLorPjiljsvsAwDT4X60Go/b7WZ32gfiAQWsuJHd6/V6PB58qIzH4/F4PPj65s2b8aZb\nxMskqKqqWrp0aXJyMuEGDTGQd6XKDMTr0aLtxZOb0aqKafGEZsaLNxcrzkjy8/OnTZuGy0lJSQih\nkpISh8Nx7NixGzdu4Cd9zZgx4+rVqxEREdevX09NTV22bJnm5iDvSpUZiNejRduLJzejVRXT4gnN\njBdvLlZ0V8OGDRs2bJj4+oIFC7gy/wgDu9LlHbdvEawtAAAAIGTNxUAAAAAAEADuio3YldxWQ6bD\nP0yLp9ui7cWTm9GqimnxhGbBFruyv7tyOp2ub8dyZf51XMjIyND7hE2+s3E6ndxSsrhdbEl+ai2u\nyq+9Qh/l+huIDWFZYRz4ZYp9JBwHwtOESfRT7CPhOBCW+QT4WdP9HEnKFD9run0kKRv/fSbpIxuY\n/RQofcEPwnKN7GDwMwNpvVHSmAmC/JmBJgLizYJp8Uw8M9D+syt7A/nIAAAECeCu2IhdGdMixK5U\nAbErPcxoVcW0eEIziF0FHUznXdFtEfKuVAF5V3qY0aqKafGEZsGWdxV07orbdgEAAAAwRNC5Kxsg\nGa+CIBYAAPYG3JV1I0nGtwixK1VA7EoPM1pVMS2e0AxiV0GHZSNJVMCzLmtGUCB2pUdVhNg+/MO0\neEIziF0BAAAAgOWwv7uilfkvQBwrIowe0XoqgUKZU2JAW7YsEz7VAspQtk2ZDczOU9YX8VMtXCM7\nuEZ2wGX8CImioiKf+odTiJ89QW4saFHzIzAElSgY4Bb9QmJGWBVJhj9FVYRmFMXTbdH24gnNQLwq\nM4pVwVMt2MDesSt+i37nfxC7UgXErvQwo1UV0+IJzSB2ZSs+Tr6m9i3s7gjXppzd/gIAEFTY3F1x\n2DI7OBBPA14KAAC2CBZ3pYBls6CMbxHyrlQBeVd6mNGqimnxhGaQdxV0WDl2RWsOZM0ICsSu9KiK\nENuHf5gWT2gGsSsAAAAAsBxB7a5sGdCiBQS3AACwFPZ3V35T5JxOp5Fru/zmDEsDFPRRzhWRXM/L\ny6OljXAc+P9XsBf0UbMNIk4TJtFPsY+E40BYJvkcCT9rup8jSZniZ023jyRl47/PJH1kA7MTv/QF\nZwTzU4MFZc3ZvoGkCRvwRrGxwhuVjQVlQpg+CBzEmwWINwtIEwYAAAAAOoC7EgIBLQAAAAsC7sq6\nWVDWbJGrKvC9GISqKJrRHXmKLdpePLkZraqYFk9oZrx4cwF39XNSgpGTKitnemmrSpvrgrwrPaoi\nxPapS0yLJzSDvCsg6IA96wAAWB9wV4AKwLEBAGAW4K6UFnZ1WiFkOnalymMpgjHZVQAAIABJREFU\nG0PsSo+qCLF9+Idp8YRmELuyG37T6DIyMoz8OJ1OJ7eUTCvl06w+BqiNcBywmdo+arZBxGnCJPop\n9pFwHAjLfAKxof45kpQpftZ0+0hSNv77TNJHNjA78UtfyNOE5S4GQ5ow3fxiDE6Z1JBfbAWYzvcE\n8WbBtHhIEwaCEckFQAh6AQAQIOCujJ4RN5leyHTsiiIQu9KjKkJsH/5hWjyhGcSu2MDr9R49enTX\nrl3bt28PsCpuVdqwGQDTeVcUycjIIBlzSF1SVRUhtk9dYlo8oRnkXbHBnDlzXnjhhQ0bNsydO5fw\nLUG1HmXlzlpZGwAAloVVdzV37txjx45NmjSJSm1lObFU6rEymp0EeBcAAKwAq+7K4XDQqioYIknW\nXLwWqJLzixD+UVUVIbYP/zAtntAMYldBB1uRJF2fzqcrYuUkqppML4Twj6qqCLF9+Idp8YRmwRa7\nCjNbgO4kJiZ+nPxzmVv0c7vdnIFkmSu4vh3rdhcRGqst62dsTCsKBuXl5V3ecWtupcs77vOT48UG\nxlBeXm5wixQB8WbBovisrCyzJajA/u7K5XKJQ1Px8fEIobJb5TLe9TIpA4QQQm65sl8DKm80thWi\nN0oZU3njLyrhXTQOUxqlBYg3C+bEu1wuXEhMTDRXCQmwGEi0sFuWEyv5/EBtezRsE7sKcBcG5F3p\nURUhtg//MC2e0AxiV2zg9Xo9Hk9DQwNCyOPxeDwezVWxFbtipUUSAll51xYJYzr8w7R4cjNaVTEt\nntAs2GJXrLqr/Pz8pKSkiRMnejyepKSkpKSkQDwWYCnkJm1qrwMAYCdYjV0NGzZs2LBhOlXu+nYs\nQhU6VY6BOywAAIAqWJ1dUUTDwq5kHIv8cCzbxK4CBGJXelRFiO3DP0yLJzSD2FXQEQyRJGsuXusa\nXZCcvzId/mFaPLkZraqYFk9oBrErALAbsPQKADbA/u6K5FRNs6BysiphWW/004+dDYk9iVsiqYfw\nNGEoQ9k2ZSawv7uSmwhLnmZtDFyLfG26lg3oowZtdMeB79vkcPJOMUcI9dt+U7LO1NRUWvoJ+4jN\nArchL/MJxIZuHwnLgs9Rsw31PpKUCcfB4D4ygf3dlV+CIZJkze+lZYMxxrdoe/HkZrSqYlo8oRnE\nrgA/BMNpI8EGBLcAwPqAuwJ+ge1v3LbvIADYlWB3V65vx7KSBRXIfdaaMVXLJhJx+I2E0WoR8q5U\nYXvxhGaQdxV0UF/Y9btaCLEr7Ab4qgKf9NBa6ydXArErPcxoVcW0eEIziF3ZDc55kD91Qm8lAEVg\ncQ8AggT7uysAIEHZ7cHTdQHAdILXXXFzHbmFXf0mQ06nE0/1DEsDNGDxWm1Kr9PplJSnOb1R7Vq/\ngqfxmybcZHohoX7CPvL/H4gNeTkvL4+KDd0+Epb5jQZiQ72PJGXCcTC4j0wQvO6KI8CFXd+iQeIb\nn/LCI9eiYWnCBixea9BGZRwULmo7EwsRpwmT6CfsI///gdiQl8ePH0/Fhm4fCcv8RgOxod5HkjLh\nOBjcRyYAd6U7psfMACrAuh8AmAu4K//AfYppDPv4+A1xZfjyAAAtwF3pu4ArGQAzfslYbYvG3GRV\nqfIriWKeCqESayYSGd+i8dk/thdPaGa8eHMBdxUUWVDWXKSmq4pingrGr4MMRL+gcuriqdRGN62H\non7biyc0M168uYC7AgylyfRC36JBZquwFrBgCAAkgLsyGcgdZh1yZwNuCQACAdxVoAu7GvyN2hYD\nd2nWXLy2bDAmwBYV3JLcS+KqxJaWHS5rhn+YFk9oBrGroMPESJJCip9kFrPmsgF91KCNZBz8tsVB\nca2f8DRhQv0kqdO4KuX+im0CKYvr12ZDPg4U9WfwzksMxIZ6H0nKhONgcB+ZANyVmdBK+SQp640B\n+rmMbL37RZgmrEdZcvold+oxlKHMyv2BCuCuLAT1OJZ+wRLxjgkIzASO2jGEMQeCCnBXdGbEqjyN\n9fOujMGywZgu77hptUgIuXhaLsr24R+mxROaQewq6GAoC0rz9Muas37LJhJRbJFiVQo2GnyY7VOX\nmBZPaAZ5VwBgW9hdPSNXzp8asttfABAD7sp8KD4DVzz9knxgPGAF4HMBAFWAu7JiJIlk0c/i0TIS\njI9d0cVSwRi1zs/24R+mxROaQewq6GAodiWHX9dlzcVry0aSNLSoeaqksDVfw3MF/TbElXFtyrKZ\nDv8wLZ7QDGJXgE0gmX4x+gQ/WEZDMAhA8AHuykw0P8GBsE4OBdelx0IBrScvGL+IIfbfevRF734F\nro3kSRxQtlOZCcBdqfjMAtkTwXcYXJkkw9zvRF7gitRO/PVYKCDMqOc/S8avvV90+rdHS79yvygG\nKgT1yz0RQ2wvtywpvq5tHPiPICKxVyg7eU9FCsRGgN6ftapxMLiPTADuir3PTBnyLRiSHtRI7BS7\nMqAqWs2RP8iKMNPLmuEfiF2pNbM+4K7YgFaEifqeQ4B1IAYGsEKY2QKkOXv27Pr162trawcPHjx4\n8GCxwdGjRy9cuMD9ee+998bFxRkmjwo4I8pltgyAFvgDZXHrCjmSHZTrte1HAzAYK86uXC7XyJEj\nY2Ji0tLS5s2bt3btWrHNrl27Vq5c+eUtrly5ork55uKNAaJttVDyOj+YpyGwR/d3Pf9zFNwodcqV\nNiuRiEpfNETLSNpVsIG8K+otQt6V+SxatGjMmDETJ0588skn58+f/+abbzY2NorN+vTps+AW/BMf\n1GLkwq6Ge7rAZ1B8BEYgMgJEj15A7IqDxK9YM1pGCMSuqFfFBFZ0V0VFRX379sXlAQMG1NXVSf46\nqKurO3jwYGlpqbHqSIGnH2Ekp2JA0CLIVjZRCcAclnNXtbW1DQ0NnTt3xn+GhIRERkZWV1eLLQsK\nClasWDF69OjBgwe73UQnPvgFdhkYDIlTV+XnOGNVHyW4UjF0fQl4JiBwLLfVwufzIYTat2/PXXE4\nHOLFwNzc3Pnz5yOEPB7PlClTJk+evGfPHvJW+O5tx44dSSIDuZsd/42SPpK76NeDCgyU7VUZq7JU\nFqzQX/ynpIHb7S4vL5drgv//spxY1Hm9WIPr27Fud5GgBjljXC4uLkaolaQkgVTlVuSMubLr27GJ\nndeLW5R8I2csKAsq37FjR0pKCnfRM6Mf6rxe0LrARk6eQllVbXwbydqaTC88PzmeXIZybV3ecZ+f\nHI+/NrjMLwgoLi7mxMtBYkPXjP+dN6bFwKvKysry+3brYDl35XA4EEKnTp2699578ZWbN29GREQI\nzNq1a8fZT5w4MTs7u7a2VmwmR3x8fNmtcnZ2dtmGyRreyC+LDeLj4xFSchUCA8nayI3lpinkMiQt\n+RclNcTH47vJL2rwzOjXcWGR+I24wHvjzzfx+PgK7o0J2yp4xsKaRbUhhNy4hvj4eFRUKGnM/79U\nzYj/kWFXITB2fdsPIRQfX8E3jo+Pd22ITey8Pj4+3vVtP9wLwRt5Nf+3jP0uX2p2djZuBddWdsu3\n8Y2xjWBABLX9sub/GuPvBv+iqDZ0Szzi24jfKFnmvxGXFVoUG/PLvA9LYPBfJC9qsDHFzFJVuVw/\n701OTEwkqcpc6C8GVlVVnTp1atKkSaNGjVqzZo3b7fZ6veRvdzgcsbGxFRUV+M/Kysra2tquXbsq\nvKW+vh4hFBZmOddrA/wuqem6fEprJzRfpAHrfpJNUBko/oD47Yhki5K54WU5seQftKqO6PT1gABY\ncELTXX388cd9+/a9//77R4wYcejQoXPnzi1ZsiQrK6t79+7Z2dmVlZWE9YwYMWL16tV1dXUIoXff\nfTclJQX/Oti8eTO3qZ3bfFFVVbV06dLk5GQ8LSPEsjEqCKIQ4vp2rN4fojGfhbbInH6oetaJnEf0\n6ykVasPuR2ysyi2BD7MldNxVZWVlWlran//85//3//7f559/7nK5jh07duTIkRMnTpw4cWLv3r1R\nUVH9+/d//vnnSWqbMGFCx44d09PT+/Xr53Q6Fy5ciK8fO3bsyy+/xOUZM2YkJSWlp6dnZGTU1dUt\nW7ZMs3jbJCWYhUWyQdV+jiS3fmUbayYSYQReR9wR6l97Vfl85HM1v8l/SGbmB3lXas2sDzV39d57\n7x06dOi5557jokqY8PDwuLi4tWvXnjp1qlu3biS1ORyOZcuWnThxoqioKD8/v1OnTvj6ggULOLdU\nVFRUUlJy5MiR06dPf/jhhzExMZrF2yYpgS7Ge6AApw4knyPd2YmgRVWzPUEWM/9BtCTNacjGE0Q3\nSU7YkmvO4MUJZd/G9V0gPnpxDnddsgZchrwrhqDjrnr06OF3g0poaOi0adOoNAcQYvrSInUBpveI\nEG2BpeDELOcneVFu5ud3jRQ+UAOwXN4VYC60DkkB1GKd0dOgxDoPPdGgRFWgiyS259fhAdrQxV1V\nV1e73e6zv0SPhqigOSRgG6wW7deGTgv0xjwEjwS6zdkmnuEXXQOQGmJvTqfTb+CtLCcWYldi6Lur\nF198MT09PSsr6xEew4cPp94QLbiFXbo3YiqnABsArVOM1bZF/YRT5QV6Kh+uQCcXHZHTT+vHDRbP\nP9wvcDjxAeL3s9M71UHhVbXfMb1P3eUfzIjLeHAk2+KfCRmIjaQGZRvLQtldXbhwYe/eve+9915J\nSYmLx+nTp+k2ZH3U3g7Ed1tj5jGaTzEmh9s6qO2kXetgrn5VWwMMg0rfddp3rvAZcS1mZGRw/9D0\n/kypnKSsU5kJKLurmpqa2NjYAQMGqMqCsismrptZdsnOSGAQJNFvWBRSGhj9LKywDQTgoOyuYmJi\nJB9Ha2U0z4jhuxUI1EfP6XQyek8EyMFfG4oftJH5Z+KlPwUgdiWGsrtq27bt6NGjZ86cWV1d3fhL\n6DZEEeZmxDrBuveFz5EunEuwQg64TpTlxFr2mDTIuxJDf6tFRkbG7t2709PTe/BIShI/9Bz4Bax7\nCxKs3EcbzMwoDi+5i7LyZ2pNCEcMniMlhvJjYaurq3/7298OHDhwxowZISFBl9Tl+nYsQhV6t4If\nT+CiUZVAsG/RoDI6+8UAFXCPUVcLxW+CReBOTgkQG/z4AMRQ9ijl5eXR0dGrVq26++677/oldBui\nCIsLuzr9a2RxKDiMFA93Q9tA8rXRYwZJZfIEsauAiIyMDA0NpVun3thmYVcO8nurBYeC/E5htV3U\nppCwrYJRV+r6dmzCNt1XJsRQ/M57ZvSjVRUiW4+F2FVAxMXFtW7d+uDBg3SrBfRD79TIQDSQQPfu\nTP4EBLPGTVVbuu6SMMvBkCMYK+53iRU+Rz6S7QoebKF3mQkou6vKysrr168/99xz/fr1e4jH4MGD\n6TYE0MIKKYSWSl1UmKCoSr0kmbSV5cQmqgzVGDlW7M7VMKyk0EKaMCH0d0O0aNGiW7du0dHRLXg0\nb96cekO0sNrWJhP3DVP5tWXWw8iN+amowcGwizV/fRMu3iZ2Xk/yTSPso3JtenwrLk/b5tcm2GJX\nlHcGtmvXbvfu3XTrBEjQvLuMT0ZGRtli2VfxMwsM2zqotkfK4s0isfN6Y/aLkoCX71SFAwMZUiY2\nmlrza4Mg70qKoNtrDuiBrjNUwp/JAADYGzruyu12//TTT37N9uzZQ6U5AACYQzCxU3jAoCSub8fy\nF9xUTRP91mYpgmrBWRV03FVtbW1GRkZWVlZpaanH4xG8+tNPPy1dujQtLe1vf/sbleYMhulosyrI\n17gtOCa2WaBXwOCJZjCcsGUFVZLROHhmoBg67qpHjx6nT5/OyMjIyclJSkpKT08fNGjQQw899MAD\nDyQlJWVkZGzatGnJkiX//Oc/qTQXnBjgIdhd4zbl4W8GoNMPbcJpDd0hteYHZE1VCGJXUlCLXYWE\nhMyZM+f06dN79+4dOXJkeHi4w+Fo06bNm2++6XQ6i4qKBgwYQKstAOCDb+jk2b52ygsGWOTWBhxA\nHfS3WsTFxf3xj3/Mz8//5JNPdu3aNXjw4LZt21JvBaCFKSmEghmDru26vh1LsieYBKulXgbYlt+0\nKlX1i2vzOzUMRL/a3xyspNAamRrM3CIh7AwMdshPE9bv96DeqYtmnQarNxZZ5Cka2YzETHxztEIq\nrpxN9OIc/Z7ZobAYy2lwOp1+v2MkNnL1K9hYFnBXNkRzojFzX18N6P2kBjtt66IbGrHmtwuroruH\nhUptELsSA+5KI7D0DOgBPhPErNat/xhA66DqY7LTjxgTAXdFGfIvMf/3lzbnBy7TdOjO1SAhmjqm\n+H4bn79sLuCugP/Cjy5Y+b4JfppD1/kQ3bQevQP7Cpm/Cr8DLLvdAPKuxOjirhobGysrKyt+iR4N\nAarw+6MvkDVudtc67LoCpnauJv56iL8Pkl8hPSIohk00LRvXgdiVGMqPuEUILV68eOXKlYKLISEh\np0+fpt4WoI0AAyQkz7+hdYp5IBoAALATlN3VDz/8sHLlytmzZw8bNiw8PJxu5UCQgFd1YMUPCB7w\nXgz4zitDeTHw8uXLMTExY8eObd26deQvodtQkKDT11ecKohXw/Ly8oyUoaDHIvULdlJQSbekuElM\n77GS+z5okyE2I9dPuBhAuMdBro+0xpNK+rPT6fT7fSOxQfJ9ZC6mRdldNW/evLGxkW6dbMHE7yO5\nVMHx48ebIecXGvptvxlIVZIxD1qpkXxtVki31DulWvx9kPS1GiIo+HdA9OIcU0Kecn3UI508kEr8\npgCT2CD5PjIX06LsruLi4jp16rR//3661QIkqPKUTLhVwDronV7NEPgnkWCHDmxeNwDKsavKysrK\nysqJEydGR0c3a/bfR7OEhIQUFBTQbQtQhXXOtBVD5Shk/RDEFax5SK41VamiyfRCl9kajCdhW4UL\nNg2RQX8je4sWLbp16xYdHd2CR/Pmzak3ZB2Y+9Wpdw6KhhUeu+4m52PKBCWQgaWb1mPNB+ETiid8\nLqICaqdfkHclhvLsql27drt376Zbpxxnz55dv359bW3t4MGDBw8ebEyjDKFhq3dGRkbZ4v/+CVuV\nmIPuBEXwfQjQzJpg8X7npsb3EfKuxOj1VIvKyspz58653e7q6mo96ne5XCNHjoyJiUlLS5s3b97a\ntWv1aAXg0PwL3e+Ugt38YnOx4FNZSQjkwHtjgEdhWRb6acJVVVVDhw69cuUKdyU5OXnbNjoHDnEs\nWrRozJgxEydORAjdfvvtubm5Y8eODQ0NpduKdYB/Pxy2Hwq8fGf7blKEcLj8DmxwBs8Ygv7sasCA\nAa1atSosLCwuLi4uLt60aVNpaWl2djbdVoqKivr27cu1WFdXx+L6rNWW2qw/hjAVM5JgCI2o6qOR\n4TeIXYmh7K7OnTsXFhaWn59/xx134OzglJSUw4cPl5aW1tTU0Gqltra2oaGhc+fO+M+QkJDIyEid\nVh2pYDW3JIdt1rgBKrASGgnEixjQR23yIHYlhvJi4M2bN9u0aSO42LJlS4fDcePGDVrPtvD5fAih\n9u3bc1ccDgcr6cl+XZfb7aZiwP+/whsljcUixWZ+ZYjfSGJJCGG1crK5NyobEI6YKhnkvVNlrPAW\nKr0zayhUoWoo/L4lwLESvEXVwCr3giJZWVkGtEILyu6qZcuWFy9e9Hq9ISH/nbeVlpZ6PJ5WrVrR\nasXhcCCETp06de+99+IrN2/ejIiIoFW/ucTHx5fRMOD/X+GNcsZ+zfzKEL+RxJIQwmrlZHtm9MP7\nHuUMFOrXMBSqRgzRHgrl3gm2MwQ4FNzAalYrwLdoUJPppDtUVQ2F37cQdr/LO24XwYipGlhsT9Ll\nAHG5fo7WJSYmGtBcgFBeDOzUqVPHjh179eq1bdu277///vvvv1+8eHF2dvbw4cOxj6GCw+GIjY3l\nDiWprKysra3t2rUrrfqDFkbXuC3yQAGFPZDWTDnyi8L3gT/mYjNzH4GharT54hVy1Iz/pwGxKzH0\nt1rs37//vvvumzVrVmZmZmZm5urVq5999tmFCxfSbWXEiBGrV6+uq6tDCL377rspKSnG/BixNwpr\n3IzecBmlyfRCKzz0KBhCI/r1ET8MRdvn6Fs0CGJXYuhvZEcIrVmzBiFUVVUVGhrasmVLPZqYMGGC\ny+VKT0+Piopq2bLlqlWr9GgFCAR+6qXB27JJ9oL/vMpEQxi+MfkCrwgwBLxh3QaPrQo2dHFXmNat\nW+tXucPhWLZsmX71A+SYPgkADCaxMx03H+RAmpda6CwGnjp1atCgQRcuXKisrBwkxUMPPUSloSCH\nv7CukIGkeSpjmzVuq0H36QyGPV/R+qGRwIfCsn2E2JUYarGr6OhoriCmXbt2tBoC9EPVGrdfpyg4\nGwmepsMceoRGKJ5OqYzcBhzBxsUA+6hfiBFiV2LoLAb26NFjy5YtuMwVAKvB3Slg+S7Y0Gn5Dpaz\nlFHes4rn3BDyJIfyzsDKysqpU6cKLlZXVz/zzDN0GwpO9DvsHG/802MuYtZCBJWD6gnr1wO19evR\nLysclG61cUY8J0QyVoHoV1t/IGUmoL+R/csvvyS8CEii4DP0PuxcDwxbiBCEMUgOBQ8ES41bYuf1\n0YtV73KTTE4Qt4sn5Qp6dL3raR5nuqokayP5jgWi32+dJDaEmplArwNE+NTV1fEfcgEAgK7Q2ovB\nYmiEC01hVX6HQixecgWPsDaKQOxKDLWN7DNmzDh27BhC6PLly4MG/eLzvnTpUu/evWk1BIjxe46i\nZBKSKadUsBs2g0PKdcLcgYXwG0NQm/TcdtttrVq1wg8GbMWjTZs2r7/++qZNm2g1BNDFUnvnDD6/\n2ICnN1nkAVGAZug+z0W5Nis8ysTKUJtdvfrqqwihysrKuXPnLl++nFa1QIBgb+R0OqPNVhKcyI08\nE3M1wq+N2m+Xb9EgA+Y0dL/zxv8LEsSlNNuQm1kfyiGldu3aga8yEsLfYvb4sloKwvwhpkeeUHz0\n4hyFoVA7YzY48KaMIBLmF4pTMYhdiaEzuzp16tTkyZPXrFnTvHnzJ598UmzQpEmT/fv3U2krSFD4\nR26p5TvlAJjgd3Qg0TLJ0YDAgw2w04cYSCoVrBv7hdpiIP+pFuJXYWcgQ5iyBYOEhG0ViJEHw5u+\n1hfIz3x4JKBOmP6tYB06XgQ/1SIuLq5du3ZbpICtFlTQnEpikXxAI1M+A0yTpPWsIMIui1fABG90\nfTv28rRtgetB/p5NxW83Ly9PTo+uyH1GRSObaauEoiTD0nKdTqffekhskPznaJHbAjm6T3o8Ho/e\nTTCNqmV6zWvQgSxeK6xRqL2hG5laq7as0wMGKX5kclU1mV5IcR2J38r48eP9tq4HCp8XXiQn+VdD\nKFhtbYF8x1SNYUZGht96SGyQ/OfIXEyLsrvyer2PPfYYd87vmDFjkpKSkpKSDh8+TLchW2JYBiI5\nuj4lFvbsmojxTxy2IDAIbEH5vKuLFy9WVFTExsYihP75z38eO3Zs1apVx48fnzBhwldffUW3LQAA\nAKsBOyb0g/Ls6tq1a1FRUbi8YMGClJSUgQMH5ubm3rhxo6qqim5bADk6LVLjjX/8KzBhEkBx5K15\n5JIpqD0LKnAXYs3Bh/OuAsXn+3kb5w8//DBnzhxcDg8Pb2hooN4WWxiw8iDXBHOL1MguGf4kI+93\nERgbGPMhcr8/BI+1ldtqSPitJlzodiwsIjGjmJNEQsK2CpLaCMUTAnlXYuinCVdUVJw4cWLnzp03\nbtyIi4tDCHk8nrq6uhYtWtBtyx7QfcSL3hi/wV3hNqfHE5vkbr6wwuMXtr7JhOjxuUPATDP03dVT\nTz31xBNP/PGPf/ztb38bGRmJEDpw4EDz5s0jIiLotgXYAAvuLjES63cf7q2AdaC/GDhnzpzjx4+f\nOHHij3/8I76SmppaUFBAvSGAHNssXgeCKbnPELsiBHvuLu+4SYz9CkvYVuF0Oin6Wrq1cSjM3iB2\nJUaXvKuIiIjq6upz58653e7q6uq2bdu2bdtWj4aCDc1fO50Wr/U78Va8iUNt/ZLlxM7rjUyTVD7V\nUBmxTitHIAgXzRTu+KrSbA2OXfEfEWtAmjAepYyMDJJxUNsWpAn/gqqqqr59+/bv33/o0KFZWVnp\n6ek5OarPOWUFg5dKDLhhqeqRWj0mpglz/WIlTVJOp3IwDykG/K2W1i3XFmGarWH/+vzqEfwuoXWy\nMK2UZL1P1jYM+u5qwIABrVq1KiwsLC4uLi4u3rRpU2lpaXZ2NvWGAFMIhmCG3z7qEXOi+3AKAyDc\n0GhNTNdmugAWoeyuzp07FxYWlp+ff8cdd0RGRkZGRqakpBw+fLi0tLSmpoZuW0yj+aav7Y3MzfoZ\nQvkTsXjsSvmmqbZFurdgwcDyK6cY15FrLsDaAh8KiF2Joeyubt682aZNG8HFli1bOhyOGzdu0G3L\nalh52sHcrN82UBx5yap0/dYJWkzYVkF+F+YLU3iXeEKp7XmAiZ3XC94Y4Mqhcm0C/E6LNXxMkHcl\nhrK7atmy5cWLF71eL/9iaWmpx+Np1aoV3bYAwIIo320NXgISNydecuTb4DLJmqTmjnBv5DsYudr4\nxpIiJd/IN5b0mnIGXJnEgNxYIFJyePlvZGtN2Egou6tOnTp17NixV69e27Zt+/7777///vvFixdn\nZ2cPHz7c4XDQbcs22GwV2/Q7sq4YOYdO2FZB3pzY6xAa+61BVW24TBKHM+amLG5Fc7tUBEu6T1Wj\nHczQ32qxf//+++67b9asWZmZmZmZmatXr3722WcXLlxIvSEmsMgKoW0Wr5nD6XSqvQHJfWfEO/tJ\napb7jU+CZUMjesSuqNcm6d58iwZJViU2htiVGF3yrtasWeNyuQ4dOnTkyJHS0tL/+Z//0aMVgBzJ\nxesm0wtVeVOLuF5JFLSZKxuPvLZJktx1UzZwB25GEZ3yrhTmT+LaVE22+MYUR9WyH5BO6Hg8Y+vW\nrVu2bKlf/UEC/96k+VeSTs9z0y9NWFudTqeTuy9IpkMKftsGokfOZ/iygxrwAAAeq0lEQVStUyHm\nQSsFldY4B/L5GlmmBXm7+GtGft36ZSag6a7q6ur27NkzYcKEs2fPUqwW4LDaryTxzjFV9tQ1aCvT\nnazwU0c1v5diORCsnNZKpY+CGZI2DbgSs/pu+hgaCTV3VVlZ2atXr2nTphUWFj7yyCOPPvoorZpZ\nxMQFKHHTpoix5sqhtp+TkvMhPpJbzgJpkSLaHmeFsU5oROBX9I5d8ZvjT5io7LagOKrW+YCMgZq7\nevrpp5s3b+50Ok+fPv3ee++VlZW5XC5alQNBhX6boyTzdeg2IdgUx9wPWD6WDY1QjOsQeiC98+e0\nmVn2A9IJmrOrjRs3tm3bNiQkZMCAAV26dHn99ddpVQ4A5Pj1QH4nSapqsyuBTMsMRuBytO2hh2wn\n60PNXQkOYMQPZddc29mzZ+fOnTtz5ky5k0eOHj26jceFCxc0t2V9bHbH1DB5orUXXL8WDcawxCYN\nrZib2KRrhYC56LgzUDMul2vkyJExMTFpaWnz5s1bu3at2GbXrl0rV6788hZXrlwxXqdFIE/T0Ywe\n925T/AH/JBGDW7QxZuVdKXsjgSo5Y+PDPxC70kwYxbqGDBnClT0eD0IoKSkJ/xkSEnLixAnCehYt\nWjRmzJiJEycihG6//fbc3NyxY8eGhoYKzPr06TN//nwKugH7krCtAt1y24md1/vkF/F1WvuyTdhA\nAUvFrji3RDeuA7ErK0DNXd1333319fVyr4qdjQJFRUVPPfUULg8YMKCurs7pdA4YMEBgVldXd/Dg\nwTZt2vTs2VODYIA6TaYXmri7JrHzep95rQMkaE6tpWsMMAo1d7VmzRoq9dTW1jY0NHTu3Bn/GRIS\nEhkZKRkGKygouHjxYklJSUxMzKpVq+Lj46kI0Mz5yfGeGVre6Ha7BQW/lroiboV/5fzkeM4tSerh\nX+TKjoVF6B03/lP8f4Gx+KIqY4U3KhtrfqMxrRB2vLy8XHxRVStUZARofH5yPHdRrmw1uJFniKys\nLLMlqIDmYqA2vF5vY2MjLjscDp/PhxBq3749Z+BwODgDjtzcXLwS6PF4pkyZMnny5D179hglWZr4\n+PiywN7ot4b4+HiEdP+3KpbBb/fSpUvKevgXuTIu4JrFZSljhJA7sfN63y3jJrfKCE+k4uN5Bkjw\nRn6L3P+x+IyMDBljpTdqNlZsUVUrv3ijgkHHjh151SLyN3Jl/htJjG+dDe//jVIGCmVZ+KfRB2JD\n3YyWeEKzwKvi0o0SExP91mM65m+12LdvX8otPB4PfnD7qVOnOIObN29GREQI3tWuXTtccDgcEydO\n/Oabb2pra3XVafGtYpLosaFAj3Vwqz0ETxXKQS9WwgaBLKaRRJIUrvPL5ycTrZEYH9eB2JUVMH92\nlZWVJZiQxsbGVlT87BsqKytra2u7du2qUAOOmYWFmd8XwC+BR5j4uycAAAgezJ9diRkxYsTq1avr\n6uoQQu+++25KSgqeYm/evJnb1M5tzayqqlq6dGlycrIe52lZNuHJssICx2/XSFKOAAEBjpjaPFxJ\nY/jUgACxoruaMGFCx44d09PT+/Xr53Q6ubOyjh079uWXX+LyjBkzkpKS0tPTMzIy6urqli1bRqt1\nzU/aDhKYHgdLnczE0GMjlBH3UdIzGZy6RDcnCfKurIAVF9AcDoek+1mwYAFXLioqMlARQghlZGSU\nLUYIocTO613fjg2kKqaXszIyMtB2OuL5C4PGhAbJF/FpzeH8tqihFc1+jrqDFD+MXAGDwz8Qu1Jr\nZn2sOLsCDEDuvG3yOxq5g7Hx0qW9gcwnwFKAu/oZK2z8s87SkB6nDNtm7SvIAbcEmAW4K1Ly8vJo\nVaV5KdmsNWh+u/xx4K5rcPaB9EXzqalOp9Pgk1X9Vmvi6bF+HQ//sy4a2UyuTv7wKrTL/3+A+kk+\nR8LPWvL7rLmPJGXCcTC4j0wA7oqU8ePH06pK81KyWWvQ/HbxOOAZlbIeZR8WSF80n5qakZFh8Mmq\nfqulezIs+cm2JJMk/ndeoU7+8CrY8/+vtl+CMsnnSPhZ0+0jSZlwHAzuIxOAu2KJQOL/xqx2WmFN\n1cbotBAnl7QLAJYC3BXwM7AhAgAAKwPuihTjEy+sCYi3VIuaJ0O2T11iWjyhWbDlXYG7IsX4xAvj\n8buUl7CtwnjxFLcUGi/erM9abn1P1XFQTKcuMS2e0AzyrgBmsFSgCGIeAADoCrgrQ7GUg1HAt2gQ\nK1ItAt2H8qmqFn4oAEFCULsrVXfkYIhdSQoTjBLT4R/LilebESVpoMpv2T78w7R4QjOIXQHSBEPs\niuJyuV/IN+VrbhEHvfitEFZFeOqSZIuCK4FHF3SaPNk+/MO0eEIziF0B0pD/QrHmMpqc/svTtkna\nBPjEBF3R+4kPX331FT2xEig8IUKhjP2WWU/BgLK9y0wA7ooUir9QTHl0npx+uk9V0ICG2YMemvnl\n1NRUtZKUEfSRa8u3aJDefYEylKn/mzULcFek6P1LBGfpmjszo7hcTqs5U1qkiIY+yvlvuuJtH/5h\nWjyhGcSuAGmM/CUS+AMmOLenyv8ZGbsir4rpqCF5H/3OMumKt334h2nxhGYQuwJYhdwzWTO6phbj\nN3CrbRG2mAMARcBdCdF2i7HUA/fkvJE9vFTgGONFAkmlAgBADLgrUuTWf+3kAyB2Rd2XcC0q12zw\nyNNt0YKfI3lz1hRPaAaxK0Aa26z/KgCxK1WQ+DbjAxWE2D78w7R4QjOIXQUjapfy7DSjAgAAYIKg\nc1d+PY3cxFlzqp1FfJudUg6NTBOWS+nlX5crWyH9E8pQtsK/WSoEnbvyi9o5uNg+kDN/+TgWFgVe\nCYfaVEFaNgpoXnnXnA7pdDrVpgkr2HPCArTh+hi4DXmZTyA2dPtIWCb5HAk/a7p9JCkTjoPBfWQC\ncFekWPyjpTKHs3LsisovAKbDP0yLJzejVRXT4gnNIHYFGIHamy/fG1lkdZEt4IQOAGAdcFfMIOml\nwHXpBGRNAYDVAHdFijXDklbIxdF8B7dIyouu+o1PsiHE9qlLTIsnNIO8K0Aaa67/KqtSeyM2LIKC\nhWlYeSfskZwZ0+EfpsWTm9GqimnxhGYQuwKAIAUW+gDAyoC7MggTHyoYDHdhCDUBgO0BdyXE6XRK\n7l/Iy8vj20iWjUe5dbWpgiR9pDgOqlbeqaRDOp1OXFY+mZefJqxQJ7/aQGwI+0hxHBC9z5puHwnL\n/EYDsaHeR5Iy4TgY3Ec28Nka18gOgv/4F7kymrafK/PfiKbtVzbmypIXBQbKxgIDNG3/+fPn+dd9\nPh+/QFLmxkHbG41phfCNYmOFSnxSyF0Xc/78eUJLCwLizYJp8QkJCWZL8A/MrkyDZA96sO1T13Xt\nDhYGAYBpwF0BssD9HQAA6wDuihRd13k1P7TCmrk4hH4u8KwRyYYUWrfmcNGtihDj03oo6re9eEIz\n48WbC8Puyuv1Hj16dNeuXdu3bzegOWvmLlg2F0cOvi+hmzVC8fQpQiDvSg8zWlUxLZ7QLNjyrsLM\nFqCdOXPm5Ofn33XXXadOnRo5cqQBLSZsqyjLiTWgIQADq5EAAHAwPLuaO3fusWPHJk2aZLYQ7ZiY\njAUAAMAWDLsrh8NhZHPWXP+1bDyDVnN0zSw7XLYXT25GqyqmxROaBVvsiuHFQG243W7lMv8inw4d\nOkjaKF8MxLi8vNzvGzt06ICQm4okfh8VjMUtau648pCenxyPCx06dDg/+b8G3HUB/NrkILFBvxz5\nAGsjbNH24gnNQLwqs8CrysrK8vt268CMu/J6vY2NjbgcyLwqPj6+TLEcHx+PECqTeiN3Xe6NyjWL\nW8H3fblyx44d4+PjBddJ3igy1veNUsayb/RnLCibhhU0aAbEmwVz4l0uFy4kJiaaq4QEZhYD9+3b\nl3ILj8djthztBFvmLwAAABWYcVdZWVkltzA4aoUQSuy8nmT99/K0bQaI4UNxVdq3aJBZsSvlHYAQ\n/lFVFSG2D/8wLZ7QLNhiV8y4KzFer9fj8TQ0NCCEPB6P3lMuDWkQBmz8s2wuDoeCK5KsSmwPqUuq\nqiLE9qlLTIsnNIO8K2bIz8+fNm0aLiclJSGETJl4AdqAnCoAAFTBsLsaNmzYsGHDzFYhS8K2CjS9\nUFwOQsAzAQAQOAwvBhoA/z4rXv9N7LxecCPWaY2YxYfg0WqOrpllh8v24snNaFXFtHhCM4hdAdIY\nH10gwax4BpUJk/HRBabDP0yLJzejVRXT4gnNgi12Be5KSOCnhZKX9caAvnB+i8S+aGQzI7VpLhOe\nJgxlKNumzAZmnw+pL5KHBfMPCPaJjveVvMgdRCs+QVihHOBxvfhwUkaP6yU/uteCMH0sLIg3C6bF\nw2nCtsL46AIJlo1nnJ/sP70fYld6VEWI7cM/TIsnNIPYFSAN4VKywQ9Zt2w8g1ZzdM0sO1y2F09u\nRqsqpsUTmkHsClANPFcJAxvWAQDQD3BXwQV4FAAAGAXclR+4mZMe0YXAnYdl4xm0mqNrZtnhsr14\ncjNaVTEtntAMYleANJbKu+L8nN7xDMkn+NGaokHsSo+qCLF9+Idp8YRmELsCAGlgIREAABMBdyWE\nML1OcpHQCml6qlJxBYeGQBkDacJQDrYyG5id+KUvgacJc+WioiKfv9TaoqIiigm8CmnCCi0qGItb\nFNfm85fhi238QpIySVgVRTOK4um2aHvxhGYgXpUZxaogTdhWWCp2pVOLkHelCohd6WFGqyqmxROa\nQewKAAAAACwHuKvghb91Qq4MAABgEcBdkWJ8ZgwJls3FodUcXTPLDpftxZOb0aqKafGEZpB3BUgD\nsSvqLULsSo+qCLF9+Idp8YRmELsC9AIW2QAAADQD7iooAE8JAADrgLv6BYmd18ul0eXl5UleV2tD\nN03PyTvPV6F+p9NJkhqsdx8FZb48BXv+/wn7GIgNIk4TJtFPsY+E40BYpvhZ0/0cScoUP2u6fSQp\nG/99JukjG5id+KUvCmnC5AcEG3Ncr7Y0YXElkgbGw/TJqiDeLEC8WUCaMGAmsAAIAICdAHeFEPt3\ndtb1AwAA+CUo3BWVuznJOq/xa8GELVI0o9hHpsXTbdH24snNaFXFtHhCM+PFm0tQuCsBJGfVJ3Ze\nL7hicGYMoYtlOnWJafF0W7S9eHIzWlUxLZ7QDPKughcSN2ZxYFUQAAC7Au7K6oAHAgAAQOCuyNE7\nuqDNLTEd/mFaPN0WbS+e3IxWVUyLJzSD2BUgjdWeGYjdG9PhH6bF023R9uLJzWhVxbR4QjOIXQUd\nCk+yMKvMh+Q6lOmWCZ9qAWUo26bMBmbnKeuL4Dx7ybPt1T6cQqenWvhED6dgOkkexJsFiDcLpsXD\nUy1shfHRBRKYDv8wLZ5ui7YXT25GqyqmxROaGS/eXMBdkUK4lGzwRj6mwz9Mi6fbou3Fk5vRqopp\n8YRmwRa7CjNbgARer/err74qLy9vbGwcOXKkpM3Ro0cvXLjA/XnvvffGxcUZIw8AAAAwHiu6qzlz\n5uTn5991112nTp2Sc1e7du06fPhwSkoK/rNLly503ZVv0aAm0wspVqjQkAGtAAAAsI4VFwPnzp17\n7NixSZMmKZv16dNnwS1SU1P1VmVkdIHchzEd/mFaPN0WbS+e3IxWVUyLJzSD2JX5OBwOErO6urqD\nBw+WlpbqrQfDrf8q+JIA14g1zLSYDv8wLZ5ui7YXT25GqyqmxROaQeyKGQoKCi5evFhSUhITE7Nq\n1ar4+HizFUkAa30AAABUMN9deb3exsZGXCacVyGEcnNz58+fjxDyeDxTpkyZPHnynj175IwTExMR\nQh8nI4SQ2+3mrnNlyYuGGctRXl7u18aygHizAPFmwaL4rKwssyWowezEL19+fn7PW9TX13PXDxw4\n0LNnT5IaSkpKEhISampqxC9pThP2iTJ/i4qK/BoLbEhSgxUgyTrELRppRliV7cXTbdH24gnNQLwq\nM4pVMZEmbP7sKisrK0APX19fjxAKC9O3L1Z7ZqCqFq0ZQWFaPN0WbS+e3IxWVUyLJzQLttiVFbda\neL1ej8fT0NCAEPJ4PB6PB1/fvHnz2rVrcZnb61JVVbV06dLk5GS/C4niExeNhB/EgoAWAACAWsyf\nXYnJz8+fNm0aLiclJSGESkpKHA7HsWPHbty48cwzzyCEZsyYcfXq1YiIiOvXr6empi5btkxVE8a4\nLnBLAAAAtLDi7GrYsGGuX4JnTgsWLODcUlFRUUlJyZEjR06fPv3hhx/GxMTorcqaaRBMpy4xLZ5u\ni7YXT25GqyqmxROaQd4VII01l5KZDv8wLZ5ui7YXT25GqyqmxROaQewKAAAAACwHuCs6QJgKAABA\nV4LLXcntsCA5bTMvL4+KjdPp5P4M/CRQ/v8V7PmNGtNHkjLhOFDsI+E4EJ4mTKKfYh8Jx4GwbPz3\nmaJ+ip810/9mKfaRDcxO/NIXfmowxXN+xWU9YPpwUhBvFiDeLJgWz0SacHDNrgAAAABGAXelHYhX\nAQAAGAa4KyJ8iwaRrPMSrgVTXDKm2yLFPtJqjq4Z3cV6g78STIsnN6NVFdPiCc2MF28uTXw+n9ka\ndKQsJxYhlLCtAiGETwfGUyJBmZsnyZU5JC/qhNvttuapKCSAeLMA8WbBtPjExESXy2W2Cj/A7AoA\nAABgAPu7K7nN6xB5AgAAYAj7uyta4PVfZSfH9Do4oRmI16NF24snN6NVFdPiCc2CLXYVvO6K+wgF\n2yhI0ugCsXE6ndwjvAJPOcRV+bXPyMgwuI8kZcJxoNhHwnEgTBMm0U+xj4TjQFjmY8z3maJ+ip81\n3T6SlI3/PpP0kQnsv9UisfN65e0VfPxutTASpiO3IN4sQLxZMC0etlpYCwhWAQAAsEsQuSs55GbE\nAvdmzaVkpsM/TIun26LtxZOb0aqKafGEZsEWuwqixUDEW98jWQw0HabXFkC8WYB4s2BaPCwGAgAA\nAAAdwF39AutMrQAAAAA+4K7YXkpmOvzDtHi6LdpePLkZraqYFk9oBrErW6E2dmUpmF4KB/FmAeLN\ngmnxELuyNLTSLaFsvzJhmjCUoWybMhPYf3aFH8eOgdmVYYB4swDxZsG0eJhdWRe+r2J6Kdn4BXoQ\nr0eLthdPbkarKqbFE5oZL95cgnR2xQRM/1gD8WYB4s2CafEwu7IcDPkqAAAAgE9wuSsAAACAUWzu\nrh49GeXXhumlZKbDP0yLp9ui7cWTm9GqimnxhGYQu7IVTCzIysH0UjiINwsQbxZMi2fiVmnz2RUA\nAABgD+zvrqyQggdltsqQJgzlYCuzgc/WJCQk+LUpKioiqYrEjGJVPp/v/PnzBrdIsY+2F0+3RduL\nJzQD8arMKFZFcqs0HYhdWReml8JBvFmAeLNgWjwTt0r7LwYCAAAANgDcFQAAAMAADLurs2fPLl++\nfMaMGXPnzuXHxtXCdBqE8cklIF6PFm0vntyMVlVMiyc0M168uTAcu0pPT3/wwQf79Olz5syZjRs3\nzp8/f8SIEQIbJhZk5QDxZgHizQLEmwUT4sPMFqCd/fv3t2zZEpdvu+225cuXi90VAAAAYA8YXgzk\nfBVCqH379g0NDSaKAQAAAHSF4cVADo/H8/jjj2dlZf3hD38QvJSYmGiKJAAAALaw/mIgM+7K6/U2\nNjbissPh4L/04osvXr16dc2aNaGhoWZIAwAAAHSHmcXAffv2pdzC4/Fw11966aUff/xxxYoV4KsA\nAABsDDOzK0lmzpxZVla2bt06fhwLAAAAsB8M7wx89dVXT548uW7duoiICDzfEiwSAgAAALaB4dmV\nYBuFw+EoKSkxSwwAAACgKwy7KwAAACB4YGarBQAAABDMgLsCAAAAGIDhrRbKnD17dv369bW1tYMH\nDx48eLDZcpTwK/Xo0aMXLlzg/rz33nvj4uIMk6cKr9f71VdflZeXNzY2jhw50mw5SpBIZWjkEUJn\nz54tKChwu92RkZGPPfZYamqq2YpkIZHK0OAfP368sLDw0qVLoaGhDz74YFZWltmKZCGRatmRD/3T\nn/5ktgb6uFyunJycBx988K677lq4cGFoaOg999xjtihpSKQuX7589+7dXq/34sWLFy9evOuuuzp0\n6GCKWr/Mnj37zTffvHjx4tatWydOnGi2HCVIpDI08gihrKysVq1a9enTp6qqat68ebGxsd27dzdb\nlDQkUhka/A0bNlRVVaWlpTU0NCxfvvzSpUsDBw40W5Q0JFKtO/ImnmSsH+PHj3/99ddx+cCBA717\n925oaDBXkhwkUmfNmjVr1izDpWmhvr7e5/MdOHCgZ8+eZmvxA4lUhkbe5/NdvXqVK7/99tuZmZkm\nilGGRCpbg8/x0Ucf9ejRw2wVRMhJtezI2zN2VVRU1LdvX1weMGBAXV2dZU98IZRaV1d38ODB0tJS\nY9WphqHUN0KprIw8Yuqhz4RSGRp8jpqamvbt25utgggFqdYceRvGrmpraxsaGjp37oz/DAkJiYyM\nrK6uNleVJORSCwoKLl68WFJSEhMTs2rVqvj4eGOVBi8sjrzH41m3bl12drbZQvyjLJWhwT9x4sSW\nLVuuXbv2/fffL1682Gw5SpBItebI23B25fP5EEL8Xw0Oh4N7PK6lIJSam5t7/PjxDRs2HDt2LCEh\nYfLkyYaqDGIYHfnp06dHR0dbPHaIUZDK1uC3bt36nnvuadeu3Q8//HDixAmz5SjhV6p1R97s1Uj6\n1NfXJyQkHDlyhLvSu3fvffv2mShJDg1SS0pKEhISampq9FenHSZiVxhyqUyMvM/nmz59+pNPPnnj\nxg2zhfiHXCorg+/z+b7++uuEhIQff/zRbCH+IZFqqZG34ezK4XDExsZWVFTgPysrK2tra7t27Wqu\nKkk0SK2vr0cIhYXZcBXX4jAx8jNnzvzmm29WrVoVGRlpthY/qJLKxOBj7r77boTQ+fPnzRbiHxKp\nlhp5G7orhNCIESNWr15dV1eHEHr33XdTUlIssvYqRk7q5s2b165di224zRdVVVVLly5NTk627I4G\nr9fr8Xhw5Nzj8fCPerEaclIZHXl066HPeXl5+KHPVh58OamMDj4ntbGxceHChdHR0enp6eZKkkNO\nKhMjb89nBno8nilTphw8eDAqKqply5arVq3q1KmT2aKkkZM6c+bMGzduLFu2DCHUr1+/q1evRkRE\nXL9+PTU1dfHixTExMWYLl2bPnj3Tpk3jXykpKbHId12AnFRGRx4x9dBnOamMDv7gwYMvXboUERFR\nU1MTFxf3+uuv9+rVy2xR0shJZWLk7emuMNXV1VevXrWso+LjV6rH4ykpKUlKSrLmrd/GwMibCEOD\n7/F4ysrKunbtGh4ebrYWP5BItebI29ldAQAAALbBnrErAAAAwGaAuwIAAAAYANwVAAAAwADgrgAA\nAAAGAHcFAAAAMAC4KwAAAIABwF0BgBEUFxe7XC6/NmfPnjVGDwAwB7grAKDJrl27PvzwQ8HFqqqq\n559/PiTEzz+3y5cvT5w40ev16qYOABgG3BUA0MTpdP7rX/8SXFy9evU999yDnyiqwODBg0NDQzdv\n3qybOgBgGEs8ZxcA7MG5c+d+/PHHmpoa/JDQDh06xMfHNzY2bty48c9//jPf8uzZs2VlZQihFi1a\n9O/fPzQ0FF8fMWLEhg0bRo8ebbx4ALA44K4AgBo7duz46quvEEIvv/wyQignJyc3N/fAgQM1NTUP\nPfQQZ7Z48eK8vLzU1NTmzZufOXPm/vvvX7BgAX5pwIABixcvdrvdlj1DAADMAtwVAFBjxowZlZWV\n//nPf1atWsVd/Prrrzt27Mg9KtTj8axcuXLx4sXDhg3DVyorKznjbt26hYSElJSUgLsCAAEQuwIA\nfamoqBAfuXn58mWu3K5dO64cEhLSokWLM2fOGCQOANgB3BUA6Et9fT1/T6DD4cjNzf3LX/7St2/f\nqVOn7ty5U7AVMDQ09Pr164bLBACrA+4KAPSlWbNmtbW1/CsTJ0787LPPXnnlFYfDMXv27MmTJ/Nf\nvXnzJhOHtAGAwYC7AgCaiJOrevXqhfdf8ImJiXn00UcXLFjw2muv7d+/v7GxEV/3eDy1tbV33XWX\nEVoBgCnAXQEATbp3737ixInS0tKampq6ujqEUJ8+fWpra8+dO4cNzp07t337djzf8nq9J0+ebN26\nNbeR/cCBA2FhYf379zdLPwBYFnBXAECTESNGpKSkPPnkkykpKX/9618RQnfffXf37t0//fRTbOD1\net9666177rnnnnvu6dWr1z//+c8lS5Zwby8sLBw6dKilThwHAIvQxOfzma0BAGzOJ598smDBgsLC\nQu5KXV1dcXHxHXfcwQ9TVVdX9+/ff+vWrYmJiWbIBABLA7MrANCdrKys6Ojo7du3c1fCw8Pvv/9+\nwZaKtWvXDhkyBHwVAEgCsysAAACAAWB2BQAAADAAuCsAAACAAcBdAQAAAAwA7goAAABgAHBXAAAA\nAAOAuwIAAAAYANwVAAAAwADgrgAAAAAGAHcFAAAAMMD/B/CSya4PgdtjAAAAAElFTkSuQmCC\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x = Accel./-(wd^2);\n", "ax=max(x); % max amplitude\n", "initialvel=-ax*wne;\n", "vtb1_1(m,c,k2,x(1),-ax*wne,3.2);\n", "hold\n", "plot(Time_domain, x)\n", "grid on\n", "grid minor\n", "title('Case 2: position vs time')\n", "xlabel('t(s)')\n", "ylabel('Position (m)')\n", "legend('Estimated Displacement','Experimental Displacement')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3 Mode shapes\n" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA\nB3RJTUUH4QoZBDkZG4vvlQAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJ\nbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAyNS1PY3QtMjAxNyAwMDo1NzoyNSt0R0gAACAA\nSURBVHic7d17WFTVwj/w5cAOBwVFQQotQUssoQLxQVPzSlGeNLw8ebxkF61EE2+Ybz/FXuOcTPNS\npimmR3nJUo95Oz6iBOIFypNoCWKICicHPDoIgeAwDDPz+2PXPnNmhmEPs2bf5vv5a8+ezd5rFnvm\nO2vtNWu3M5vNBAAAQGwqsQsAAABACAIJAAAkAoEEAACSgEACAABJQCABAIAkIJAAAEASEEgAACAJ\nCCQAAJAEBBIAAEgCAgkAACQBgQQAAJLgLXYBJC0vL0+j0fj4+Lz88suec2gBNDU1HThwgBAyYMCA\nXr16UdyzkurNZDLl5OTcvXuXEPLII48MGjTITQeSS6XxPG0EqzdLcqlDqTN7trfeesvPRmBgIPts\nQkICIYR76Kzc3Nw5c+bMmTNHp9M5u5mLh5a4qqoq9vT78ssv6e5ZMfVmMBiioqK49+nkyZPddyy5\nVBqf00bIerMklzqUOE9vITU2Nt67d89qZVNTE7sQGxtLCPHz82vbzi9fvrxp0yZCyKpVq1zfDDzK\nP/7xj4sXLxJCZs+e/dxzz3Xv3l3sEskD6k3WPD2QOKtWrerfvz+77OXlxS688MIL/fv39/HxYR9e\nunTpzp07Pj4+Q4cOLSgouHr16rPPPtu9e/erV6/+/PPPzc3NarW6e/fu/fv3V6lUP/3009WrV9k/\nzM3Nbd++fUBAAHcIDp/Nfvrpp5KSkoCAgOeee85xSdhnf/nll8LCwubm5kceeWTw4MFWR3T8rOVB\nq6qq2EOw3RFPPfVU37592aeuXLny4IMPjhgxwuqviouLL1++3Nzc/OCDDw4bNkyl+s9FysbGxszM\nTL1eP3DgwI4dO9o9KM+yEULs1rntS7CqN0KIyWT6/vvvtVqtTqfz9vYODQ0dMGAA96xtxdp9mQ7K\nyadgrVbXpUuXfvjhB3a5Z8+eHTt2DAkJsfpbxy+kJTdv3iwoKNDr9T4+PiNHjvT397fdxvGenTox\nXK9PPqcNp6V6c/B+ofJmsdKGV0qxzuVN7CaayF599VW2HrKysmyftWqGcw/fffdd9q8OHDiQlJRk\nVaVTpkwxm81jxoyxWj969GjbQ7S0GXes2bNnc089//zzDkpiNpurq6utdvjUU0/9+uuv7F85ftZu\nwQIDA6dMmcJtv2PHjgULFnAPx40bx21fVVX1/PPPW+78kUceOX36NPtseXl5jx492PUMw6xfv55d\n5vpenCpbS3Xear1Z/sc5UVFR//73v63+/I033rD8j+j1ej7ldFAwK46riy2GJdtOKscvxFZtbe2E\nCRMst+/WrZvVq+azZ6dODBfrs9XTxkpL9Wb3/ULxzeL6K6VY57KGQPr9PBg4cOCYP6SmprLP2g0k\n1vDhw0eOHLl161b24aJFi/bu3Zuenv7qq6/Onj3bbDZbNrkmTJgwefJkbreWWtqMO1a3bt1ee+21\nxx57jH3IBqdtSQ4dOmT+48Tt0KHDmjVr0tPTH3nkEUJIv3792GM5ftYK9+YJDQ197bXXAgICCCE+\nPj5qtXratGn9+vVjn83Pz2e35z5e33rrrZUrV3br1o0QEhAQwL6pBg4cyD47ZcqUd999l+sF5T5Z\n+Jft+vXrLdV5q/VmNpuTk5MXLFiwe/fuvXv3pqSkMAzDltnqz8PDw5OSkp599ln24QcffNBqOR0X\nzIrj6lq7di1XY8OHD588eXJOTo7VHhy/EFsvvfQSu8MhQ4Z8+umn7777Lle9Vue54z07dWK4Up98\nThsrLdWb3fcLxTeL66+UYp3LGgLJ+osJISQhIYF9tqVA2rt3L7uG6x/44osvKisrrXbOXhkihNy7\nd89BGexuxh3r8uXLZrM5Pz+ffbhx40a7JTGbzVeuXGFX/vWvf2XX7N+/n12Tk5Pj+FnbUnHvgevX\nr5vNZu67GFuA48ePsw8zMjIsD/3qq6+yf/7NN9+wa1atWlVaWsouT5s2jX2WC3L2k8Wpsjmuc8f1\nxjIajYWFhceOHTt69Cj7+RUSEmL159xXV3Y0V48ePVotp+OCWXJcXeya9PR0do3dtnurL6SlI1o2\nFrkv77YX5B3s2akTw5X6bPW0uXPnzjcWysvLW6o32/cL3TeLi6+Ubp3LGgLp90BatmwZd2afPXuW\nfdZuIHXo0IH7c6PRGB4eTv7g4+Pz/PPPcx0vrgeSn58f+5AbeTFnzhy7JTFbnOK20tPTHT9rWyr2\nPRAQEMA+/PLLL9mN2W9hp0+ftvxbbue7d+9mt9fpdOyaadOmHTp0yOpAP//8M7uG/WRxqmyO69xx\nvZnN5i+++MJ2lAq3Pfvn3Ks2m82TJk1itzEajY7L6bhgdv9TdquLXdNqIDl+IS0d0e5nltV57njP\nTp0YrtRnq6cNdyzLzRwEkuX7he6bxcVXSrfOZQ2DGn43bNiw0aNH89lSrVZzyyqV6vvvv9++fXt+\nfv4vv/xy5cqV48ePnz17trq6+oEHHnC9VNx4ilZLwhaGXRg+fDjbIcB55JFHampqHDzb0iG48R2c\nPn362G7GHbq5uZldaGhoYBe8vb25Z00mE7tgNbLRccltj9VqnbdUb8XFxey1peHDhycmJjIM8/HH\nH3MtGw73KqyO67ic/E8Gx9Vlt+RteyG2R+QO5OKeeZ4YrLbVJ3eStHTaBAQEWF6V4a42tcTqnevg\n0G17s5C2vlJ31LlMIZBcYjKZGIZZvHgx+3DPnj2TJ09uaGjIy8sbMWIEdxZWVlY6OHV4btYqblhO\nZGTkZ599xq3/8ccf+/fvf+vWLQfPtvmgLHZ8PCFk586d06dPJ4R89dVX7JpBgwZxBdu5c+eMGTMI\nITk5OfxLbnUsx3XuuJy//PILu7BgwYKxY8c2NjZaXhbm3Lt37/vvvx80aFBNTc2pU6fYgrVaTv4F\nc1xdjl+CUy/E8oheXl5Go3Hz5s2vvvpq+/btCSE3b958+OGHXdwzH22rT61Wyz5s6bSJiIj4xz/+\n0bYiuenN0rZXevDgQfYhxTqXKQSSS5qbm7t06RIfH//QQw95eXmdPXuWEOLj48OetdxXoSeffNLP\nz2/69Onr1q2z3QnPzVrVvXv32bNnf/HFFxs3brx69WqvXr0aGhouXrxYWFio1+sdP+tie+6hhx5i\nd56TkzNo0KCHHnro8OHDhJBevXqxH3+TJk3at29fbm5uTExMYGAgW1E8S25VNsd17hiX90uWLMnM\nzPzhhx8qKipsN/Py8oqLixs9evTFixfv3LlDCJk7d26r5SSE8CyY4+riU+E8X4jlEefOnfvpp5/+\n/PPP4eHhsbGx9fX1BQUFt2/fdnHPfLStPoODgx2fNq5w05ulba/UHXUuV2L3GYqsbcO+uQ0MBsPj\njz9uWZ89evRgR2CzFixYwHUNz5w5s6Vi2G5mdSy715Ds/iw8JSWlU6dOXHnUavWkSZOMRiOfZy1x\nI03Zh1y3dVVVldlet7XRaHz//fcte0VeeOEF7tp+dXU11yMaEhJy5MgRdtlyuBTPsjmuc8f1Zjab\n16xZw3V9JCUlWV1z4v5848aNbL+fl5fX+++/z6eGWz0ZLDmuLjOPa0iOX4hdK1euZAdosZ566im7\nleZ4z06dGK7Up5nfaWPFwTUk2/cLrTeL66+UYp3LmqcHEhVGo/HKlSs5OTkOfgUisNLS0pycnJbG\nejl+1hVsVZw9e5YbwWXp119/vXjxYktvaafK5kqdNzQ05Obm2h1pYvnJZTAYTp8+3dKAlJbK6VTB\nHFdXqxy8EAeuX7+ek5NTW1tLfc+2XK9PM7/TxhXU3yxte6W06lzW2pnNZgIAhBBCxo8ff+DAgcDA\nQO4CBrgC9QlOkdw1pNLS0qysrLKyMl9f33HjxkVHR9tuc/78+fLycu5hTExMaGioYCUEAAB3kFwL\nacCAAcOHD4+Njf3ll1++/vrr1NRU2+lAli1bdu7cOW5O38mTJ9vNLQBnNTY2siN3W502DfhAfYJT\nJNdCys7O5iZ87NSp0+bNm20DiRASGxubmpoqbNFA+djx0EAL6hOcIrk7xlpOP9ytWze7PzQjhOj1\n+jNnzly+fFmocgEAgHtJroXEMRgM6enp48ePt/tsVlZWRUVFUVFRcHBwWlpaWFiYwMUDAAC6JHcN\niTNv3rza2todO3bYTpih1WqDgoIIIQaDYf78+eXl5UePHrXdg+XEYpmZmW4tLQCA9En8u7tEW0iL\nFy++c+eO3TQihLBpRAhhGCYxMXH8+PE6nc5qYjdWSUmJewsqbWVlZRI//wSASiCoBEIIKoGQsrIy\nsYvQCikG0pIlS65du5aenu7r69vqxuztxnlOSQkAAJIluUENy5YtKyws3LZtm1qtNhgMBoOBXb9n\nz55du3axy9xNbmpqajZu3BgZGcnezwoAAORLcg2Lffv2EUKGDBnCPmQYpqioiBBSUFDQ0NDATvqb\nnJxcW1urVqvr6+ujo6O5+wkBAIB8SS6QWrrqs3r1am45Ly9PqOIAAIBAJNdlBwAAngmBBAAAkoBA\nAgAASUAgAQCAJCCQAABAEhBIAAAgCQgkAACQBAQSAABIAgIJAAAkAYEEAACSgEACAABJQCABAIAk\nIJAAAEASEEgAACAJCCQAAJAEBBIAAEgCAgkAACQBgQQAAJKAQAIAAElAIAEAgCQgkAAAQBIQSAAA\nIAkIJAAAkARvsQvQFqWlpVlZWWVlZb6+vuPGjYuOjha7RAAA4CpZtpCmTJlSVlYWGxvLMMz06dMP\nHDggdokAAMBVsmwhZWdn+/v7s8udOnXavHlzQkKCuEUCAAAXybKFxKURIaRbt27Nzc0iFgYAAKiQ\nZSBxDAZDenr6+PHjxS4IAIDUdSocKXYRWiHLLjvOokWLAgMDExMTW9ogPDycXcjMzBSqUBKi0WjE\nLoL4UAkElUAI8fhK6PV52d1hYheiNTIOpMWLF9+5c2fHjh1eXl4tbVNSUiJkkSQoLCxM7CKID5VA\nUAmEEE+thHaLcgghd4e9WRuZ00Xswjgm10BasmTJtWvX0tPTfX19xS4LAIDksDnEMq8dWX1YxLLw\nJctAWrZsWWFhYXp6ulqtNhgMhBCGYcQuFACAJHBRZF77+0Wj6sNhXcaW1ZaViVcoXmQZSPv27SOE\nDBkyhH3IMExRUZGoJQIAEJ9tFMmLLAMJV4YAADhWvXNWz7LNI2FL1EayDCQAACDybxJZQSABAMiP\nwqKIhUACAJATRUYRC4EEACADbcshGV1AIggkAACJU3CTyAoCCQBAolyMInk1jwgCCQBAahwP41Yw\nBBIAgFRQ7J2TXfOIIJAAAKTAcy4UOYBAAgAQjcf2ztmFQAIAEIFbm0Ry7K8jCCQAAIGhd64lCCQA\nAIEIE0UybR4RBBIAgLuhScQTAgkAwF0QRU5BIAEA0CdWFMm3v44gkAAAKMIwblcgkAAAKEDvnOsQ\nSAAALpFOFMm6v44gkAAA2ga9c9QhkAAAnCOdJpEluTePCAIJAIA/aUaRYiCQAABaJ/EoUkDziCCQ\nAAAckHgOKYzkAslkMl24cEGj0RiNxgkTJtjd5vz58+Xl5dzDmJiY0NBQYYoHAB5CRlGkjOYRkWAg\npaSkHDt2rHfv3sXFxS0F0sGDB8+dOxcVFcU+7NWrFwIJAGiRURQpjOQCacWKFampqadOnZozZ46D\nzWJjY1NTUwUrFQAoHnJIdJILJIZh+Gym1+vPnDnTpUuXfv36ubtI0Krqw2GON1BGfwIoFaJIIiQX\nSDxlZWVVVFQUFRUFBwenpaWFhdn/QAwPD2cXMjMzBSydVGg0GjftuVPhf71vayNzWtryd/+dWK1v\nT4/7KkFGUAmkhUro9fnvX5VuzA0jhJSVye+bU6fCkbWRObU8Sq7RaFr6qJQIWQZSUlIS219nMBjm\nz58/d+7co0eP2t2ypKRE2KJJDt3zj2sJWbV4urRejv9sX304jMszYVpOEn8TCgOVQCwqQUlNoupC\n5fxzZRlIQUFB7ALDMImJiePHj9fpdGq1WtxSKVtLUdQG3B6qD4exu0WHHghGSVFEFDS+jiXLQLLU\n1NRECPH2lv0LkSaKOWSL3SdiCYTR6/MyQsqIUqJIkST3OW4ymYxGY3NzMyHEYDCQP4Y57Nmzp7Gx\nccaMGYSQ/Pz8Z555hhBSU1OzcePGyMhInkMhwCnC5IRlg0mAw7XZ1Ykh/Dfu8/dK95UEnMW2im7M\nDVNM1xZLYc0jIsFAOnbs2MKFC9nliIgIQkhRURHDMAUFBQ0NDWwgJScn19bWqtXq+vr66OjoTZs2\niVliJRIlG7qMLZNCa6ml4OGfMVcnhri+E3CdVe+cHAcseJp2ZrNZ7DK4S3h4uIcPaigrK3P2K6Ho\neUC9DK1WglV4uC8zbFNKsHxqw5kga3YvFCmvEpxtIUm/BiTXQgKxuPVykVO4a0tuKoxgCWTF6kBW\nDSk0nqhgo8gTrhIpr7+OIJCAJYWGkRWuB49WqaT26W9ZBstwkkLZZEdhY+c8FgIJpJhGLCqZZEge\nfJUQIu0Peqtwsl0JLfHMKFJk84ggkEDiZ3abu++4j3VmTZ7E+82tcDmEZHLMc3rnPAcCyXNJtmFk\ni39TyfZDXL5jq9iXgN48K57ZJPIQCCQPJaM0YjnOJKldH6LIssHEvkyFvUD+EEUsifdquAKB5Ilk\nl0YsrvvOsuSe03qwajAp/vVykEOeA4HkWWQaRVbYTPK0z2WWR8USosjTIJA8jtzTqCpdTwghJEzZ\nn8WOKT6WEEUtUXB/HUEgeQ4FtI2seucU8IpcZDUeTwGxhBzycAgkjyD3z267F4rYYQ4ilUha+vy9\nUu6tJUQRH8puHhEEkieQdRo5HrPAZpJMXxpd8m0tIYqAg0BSOPmmEc/hc8gkK3JpLSGHnOUJ5zkC\nScnYO4XL7iR2diQ3MsmK5ZAHCWYSoghagkBSrOrDYbWROfKaNceVHxUhk6xYNpWYNXliF4cQTPbj\nAg85vRFIyiTHq/2u9DLZ/c0scE0lQ/JgIl5TCU0i4AmBpEDcdaNamUzjRuuCB/ru7Orz98qyMnF+\nR4woAqcgkJRGXqMY3PEpiUyyi+vBEyCTkEPQNggkRZFRGrlpDjrqt/VTEq4Hj7itqYQocgfPOZ8R\nSMohuzRy02ciMskxdzSVuBwiiCJwAQJJIeSSRsJcycAkDo5R/K0SmkTu5lFfrRBISoA0ssuj3snO\ncv23SogioA6BpBAS/+QVfogX13FHInNa39pTta37Dj8nEoynfamSZSCZTKYLFy5oNBqj0ThhwgSx\niyMy6Z+yYs0XwGZSp8KRJEzS9SMu/iMd0CQCd5NlIKWkpBw7dqx3797FxcUeHkgSTyPRZ1TDAAee\nHDSVkEMgGJXYBWiLFStWFBQUzJkzR+yCiEzin7Psp5voc6nVosuOH8umEqvdohyudw5pJDyJv8Hd\nQZYtJIZhxC6C+KQ8ikz0hpEtD3xvt4FlJoX3zCBoFYGwZBlI/IWHh7MLmZmZ4paELnYa79rIHMeT\nA2k0GqFK9DtD8mDyx1SeZdKYuEij0ZDInE6FI9nZZsUujjh4ngm9Pi8jhJCeGSX/mlbyr2nMmjyJ\n/BOpEP7t4IpOhSNbfYM7S6PRSHy2ZYUHUklJidhFoI8b5N2Fx8ZCnn8SbBixwsLCSNjvAxw8tp3k\n+EywGThXyU7JKsH/pisk/nFsqbpQTqWlReGBpFRS+1SVbBRZwq9lbTkYsCDk3HdgxWN7mGU5qMGT\nSfBMlUUacZBJLD4DFmyHOQC4lSxbSCaTyWg0Njc3E0IMBgPxmGEOUksjeUURwShw58dwc5kko/8y\nyJcsA+nYsWMLFy5klyMiIgghRUVFis8kqX2Myi6NWB7bcefKz4nY7jvZ/a9lSmrvdCHJMpDGjBkz\nZswYsUshKEl9hrrpzhGC8aj7+NH6WSubSUS2/3SQBcqBVFNTc+vWrU2bNt29e/e5554bMWJEz549\nVSpcqaJAIh+givlUUnwmWUYRldHbGOYA7kYtkI4cOfLXv/61urqaENKxY0eVSrVhw4aPP/6YENKv\nX7+tW7cGBQXROpankc5Hp2LSSMEdd26d6QeZ5G7SebOLgkIgabXa+Ph4hmFmzpw5btw4y+DR6/W3\nbt1asWLFkCFDhg0blpaW5vrhPI1ETlDFRBFHeR13wkw6h0wC96ETSF9++WVUVJTtUz4+PqGhobt2\n7TIajZ9++qnrx/I0EvnEVF4asZSRScJPfopMchMFnI0uohBITzzxRKvbeHl5cePigCcp9CnJffwC\nH/L9FBBxHm4MBwd3EGKUnclkwrgGZ0nhJrBKbRhZkuPFJEndDwKZBBTRzAmtVvvaa69Zrbxw4UK/\nfv0oHsUTII0EJotMYudWkNT9INjbi2AqByrk21KniGYLKSgoqLCwMCIi4tChQ7179yaErFmz5ssv\nv5w+fTrFo3gI0dPIQ6KIyKGRJKkmkS38bBZoodxlV1BQsGTJkhdffPG9997bt2/fjRs3/va3vz3z\nzDN0j6Js4n44elTDiCPN0Q0SzyFLyCQXSfD0EwX9a0irV68ODQ39+OOP/f39i4uLvby8qB9CwUTs\nrPOE8QsOSCqTZBRFHGQSuI5+IH344YcZGRmvv/56RkbG4MGDs7Ky/Pz8qB9FkURPI3yaiEuOOWQJ\nmdQ20vkmJDrKg9/i4uIyMjL27NmzdOnSoqKiJ554IiYm5siRI3SPomBIIxGJdTFJgqMV2gxjHMAV\nNFtIWq22Q4cOV65c4QZ579ixIyMjY8mSJS+99BLFAymSWN+S8JXWipD/CLk3iexCO8kpaB5ZojzK\n7uDBg1Yrp02bNnKkct5sbiLKSYmGkS1hGkmKzCFLyCRoGyF+GBsSgia8I6J0EyGNWuK+0Q2KzyFL\nyCRoAwrXkMrKyu7evdvqZkePHnX9WMojykAGpJHAFHOJyCm4ntQq9NdZodBC0ul08fHxYWFha9eu\n7dOnj9WdW+/evbt79+6dO3f6+/t72l31eBLyjEQU8UGrkeRRTSK70E4Cp9CZXPXKlSupqakTJ040\nmUz+/v5+fn7t2rUzGAzV1dUGgyEwMHDDhg1Dhw51/VgKI/D3I6QRf65kEnLIEjKpJWge2aJzDUml\nUqWkpKSkpJSXl3/zzTenTp0ym81+fn7Lly+Pjo7u2rUrlaMojJCnI6KobZz6HyGHHEAmAR+UBzWE\nhoYuXbp06dKldHerPEgj6eM/4g5R5Bjun2QLzSO7hBhlByJCGrkPcog/3D8J+EAgiUCYL0eIItfZ\nvZKEHHIFMgkcQCAJTcg0YtbkuftAisdlEpdDBFHUVui7A8ekGEilpaUZGRk6nS4uLi4uLs52g/Pn\nz5eXl3MPY2JiQkNDBSueKwRuG5WVoZPaVe0W5dwdxv5cbDtyyHXIJIILSC2jfMfYBQsWWK2sq6ub\nMWMG/52UlJRMmDAhODi4f//+K1eu3LVrl+02Bw8e3Lp16z//UF1d7VK5hSLMjAzsW92T3+1UcBOe\nmteOZD87kEa0cNeTxC4ISA7lFtI///lPnitbsnbt2ilTpiQmJhJCHnzwwaSkpGnTptneVCk2NjY1\nNdWVogpMgBkZcNHIdS31y0nqbkkK4MmTOOBEcsDtXXZ6vZ6b/JuPvLy8qVOnsstDhw7V6/X5+fm2\nP6rV6/Vnzpzp0qVLv379qJXVzdx3FiKKXITrQ8LDD2bBFp1ASk5OLigoIIRUVVVZze1969atp556\niud+dDpdc3Nzz5492YcqlcrX17eurs52y6ysrIqKiqKiouDg4LS0tLAwMW/73Sq3fidCGrWZUzmE\nRhJ1HphJOIUcoxNInTp16ty5MyGkoqKCXWCpVKp58+a9/PLLPPdjNpsJId26dePWMAxjNBqtNktK\nSmL76wwGw/z58+fOndvSzK3h4eHsQmZmJt8XQ1unwpGEEDeNLzAkDyaEMGvy7O5fo9G446DyYlsJ\nvT7/T13dmPv7Vxk+/6BOhFQfDquNzGl1S6mR7JnArMm7OjFEmOGgUqiETm77KOBDo9FI/Ls7nUBa\ntmwZIUSr1a5YsWLz5s1t3g87MWtxcXFMTAy7prGxUa1WW20WFBTEbZ+YmDh+/HidTme7GSGkpKSk\nzYWhgrt01IX2nnk2jCR+/gmDrQQK/XJhZdWHw2RapZIt9lVCDMmDhWkniV4J1YXil0HKKN+gz5U0\nIoQwDBMSElJZ+fupqdVqdTrdo48+6uBPmpqaCCHe3lIcv+6+gQzopuOv1+dlhPz+L3D9+hA67qjz\nnIHgOHNaRf9z3Gg0spN8W67kf4++hISE7du3P//88z4+Plu2bImKimK/UOzZs6exsZEdQZ6fn//M\nM88QQmpqajZu3BgZGWl1zwspcGsaKf6t6zqMU5ARz8kkcIxyIK1bt27r1q1WK1Uq1ZUrV3juYfbs\n2SUlJQMGDPDz8/P3909LS2PXFxQUNDQ0sIGUnJxcW1urVqvr6+ujo6M3bdpE8SVQRD2N0DBqlW0O\nUe+yRyPJHRQ/EBznDB80A+n27dtbt25dvnz5mDFjfHx82rYThmHsBszq1au55bw8qc+I446TD2nk\nANpDCuCBg+7ACs1AqqqqCg4OnjZtGsV9yhH1NEIUtUSseU7RSHITpWYSzhaeaAZShw4dbIdoexqk\nkbtJpDGETHITpWYS8EEzkEJDQx9++OHs7OxRo0ZR3K2M0J2tDlFkSSI5BMJQUibhiwt/NANJq9Vq\ntdrExMTAwMD27dtz61UqVVZWFsUDSRPdYXVII5bEcwifNe6AQXcei/Iou44dO/bt29dqZbt27ege\nRYIophGiSOIhxOF/j3NwFjLJM1H+YeyhQ4co7lBekEYukksOWcKVJPdRRibh9HCKWyY4KC0tbWho\nePrppwkh9+/fJ4T4+vq640ASQeWc88wokmMI2cKHjpso/sdJYIVyIGVnZ7O3MiKEnD17Nigo6Lvv\nvvvkk09Onz5N90ASQaWnjnvLeU4aKSOHWOi4cytZD7rDNxVn0Qwko9GYmJi4zFD19wAAIABJREFU\nfPnyadOmDR48mF0ZFxeXnJxsMBgkOLuPiyimkUzfb05RUghZQcedW8k6k8ApNAPpxo0b3bp1s/ph\nrFqtZhimrq6ua9euFI8lEa58BnlCw0jBIQRCkmMm4TtKG9AMJJPJZLcZZDQanbpprCy4eLYpuGFk\nGULEY3IIjSQA19EMpC5dulRUVFj1zmVnZ6tUqoCAAIoHEpeLPXWKjCLPDCEQkhwbSeAsysO+Bw0a\n1L9//48++ogQotFotmzZkpGRkZycTPEoUtC2NFJYFCGErKCR5G4yyiScCW1DeZTdzp07ly1btmTJ\nkubm5smTJ6vV6vfee++NN96gexQRte08U8zlIoSQY8gkd5NRJkEb0P8dUmpqampq6r179wghfn5+\n1Pcvljb31Mm9YYQQAqmReCbhS0mbuevO30qKIo6zJ5l8owij41yBzyO3UsYMDmAX/UCqq6u7e/du\nc3Oz5crHHnuM+oEE04a2keyiCM0gWvA7WQFIeQYHfB1xBeVAmjdv3vHjx61WOnULc6lxNo3kEkVI\nIPfBlSQB4GKSItEMpPLy8uPHj3/55ZcDBw5UxrwMTqWRxEcuWCUQQQi5GTLJ3SSYSfinu4hmIN2/\nfz8kJGTo0KEU9ymWNkSRpN4Y5D8J9J+XgAQSDDruBCO1TAJX0Ayk4ODguro6ijsUC880klqTyLYN\ndGNuWFgYPhZBsSQ1wAHNI9fRDKSuXbv++c9/XrJkybJlyzp06GD5lJeXF8UDuRWfNJJCFPHpgisr\nw9tDNLiSJAwpD3AAZ1Ee1PDMM8+8/vrrVrfpk8ughlajSMQcso0fgi44yUMmCUMKF5Pwj6aCZiDV\n1dW9/vrrw4YNS05OdvdsqqWlpRkZGTqdLi4uLi4uzvUdOk4jIaMI2QPQBqJnEriOZiBpNJrAwMC0\ntDSK+7SrpKRk0qRJ77zzTteuXVeuXFlZWTljxow2781BFLk7h+xmD0H8KAgaScKQ1MUkaDOageTr\n6yvMtaK1a9dOmTKFvTXtgw8+mJSUNG3atDYc2m4UWfZHUzm5W0odguDxDMgkYYh4MQn/X1poBlJo\naGhAQMCZM2fcPfI7Ly9v6tSp7PLQoUP1en1+fj7/g3LjcblzyOo8djaHHEQOC8EDIAApXEwCV9AM\nJK1WW19fP3PmzMDAwPbt23PrVSpVVlYWraPodLrm5uaePXtyO/f19eUz3NzydyFV6XpCSFV6iz1y\nrWaMJeQNgHQInEloHlFEeZRdx44d+/bta7WyXbt2FA9hNpsJId26dePWMAxjNBptt/x+TZNtCBFC\nwntmkJ7/vam9X/DwL5JkR1drNBqxiyA+SVRCZA45HFYb6cS3HLokUQmCYNbkGZIHX50YwqzJs3rK\nTZXQScKfAFY0Go3Ef5hI+QZ9VgO+3YGdlKi4uDgmJoZd09jYqFarbbcclPxASUkJ97DL2N8XzO4u\nopRI/PwThiQqIayMiPpVWhKVIIy/V16dGGL39VKvBLZ51IXuTj2YWwZnl5aW/vTTT+zy/fv379+/\nT3HnDMOEhIRUVv7eJNdqtTqd7tFHH6V4CACQNfxaVqYoB1J2dnZ4ePif/vSnV155RavVEkK+++67\n+Ph4ukdJSEjYvn27Xq8nhGzZsiUqKsqDvv2BbGGCO4G5O5Nw9Yg6moFkNBoTExOXL19eUlISGBjI\nroyLi7t9+7bBYKB4oNmzZ/fo0WPAgAGDBw/Oz89fs2YNxZ0DgAKw4xrQTpIXmteQbty40a1bt2nT\nplmuVKvVDMPU1dV17dqV1oEYhtm0aROtvQEIBr9JEpJbO+7wf3QHmi0kk8lk9zZIRqPR3TMJAcgI\nOu4Eg4tJ8kIzJ7p06VJRUWHVO5edna1SqQICAigeCEC+8LVaYO7IJDSP3IRmIAUFBQ0aNKh///5H\njx4lhGg0mg8//DAxMXHBggUUjwIgdxjdIDy0k2SBck/azp07x44du2TJkqqqqsmTJ+/fv/+9996b\nOXMm3aMAAPBHd+IGNI/ch/JMDYSQ1NTU1NTUe/fuEUL8/Pyo7x9AATC6QWC/d9xhmjtpo9xCun79\nOrvg5+fHplFNTQ23EgBALMyaPHTcSRzNQNJqtS+++KLVyubmZtuVAIArSXKEdq1bUQ6k7t27W60M\nCgry9va+e/cuxQMBKAMySWAujrjDP8vd8PMgAPAgLmYSmkduRXnYd0VFhdWtibKzs5ubmylO0wCg\nMPjeLbC2ZRI66wRAOZAGDhw4YMCAHTt23Lx58+bNmytWrEhMTJw1axbFowAoCT7jADiUu+x27dr1\n0ksvrVmzZvTo0aNHjz5w4EBSUtLixYvpHgVASXAlSXiYUkia6P8O6ZNPPvnkk0/q6uratWuH3yEB\n8IQeIYGxmcTzN7P47wjDXYMa/P39kUYAPOHDTixoJ0kK/RbSRx99dOzYMaPRyK1RqVRnzpyhfiAA\nJcHcDcJjG0mttpPwfxEM5RbSs88+++233z788MNVVVU9e/Zs3759VVVVr1696B4FAICKVrvskEZC\nohlIFRUVt2/f/vHHHz/77LOAgIDdu3dnZ2cvX75cp9NRPAqAUmF0g1jQcScRNAPpt99+Y2dqUKlU\n3F2Rpk2bdvny5fv371M8EICCIZME5uBm52geCYxmID3wwAPsQvv27evr600mE/uwubm5oaGB4oEA\nlAoff6Kge38KaDOagRQcHFxRUUEIUavVXbp0+d///d+KiooVK1YQQoKCgigeCEDB0HEnCttfJuG/\nIDyao+z8/f2/+uordnnv3r0vv/zyN998wzBMWloaxaMAALiJ1Yg7NFgFRnnYd0xMDLvw8MMPFxQU\n6PV6Hx8fuocAUDwMAReFZSMJ9S8K9872jTQCABnBlELiotBCKi4unjt3roMN2rVrl52d7fqBADwH\nGkli6fP3StS8WOh02QUGBjp4VqVyrh1WWlqakZGh0+ni4uLi4uJsNzh//nx5eTn3MCYmJjQ01KlD\nAEgfMgk8DYVAeuKJJ/bu3ev6flglJSWTJk165513unbtunLlysrKyhkzZlhtc/DgwXPnzkVFRbEP\ne/XqhUACACrYLwH8510FiujPZeeitWvXTpkyJTExkRDy4IMPJiUlTZs2zcvLy2qz2NjY1NRUMQoI\nIBw0ksTi1FzgQAv9QQ1arTY5OXnMmDEvvfTSunXrnJ03KC8vb9CgQezy0KFD9Xp9fn6+7WZ6vf7M\nmTOXL1+mUGIAacMPYgRjFf8Y4CAwyi2kjRs3fv755x07duzcuTMhZMeOHVu3bk1LSxs2bBifP9fp\ndM3NzT179mQfqlQqX19fq3uis7KysioqKoqKioKDg9PS0sLC7L9jw8PD2YXMzMy2vB6Z02g0YhdB\nfLKvhMicToUjy8pcaiTJvhJoaLUSOhWOJIRwVc2syTMkD3ax5iVFo9G09FEpFWZ6qqqq+vTpc+jQ\nIcuVH330UZ8+fVr6E6PR2PQHs9nc0NDQp0+fhoYGboPY2FirHZrN5jt37rALTU1NiYmJL774ot2d\nOziuh7hx44bYRRCfMirh7qFQV/5cGZXgolYrwW4ll0x4yD3FEYH0TwOaXXa3b9/u3r372LFjLVcu\nXbrU29tbq9Xa/ZMTJ05E/cFgMDAMQwgpLi7mNmhsbFSr1VZ/xU1ExDBMYmLitWvXMKE4KB467sSC\njjvB0AwktVptNptt15tMpo4dO9r9k/j4+KI/MAzDMExISEhl5e8XErVarU6ne/TRRx0ctKmpiRDi\n7S250RkAFGFcg7u1NHgE4xqERDOQwsLC2rdvv3//fsuVK1eu/NOf/mTbymlJQkLC9u3b9Xo9IWTL\nli1RUVFsp+eePXt27drFbsMNc6ipqdm4cWNkZCTbtAJQNjSS3MTxUEZM3yAYmg0LrVbb1NT0/vvv\n/+Uvf2EHNdy6dctkMvXo0WPUqFHsNkeOHPH19XWwk9mzZ5eUlAwYMMDPz8/f35+bmLWgoKChoYH9\nTVJycnJtba1ara6vr4+Ojt60aRPFVwEgTZgFXEQYBS4Myj1dHTt27Nu3L/fQz8/P2T0wDGM3YFav\nXs0t5+Xlta14ALKGnyW5A6pUOmgGUlBQ0KFDhyjuEACsIJPEgkaSANw72zd3I3MAAAlyKt1xMcnd\naAaSyWQaN24cN0ZuypQpERERERER586do3gUAA+Hi0m0tK0akUnuQzOQKioqKisrQ0JCCCEnT54s\nKChIS0ubNWvW7NmzKR4FAAhG3FHibOcnuuzcimYg3bt3jxvFsHr16qioqGHDhiUlJTU0NNTU1FA8\nEICHwzUk17X5Uhw67tyH8jUk7oext2/fTklJYZd9fHyam5vpHgjAw6HjzhWuVx0yyR1oBlJQUFBl\nZeWlS5cOHDjQ0NDA3qPIYDDo9fqWZmoAAFcgk9rMlVYmOu7chPKw76lTp06aNIkQ8vrrr7M/gM3N\nze3QoQP/mRoAgCc0ktqGSqVhFLg7UP5hbEpKSnJyskql8vHxYddER0dnZWXRPQoAsPCzpLZBjUkT\n/d8hqdVqLo0IIV27du3atSv1owAAC+0kp1DMb4xuoI5CC6m4uHju3Lk7duzo0KHDK6+8YrtBu3bt\nsrOzXT8QAIArqLcm0XFHF50uu8DAQKsFSyqVe+eDAPBw6LgTETKJIgqB9MQTT+zdu5dd5hYAQEjI\npFZ1KhzpvvpBJlGBtguAcuBiUkvcWjOIIlrodNndvXvX8QYY1wDgbmwjCe2kltRG5nRx5/7RSHId\nnUENCQkJDjZQqVRXrlxx/UAA4BhG3NnFhnRtmRtzGiPuqKDQZffEE0/k/2Hu3Ln+/v7ffvvt2bNn\nz549++mnnxJCtmzZ4vpRAIAnZJIokEmuo3MNif2xEcMw27Zt+/HHH/v16xcUFBQUFBQfH3/27Nk5\nc+ZQOQoAtIrtr0MmcYTsw0QmuYjmoIZbt27ZDvsOCgry9vbGbN8AgkEmcXBFTV5oBpK3t/etW7es\nVt68eVOn03l5eVE8EAA4hk9hIlIaoZHkCpqB1Lt3786dO/fv3//48eMVFRUVFRXbtm0bPXr0oEGD\n/P39KR4IAPjw5EaSiK8dmdRmlCdX/f77799444158+b9vndv76lTp3I3RgIAwWAUuLgvHKPA24By\nIBFCduzYQX2fANAGXCaRyByxyyIo0WMYjaS2ketMDSaT6fz58wcPHty/f7/YZQGQLg9sHkmnoxKZ\n5Cy5BlJKSsrbb7+9e/fuFStWiF0WAEnrMrasU+FIsUshEDaNpBDDbH8dMskpcg2kFStWFBQU4BdO\nAHzURuZIp93gblJIIxauITlLroHEMIzYRQCQE0+YVUj0S0e2cDHJKXINJABwlrIzSYJpBM6iP8rO\nHUwmk9FoZJedahuFh4ezC5mZmfSLJXkajUbsIogPlUAsKqETIdWHw2oVN+iuU+HI2sgcx9OninUm\nMGvyrk4MYdbkiXJ0SxqNJixM0t9I5BFIJ06cWLx4Mbt88eJF/plUUlLitkLJg8TPP2GgEghXCWFl\n1YfD3HqrOuGxzT4+/2XRzgTcVZYfeXTZxcfHF/0BV48AXKGwme6kM6yuVbiY1Cp5BJItk8lkMBia\nm5sJIQaDwWAwiF0iANmQxcc3HzJKIzSP+JBrIB07diwiIiIxMdFgMERERERERCCTAPhTzAAHWaQR\nCyPuWiWPa0i2xowZM2bMGLFLASBjXCbJ6DPdEobVKY9cW0gA4DqZfqDLd9JYNJIcQyABeLT/TMAq\nE7Ju1RFkkkMIJABPJ7txd/JNIxYyqSUIJACQRybJt6cOeEIgAQAhks8kuffUWUEjyS4EEgD8rsvY\nMmleUlJYGnGQSVYQSADwX6TWVFJqGuGnsrYQSABgTTqZxF40Ul4asdBxZ0WuP4wFALfiMkmsMFBq\nwwgcQCABQItEmc3Bo6KoDyYCt4AuOwBwROCRDh6VRix03HEQSADQOq4Hz32xxP3MyKPSiINMIggk\nAOCJiwrqscTt0DOjiGDE3R9wDQkAnGCVSa5ECJdqHptDlnAxiSCQAKANuAixbCrxyRWrphWiyIqH\nZxICCQDazjJRePbjIYRagtENCCQAoANJ4zoP77jDoAYAAAnx5HYSAgkAACQBgQQAIC0e20hCIAEA\nSJEHZhICCQBAcjxzXAMCCQBAijyw405yw75NJtOFCxc0Go3RaJwwYYLdbc6fP19eXs49jImJCQ0N\nFaZ4AADgJpILpJSUlGPHjvXu3bu4uLilQDp48OC5c+eioqLYh7169UIgAYDyeNrPkiQXSCtWrEhN\nTT116tScOXMcbBYbG5uamipYqQAAROFRmSS5a0gMw/DZTK/Xnzlz5vLly+4uDwCA6DzkYpLkWkg8\nZWVlVVRUFBUVBQcHp6WlhYUJceswAADhec7oBpEDyWQyGY1Gdpln24gQkpSUxPbXGQyG+fPnz507\n9+jRo3a3DA8PZxcyMzNdLqz8aDQasYsgPlQCQSUQQuRfCVcnhjBr8lzZg0ajkfh3d5ED6cSJE4sX\nL2aXL168yDOTgoKC2AWGYRITE8ePH6/T6dRqte2WJSUltIoqUxI//4SBSiCoBEKIrCvh75VXJ4bI\nuPz8iBxI8fHx8fHxruyhqamJEOLtLde+RwAAPjxhdIPkBjWYTCaDwdDc3EwIMRgMBoOBXb9nz55d\nu3axy/n5+exCTU3Nxo0bIyMj+Xf3AQDIlOIvJkmuYXHs2LGFCxeyyxEREYSQoqIihmEKCgoaGhpm\nzJhBCElOTq6trVWr1fX19dHR0Zs2bRKzxAAAQlF2O0lygTRmzJgxY8bYrl+9ejW3nJfn0pU9AACQ\nIMl12QEAgAMK7rhDIAEAyI8iMwmBBAAgM0q9hoRAAgCQH0V23CGQAABAEhBIAACypLxGEgIJAECu\nFJZJCCQAAHlTTCYhkAAAZExJI+4QSAAAsqeMRhICCQBA3hTTSEIgAQDInjJGNyCQAACUQAGZhEAC\nAFAOWWcSAgkAQCHkfjEJgQQAoCjybSQhkAAAlEPWjSQEEgCAosh3dAMCCQBAaWSaSQgkAACQBAQS\nAIACybGRhEACAFAm2WUSAgkAQMlklEkIJAAAxZLXKHBvsQvQRqWlpVlZWWVlZb6+vuPGjYuOjha7\nRAAAEnV1YogskkmuLaQpU6aUlZXFxsYyDDN9+vQDBw6IXSIAACmSRRSx5NpCys7O9vf3Z5c7deq0\nefPmhIQEcYsEACBN7OgGZk2e2AVphVxbSFwaEUK6devW3NwsYmEAACSuz98rDcmDxS5FK+TaQuIY\nDIb09PTx48fbfTY8PJxdyMzMFLBQUqHRaMQugvhQCQSVQAhBJRByO2lPH7HL4Jg8AslkMhmNRnaZ\nYRjLpxYtWhQYGJiYmGj3D0tKStxeOGkLCwsTuwjiQyUQVAIhBJUgefLosjtx4kTUHwwGA7d+8eLF\nd+7c+eKLL7y8vEQsHgAAuE4eLaT4+Pj4+HirlUuWLLl27Vp6erqvr68opQIAAIrkEUi2li1bVlhY\nmJ6erlar2TaTVVceAADIi1wDad++fYSQIUOGsA8ZhikqKhK1RAAA4BK5BhJGKwAAKIw8BjUAAIDi\nIZAAAEASEEgAACAJCCQAAJAEBBIAAEgCAgkAACQBgQQAAJKAQAIAAElAIAEAgCQgkAAAQBIQSAAA\nIAkIJAAAkAQEEgAASAICCQAAJAGBBAAAkoBAAgAASUAgAQCAJCCQAABAEhBIAAAgCQgkAACQBAQS\nAABIAgIJAAAkAYEEAACS4C12AayVlpZmZWWVlZX5+vqOGzcuOjradpvz58+Xl5dzD2NiYkJDQwUr\nIQAAuIPkWkhTpkwpKyuLjY1lGGb69OkHDhyw3ebgwYNbt2795x+qq6uFL6csxMfHi10E8aESCCqB\nEIJKkEMNSK6FlJ2d7e/vzy536tRp8+bNCQkJtpvFxsampqYKWzQAAHAjybWQuDQihHTr1q25udnu\nZnq9/syZM5cvXxaqXAAA4F6SayFxDAZDenr6+PHj7T6blZVVUVFRVFQUHByclpYWFhZmd7Pw8HB3\nllEGUAMElUAIQSUQQlAJktfObDaLeHiTyWQ0GtllhmEsn5o3b15tbe2OHTu8vLys/kqr1QYFBRFC\nDAbD/Pnzy8vLjx49KkyBAQDATUTusjtx4kTUHwwGA7d+8eLFd+7c+eKLL2zTiBDCphEhhGGYxMTE\na9eu6XQ6gUoMAADuIXKXXXx8vO3AjyVLlly7di09Pd3X17fVPTQ1NRFCvL2l2/cIAAB8SG5Qw7Jl\nywoLC7dt26ZWqw0GA9ds2rNnz65du9jl/Px8dqGmpmbjxo2RkZFW3X0AACA7Il9DsmV11ZFhmKKi\nIkLIkiVLGhoaNm3aRAgZPHhwbW2tWq2ur6+Pjo5et25dcHCwOMUFAABKJBdIAADgmSTXZQcAAJ4J\ngQQAAJKgzMFppaWlGRkZOp0uLi4uLi5O7OIIodWXzGfWWrnj/3+/cOHCjRs3hg0bxv2EQBn41IDR\naNy7d+/FixcfeOCBESNGjBo1SuBCuhufSjh58uTx48eNRmNERMTkyZN9fHwELqTATCbThQsXNBqN\n0WicMGGC2MVpkQJbSCUlJRMmTAgODu7fv//KlSu5sXkKxucl85m1Vtb4/9+1Wm1ycvL/+3//71//\n+peQJXQ3PjVgMBimTp26f//+J598smfPnocPHxa+nG7FpxK2bt36/vvvR0REDB069Ntvv33zzTeF\nL6fAUlJS3n777d27d69YsULssjhkVpxZs2Z99NFH7HJubu5TTz3V3NwsbpHcjc9Lrq2t5ZY/++yz\n0aNHC1c+QfD/v8+aNevAgQN9+vT58ccfBSyg2/GpgU2bNiUkJBiNRsFLJxA+lTBixIiMjAx2+dq1\na3369GloaBC0lIJramoym825ubn9+vUTuyyOKLCFlJeXN2jQIHZ56NCher2e+92SUvF5yTxnrZUv\nnv93tk0wZswYQQsnCD41sH///unTp2u12jNnztTU1AheRrfjUwkPPfTQ/fv32eXGxkZvb2/Fd9nJ\n5ZeaSruGpNPpmpube/bsyT5UqVS+vr51dXXilsqtnH3JjmetlSmelXD37t0NGzZ8/fXXghfQ7fjU\ngNFo1Gg0J06cWL9+fa9evX788ccFCxbMnDlTjPK6Bc/T4IMPPvif//mf69evMwxTWFi4atUqu1OU\ngfCUFkhms5kQ0q1bN24NwzDc/K2K5OxLXrRoUWBgYGJiohCFEwrPSli5cuWbb74ZHBxsOXGiMvCp\nAZPJRAi5detWdnY2wzDnz5+fOnXqiBEjevfuLXBp3YTnafDvf/+7traWENKhQ4fGxsbKykohCwkO\nKK3Ljm2ZFhcXc2saGxvVarV4JXI7p16y41lr5YtPJZw7d+78+fM9evQ4derU6dOnCSEXL14sLS0V\nuKhuwqcGvLy8VCrVxIkT2Y1jYmL8/f2VdFMxPpVgMpmSkpLmzJmzatWqpUuX/t///d+GDRuUVAmy\nprQWEsMwISEh3FcerVar0+keffRRcUvlVvxfslOz1soLn0rw8vLq16/fV199Rf74Kp2dne3r6/vY\nY48JX2Dq+NSASqV69NFHLVsMbJtJMfhUgl6vb2hoeOihh9iHQUFBPj4+N2/e7Nevn9DFBRtKayER\nQhISErZv367X6wkhW7ZsiYqKaun2fYrR0ku2nJG2pVlrFaPVSoiJiUn7w+bNmwkhixcvnjp1qrjF\npojPaZCQkLB37172kv7Jkyfv378fFRUlYpmpa7US1Gp1cHDwiRMn2O1PnTql0+mU8aXEAZPJZDAY\n2KFMUn77K62FRAiZPXt2SUnJgAED/Pz8/P3909LSxC6R27X0kgsKChoaGmbMmEEI2bdvHyFkyJAh\n7FPcrLWKwacSlI1PDbzxxhtXr14dOHBg586d6+vrP/nkk4cffljUUlPGpxI2bNiwaNGiAwcOdO7c\nuaqqKiUlRTFX0Vpy7NixhQsXsssRERGEkKKiIgkOvVPs5Kp1dXW1tbUKe7M55oEv2RYqgU8NGAyG\n8vLy3r17q1QK7CMh/CpBq9XW1dWFhYUptRLkSLGBBAAA8oKvBgAAIAkIJAAAkAQEEgAASAICCQAA\nJAGBBAAAkoBAAgAASUAggQe5dOnS+fPnhTziTz/9dOHCBf7bX7x4saSkxNmjWL4uJU3QB54GgQQe\nZO/evW6duePgwYPsXHmc3bt3b9++neef19TUvPXWW/x/p3ny5En2Ft2Wr6uqqioxMVFhk9SBh0Ag\nAVCTn59/6tSpNv/59u3bn376af7zqp04ceKZZ56xWhkXF+fl5bVnz542FwNALAgk8GjFxcVHjhzJ\nysqyvH1qaWkp22+Wn59/5MgRq3sTGAyGrKyso0ePlpWV3bt3Lz8/n52q8vr163fu3Pntt9/y8/Pz\n8/PLysos/8ruriwZjcavv/56woQJ7MO6urr8/Py7d+9yG2i12vz8fJ1Ox63Jzc0dPXq07a4SEhJ2\n797tVD0ASIECJ1cF4EOv18+fP//06dOxsbHV1dW//vrr+vXrhw0bRgjZtm3bb7/9ptfra2trVSrV\n5cuXFy9ePGvWLEKIVqudMmVKfX19ZGTkqlWrhg0btm/fvh9++CEgIODbb79lLxe99957hJCJEycm\nJSURQgwGw4wZM2x3ZSU3N/f+/fujRo1iH/r7+2/evLmpqembb75RqVR6vf61117r2bMn1yQ6d+5c\nQ0PD4MGDbXc1dOjQdevWlZWVKX6ee1AYBBJ4qHXr1l27di03NzcoKIgQsnnz5oULF545c4a9WdSp\nU6fWrVs3ZswY9qkNGza88cYbXl5eH3zwga+v76FDh3x9fXU63Z///Gduh8nJyVqt9rfffrO6TNXS\nrqzK8/PPP/fo0cNyAub169ePHTt29erVS5cu/fDDD+/fv79q1Sru2ZMnT44aNcruBae+ffuqVKqi\noiIEEsgLuuzAE5lMpq+++mrmzJlsGhFC3nnnnfv37+fl5bEP+/Tpw0btKzewAAADYElEQVQIIWTc\nuHHNzc0XL140Go05OTlvvvkmG1pqtfqdd95p9Vh2d2W7WWVlpdWt5IKCgtavX/+3v/1t5cqV+/fv\nX79+vb+/P/dsbm7uyJEj7R5RpVJ17Njxl19+abVsAJKCFhJ4opqaGoPB8PXXXx87doxbqVKpGhoa\n2OVHHnmEWx8QEEAIaWpqqq6uNplM7ENW586dWz2W3V3ZbtbU1GTb3Bk4cOBrr722c+fOhQsXPv30\n09z669ev/+tf/xo+fHhLB/Xy8qqvr2+1bACSgkACz/Xcc89Zfsq/9dZbjke4dejQgfz3bb8p3nmz\nffv2VVVVVivr6upOnDihVqtPnTr19ttvc+tPnjzZv39/Pz+/lvbW2NjoyTeFAplClx14oq5duwYE\nBBiNxmf+G9eDZ5evr29AQMClS5e4NVbdYq7c6u3JJ5+0/Qnt0qVLfX19Dxw4cOPGjc8++4xbf+LE\niRdeeKGlXRkMBp1Op/i7oILyIJDAQ82ePXvr1q1HjhxhWzx1dXV79uzR6/WO/2rmzJnbtm3bv39/\nXV3dwYMHd+7cafns448/funSpcuXL9+/f7/VXVmJjY3V6XTXr1/n1uzatSs/P//zzz8PCwv7y1/+\nsmnTpvz8fEKIVqv9+eefR4wY0dKucnNzvb29udvVA8gFAgk81IwZM+bPn798+fLIyMinn356wIAB\n+/bta/WvZs6c+c4776xZs2bQoEGHDx9OTk4mhPj4+LDPJiQkREVFvfLKK1FRUR9//LFT5Xnsscce\nf/zx7777jn1YXFy8evXqlStXsiPlRo0aNX369EWLFt29ezcnJ6dv374hISEt7SonJ+fFF1+0HLAH\nIAu4hTl4NJPJdP369fr6+oiIiDZ8gu/fvz8lJaWwsNCVzjpOZmbm6tWrc3JyHG82e/bsxx9/fN68\neXafraurGzJkyL59+8LDw10vEoCQ0EICj6ZSqR577LGoqCieaVRWVpadnW0wGEwm0/nz5zds2DBx\n4kQqaUQIiY+PDwwM3L9/v+PNXn755VdeeaWlZ3ft2vX8888jjUCO0EICcML169dnzZpVUVFBCPH2\n9h47duzy5cvZnyUBgIsQSABOMxqNer2+ffv2tNpGAEAQSAAAIBH4fgcAAJKAQAIAAElAIAEAgCT8\nfwFB/N4/IFV8AAAAAElFTkSuQmCC\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "% Plot first 3 mode shapes\n", "beta_1 = 1.87510407;\n", "beta_2 = 4.69409113;\n", "beta_3 = 7.85475744;\n", "sigma_1 = 0.7341;\n", "sigma_2 = 1.0185;\n", "sigma_3 = 0.9992;\n", "close all\n", "figure\n", "hold on\n", "x = [0:1/1000:1];\n", "Xn_1 = cosh(beta_1.*x)-cos(beta_1.*x)-sigma_1.*(sinh(beta_1.*x)-sin(beta_1.*x));\n", "Xn_2 = cosh(beta_2.*x)-cos(beta_2.*x)-sigma_2.*(sinh(beta_2.*x)-sin(beta_2.*x));\n", "Xn_3 = cosh(beta_3.*x)-cos(beta_3.*x)-sigma_3.*(sinh(beta_3.*x)-sin(beta_3.*x));\n", "plot(x,Xn_1)\n", "plot(x,Xn_2)\n", "plot(x,Xn_3)\n", "title('First three mode shapes of a clamped-free beam')\n", "xlabel('length (x/l)')\n", "ylabel('displacement (x)')\n", "grid on\n", "hold off" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "hide_input": false, "kernelspec": { "display_name": "MATLAB", "language": "matlab", "name": "imatlab" }, "language_info": { "codemirror_mode": "octave", "file_extension": ".m", "mimetype": "text/x-matlab", "name": "matlab", "nbconvert_exporter": "imatlab._exporter.MatlabExporter", "pygments_lexer": "matlab", "version": "9.2.0.556344 (R2017a)" }, "toc": { "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "toc_cell": false, "toc_position": {}, "toc_section_display": "block", "toc_window_display": false }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }