{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Time Series Recitation Two Code" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Generating iid data:\n", "\n", "$$X_t \\sim N(0, 4)$$\n", "\n", "$$ (\\sum X_t)/T \\rightarrow E[X_t] $$" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXmYHFW5/79vV+/T07NPZrInJCQkkEAIMQiyKQhRQBQVr3oVF1wubtcN9F5Erz+vet19UMGNRWURRUBB1rAnQCJk3yfbZCaz79N7n98fVaf6VHX1TM3ak+T9PE+edNdUd51a+nzPu5z3kBACDMMwDOMGT7EbwDAMwxw7sGgwDMMwrmHRYBiGYVzDosEwDMO4hkWDYRiGcQ2LBsMwDOMaFg2GYRjGNSwaDMMwjGtYNBiGYRjXeIvdgPGmurpazJ07t9jNYBiGOabYuHFjuxCiZrj9jjvRmDt3LjZs2FDsZjAMwxxTENFBN/uxe4phGIZxDYsGwzAM4xoWDYZhGMY1LBoMwzCMa1g0GIZhGNewaDAMwzCuYdFgGIZhXMOi4ZK1O1txpDtW7GYwDMMUFRYNl1x7+6u49MfPFbsZDMMwRYVFYwT0JdLFbgLDMExRYdFgGIZhXMOiMUJePdCJJo5tMAxzgsKiMUI+9YeN+PXzDcVuBsMwTFFg0XBBNivM133xNOKpbBFbwzAMUzxYNFyQFTnRSKSzSGdYNBiGOTFh0XCBYmgAADL2DQzDMCcILBouUC0NAEixaDAMc4LCouECu2iwe4phmBMVFg0X2A2LNFsaDMOcoLBouIAtDYZhGB0WDRdkbZYFWxoMw5yosGi4IM89lWHRYBjmxIRFwwV29xSn3DIMc6LCouECu3sqleWYBsMwJyYsGi5g9xTDMIwOi4YL8rKn2D3FMMwJCouGCzjllmEYRodFwwX2EAYHwhmGOVFh0TDoHkziPb9ah39ubc77W37tKbY0GIY5MWHRUHjlQCeae+J52/NSbjkQzjDMCQqLhoFP0y9FyiFewVVuGYZhdFg0DLwaAQBSDlZEfsotu6cYhjkxYdEw8Hn0S5FMD29pcMotwzAnKiwaBh4PweshpB2C3PZNPLmPYZgTFRYNBZ/mKeCe4tpTDMMwAIuGBZ9GrtxTnHLLMMyJCouGgt/rKZA9ZX0vRH4RQ4ZhmBOBoooGEV1KRLuIaC8R3eDw9w8TURsRvW78+9hEtkd3T+WLhpM7iq0NhmFORLzFOjARaQBuAXAxgEYArxLRQ0KI7bZd7xVCXD8ZbfJq5BjTECJ/G8c1GIY5ESmmpbEKwF4hRIMQIgngHgBXFrE98GkeJF24pwDn+RwMwzDHO8UUjRkADivvG41tdt5FRJuJ6H4imjWRDfJrHseJe/ZAOMAT/BiGOTEppmiQwzZ77/wwgLlCiGUAngRwh+MXEV1HRBuIaENbW9uoG1Qw5dbB1GD3FMMwJyLFFI1GAKrlMBNAk7qDEKJDCJEw3v4awJlOXySEuE0IsVIIsbKmpmbUDfJpNGz2lF/WqGLRYBjmBKSYovEqgIVENI+I/ACuAfCQugMR1StvrwCwYyIb5NM8w87TmFsdBsCVbhmGOTEpmmgIIdIArgfwGHQxuE8IsY2IvkVEVxi7fZaIthHRJgCfBfDhiWxToXkaGUM07r1uNT51wUkAOOWWYZgTk6Kl3AKAEOIRAI/Ytt2kvL4RwI2T1R6vZ+iU26BPg1/TAHBMg2GYExOeEa5QaHKfNCo8REoJdbY0GIY58WDRUPAVLCOiWxVEujUC6JVu7331ENY3dExqGxmGYYpJUd1TUw3/MFVudUtD19l0VuCrf9kCADjw3bdNXiMZhmGKCFsaCsOl3GrGmhsA0BtLTWbTGIZhpgQsGgoFYxqmpZFzTx3oGJjUtjEMw0wFWDQUCs3TkJlSpATCD3YMAgCqI/7JayDDMEyRYdFQ8BWscqv/r7un9EsmLY1p0SAA4M51B/D3zU15ny1E10ASDW39Y2swwzDMJMOioeDKPWWzNMrDPgDAH9YfxF82Nro+1s+f3osP//7VsTZ5SLoHkzj9W49j48GuCT0OwzAnDiwaCj7Ng3RW5K2fIQPhHsq3NKTrKp0RSDi4tgrRE0uhezA5Dq0uTGtfAt2DKRzk+AvDMOMEi4aC32sUI7S5qLJZZZ6GYWlIXZGikcxkEU9lXB8rnc2OSGRGg4zPpHn2OsMw4wSLhoKvwGxv6Z5SU24lUmBSmSziKfcikMroouG0KuB4Ic+DS54wDDNesGgo+GTZ8zzR0P9XJ/dJVPdUPO3e0siJzcR16NLCYEuDYZjxgkVDQQqCfcnXjEMZEYkUmGQmi8QILQ0ASChC0z2YHNegdcpwT2W4ThbDMOMEi4aC33RPWUfm0oWkUb57SloaurtpBDEN4xhqXOOudQfxvl+vd1wpcDRI8WNLwz13rTuAhze5T51mmBMNFg0F6Z6yr/8tO3E1ewoASvxazgWUESOKaSRNSyP3mf5kGsn0+AXIpTCNRjS++fC2Ec07OV747we34TN3v1bsZjDMlIVFQ8FdTCNnadRGg0hns8hmBdJZMbLsKenWUgQildYPNJhMj6r9dsYSCH/w9Sas3Tn69dYZhjk+YdFQkKKRTNvnaRgxDQ8solETCSCdEeYqfumsyLNSCpEy3VMZZZv+2dgIxGcoTPfUKILtiVRmRIH94wG3945hTmRYNBT83mFSbm3uqZpowBCKXKccd+laMgPhikvLFI3k+HTWsl2ZUSxNm8xkER+ndhwr9HDlYoYZFhYNBTfuKU0JhFeEfUhnspb9Ey6thJRDTENaH4Pj1FmnRhkIz2YFUiNMIT4e6JrgGfoMczzAoqEgrYi8lFtlRriKLDui7u/W0pAduSWmMc7uqVRWWhojE43kOFs8djJZgf7E+MRtxpPOAbY0GGY4WDQUcu6pAim3tnRbr4eQztjcU24tjXT+PI3xdk+lRllGRFo/I8kGGwn3vHoI531/7ZSbqd45wJYGwwwHi4ZCwZRbxT2l4tU8yGSFxT3lWjSy+fM0pFiNm6VhBsJH1vlL66eQeyqbFWPq8Pe1DqBzIDluWWLjBbunGGZ4WDQUCsc0cqXRVbweQiqbtYnGCAPhDpbGeMU0RltGRLqnCgXC3/T9tXjLj54ddbtkdd/xEsfxgi0NhhkeFg0FM+W2YJXbnGosn1kGr8cDIaxC4ToQnnaYp2G6p8ZnBC6/e6RWgTyHQvGZI90x7G8ffbl1OaKfqJjJaOlSRGMy02/jqQx++cw+TvlljglYNBT80tJI57un1HjGppsuwb2fONucs6G6pNzO5nZ2T41zIHyU2VNuA+EjLXeybl8H/vtvW9E1qAecx8uicsMdLx3A7pa+IffpVNxTE122XuUnT+7B9/65Ew++XngGvhACv39x/4SvwXI8c/ND2/Aff/pXsZtxzMOioeAbYp6G6poqC/sQ9GlmHSq1k3cd03CcpzG+Kbfp0WZPKTGNoUq3t/TFR/S97/v1ety1/iCae2IAJs89FU9l8I2HtuGudQeH3E+1NOKpDNbt68CrBzonunlo7dWv41B3qbErhm8+vB3/3Hp0wttzvLK9uRdbj/QUuxnHPCwaCjLl1i4aGSEsrimJtD7UEbmbuQ2ZrDAXcVLTdcd9RvgYs6eEyB9xqyJyuDM2ou8t8WsAgJbeBIDRu6d64yl89PZXcbhz0NX+7f368Rrah16TvXMwl3K7p7Uf7/v1erz7V+tG1caRIAcJ8vqo2ONcA1PMpXcs0R9Po7OfLbWxwqKh4C8Q0xBCnw1uR8ZArJbG8G6NQpMBxz3l1qw9NbrsKQB55d5VEXHbaUuiIZ/lvdvzFEJg19Gca+nlhk48tbMVN/51S95+tz67Dy29VguorU8XjX2tQ8dh+uI50fjmw9sBAKUBr6s2joUBI4blsWVatPTGsfQbj2HjwS7zGRtJfbOJ5qFNTTj//9ZaYjEv7W3Hqv/3JHrjU2/OS18ihb5E2vJ8MyOHRUNBVgjJWyM8K/Iyp4AClkaBH/VAIo0eYyRrEQ3lAZbzPca7jMhIak/1xVMYUCbe2a0e9W+HRigaZTbRGHTZAW482IW3/uS5PNfCzqO9lvdNPXH876M78eiWZst2KRpHe+MYSKQRT2UcXXaxZAbRoC4Su4zvjgQnXjSkFWHvzI50x5BMZ7G/fcB8rqZSmvKmw9042DGIvniuTdube9Hal8CBMSRKjAcbDnTmPS/9Rjs5LjQ2WDQU5DwMe4eSFflzNIDc8rAxF4Hwb/9jOz5256sArJMH1f2lq2q4zjSbFdjWNLxvdjSB8NNufhyf/MNG873srO546QD+ubXZMpP7cNcILY2gVTQefO0Izvv+2mFHfs09uuXQaBxP1ohq709a7lWvsT1ms45aDdEAgP3tA3jbz57Hr57dl3ecwWQGFSV+ALm5OWqHOFEUEg05eJBCp+5bTF7Z34mbH9pmXlf1Gsl709Sds/aEEGYca7L4yv2b8f3HdlnaIJ/dThaNMcGioSAth4zd0hAir4SIvr+Te8r5R93UHTdHvOkClsZQ7in1M1/88ya87WcvYG/r0D76lMM8DT2e4iwicgSmaow8t9+9uB/3bWi0iEbjCGMawhbqffVAJw51DlrcQk5IV0eHEajuVQoLbm/KWRs50bBevzZFNPa19eNQ5yAauwbxwd++jFsV8YilMigP+y2f7U+kJ3zmukyxtpevkQLRn0ibbs/JSFNu7olhzxCZZo9sacbtLx1AQ5v+/KmuKHkPjioisW5fB87+36dx57oDI2pHOpPFu375En78xO4RZeol0hkc6BhA50ACzT0xHO2JI5HOmoM1jmuMDRYNBWlN2PvUrBB5JUSAnKURTw4f0xhMps2HNmkRDTWmkXNPdQ0k8fg2PVPmty/sx4KvP4qewRR2t/ThgdeOAMCwLgBzuddsrsNZ8T9P4PHtLY77O7mbpAj2xdPojaVME7+mNIDm3pGJhn2U3Gt8l5wPkspksbmxO+9zvTF9vw7jx65Wo1WD2+b32UWjP4HSoBdEwK6jfUhlBGLJDLYc6cE2Q3QyWYFkOouKcM4akoFpVSh7Yqm8mMloUIVbBrdlfKujP4Erb3nRTBEeSKTN52QyLI3vPLIT1/+p8EJU8vxlnMnJ0mhWrpE8j5se3IbWEWTcdQ2msPFgF3761B7cu+Gw68/tbx9AVgBdAyl86c+b8JW/bLa08Xi0NH765B78/Kk9k3IsFg0FqQv57inh6J7yjsDSGEhkHNe3cLI0BlMZfPj3r+C6uzaiL57Cz4yHoTeewiv7cymgTcOY/Cnb8dr7E+iJpdDQ5iw2BzvyRSOW0tNu++Ip9MXTZtB2YW0ER3vi5gjwX4e6sHZn65DtKTRKltsf2dKMK295Ma9jkSNZOWNbFQ21M5CjXLvfv60vgRnlIUSDPhzp1q9ZPJVFLJkxYzTyHlYolsaMipDlewHgpge34g3feQov7m0f8lyH4sntLZh34yM4ZFxvef7y+djd0o9Nh7uxwUj3Vd1T45mmvPVIj8VSkzR2DZpWnRNSNKQFq1oaPaalEceO5l687WfPY49iEe8+OrR1rKJ+787m/HYCQGtfPO+52tOiH6N7MIkjXTEctlmzXS5m/h/sGMCfRyBUbkims7jxr5tHnEBSCCGEWdLnx0/uxg+f2D0u3zscLBoKpnvKJhqZLIZMuVVHf4XWCdctjdyI2tzfYT2NeDKDTY16zCKWypg/xHRWmDn9modwpGsY0bDN05Ad7GAyDSEEfvLkbsuENydLI5HKmqZ9XzyF/oR+fgtrI0hlhNm5vPMXL+Ha218dsj2FRsmyQ2ztTUAIayett1vGMHQ3U288hSoj9mARDWO/WDI/plFTGkAk4EWrke47mMogkc6iL5G7JgBQrlgaM8pDeceQnezXHrBmbo2ExwwLcu2uVgghTCGWMQ1p2ciYQX8iY3aMffEUPv3Hja5iWsPxpT9vwjcf3pa3vbU3YUl4sCNTpiWOlkZ3HE9ub8G2pl6LZTuSUb76vU09uYHE3zc34UD7ABLpDNb89AV8+f5Nls9JkRpIZnC0N46W3rjFWnRTzfjuVw7jy/dvHrIac+dAEj99cg8S6YxlMAforrXfPN9gGUTubunD3a8cxhMFLH0A+PlTexznknzk9lfz4nCfuGsjbvzrlkmfe8KioUCmeyq/yq3mcKXsgfBIwJuXoioZTGZMd5EaCLfO0zAm96VyD2pc6QCT6SxaevUOcFZFyBw1F8Je5VZ2vgOJDLoGU/jJk3vw9825TKNDDpZGPJUxO/HeeNp0Ty2YVgoAWPOz53Ht718x9x/K/y+vU1WJH5UluRF9zgWWMt5br6F0T3UqMY3aaBCahywjSLlfPJVBOpPFNbetw3O729BuiEZJQDOtGNm5yc5RXmcnS0M9ht+rPwgHOwZHHV+oLwsCAPa06q5G+bglTNHQjydjMQOJtOnC29c2gEe2HMWzu8e2FO9gMo3dLX2WJAFAT7Jo7YsjZlxDO0KIPEvQcg+M56O5N4YtRmfW1pfAnKowAN1iuPznL6BpmGcXyA0eSgNeM5CeyQp8/p7XcetzDXhiewva+xP4x5ZmM74CAHtbcwOheCqLwWQGRxXR6RywnrMT0hoZygX8wGtH8OMnd+Oa29bjPbeuw2uHusy//etQN779jx14cW874qkMvvvoTmw2BoKFfrfxVAY/fGK36X5WWd/QgQ0HuizbdrX0YeOhLrygWL32OWYTAYuGDc1DjoFwJ/eUtDRkp1ca9Bac3DeYzJii4DRPQ60cq46U1ZFZKpNFS18c06IBTC8PDfvDs68R3q+MquWPqF8ZzR3szP+BxFIZsyPoT6TN0fzC2ggAvUNYuyvXgR0dwt8fS2Zw3XnzsfbLF6AkkJvIJsVEHsdurZmBcCWmURbyojTotYxGZecVS2XQMZDE+oZObDzYhba+BGoihqVhdJIy6C9FQwp1hcXSCBvfmztG92AKAUM4hpssKHlyewsu/tGzpshIi+uPLx/Cf96XGyWbloZxPGnFDSRz7ilpbbXZOvtfPbvPMpelEDf8ZTNuWbsX25t6kRX539M5mDSf04FE/rPcNZjKWzrAydJo6UmYogEAi+v0Qcazu9uw5UgPXj+cH7uyI+/7yXWlaDaysdr6EkhnBfa19eO+DY2oLQ3Ap3nwh/WHzM/taenPi0HuVURFncRZCFkfrWEI0dhixN9eO6T/73QdeuMprN3Zil89uw+/e3E/ABT0EMh7Ie+xJJ7KYDCZydveOZDEoc5BPL+nzbJtomHRsKERwT5YLpxya8Q0kjnRcBp9ShdEMpOFEAJpIzBNlBtdppQJeKpPXvU5JwxLoy4axIxyF5aGLXtKPtS62R4ztuV+QM6WRtayjxSbk2oijscs5K9NZ7JIZrKIBLyIBn0I+3LzH+Q1kyPLfEvDmj2li4bPEA11lCvdUxlTYNr6E0hmsqgo8aMkkBOZbqPjkEIq21DuZGkkrNlBZ86pAAD84pl9+NHjubTOQvz2hf3Y09pvupRkh2JPuJDPgnSZqWIvr4n8TLuSAdTWl8B3H92Jj94xtHswkc7gr68dwRPbW8xRb38ibXne1CB/v8OcEKckANWd2BNLIeTTkMxkzVRpAJhVEUZp0Gu6jlpdJBNIy3FRXSk6BpKIpzKmxbG7pQ8v7W3HVWfMwEk1ERzs0Dv3VEaf17J0etTyXXJiZyTgxcObmvBff9uCdfs6Ch672zin/QXifwDMa2i21yG20xtLm+4ome3Y2J37jQghcNtz+3CgfcAUBbs4yOde3Z7KZNEX1ycqvnqgC9WRAID8QcBEwKJhgyi/EF82Wyjl1uqeKgv5HAOViXTW/LGnswLJtP6mxO/NiYYxegv6PBbf/1bFd53KZNHaG0dtNIjp5SG09iWGnONgz56SndFAIo2jPdJfrm/rGUxZ/MaSWCpjGUE198QQCXhRHfGbrhqVxgKjKDn3JGxkJIWUkhnS9dJruqecs6y6BpPIZgV6Y2ldNAI+WyBcWg0Zc8QlRS4S8CKizO6WP2q7aKhusxnluhtJHiOd0WMgy2eVw0PAPzY342dP7807118/14C5N/wDu1v6cPND27CuQe+c5Mi70FrkCZulIVED4ZI2xUVkn+QohMCGA515nc/2pl4k01kc6hy0WAHtfTkBalXiFfZ2ADnRmFddAkBfHqBPyVpLprN4xxnTzf1rSvXOrL48hMoSv/m82t1iQP48Ffk8SCvlaE/cPH73YArprMDKuZWoLQ2Y33ewY0DfPqfS8l3SKlw2swwAcP/GRnzunsIZYtISPdBhFY22vgR+83wDeuMpNLQPWEq/dA9aBxeA/sw+vcuaIKJaGke6Y/jOIzvxy2f2mQMB9X4AejYdoIvGQ5ua8PCmJsuxMlmBCxbV6O3rP85Fg4guJaJdRLSXiG5w+HuAiO41/v4yEc2d6DZpHnLMnhoq5TaWysJDQEnA2dJQg4rJdNa0NEoCGpKGK0Z28KdOL7N8Vh39DyTS6BhIYlppEDMqQhACFl+tHXv2VC6mkTbz6OUPfvMR3cSWo2jZucfzRCOOkoAGIjJ980DO1y8tjbvWH7T8KOV1kWIR8imikbS6pwpZGpmsQE8sZbM08gPh8WQGHYbfWo52S4NW0TCPncoincmaYl8a9Jr32u6eku2rLQ1gVmXY/A67H/mXRsDyGw9uw+0vHdCfDb+GrUf0zr0nlsIZs8txx0dW4YWvXoi3LatHediXFwiXDCQyeS471dLYYWQW1UX1+/Hl+zfj6l+tw/ce3Wn5zL8MN0rnQBLrGzrMe6x2NBZLI5EvblJUzllQhRK/hpkVIdMSk/dpyfQyPPq5N+ET583HO8+YAUCP46hWnD2YvqelD8u++ZiZZg7oz6vXQ6ZV29QTs1gvAHDG7HJDNOLG9+jisGpehWW/fcYo/4fvWY4n//M8fOWti9Hal8DGg51mlpqKrMQs3VN/evkQLvvp87j7lUP49j924J9b9HZ++sIFmFWpW6TqYEA+ixsOdKF7MGURl67BlGndSU/CUztbzHMoZGnEU1l89u7X8Jm7X8P/PrrDso8Ujfbj2dIgIg3ALQAuA7AEwPuIaIltt48C6BJCLADwYwDfm+h2jcQ9ZU7uS6bh0zwI+zVLQbl0JotEOmOxHFKZ3KJNJQFvnnvq0lPrLMdQHyAZw5gWDaDS+AEWGrUC+Ysw9ZvZUxkz9iCFZJPhYz57fhUAfWROJEUjd4ym7rjZ+crRJpAbJW5q7EZj1yAe2dyMhzc15c1klh1V2O8Q0zDOxSmmITvE5h49SBsN+lAa9FknlqkxDaNTle6M0qAXJQXqSA0kcvco5NPMmEWVYU3J75Wjz/Kwz+Kea+1LmCViAGCW4dZa19CBRdNK8ffPvAmr5lWaWS49sRSqIwGcf3INZlaEccu/rUBl2G8mRdhH+P2JdN5gRHVD7GjWYxkC+oDg/o2NAICNh6yB038p75t74jh3QXXed6md+R0vHcQ3Htxq+Q4pKl+5dDH++fnzUFHiN0VVPotlIR9OqY/ixjWn4GQjYaK+LIhKJV4kO8hUJoufP7UHX//bVsRTWbPtgG45lga95uCkuTtuGSTNrgyjOhJAbTRgVgeQLqAzFUvDpxF642kEvB7Ul4WwoLYUp87QB2fv+uU6XG0rSimEMO91Q1s//nWoC197YAt2NPfi8e26WGw8qF/Ld66Ygee/chGCPg/a+hK4c90BpDNZ0+qVAfqz5untkbXMpLUh7117f9J0Y3UOJs0khEMdgzjoEFd5eFOujL7mody9PM4tjVUA9gohGoQQSQD3ALjSts+VAO4wXt8P4M3klPs6jhDlVuqTZArMCFdLo/s0D0I+q6Xxibs2YvF//9MyckxmcjNTSwP57qlo0Idp0YC5vzqiPNwlRSNoruWRcihGKITAA6/lZm/bU24HkmkcNToH6bJ6/XAP5teUoM74gfq9HoR8Wp6l0d6fMEXj/65ejv+48CTLsZ/Z1YarfvESdrX0IStyflw5spIWhsU9pUwg1N9bJz/GU1mcVKsLlHTFlIV9iNotjVhuzoV0T0kzvjToc7Q0AN13L4Ur5NcQ9GkI+zX4NA+iQa/5vdLPXR7y4yPnzMNZc/XR7Cfv2ojl33rcFLsKxcV1/qIaLJkexWkzyrCntQ+xZMa0lFT8Xo+ZFNGXZ2mk86yvnljKPJ60NLoGkuYcnMV1pdjfPoD33roOtxsB2E2Huy1C/+ZTagHYLA3F7fXw5ibc/cphSxbV5iM9mFWpz3mZVRnWhdu4LqpoSNacVo//ecepWD6z3JKZ9tqhbpzz3adx23MN+OETu/HK/k5UlfjxzO429MRSuPmhbXjtcBeiIR+mG6nPR7p1S2NGeQg+jXDG7HIAQG1pEJmsQOdAEnta+zGzIqS7T42BnPx8qVJHbOn0qOU3rXoDBoykleqILoj/UDIMpbW4y0hVryrRf6vlIT8efP0IbnpwG17e32kONKTL96y5umicaTwzjaZo9KLOyAR8xkgoEUIXjmxW4KpfvIibjeKZkkjAa/YXmoewsDaC8rAfkYB3asU0iOhcIrrWeF1DRPPGeOwZANTZM43GNsd9hBBpAD0AqsZ43CHRPJQnGkIIxyq3suOOJbPwaYSwX7MEFZ/a2QohgHteyWV2dPQnTRdOZYkfPbEUhBCme8rnJTx0/bn4+ppTAFjTA2XtpdpooOCCUYA+ovzCvZtyZUuyVrfHYCJjcU8JIbCpsRunzyw3O9aA14OgTzNiGlZrRhbxqykNYG5VrhOaa6RVtvUlzE5bzgPJuaf0z6ruqVz2VH5MQ4rCuQtq4Nc85g942EC4LYskEvAWLD7YH8+N5KWlIetklQZ95jFkpxgN+XDuwmp84eKTAeRiFc8aP3o1MHz6LL1jm11VYmYrOYlGwOspaGmkDbecnQ7b6LpzMIm9bfr1vvrMmRACeHl/J25+eDt6BlNo7IrhMsWSPf/kWhDlLA0ZC5FWnRD6IEd2cOlMFusbOnDOSdXmd6jCLa9/VLnOIb+GD66eA4+HLGLan0jjSHcMPzUmrv7+w2fhF+9fgWQ6iztfOoDbXzqArUd6EQ3qa9fMKA9hT2s/jvbEMaPEFC9qAAAgAElEQVQihJ9ecwY+++aFAHR3IaBbL1uberCwNgIiQlnYh6qI3+JGlZQEvJivCKhaR01aGQtrdStpe1MvokGvOUgE9MFQWchnumXLwz7TpdXYNZg31+gcwxK4cJEu1DuMwc+Oo71YMafczEaUtPclsb9jwHGSpbQqAOBra07BFy9ZBACojvgtg8yJwpVoENE3AHwVwI3GJh+AP4zx2E4Wgz3J380+IKLriGgDEW1oaxtb/rqHHGIaWWf3lJwRHjcsjXBAs7ii5EjoDmXxn8t++jy+/Q/dHzmzIoxkOovuwZTpsvJpHkyLBvGes2bpxxa5meoyQFYe9sNnPKxOxQjtNalyloYa08il3B7qHERbXwJnzKkw3UZ+rz7a7o+n0RtPW0ZlddGQ+VodvV9/0UI89vnzLMeWI7Lh3FPZbK6gnJq2LH989WVBvGF+pTk/obLEj9KgD/2JNPoTabz/N+tNq0J3T1lHXEO5p/QMItk+L4I+zRyVqnET6YKSEwDVoDkA/O11Pb++20E0pEi09ScwmMyg3MHSKBTTAGDGaFTa+hLoi+sB4bKQD92DKew82gevh3DF8lwwujzsw3bDGnnD/CpUhH2oLwuiriyoj+53tWJ7Uy+e29OO3S39uO68+Zbj7DNcLFubetEXT+ONSqdVGvShoX0Af3vtiKOloSLTmWUsENDdmqvmVuLCxbU4fXY5iIBnlDko0ZB+HxbXlWLX0V4c7Y2jLhrEmtPqTRdhrWGZP7WjFQ1tA3jLkmnm8aojAXzuzbq42zvUVfNyLqwfP7Ebn7lbj8HJ5+jkafr3b2/uxfTyEOYqItOfSKMqkrv/atn/xq6YxW3q0wjLZ5bh7585Fx9YPQenzSjD49ta0BdP4WDHIE6pi5oBenlt2vsTlpI6akj1nAW5cfO/rZqNi43zrSkNWBIkJgq3lsZVAK4AMAAAQogmAKVjPHYjgFnK+5kA7OtdmvsQkRdAGYC8qJUQ4jYhxEohxMqampoxNcrjYGkUKlhod0+Ffbq7SXbSwxW6m2n4vo/2xk1zUwqROhKXZr18EINej3lse5E7AHllQuwpt30JXQjCft2SkKmHq+ZWWiyNumgQLb0J9MXTqInkXGanz8oF69XReySgGaay/gOqiwbNuQNqzADIWRyAPuu8L5HOTXIzXDE9gyn8yCiNEA15ccGiWqSzAstmluHs+VUoDXqRFcDana14ca9+Dn7Ng0xW5KWGlgZ8iATyFzkCjJiBYd0EvB7d0ghJSyNnzZgxDeNvlbbihk/uaEVHvx7fOGdBFW56+xJzlCs7UmllloUd3FNDiYbS4akdi7yn0srbcKALc6rCqI0GccGiGlSEdTF5eb9+fZbUR3HG7Aqct1D/nfTEUtjc2IP3/Xo97nnlEKojAbxv1WzLseUgRD4nbzwp12nJa/r5e183g7p2MZVIS0PGOSSr5+udd8CroT4aNONrQK4y8qK6UjS0DaC5J2a6myS1pfo1/vnTexD0eXC5IZiXnVqPty6tw9knVeE3/74Sv/3QSsvnbnr7Ujz9xfMBAI9tazGykpLmHA05gbUnlkJNaSDPGqhWfhPldtGI5e5hRdgPIsKpM8qgeQiXnlqH1w934/Ftegxj+axynDZTH1zI30d7fwKbDufP9J5eFsRJRjuCPo/FzXv+yTVmIstE4lY0kkKfJi0AgIhKhtnfDa8CWEhE84jID+AaAA/Z9nkIwIeM11cDeFoMtf7oOKARwR4mKFh7SsuVHZHuKSDnvx+u5LfMwNFFQ9/Xbyw569PIzOKpsJXMCPo0c46Ik3tqX5vN0shYJ/dJFhgP39M7W1EW8mFhbQRhQzT8Xg/qyoJo7tFHTao/epnxgAOwjN7Dfi88HsLqeVWoi+qWwY7mXgghEDMmz4Vt2VNBnwexZMZizstlZr/459fNGetlIR/evqwe5yyowo/esxxezYNSo0NRZ0cvNEaH9jkskaAXJf5CgfA0Ysk0Qj4NHg9hblUJFhijWDWtt8foCKSgqNlAa06rQzKdxZ9ePoTuWAqnzyrHR86dZ1YZkEIqS7Xku6c083lR4zRydN4xkDRHm9Il2NwTN0f3s41tGw92mff19mtX4XvvWgZAD5zWlgZQUxrAbz+0Ev/7ztMAAB9cPdc4txReO9SNs0+qQsivWQYtd647iN+9sB+NXYOoKvFbOkt1zPLMrjaUhXx5lYJz56Jvl+6v1fMrsbiuFGuW1Zv7zKoMW6znsHHPFtWVIp0VSGVy6aUSmdabyghcsXy6KTRfuPhkfPJ8Peb2liXT8OZTplk+F/JrmFddYsls2tTYgwNGxuLJikiooiG7gmrF0lDLzzR2DVosDbuIvnWp7iKUA6Lls8qxzAjMp5Q6cZsbu81yOVmhH2NudQlmVej9hn3Qcv1FC/Hlty7GRON2hZn7iOhWAOVE9HEAHwHwm7EcWAiRJqLrATwGQAPwOyHENiL6FoANQoiHAPwWwF1EtBe6hXHNWI7pBg85lUaHY8qttAoA3a0kVT+WzKA06Bt2Sr+0NFp64mZWhRQDIkLIp6E/kTY7DvkgBrwe05fq5J4aztKQLKiNYHNjD57e2YoLFtXC4yFz5BjwejC9PIQntrdgWjRougkAYHF9bqSouqfkLO9vXrkUXYNJbDzYhQdfb8K+tgGL+wfQg7DdsSSe3tmaNxckkcri4c3NeHJHKz51wUmoKvHj9FkV0DyEP35stbmfdCE9s6sNp80ow60fPBNP7WzFf/9tq8UVEfZr0DxUOKZhWBry/v3i/SvMv0VDyoTAWBKRgNe8R36vB6UBL/oSabxpYQ364mnc9lwDMlmB8pD1By1FQhaFtK9i6NdU91Suw6kpDZi+8oqwHx0DSSyqK0Uyk8Wd6w5gbtVSADlLAwAW1+UmtkkB2dc2YHa2RGR2fDddvgTLZpbh8/e+jqO9cZxi3NuSgNe0vo50x/Ctv2/H25fV54ndJ8+fj2jIi588uQd7WvuxfFY5CnHytFKE/Rquv2gBskLgO+88DfVlVqthdmUYLyt1nHpi+n1cZMzVmF4WxKq51jkYQUXgZJzDLUSEWZVh7DQs4s/d85rpnppXUwKvh5DOCtSUBnDF6TPQE0thfUMndrX0mUFwwDqAaOyKWQLrqhsL0O/J0ulRbGvqxYLaCMpCPvM39cGz5+CeVw7hqR2t2NzYg2vPmYtfP68nMrx9WT0W10VRVxaEh1BQnCcaV6IhhPgBEV0MoBfAIgA3CSGeGOvBhRCPAHjEtu0m5XUcwLvHepyR4PFQ/uQ+4bxyn1fxzcqUWyDnikllBKpK/AUrhkozu6U3gdnGj96nFLkKGqIhH47+RBpeD8Gr5dxTUpja+xP44n2bcOOaxXmFB9VAeMgIbgO5QF86K7Bijv5jl526np4YRCKdxeHOQSyuz3VEAW/uR2q3NAA9u2taNGiO7J/d3WbW85Id86kzynDqjDKs29eh17dSRmbxVAa3PrsPJ9WU4MuXLMpbBlUiRaO9P4G3L6vH9PKQxQUlA4Nyv4LZU3E9piFH1+rx7IFwe6dZUeJHXyKNurIgLlxUi+f36HWAym3uJ/m5Q0apFqeYRn8ijdtf3I94KmvGUmpKA9htzD2oMJ6lyhI/brxsMT75h3/hbiPJYo6SkKC6j2Yr80n+/ew5jucvO2QAOMW4z6VBL9r7Ezh9VrlZ8qNzIIlSW7tro0Fcf+EC3LJ2L1IZgXmKeNlZUBvB9m9dCgD4/bWrHPeRNaqCPg/iqayZAjy/OoJo0It3nTnT8Xn4wOrZmFNZgpkVhY9fCCka1ZGAJcW9IqxbVUd746iJBLCgNoJvXnkqPvCbl7Grpc9iccn769PInEsi41SVirhI3rViJrY1bccZhsgGvBq2fvOtCPn0Ac4vn9kHzUP40BvnYk5VCRbXlWKlIpb1ZaGCbsCJxpVoENH3hBBfBfCEw7bjCqfsqUxWOFa5VbMprO4pKRpZzKoMFxSNEr8XVSV+HO2NmxPw1CBhyG9kZhgPpBBA0DiGFJdkWp+ctvVID57d3ZZXyE4zRkqyvPnMijD2G3nfqo/2FGN0WqK4p6Q/vqknjnMXVuMX719hSVsErB2xvVOeVRnG/JoSPLe7DSsNX6vq9tDPURcx1T31wt52NHbF8N13nlZQMACY7ikgl1ESUsqTLKiNoL2/02yXPRDuId2KlDOuVf9w7hheDCT1mc7/OthlTuSSVJb4cahzEPVlQcvzYB8FBo2sLGlpOKXcHumOmemVddEg+uL9pr8eyLkjykI+XLKkDpqHzPkCc5TO+ozZOb+2V/Pg6jNnYnZlGBcttrpnJCfVRMwR9VJDNOQ1+8G7l+PpnS34ziM7cahzEPMdysd4NQ9mVYbR0DZgEa/RIF22b1pYgye2t+AdxuRAv9eDp790QZ7YSr79jtNGfcw1p9WhMuxHLJXBQ5ua8MWLT8bBzkH4NA9qSg3RKM11/NJdrFoQ8n6eNqPMnEQ5szyEhvYB08WkcuXp03HL2r24aHGtuU1e84+eOw93vnQAl55aj5kVYXxgdb7Yf/KCk/LcU5OFW/fUxdCzp1Quc9h2zOMhgn1JbVHIPaVZ3VNypC1jGqlMFrMrwwWLs3k1wrRoEC1KTEO1NGQHWxLwwqcRUhlhTjyT+21r6sXXHtiCryi+zLPmViDk9+K53W0I+zT0GXn+qYxAbWkA+9sHUBbyWX4IcrQp/bsBr2ZxHZw5pwJrTsv5niWqPzjs0OmePb8KD21qwpLpUXg9lFd6JOjV54LIzjTg9Zgpnm86eeikBjmn5NKldeacA7XjXz6rHOsbOk1xkS5AOYov8XuRyGTRn9QtDaf2y8/+eeNhHOgYxFcvtfqM5WivPhqyuKTslgagdyyy3MW0qDUNNGC7LnVlQexp7cfJ00pNcZPHKgv54PEQKo0BBwDMUSwK+zX+wbuX57VFxe/1YF51CToHkuYzIV2N06IBs61N3TEzG8zO3KoSNLQNWOaBjAYpOkvqo/jF+1dYhFgd2Y8nV50xE1edMRPbm3rxxpOqcI2SCCCvh5oIIkXAKaZx7sIaUzRmVOii4WQRVEUC2PBfb3EcjFZHAvjn588b8nw/6CAkk8WQokFEnwLwaQDziWiz8qdSAC9OZMOKhcdhcl9B95RnaPeUbpoWHg14PYS6siCO9gwtGpGAF37Ng1QmY/pvpUXS0D6AVEaYE7zu+8TZOGtuhbneQ8AQDZkRIjvaumjQYjVIq8Kr6dlDfs2D+vJcx7ZqnvP0GK/mMV0JTimt9WVB9MXTaOmNO3akIb+Gp3d2YOPBLiybWYZ0RpjpoRUO+6vMKA/h5a+9GbWlAfPHp1oyy42AvTxP2b7qSAB98TSCfg1ahnDrsw0oDXpN14yK/Oxd6w5iRnkIlyy1ztivCPsR9HkQDXktcR+nEbEUjUqjeKKKvaOXx60I+zC7MowDHYOYXRXGydMiZiJCVYnfnGNRFQngvStn4bLTrO1zy3vPmoXeWMq8jpGADyGfhkjAa1pNWVE4nVYG5+eOUTQW1EYwuzKMVfMqLb+FyWDJ9CiW2AodSmFQB1iVpmjktp05pwIXLKrBR86Zi1+s3WvGQdT97Qw1T1ktUzPVGO6u/AnA5dCzmC5X/p0phPjABLetKHgoP6bhxj3l1cgc5aoxDaeifhIi1dKQ7qnc/gHV0jC+J+AzLA3jvQy4yRHsvOoSEJE5QpbxDJmCKjMvppUFLe4d9fxmVIQwrSyI6pIAfBqhOhKwBFrtRAJeeCh/tAzA9OfutvmAJWYcgQh/+vhqc4Tr1zx5riwnpkWDlrarn5HnKjvgsF8DUe5HHPJpeIORq98XT1tWUZTIiWoNbQNYOC2SZ3G+f/Vs/NfblhjB5dzf7Cm1QG40KhMgVOzPyVKjBlnIr5mj77KQD49/4XxzfoHslCIBvV7W965ehgsW1WI0fOxN8/GfxiQx2caF0/RJcqoA2gP4ktNnl6M04MX8mrGJRiTgxXNfudCcDFdsTEtDEQ3pllKf5/qyEG6/dhXKw3786L2nA4AZr6iOFMeNNFEMaWkIIXqgz8J+HwAQUS2AIIAIEUWEEIeG+vyxiFPBQiEAj0Pfr3YgfsU9JdNLU5msJUbhRGWJD12DSbPDssQ0TNHQzBngQSMI7fNYRUOKgux0pStGCpj8uwyM1iuWht2KeuBT5yDo95jpp0umR4ccFUWMcihO+8gf2J6WfrOUgoq0nGZXhREJeM335WHfkMcshOqeitgC4ER68bsl9VFsPNiFoM+DWz+4EtuaevC2n72APodZ1zJ9M5nJOoreitkVWDE7Pzfenj0F5EbpsxyCtTK5oDrix9ovXYDOgSR+8Pgu1JeFzHiFPe1bukmiBbLCxsINly02rV/VQixkaVy+rB4XnzLNMS50LPOO02fAr2mW8758+XT4NY8ljqRyxfLpuGhxLcLGszxaIZ+quA2EXw7gRwCmA2gFMAfADgBLJ65pxcHjWLBQmJ20ChGZAUSvLRCezQqkswI+zYN7r1uNHc29eTVkAL1zyYrcxDHHmIY/l+ZpWhqGuPQrokGUExUpCDKNU84Alw/6tLIggj4NX37rIjMeIFFHyXd+dJUphoUoCXjzaiNJ5Cgrkc46jriCxvnIkg6y83RyZblBdlqah3ICqlhUj3/+PLT0xXHX+oPm9V06vQy3X3uWJegsUT+rjjYL8aePvwGPbjnqaGGWGULiZGlIK63UKMRYGvThpRsuQl00aM4Mti+6JUVMbeN4EfRppoCrHWa0wLGI6LgTDABYOK0Un7NNRowGfXj3ylkFPqEjByofPHvuRDWtaLgdonwbwGoATwohziCiC2FYH8cbHk+hmIbzqFdmJ1VHAjn3VCJjFhL0aR68YX4V5lWXOIqG/EHKeQUW0fDn3FOyU5GioHn0XPtBpax4iV8zs43s2TuyQOHJ00rx7jNn4hKj9MB/XLhgyOthz6N3Qs3pt6Pmslc5jNRl5yqzcqQoOo3U3SAD85cvq0dpQL+2aqfn8ZApFgHFlVVoNKjGfdwEYt94UjXeeJKza0W2w9E9Zdx3NQNNXns598EuWvJ6qrGUiaA06AORbnEXsjSYEwe3T1tKCNFBRB4i8ggh1hLRhJcpLwaaQ+2pjCgctJJupXnVJaY5qi7tKjuDQkE9OaqX+eGqeyqouKfk5+XInIjg83gss7zDSoez5rR6vH5Yn1H6wyd2o6U3Dr/mQXnYh/8bJptmpMyqCDvGMwBrWqJTpyvXYpaWRnCMlkZ52I9HPvsmLKiNwO/14OfvO8NSYwjIXVc3MROraIzNN22KhkOQU4pliUOpk7PmVuLPnzzbDOzb21No9D9eaB5CNOhznKfCnHi4FY1uIooAeA7AH4moFUB+gZzjAKJCVW6H/tzcqhJ4NX2m9mAqnataa3zQp3SqH1w9xyyqZhay65Oi4ZxyK0fk6uxXn0aW9TvUUarf68HNVyzFfa/qhYSbe2KoivhHFScYjv931akF62xFjLYnC7inZOlomS0iRXG0ogHAkgFzuVK4T2Jabb7hs3Ms7qkxpnyagfDyoSwN5/N2igfl3FMTa2kAett7YqkJt2qYqY/bnLYrAcQAfAHAPwHsg55FddzhNLlvKPeUZJ6RNRL2a4glM7kUWq81BgEAp86I4nxjDkJ5nqWhuqdyLgv5eXVE77ON7p3mGchZ6y29iQnLcw/6tIIVZIkI1TJF0SEmIJfylFk3uUD4xGWcEBECXnfZWX6jgCHgLqYxFJeeWocvXnyy4/rqMqliJAJQbbqnJn70LzOo2NJg3JYRGQAAIooCeHhCW1RknAoWZrJD51QDuYyYsE/DQCJjVp+VIqAG0lVrIRfTyHdPyU4t7NccLQ2vLTjvVJBPdkZHe+JmNdHJpjLiR1NP3HGk/rU1p+Df3jDb9N/LDnosloYbwn7NdeA2GvKhrW/sojstGsRnCtRGkjEhJ/dUIaomyT0FAGXhyTsWM7Vxmz31CQDfgm5tZKGvcyEAzB/qc8ci5FCwUBSY3KciO3W9LEY6L6bh8eQyrdTaTTLg29aXgNdjy/UP+eAhvdPKxTRyn/XbfGZhhw5HCksslZkwS2M4ZDDc6fhBn2YpsGdaGqMMhLvl5iuWmrW3hqM06EXXQHJCR9myKGIhi82J6kgAddGgWdl3IpGWxmS4wpipjdsn4EsAlgoh2ieyMVMBzUN5Jc2zQjiWEXEi7PcagfD8Gd4+zYN0NmMGPQHdr+7X9FXb7CuMvXPFTJw8rRTRoM8Un6HcU04djtpuJ/fQZCBHxPZqn06MR0zDDVeebl8ksjClQX0xn6HqYI0VmfXlFO8ohN/rwfqvvXmimmShtjSAyhK/pXQOc2LiVjT2ARgcdq/jAOeYhvPKfQDwk/eebq4cBuTqGiUdJuv5NEIslcsQAnS3VzTkQ3t/Iq9uT0nAizfM18t3+M0Z4ap7ytqmEqeYxiTU7hmOFbMrcKhj0FVZCHVy31ShusTvuITkePKB1XMQ8mu4+syh8/+LxacuOAlXrXAvtMzxi1vRuBHAS0T0MgCzdrAQ4rMT0qoiQg4FC7NZ55X7AJhVOCXRoA/72vrzAuEAlLiEtfMsC+llqOdWFy7VYU+5VbdJnCbhaYpo1RbJ0vjA6jmOlTqdMGMaE+yeGgk3Xb5k2LVRxopX8+C9Z80efsciURUJOM6zYU483IrGrQCeBrAFekzjuEUjOK6n4dY9VRbyoTeeyotpAHCMSwAwlzmdV13YN+23Te5Tv0/iFERVLY26svwZz1ONFbMrcM6CqiEFdLIZa7lvhjmecCsaaSHEf05oS6YIHod5GkO5p+xEQ170xtIFYxpAfmG/zkHrBDcn7GVE9G22QPgQ2VMAMM2hTMZUY+G0UsvqfAzDTC3cRrXWEtF1RFRPRJXy34S2rEh4HAoWZoZwT9mJBn2IpTJmIUG1Y3dKmwVgLi85VFlpexkR/butt89pZTp1HzX2wjAMMxrcWhr/Zvx/o7LtuEy51YhgMzSMlFu3loYewJWr9TlZGnbRmFddgv3tA441iXKfpbzP5sc08t1TqqVhPy7DMMxIcTu5b95EN2Sq4PHkz9PICl1M3CDLLHQYk/XUaqd+s+O3dvZ3f3w1Gtr6h8wucgqi291TTim39gwrhmGYsTDcyn0XCSGeJqJ3Ov1dCPHXiWlW8XBahCkrhON6Gk7IGbNOVWtzMQ3riL+uLDhskNrps3KbrEA6nKXBMAwzVoazNM6HnjXlVGdKADjuRKNQ7Sm3hf7K8txT1iVhfRqNqiN3tjT016fURdHSG3esaSRnhPNMXoZhxoPhVu77hvHyW0KI/erfiOi4dFl5iMbonjJEQ7qnVEvD67EEskeC3yEeIgVp4bQIHvncmxw/J6Cfy1grtDIMwwDus6f+4rDt/vFsyFTB41CwMOui9pQk557Kr1rr18gyo3skDGVpFFrLAgDSxnyRU+qjBfdhGIZxy3AxjcXQl3Qts8U1otDXCj/u8FD+yn16yu1IA+GGe8pr7eTdrOHgRMCpyq0hGk5Li0qWTo/i++9ahjXL6kd1XIZhGJXhHN2LALwdQDmscY0+AB+fqEYVE81hnoYYweS+kE+D10OOMY0z51Q4zqVww6VL65HOCsxQCtrJbCy/Vth6ISK856ypWc+IYZhjj+FiGg8CeJCIzhZCrJukNhUVj4dgX4ROLyPi7vNEhLKQLycaStrVx940+mktZWEf3v8Ga/0mnwtLg2EYZjxx29tcRURRIvIR0VNE1E5EH5jQlhUJJ/eUm5X7VGQw3OuhCS2n7cY9xTAMM5647W0uEUL0QndVNQI4GcCXJ6xVRUSjfPdU1sXKfSpRI73VTSnwseB3WAKWYRhmInHb28jFDdYAuFsI0TlB7Sk6Hts8jWQ6i2Qm6zhxrhA1RmFA+4zt8ca0NHhhHIZhJgm3UdmHiWgn9OVeP01ENQDiE9es4mGfEd4T04sJVoxgUSBZQ2qiLQ2OaTAMM9m46m2EEDcAOBvASiFECvoqfldOZMOKhWYLhHcbZcvLwu4XBZKiMZBMj2vb7PjYPcUwzCQzZG9DRF9R3r5FCJEBACHEAIDjbtU+QK/jpM4I7xocjaWhLyAUT03selVsaTAMM9kM19tco7y+0fa3S8e5LVMCzeaekpbGSJYfnVVZuMT5eMKiwTDMZDNcb0MFXju9Py6wFyyUCySVj8LSmGh85uQ+Fg2GYSaH4XobUeC10/vjAiI9piEM4egyLI2KEveWhqx0O9GwpcEwzGQzXPbUciLqhW5VhIzXMN4fl7WnZDVbvbKtHtPwegglI0i5nSxYNBiGmWyGKyMy9XrKCUZO4M4KAQ2EnlgS5WH/iCb3AcD/Xb1swk0xzp5iGGayKcrKPERUCeBeAHMBHADwHiFEl8N+GQBbjLeHhBBXTHTbZNmPTFbApwFdA6kRZU5J3r1y4osELpkexVlzKxwXX2IYhpkIijVEvQHAU0KIhQCeMt47ERNCnG78m3DBAHLLo8pYeNdgckRB8MlkZkUYf/7kG1E+gjkkDMMwY6FYonElgDuM13cAeEeR2pGHdE/JuRo9sRR3ygzDMAbFEo1pQohmADD+ry2wX5CINhDReiKaFGGR1Wxl0cKuweSo3FMMwzDHIxMW0yCiJwHUOfzp6yP4mtlCiCYimg/gaSLaIoTY53Cs6wBcBwCzZ88eVXslOfeUFA22NBiGYSQTJhpCiLcU+hsRtRBRvRCimYjqAbQW+I4m4/8GInoGwBkA8kRDCHEbgNsAYOXKlWNKWlItjXgqg2Q6O2nzLhiGYaY6xXJPPQTgQ8brDwF40L4DEVUQUcB4XQ3gHADbJ7phZvaUEOiL6wUHS4NFSTJjGIaZchRLNL4L4GIi2gPgYuM9iGglEf3G2OcUABuIaBOAtQC+K4SYeNEwAuFCAH1xvYQIiwbDMIxOUXpDIUQHgDc7bN8A4GPG65cAnDbJTTNnhGeyiqURYORNG+YAAAouSURBVPcUwzAMUDxLY8oi3VNZdk8xDMPkwaJhQwbCs1mgPyHdU2xpMAzDACwaecgq4xkh0MuWBsMwjAUWDRumpcHuKYZhmDxYNGzk3FPCzJ6KBFg0GIZhABaNPDRPbj2NvngaYb8GL6+MxzAMA4BFIw+zYKFhabBrimEYJgeLhg17TIMzpxiGYXKwaNjQlHka/Yk0WxoMwzAKLBo21IKFvWxpMAzDWGDRsGGdEc4xDYZhGBUWDRsaWbOnoiwaDMMwJiwaNvKzp9g9xTAMI2HRsCHdU+mMQDyVRYmfLQ2GYRgJi4YNGQhPpDMAAJ+XitkchmGYKQWLhg05+TuRzgIAfB6+RAzDMBLuEW1ISyNpiIZXY0uDYRhGwqJhI180+BIxDMNIuEe0IWeEmzEND1saDMMwEhYNG4ahYcY02NJgGIbJwT2ijZylYQTCOabBMAxjwqJhQyOraHg5e4phGMaEe0QbxNlTDMMwBWHRsJEXCGfRYBiGMWHRsMHuKYZhmMJwj2hDZk+xe4phGCYfFg0b+dlTfIkYhmEk3CPaMAsWpvSYhpcn9zEMw5iwaNiQIYxkhi0NhmEYO9wj2jAD4SmOaTAMw9hh0bBhFizMcPYUwzCMHe4RbXhs8zQ4psEwDJODRcOGzJ6Ks3uKYRgmDxYNG3IGeCwpZ4TzJWIYhpFwj2jDb4hEfyINgN1TDMMwKiwaNogIAa8HA1I02NJgGIYx4R7RgYDXg3RWAOCChQzDMCosGg4EfJr5mlNuGYZhchSlRySidxPRNiLKEtHKIfa7lIh2EdFeIrphstoX8OYuC1saDMMwOYo1jN4K4J0Aniu0AxFpAG4BcBmAJQDeR0RLJqNxUjQ0D5mLMjEMwzCAtxgHFULsADBch7wKwF4hRIOx7z0ArgSwfaLb5/fq7inOnGIYhrEylR32MwAcVt43GtsmHGlp8BwNhmEYKxNmaRDRkwDqHP70dSHEg26+wmGbKHCs6wBcBwCzZ8923cZCSNHg2eAMwzBWJkw0hBBvGeNXNAKYpbyfCaCpwLFuA3AbAKxcudJRWEaCzJ7izCmGYRgrU7lXfBXAQiKaR0R+ANcAeGgyDpxzT7GlwTAMo1KslNuriKgRwNkA/kFEjxnbpxPRIwAghEgDuB7AYwB2ALhPCLFtMtrH7imGYRhnipU99QCABxy2NwFYo7x/BMAjk9g0AEDAyJ7ysXuKYRjGAveKDvjZ0mAYhnGERcMB0z3FlgbDMIwF7hUdCPg4EM4wDOMEi4YDMqbBZdEZhmGscK/oQM49xZYGwzCMCouGA5xyyzAM4wyLhgM8I5xhGMYZ7hUd4BnhDMMwzrBoOMAptwzDMM5wr+gAxzQYhmGcYdFwwCwjwim3DMMwFrhXdIBTbhmGYZxh0XBAzgjnyX0MwzBWuFd0IOeeYkuDYRhGhUXDAc6eYhiGcYZ7RQf8PE+DYRjGERYNB3IFC1k0GIZhVFg0HGD3FMMwjDPcKzrA62kwDMM4U5Q1wqc6Yb8XX7l0ES5dWlfspjAMw0wpWDQK8OkLFhS7CQzDMFMOdk8xDMMwrmHRYBiGYVzDosEwDMO4hkWDYRiGcQ2LBsMwDOMaFg2GYRjGNSwaDMMwjGtYNBiGYRjXkBCi2G0YV4ioDcDBUX68GkD7ODanmBwv53K8nAfA5zJV4XPRmSOEqBlup+NONMYCEW0QQqwsdjvGg+PlXI6X8wD4XKYqfC4jg91TDMMwjGtYNBiGYRjXsGhYua3YDRhHjpdzOV7OA+BzmarwuYwAjmkwDMMwrmFLg2EYhnENiwYAIrqUiHYR0V4iuqHY7XEDER0goi1E9DoRbTC2VRLRE0S0x/i/wthORPQz4/w2E9GKIrf9d0TUSkRblW0jbjsRfcjYfw8RfWgKncvNRHTEuDevE9Ea5W83Gueyi4jeqmwv6jNIRLOIaC0R7SCibUT0OWP7MXdfhjiXY/G+BInoFSLaZJzLN43t84joZeMa30tEfmN7wHi/1/j73OHOccQIIU7ofwA0APsAzAfgB7AJwJJit8tFuw8AqLZt+z6AG4zXNwD4nvF6DYBHARCA1QBeLnLbzwOwAsDW0bYdQCWABuP/CuN1xRQ5l5sBfMlh3yXG8xUAMM947rSp8AwCqAewwnhdCmC30d5j7r4McS7H4n0hABHjtQ/Ay8b1vg/ANcb2XwH4lPH60wB+Zby+BsC9Q53jaNrElgawCsBeIUSDECIJ4B4AVxa5TaPlSgB3GK/vAPAOZfudQmc9gHIiqi9GAwFACPEcgE7b5pG2/a0AnhBCdAohugA8AeDSiW+9lQLnUogrAdwjhEgIIfYD2Av9+Sv6MyiEaBZC/Mt43QdgB4AZOAbvyxDnUoipfF+EEKLfeOsz/gkAFwG439huvy/yft0P4M1ERCh8jiOGRUN/mA4r7xsx9AM2VRAAHieijUR0nbFtmhCiGdB/OABqje3HwjmOtO1T/ZyuN9w2v5MuHRwj52K4NM6APqo9pu+L7VyAY/C+EJFGRK8DaIUuwvsAdAsh0g7tMtts/L0HQBXG8VxYNHTzz86xkFJ2jhBiBYDLAPwHEZ03xL7H6jkChds+lc/plwBOAnA6gGYAPzS2T/lzIaIIgL8A+LwQoneoXR22TfVzOSbvixAiI4Q4HcBM6NbBKU67Gf9P+LmwaOiKO0t5PxNAU5Ha4hohRJPxfyuAB6A/TC3S7WT832rsfiyc40jbPmXPSQjRYvzQswB+jZwbYEqfCxH5oHeyfxRC/NXYfEzeF6dzOVbvi0QI0Q3gGegxjXIi8jq0y2yz8fcy6O7TcTsXFg3gVQALjWwEP/Tg0UNFbtOQEFEJEZXK1wAuAbAVertltsqHADxovH4IwL8bGS+rAfRIl8MUYqRtfwzAJURUYbgZLjG2FR1bvOgq6PcG0M/lGiPDZR6AhQBewRR4Bg2/928B7BBC/Ej50zF3XwqdyzF6X2qIqNx4HQLwFugxmrUArjZ2s98Xeb+uBvC00CPhhc5x5ExmJsBU/Qc9E2Q3dF/h14vdHhftnQ89E2ITgG2yzdB9l08B2GP8XylyGRi3GOe3BcDKIrf/bujugRT0EdBHR9N2AB+BHtDbC+DaKXQudxlt3Wz8WOuV/b9unMsuAJdNlWcQwLnQ3RWbAbxu/FtzLN6XIc7lWLwvywC8ZrR5K4CbjO3zoXf6ewH8GUDA2B403u81/j5/uHMc6T+eEc4wDMO4ht1TDMMwjGtYNBiGYRjXsGgwDMMwrmHRYBiGYVzDosEwDMO4xjv8LgzDjAUikmmrAFAHIAOgzXi/Suh1jRjmmIBTbhlmEiGimwH0CyF+UOy2MMxoYPcUwzAM4xoWDYZhGMY1LBoMwzCMa1g0GIZhGNewaDAMwzCuYdFgGIZhXMMptwzDMIxr2NJgGIZhXMOiwTAMw7iGRYNhGIZxDYsGwzAM4xoWDYZhGMY1LBoMwzCMa1g0GIZhGNewaDAMwzCu+f8bKv6C0dW+LwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy\n", "import statistics \n", "mean = numpy.zeros(300)\n", "T = numpy.zeros(300)\n", "Tval = 10\n", "i = 0\n", "while Tval < 3010: \n", " X = numpy.random.normal(0, 4, Tval)\n", " mean[i] = statistics.mean(X)\n", " T[i] = Tval\n", " Tval = Tval + 10\n", " i = i + 1\n", " \n", "import matplotlib.pyplot as plt\n", "plt.plot(T, mean)\n", "plt.ylabel('Estimate')\n", "plt.xlabel('T')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Generating Stationary MA(1) Process Data: $$X_t = .5e_t + e_{t-1}$$\n", "$$e_t \\sim N(0, 4)$$\n", "$$ (\\sum X_t)/T \\rightarrow E[X_t] $$" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXmcHVWZ//95qu7et/fuJJ09IWFJgLC0kR0FVGBAVMRBHQd1FMdlhhmHUfy566g4rl9HBsR91BEQZUDZN0H2JJiFkIXs6XSW7nR677vW+f1Rdc49VXfv7tu3l+f9euXVd6nUPbWdz3mW8xwSQoBhGIZhSsWodgMYhmGYqQULB8MwDFMWLBwMwzBMWbBwMAzDMGXBwsEwDMOUBQsHwzAMUxYsHAzDMExZsHAwDMMwZcHCwTAMw5SFr9oNqAQtLS1i8eLF1W4GwzDMlGHdunXdQojWUradlsKxePFirF27ttrNYBiGmTIQ0d5St2VXFcMwDFMWLBwMwzBMWbBwMAzDMGXBwsEwDMOUBQsHwzAMUxYsHAzDMExZsHAwDMMwZcHCUSF2dQ3iuR3d1W4GwzDMuMPC4eHBTQexq2twzPu57amd+Pe7N45DixiGYSYXLBwePv37jfj1C/vGvJ9Y0kI8ZY1DixiGYSYXLBwe0pZA2hp7h59MW0iNw34YhmEmGywcHgSAtBBj3k8ybSGVHvt+GIZhJhssHB6EAKxx6O8TaYFkmi0OhmGmH1UVDiL6GREdIaJX8nxPRPQDItpBRBuJ6IxKt8kSAtY4KEcyZSE9HgrEMAwzyai2xfELAJcW+P4yAMudf9cDuLXSDRKwxWOs2DEOATEO+2IYhplMVFU4hBBPA+gpsMlVAP5H2LwAoIGI2irbKGA8PEzSTZViq4NhmGlGtS2OYswDsF973+F8VjEsMT5WQsIJjHOAnGGY6cZkFw7K8VnOnpiIrieitUS0tqura9Q/OJ5ZVQCQ5JRchmGmGZNdODoALNDezwfQmWtDIcTtQoh2IUR7a2tJy+bmRAgxLllVylXFFgfDMNOMyS4c9wH4eye76iwAfUKIg5X8QUtg3LKqAPAkQIZhph2+av44Ef0WwBsAtBBRB4AvAvADgBDiNgAPALgcwA4AwwA+MBHtGo+sKo5xMAwzXamqcAgh3l3kewHg4xPUHBUUH4/5F+yqYhhmujLZXVUTitSL8YhxpDg4zjDMNIWFQ0NaHOMzAZBdVQzDTE9YODRkFz9W4RBCICEtDq5XxTDMNIOFQ0PqxVhjHPpsca5XxTDMdIOFQ0NaGmP1VOlWBqfjMgwz3WDhyIEUkFcO9OFQX6zs/59MZZQnyTEOhmGmGVVNx51seF1VV/zXMzANws6vX17WfhK6xcHCwTDMNIMtDo1crqrRxCh0VxWn4zIMM91g4dCQEjHWIodJtjgYhpnGsHBojNc8Dl040mxxMAwzzWDh0FAzx8eYQpvg4DjDMNMYFg4dreTIWBZz4nRchmGmMywcGgKZIof6JL5yRcQVHGeLg2GYaQYLh0amyKFwdf7DiXRZ++F0XIZhpjMsHBp6cFy3FAZiqbL2o/9fdlUxDDPdYOHQyBQ5zJRFB4D+WLKs/cjV/wC2OBiGmX6wcGjINFzL8locZQoHB8cZhpnGsHDo5Ilx9I+U56pKcHCcYZhpDAuHhj5zPDkWV5Ue42DhYBhmmsHCoZFxVbnX1Cg/OM6uKoZhpi8sHBoin6uqiMWxv2cYu7oG1Xuex8EwzHSGhUNDXzq2nHTc/7j/VXzq7o3qfSLFtaoYhpm+sHBoyBpVacuTjjtS2OLoH0lhMJ4RF+nmCpgGWxwMw0w7qiocRHQpEW0joh1EdFOO799PRF1EtN7596GJaJcQwpUZ1TtcWDgSacudSeVYHCG/wTEOhmGmHVVbAZCITAC3AHgTgA4Aa4joPiHEq55N7xRCfGIi2qRWABRCZUP5DELXYLzg/0ukLE9cwwIREPKbnFXFMMy0o5oWx2oAO4QQu4QQCQB3ALiqiu1xTQCUlkJbQwjdA8WFQ49rJNICftOAfwq5qtbtPYadWoCfYRgmH9UUjnkA9mvvO5zPvFxNRBuJ6G4iWlDJBuklR+SaGm314eIWR9pyCUQybSFgGvCZNGVcVVff+hwu/s5T1W4GwzBTgGoKB+X4zDs8/yOAxUKIUwE8BuCXeXdGdD0RrSWitV1dXaNqkF7kUHb4c+tDGIilEEvmr5CbSFmu+lTJtAW/SfAZVJarat/RYdz0+42uwPxEMJp11RmGmblUUzg6AOgWxHwAnfoGQoijQgg53P8xgDPz7UwIcbsQol0I0d7a2jqqBsn+M21l5nHMbQgDALoKuKviKU9wPG3BbxrwGeUFx5/Z0Y071uzHwb7YKFo/erqLWFQMwzA61RSONQCWE9ESIgoAuBbAffoGRNSmvX0rgC2VbZKtHEJkJu61SeEo0LkmUmkk0payWBIpO8bhM8uzOKRYJSbY4ujsHZnQ32MYZmpTtawqIUSKiD4B4GEAJoCfCSE2E9FXAKwVQtwH4J+J6K0AUgB6ALy/sm2y/+pZVXPrQwAKWxy2aNiWis8kJNIWAj4DPtNAsgw3kBKO1PgJx6aOPixprUE0mP9SSwsn4ONpPQzDFKdqwgEAQogHADzg+ewL2uvPAPjMRLUn1wqAbfXFXVWyo0+mBXymbYEETAN+g8qKV0grJzlOFkciZeHqW5/Dpy49AR86f2ne7aRwtNQExuV3GYaZ3vAQU0O4XFV25z3HsTj0OMCn796Ir/7Jnm6SSltKcKSLKZ6SFscoXVWOEG3q6MMJn3sQh/tHF/OQExOLTWA86LiqQgFzVL+jMxRPYXf30Jj3wzDM5IWFQ0NofXzc6bzDfhNNNQGXxXHn2v346TO7AbjjEbLDT6QsBH3OPI4yguPeGMeeo0OIpywcGmWwXFo7IwUywoCMxTEeLrKfP7sbb7vl2THvR6c/lsSbvvsUXjnQN677ZRidVw704V/vXM9ZhiXAwqFhacohhcNnEhrCfvTlqVeVTLnnbwB2BxzwGTANKusmlIIhXVayDfmC5WlL4JYnd+St3iv3UyiVGAA6+0ZUu8dK92ACfSNJVfdrPNjbPYzXjgzi1c7+cdsnw3h5ansX7vnrAfQOJ6rdlEkPC4eG2+KwO1ufQTAMcomKTjyd6ZR1iyHgs9Nxy5k5LkVIt1z0v14e3nwI33p4G7754Nac38tU4FiysCBIN1x8HIRjOGEXexzPzDAp2rFUYQEshbvXdeDJbUfGvB9m+iELlRaz0BkWjrzEk/YkPiKCQfbiTrlwlRrROvqAacBvlhsct1x/pXjl64Rlh6oH07ce6le/KeMrxTpc2e74OHTMw4m0s6/xF46RxNjbd+PvNuADP1+j3j+59QgLyQxiIJbE3qO5Y3CDzvIJxQZaDAuHC6+rymfYp8cgQjqPxeGuUeV2VflMw7WSYDG8wfFiFod0QYX9dlD7yEAMl37/L/jCfZtd+4sV6XClVZRIZeaijBbZuY9nSrESjgqMBH/45A7895M7xn2/zOTktqd24prbns/53YDj8i3m2mVYOFx4XVV+066KYhqU12fvnjGeiU0EfHY6bjmptaoD17Kz7M9z70N2pCFHOOLOSOmP6+0J+FK0ilkcslyKJVCW0OVCWhzluKqG4oUXylKuqjGOBHNdw+FEelxFrhoc7o+51oNh8nN0MIGuwXjOe0GeQxaO4rBwaOi3Ujxllw0BbOEoxeJIah1+0GeWnI774q6jeNstz6oONJ/l4UV2pEFHOGRnPeDZT7EON562YDiVw8bqYhpOlmdxPLL5EFZ+8WFs6sifMZURjrE90LlWchxJpMbVrVYN3vuTF/HdR7ZXuxlTgpFkGkIAQ4nse6GfXVUlw8Kh4XJVJTPCQUTINxDPHeNIIyhnjpcw8t50oA/r9/fiyIA7LVbGHPLtI64sDiOrLYAW4yjQ4QpnsmONM7N8rKPvERkcL3E/z+08CgB4cffRvNuMl3DkyoybFhZHX0zN9TnSH8ubAThduf3pnSo9vhjSlZrLQsvEONjiKAYLh4bXVeWTrirK7eYA8sQ4nKyqoM8oqVOSN+qgx1KIlxjjMIhc/0+2X74vFBtIWwJCALWOcIw1QJ4Jjpe2n7qwH0BmtJeL/nGKcfSOZKdZjiTSY7Y4fvX8Hpz0+Yeq0uEIITCUyCxdvPrrj2P11x6b8HZUk/s3HsSDmw6WtK28h3JZnwPxysXSphssHC7cwfGA7qqy5Kxyt4DEcywZK7Oqgj7T9X0+pGksRzxZwfE87i55g6fS2QKz7+hwJuZSwPSW20RD2RbHnzZ24oESH0jVphKC4/t7hvFvd21AImUhGrTdbAN55qIA42dxeGfQCyEwnBy7cHz+3s0YSaarUiwylrQrF+gj6HjKmlGT2IYSaTVgKUasgHCwxVE6LBwalsvisJTFoWdVeedluGMcQpUgCfgMBByLo1imkrqZ43IOhGcCYJ6ObSTpnjCob/fBX65RpT8KriXiiE5UWRyZffzsmd0luwAkQyW4qp7d0Y3fv9yBvUeHVNv7R/JbHJmsqrF18F4XTjItkLYEEmOwsg5oYjHaGf5jQZ5vb4LBloMzZ7LkcDxVspUwksztqhJCKDGJTXHX5UTAwqHhclUl0650XKGEw31TuV1VmdGrdFUBxQPOMuvJ66oqFhyXnYWc6KdbN/t7RnD/Jju7qpBwyN+KhmyXkW6dxJJW0YwnHcsSynoqZGnJ4xyIp9SEwWMFZuvmsziEEPivx1/DrhKXvO31CMfIOMw5uXNNZhHLzioIx3A8c+/oA5QXduWPGU03bIujtPtUXnOvhRtPWZksxHGYLzTdYeHQ0B+8RMqC35ftqvJmSbksjpRQ74OacBRLTZWdrfz5UoPj8uZPeSyOP3zsHADAkNOp2JkkhWM00mWU0GbCx1K5XQBpS+DXL+zFHS/tc32uj/oKWRxSOAZjKYwk7O1kYkAu8glH/0gK33l0Ox585VDe/+vajyNOsnz8cDIzy12en1TaKjk+0zUQx0/+sgtvPMFeOOygZn0IISZkgSx1LuMp17X6677eiv92Kew4MjjmuUHFGE6kynBVud3CEr1sz2R0VT246WDZbuNKwsKh4XVV+Z0cVcMgSL3wioCryGE6sxJgwGeoNFk5iu8dTuRMO/WOeLOC43mFw13eI+lxO8lRmCWyXWze31KuKs3iiCetnCO57z26HZ/7v1dw0x82AbA70E0dfa6Ht5BwDGmd3YjTeR/uj6OzdwSXfv9pV6zAsoR6qL0zx/vLnLClYhzOqZDt1Rfuesetz+GEzz1U0v7uWrsfI8k0Pn/FCjTVBHDPXw/grT98BiOJNP68vQtnf+NxHBllZeNSGdZcVbp1WMnMqoFYEveuP1B0u+2HB3DJd5/CC7t6KtaWRMpCMi1KriqQz1WlC4n0AGzY34tLv/903lpwE8lHf/MyPvabl7MGNcm0NS4VFcqFhUNDeILjah6HllVVyFWVTFvqfcA0EDSlq8q+sF9/YAuu/OEzeO3wAAC7bPoX7n0l68JnLI6MqypXRyCFw2txBH0GAqahLA4g/yRAeTwyHTfuyczKlba4zWl/nRNQ/+8/78AHf7nGdRyyLbu7h/CDx19zjToHpXsllhkpdg/GselAH7YeGsCG/ZnR8kA8pSwx7zGUGzSX20sLw9Ve57g3FphP4mXLwX7MbwxjaWsUbfUh7OoewsaOPmzu7MP+Hjs5oaPCAfOMe1OgR3P3DVRgQmD3YBx3rtmHP208iBvuWI+OY8MFt5fxn/1FthsLGXetKDhYkchr7s3i04Pl0gpet/cYth4awPZDA+PV3DHz5FZ3eZxvPLAVf/fTFye8HSwcOvnScbUihwVjHCnLHePwzK+Qo9qvP2CvgHvlD5/B/zy/F10eN03CE+P46/5enP6VR7DT48sf9MQ4EtpvB3yGy1rI17kmnMKKtTksjljSQixpYc2eHjy/M+Mzlw+fHL3JirjS9aMfw33rO/HdR7fjiFaWXj7s/bGka9S/q8sO5h/SRukyFTdgGuqBVt85I8FSYxR6jKOzL+aybOKe85MvtnO4P4Znd3QDsN0wy2fVAsgs+AUArx7sV+3uGaxspVXdyjvcb59je9Aw/sJxz8sH8Onfb8LOI/Z9eGyo8EhcVpk9WsFzoE/kG0mk8cV7X8HH//dl1zaWlhGp4oke4dAHSHKbo0P2+ayk8BXjQO+Ia32be/7qtvT29Qxh++GJFzYWDo0sV5U2ATCdz+JIuy0OaV0EcwTHI85CSU9u68L+nszNeHTI/WB5y6rv7hqEJYADx9w+dH20qbfF7xRY1DuVmNPp/uOv1rncDPJ4amVwXBvVS7H56p9exc0PZpZ7l4KRTNuTBwdiSSRSluthlB1xj/Pw6Z20y1WltVEKoy4cUnDmNoSyOveBEtInb3tqJ/7yWhcAoE9Lx/3or9fhxt9tUO+97sCDfbkthZ89sxsf+PkaxJJp7OoewvJZUQDuyaOvdvarEW1PmSW6f7+uo6zyIbpAHHaC87Prg2UJx5GBWEkl62UcSpbhL5TQAGRcg0c9sZ41e3rGpaAm4BbO4WQKL+05hjW7M66xr93/Kk78wkNIWwLxlKWs18G4W/QGcsQ4pOB19Ex8mrXkn/73ZfzdTzIWxWuH3YPHoXgaA7HUhMdlWDg0dFdVImVlalVRxuJIpNyxAn20m0iLrFG/vo2+7R83dqrXPR7hkOmh8u+A1tFKRpJpLWCf2+LQfy+WSqN3OIGHNh/CXWszmUAJj6tK7iOVzmSZHDg24vpt18OaSKvRtS6Acr/ys4NaxpEeHB9OpJSgSjFds7sHp3/lEWw7NKDcHctmRbNSLvuL1LAaSaRx84Nb8b6fvgTAPQFwf8+wy13hnetyoDd3bKJrII5E2sKLu3uQSFlY5giHPIbW2iC26BbHUOnCse3QAP7tdxvwb3etL/n/uITDEdw5daGsEXUhvvXQNvz9z17K+d1ALIn1jutQLmbW6Zwbb5aaFyUc2jnY3NmHa257ftQlUn71/B48sfWwej/kuS8PHBvGkYG4EqYf/2U3EikLj7562NW5btjfh7vXdaj3chAS9Blqu25HOKppcew5OqyegTMXNaKzb8Tl9pVehXLus/GAhUPDm/zhyzEBMGVlu6r8JiHglBfJxBlMBH2m2kb/C7jXMPdmhHgtDjVK0joD/XXSclscATMjWpJYMo2djito3d5jGfeZzKoKuedx6LnsR4cSrs55xOMekB2w7pKQ+5efuSyOhDsTqM1ZnlcKx8v7enFsOIlXDvQpK+u4WVGkLIFLv/801u09BkCvLZR7tPVKpx2v8DlJDrrweDs9ee5kza58k/nkKPuJLXbntXy27ar68ltX4nt/uwpXrZqLrYcG1IN8rIwHWg5c9AQKyxL43P9twsaO3FlSQ7qryrEIZtWFMJhIlZzNtP3wALoH4zkXMPr1C/twzW3PIZZMo2vQbT32DMax9VB+S0XGlPTsMul20d0vQoiCbR2Mp9T9dNtTu/C/L2ay+fRn50h/XN0TBx1xO395CwA7oaNDs9i3HR7Ajb/bgPs3HsRHf71OtbW1Nqjudemq0v9fMcYzgyyWTLsEoX1xI2JJC8c0y1kOwiYig0+HhUPDe8nlzHHDINV554pxyLU3EilPVpXPHRzXhaPQhLd88zf0gKc+Wk56tg+Yhmq7JJa01HyHWNJSHVFmAqBb5LydsR4vGU6kVX2s4URKG11nbl65n55cFkcsM49jJJlW8YFDngykwwMxdBwbRkPEj5aaIAC4gufK4sgT4/jrPltgZOcul/QFsgcJ0uKQJVDyCUeP89A+tsUOUh7XWgMAaI4G8fbT5+P42bWIpyxsdQKqXjdkIWQburROoGswjl+/sA+PvXo45//RR9yH+uz/N7s2BCGyBySSO9fsw/977DUAdkcnY0u51oo/1DeCZNrObJODHdm+O9d24NLv/yXvPBopsnrnJ1831QTsY06lsfrrj+O+DZ3ZO3A4+YsP431OAHg4kUKXNkDRj3+H1g45SpfndM/RIbz3J9lB5K/8aTMefMW2wptqAphdF8pyVeWyOBIpCz96aqfrOTkyEMPKLz48pjk0P3j8NXzq7g3Y0z3kmlDqNwmnzW8A4L435TX+w8sH8I0Ht4ybC7AYLBwa3lX+5EjVIKiZ415XVSKdVq4hb1aVclUls9NqC5XY8JZVl+hWhv7ASJdSImWpFQv9HuEYcSwOeUxPONkZquRIUMY48glH5v1IIo1mpyMfTqRVkLpbe6DjHlfVr17Yi3NvfgL9saQrq2okkcbsOtvi8FbJONIfx4HeEcxrCCPkuIKATFC8WDqunMug5tOkLBXL8SLnr0gR7szjqpKj8gO9I1g2K5q1v1l19nmRGUflWBwxLXYkkaPdfALkGnE7FsecersN+eIcv1vbgd+ts92VXYNxNSDJJRzdzu8OxFJKOORjIq2NTXnWgs/EODJtlxZkyG/i58/uxubOfnQNxPHi7sIpu/L7oUQa3U477nhpH36nuZt2aEFi+TvDyRQuOnEW3n/O4pyZiTKhYPvhQVywvAWRgKm5qjKuOe+CbM/u7MY3HtyKhzcfUi7CzZ39GE6k8YpzPlJpC998aCvW7XUf27cf3qbcf15+8PhruGttBz5190YVSwLs5Iv5jRGnPZnPpcXxi+f24EdP7cI//GLthKTnsnDo5HNV5QmOH+gdwZ7uYXvtDY+ryrY43OXOXRZHAeHwrgAo0QN6ejaJvnKgFKtgTlfVIJa01OCyk+fg9qd34eV9x7LncaRyz6aW9Y9kfaeWqD1i7B9JKtNeH1kmUhYsS6hRZyJl4UDvCF7YeTQrOF4X9qnf1zkyEMOBYyOY3xhWi1XZv5ly/dXbeu/6A6oDlKm18vfiKQu1oezfkd8JIZQ7K6/FoR3jWUubsr5vrbU7bSmC5VgcuuUkXR5y5KwHond3D7nmb5jOYEB2YLIN+YLse3uG0Tdsr4T3u7WZjndPDuGQWWHHhhIuF4ndRvvvq3nKm/Sq2Fcc/bEk3nnrc7jfmcS2s2sQX/7jq/jx07sA2BlqOv/027/ioVcOuZ4B+Xx1DcYhhMDtT+/Co5ol9pq2D5kGPZxII+w31TnRkYMoyYUntCLoMzHizF8aTqSxqDmCtCWwt8dtdUhBv+XJHXj91x/H41sOq/MnEzpe3N2DW/+8E1ff+rxyP/aNJPHDJ3fgTzksrKQTV/QZhJf29KiBT1NNAIuaI5jbYA+w5L0phMiyKluigSw3dSWoqnAQ0aVEtI2IdhDRTTm+DxLRnc73LxLR4kq2R3iUQ3a+hraQky4c533zCTyzo1sJR9ybjutxVekPwUAspYLvXmR9K6+ryhUcT2Syt/R5HPKmyRXj2NU1iKWtNbj56lMRDfnwmxf2qd8I+Q34DMrrqgJsq0XWd2qO2g/iYS2VWPqEa4M+NffEW2zvmR3dWoXSJIaTaUQCJurD2ZbAob4YOo6NYF5DRLnG5P8DtHRcra033LEeb/z2nwFkfOxSOBJFhCPh1BkDbNeGFzuDLHMNzlranLXNrNqQ632xzCMd/ZzLtqt4giNAI4k0rvjBX/DRX7+MRzYfwvYjg0rED/fHURMwUetYj7mEYzhhWw4D8RQ++Is1+NbD2wDYq0juPprtkpHXNJc1Itly0B7pJ9MW7l1/QD0rcqZ+Mi3wzGvdWLv3mOpw5V9prezUOv0jAzH8cUMnnth62HW+9UHIQDzlihMCtnAETANz6kLK4hhJpBEOuIVDpp5fdOIsAMBlJ8+BzyCcv7wVIb+BeDKtrKS3nTYPpkGuuAqQsWi2O1lOO44MqnMkBVyf6f2Mk8It43i5gtnyNz9w7mIAwE/+YovqndefhW9efSqaagII+gxV2sZbzPLGNx+P7/3taWogUUmqJhxEZAK4BcBlAFYAeDcRrfBs9g8AjgkhlgH4HoBvVrJN3nXFZedruooc6qNC+++R/rjjqvKUHPF7XFUei6MxEsjZDjny8Lpu9IdIBkXrw/7MPI50JoXYKxz9sRT2Hh3GsllR1If9OGlOHXZ1D7pSeINaJlauTKXhRCZ9VvqopV8dsG98g+xAeyJlZY22awKma4TYM2QLSyTgU8KxsCmi/m47NICRZBrzPBZH30gS331kGzY7nY7scPXAZM9QQo3KB+MppC2BlCXyCkciZamU5bb6EA72xbICjtL1IgcEr1+SLRxNNQHoz2058zh04djndDCyg5JzJl7YfRRDiTSe2t6F63+1Dhv297o6xZqgT2XIDcZTWLOnBy9qPve9mjhINw1gB153d2fHKmQHt6ugcNgWx/ce3Y4b7liv1nDvHUmq8+0tunjAIyBHhxI4NpTAur2ZOUOH+uMqjmXvI+OK2t8znDXJsWsgjrkNIcxvDONAr32cwwl7YKKfox9f144/3/gGXH3mfLzxhFZ8512rcP8/n4+WaBBhv+2qkvfuqfPrceWpbbjjpX3ojyURT6Xxq+f3YF+P+3zEkpZLONKWwMObD+FvTmmDaZDyFsjzn8sSlUL4usVNWLWgAceGk6gP+7F8di3mNoRBRJjXEFZWqNcVefzsWhBVXjSA6locqwHsEELsEkIkANwB4CrPNlcB+KXz+m4AF1MFz0z+4HjG9aD7n5tVgM/OrEqmLOXbD/gyAWo9ZiHTNgdiqbzCkUhZWdYG4LU47Nd1Yb97yVrnN70xjg37e5GyBFbOrQcALG2twe7uISWEmRRexzrKZXEk0mqSX7Ma5brnXESDPpUKLDudT77peFxy0mx85MLjVJA86DNUxxz2m2iI2MKxekkTnrvpIlx2yhwljouaIi7h2HZ4AD94YocaecU8VYIB4A8vd6hrNqQt1pTLJSbPnZz4tXqJ7YKS7oX+WBK7u4fUiPfjb1yG/7z61JzuD9MgZY1FAiYGtIygYugpwd96eBv292RSMeV8kKe2dSHkN/Ch85aobWuDfnVfRYM+1VkPxlK47mcv4W9vfwFr9th+dl04pFD944XHYUlLjeu7VNpC33AyIxx5AuA+g9A1EEf3YBxPbrPny6Qt2yrtG0niuFY7Xdnr089VRmfroQFcfevzuOEOOx35UN/g6FsWAAAgAElEQVSIKwlEn2uii4jO/MYIZtUFVbxNWhyztGvVEg1icUsN3rJyDn7+gdWIBHw4YY6dQBHymxhJptXck+ZoEO9evRBDiTRe3NWDJ7d24fP3bsZDrxxyDRCODMRUksGR/jj2HB1C92ACF57QimjQp6o4yAFBLku0azDjavzAOYsBZJeOmdsQVqKrV4YAbOGYKKopHPMA7Nfedzif5dxGCJEC0Acge5g3TnhT6ZSrinK7quY7o2MAWcHxoGlm1apKpDOukoFYCuGA6XLBSJJpkXM2tDs4nrE4MjEOodqsryUCAGudjuOUebZwLGmpQe9wEkecUaffWT9EuapyZGcMa+seyCwnPfOjZyiBhkhALWAls6wuPmkWfnJdu+qQAagUXMDuYKVw1IZ8mNsQxmzN5XPWcc3qXALZ6ZGyrXpnJAOULdEg0lq9K28wW4l7KlPzp31RI4gyMZLzbn4Cb/z2n1Wg+8xFjXjX6xZknR+J7KQWN9sZV6W6q+RxnLmoEX95rRt/3NipXFXHhhIQQuCp7V04e2kzPnfFCnzlqpVq/9LK0C2OIW2OzDcf3AoArpFyyhJ431mLcNNlJ2JeQxgDsZQ6T794bg9O++ojSnxlp1jrEV55Tbcc7FdWxfbDAzjty49ACOD42Y5w7OtFTcDExSfOUp9JZAf8spMFJznUF3NZHJs7M0H4fGXjFzZHnI46hZRTOy7i96E1mrmfwlqihZeQ30AsaanRf0s0gFULGhAwDazd06POnyWAy05pwwfOXYw5dSHs6xlGZ98IiOzswH2OCC9tqUE06MNALIV7/tqBl5yVLnPNppe/2VobxN+c2gYgk7WnH99ex40q45wG2X3VAq0/qjTVFI5cloN30F/KNvaGRNcT0VoiWtvV1TWqBnldQ8pVZeR2VaUtCy3RIH7/0XPgNw27yGGBGIee1ZO27E6+JuB+EA2Caz863gmAgOOqUjGOdFaMozbkg88g7Dk6jPqwH/Mb7dTXpc4NKcsVyBIphV1VadW5SovDm0JbH/bb65CkM66qFmcEvnJundpOZlIB9oNcHw447fW7vp9dF0Q06HMJrDduEkum8V+Pv+aa65D5bXu/cuTsdVVJwYqn0uqcttYGsbSlBpsO9KLjWGaioNyn/D/5kJaIPMdH+nPn2CdSFr736HaVDSUtgF9+cDVaogHsOzqMA8dGYBqElCXwmuNHP3eZPTfhdYvtTnvroQGsaLPPbchvoMZJrR6Mp9XqkNsPD0AIgT2eOIa8NvOc+0KOZl892O9KWZbWiDwmKUgyTvDUtswz95fXupUbqX1RE0J+AwPxFE5sq8NP3/86nL6g0dWGJS01MAhZs9f7YylXqRo9CJ9vpvuipgiiQb89udQ5n5GAibpw5rrr1quXsN9ELJXGgd4RGGTfhyG/iVPn12PNnh5lMQDA8llRfPHKlTh+Ti1e3nsMQtj3+HAirdoqhezIQAz/eucGZZXJ+3HNnh51/TNiFYTfNPDcTRfhzo+c7T5XzTU4NpzEj57aiZudwcC72hfgw+cvnZDYhqSawtEBQB+2zQfgTTVQ2xCRD0A9gJx5e0KI24UQ7UKI9tbW1lE2yd0hBTSLQ3ZW+mp8qbTAaQsacOaiRls4UpmSIwGfHWwmchct1DuuoN9Uo0NJTcDn3o+pB4Xd6bg+gxAOmEhqtaq8MY6AaeBtp9uGXGttUPlAl7TYo75tznwDv0moC/mx/fCgs65GHleVJ8Zx2CMcDRE/As65OObpaPXR/hzN4gj7M8FxWThRTki86MTZ6tjyEUta+M6j2/GeH2fy9PXRG6ALh7vTl21LpCwlHEG/iRVz67H10IBrvQ0ZMM/nYpS0KqG0rbt9PblnHr+87xj+3+Ov4aJvP4Vk2lJiHfIZWNgUwfr9vRiIp1RZExkfkoJxguOaWL2kCe99/UIAwJo9xzLB8VhKJWH0x1LoHkzgUF/MJXwttfaxzGtwC8cBj1WXSFtoqgkoN9wCJzV0xdw6zK4L4k6tGoE+mKgL+3CqM/9ADlq84j2rNoTGSMCVWXXuMtuxIAuC1of9rtIfW/JMPFzUHEE05MNgIqXWKgkHTJfvv5BwBP2mXTetewhz6kLqeWpf3IRNB/pUMBzInLNZtUHlVn3jCbaQrtnTY2dzRYOIhnxZ53MkmcYfXu7ANbc9j288YAtA10AcdSEfQk775jaElbDrxwcA33lkO57abovQNe0LcONbTsh7TJWgmsKxBsByIlpCRAEA1wK4z7PNfQCuc16/E8ATooLF/b171oVDTQD0VMOVKX36zHHTIJgGgYhcAWfvPIKgz1AjN0lN0OdyeUW1h0yf+yGDfn6DMhaHlo6rxzr++aLlAGwXjGRBYxg+g1TQ028Y+OB5i7HlYD/ufrlDtVnP/BpJZkbltSE7luG1OOqkxeFkvuhpyfo5lQ8dAEQCPperCgDOW9aCr161El+4ws6XOGVePf7urIW46rS5KAVpHchOXAmHJtRXnzEfV51mi2o8ZakFfMJ+E3MbQjjcH3O5RGTHVlQ4HLE6eZ5tBeQTDtmZDMZTeHyLXRLDZxB8poFFzTVqEqGc/fzIq4cR9ptY4VhuhkF49qaL8NPr2nHxSbbAvnnFbIT8Bgyy/eMjyTRWOR33rq5BDMZTKgEBgJqPIy0OOXegM0etrtl1IXV9ls2OgshOYjiprQ4DsRTa6kMI+gx1XF+6cgUuOnE2znTuOykcUY9wtNQG0RwNqODyfZ84Fx97wzIAGYt4cXPE5YrsHU4i1wB7QVMEtUEfhMjMw/A+Y95UdR3Zae84PKjmTQDA65c0IZkWWLOnR7kiZbmZ2c7cnbDfVJl2a3b3YGFTBESEaNDnWilStvuTd9m10qQrs3swkTNuprOkxbb49HMhLcyJpGrC4cQsPgHgYQBbANwlhNhMRF8horc6m/0UQDMR7QDwSQBZKbvjSZarSsUJkHMeRzItlHmoZ1XpN2bQZ7qFQ+u4gj4jy+KIBE3V6QLuYO5gPIXO3hGs/tpj2HSgD5GADz7TcNWqCngsDp9JWNgcwWOfvBCfvyKTtOYzDeUOkpMG33baPCyfFcV96zuVxSEtC0BmVdntCvt9iATs0ZlpkBrF2a4qe631oXgqyycuxetNK2arz8IBEw1ht1ViGoT3nb1Y+aN9poH/eNspOHFOnTovT/37G/DpS09ELqRQtDgPopy4p492P/nm4/GRC5aqcydFMew3MacuhGRa4NXOfnWNdx4ZRMhvFPSRA5kYR1t9GM01gfzCoXUmO7uGEEtaquPSO/fzl9sW9Ib9vVi1oN6V+DCvIYzakB+mQXjly2/BD99zhuqspDV42oIG9RtD8RRao0E1IGh1LI6WmiACpt3ppy3hil01OqI+py6ozt+72hfgvo+fh/mNEXVN3rJyDhoifjUX4bpzFiPgM3DGQikcEeca2PuTRkBrNIimmoDqDBsjAXVvbj88CNMgJWzymAEo60dnUXONEibpAvIKh1HApSPFbdvhAddvvn5pEwKmASGAd5wxH4998gKc7hyXTME+sa1WtW0okcZCxzqIBn2qD/jyW1fiu+86zfWb0pPQNRAvKhwLmiLwpgd53d0TQVXncQghHhBCHC+EOE4I8TXnsy8IIe5zXseEENcIIZYJIVYLIXZVtD35XFVajCNl6a4qS92EkYBpV4nVRv1yH3EteOtyVflMdVPLDrYm4EPKEnjHfz8HwC0clrB9u0cG4li/vxeRoGlnc2luNK/FIS2iZbOiWSIl3UN6FeCT2uqwt2dIuU2aaoKq49BdVZGAiYjTydWGfOo4GsIZV9VgLJU1urzlPWfg81eswGkLGvCOM+Y57ci2OPIhfdXN0QAWNdfk3V4KvYxxSAtEb0/Eb8JnGjCd+StKOAK2cAB2+XUZP9h6aMBVPj0fF580G9e+bgEWNUewoCmSlbopOXBsBC3RIGbXBbG7ewixVKaUy+IWu9NZ1BxRo0wgdwqwRGa0ydeywu/xc2oR9pvY2TWIoXgKNUEfGhyrSbpCDIMwtyGEjt4RdA/GXRlqUsTm1IdUp98Q9uOU+bYrbpXz9/JT2lAnv4/4lXvo3GXNeOeZ8/EGZ6VEea8vcvbbWht0iUBDxK9cmQd6R1Ab8rkGMNedswiA3dG+9NmLse5zl7jOgXxmZGwp7HSsxe4t+1ga1Ov5jW6rWCYCLGyKYNmsTAaTHCisnFuHRc0RrHKEeq5zDPoz/LrFTVjQlNnvirY6VymXVs88IC8hv4m5nnvQ+1xPBDxzXMPrqpKWg6llVbnXGBeqY26JBnF0KOEa9ct9xFMWUml7so47xpEJjssaSd7RkdxejuhlsNCe/2CqGeuybd4YhzctV0cKhy50i5oj6OyNqRzxC5a3qDiDnlUVDphq5L2wKaJe14f9TlaVvQiUN/21sSaAfzhvCYgI33rnKvzuH8/Gslm1OHlePU6ZV6/SIvMhOybZ4YVy+Kv1cyhHcMdyxDgijokfdMR9RHNVzdZiMDKoH09ZavZuIRY0RXDz1afCbxpY1BwpaHHMawxjUXMN9h4dQiyZVm69hU22WJw8r97VaX5QS8MtRE3Qp1Kf60J+LG2twY4jgxiMp1ET9KnBgN5hz2sMo7N3RGWt3freM/DkjW9QIjOnLqzuR12A37JyDv7v4+di9ZImdR83aO68SMCHb1+zSomu/L/Hz67F195+Mq4+Y55KbfcZpDr/TNzLjyZtf9eutuM5QZ+BWbWhLMtD7l92yPJ+ePyTF+Lej59b8LzNqQ8pIdDdqQCU8C30ZC/NcgYZK+fWg4jw4fOXuM6Bfq4aIn6Xq/P1S5vQNWDPhO8aiCvXaiGWzYqiRrvHvX3GRFCyVBHReQCWCyF+TkStAKJCiN2Va9rE461Vpcc4ci3klExbKmuluSaAgVhK+fUlUjikGa53XAHTUJ1XXdiPA70jWaMH+aDOqgti79FhVzA6EvDBZxiurKqgRzB8eWanA9kWB2B3emlLYFf3IAI+A5+5/CTEkmn8/uUOVyn3cMBEJJDpAGThwYaIX83jGIilCo6GTINUoHd+YwR//Kfz8m7rPR/SktCzrQyyrbL6sD+TNuw8iNLiiARMGGRfU92tl0hZquRHyJ+xOADgpLZMNph3tFeMhU0R/GnjQSS1yZmSA70jWNFWh2jQhye2HcGs2pA6nqUtNfAZhDMXNqIm6MOnLz0RFx7fmnOGfS4aIwGV3lob8qGtPqQGBNGgiYZIwM7A0jqdeQ1hPLH1iHKhLW2NYklLxqqbUx9EyG8iYBqqowdsa0W6w2RyQ2OBzDO5v4aIH+99vW09yFiLbqmcMq8ez+zoRl3Yh0bn9wI+A3UhP37/0XNcg7Bnb7oos39pcQxk5gkBdgc/q6648J86vwGPbTnsinEAdkzsUF8M7YvdWWGr5tfj05eeiCtX2fG3vzmlDb6/M3CeE5vSn4H6sN/lWppTF8KIU3l4MJ4q6qoCgK9edTL6RpK48ofPACgcs6kUJQkHEX0RQDuAEwD8HIAfwK8BFJbvKU7AtG84wyBYwp7n4RUO2Rc0OR1ZR8+wawQQ9JmIJzPBbv0mCvoN1Fj2+3rHBeMdPUj/+qxaWzj09ETb4qCMxZHOLjniM0qwODRxkaOpbYcGENLqXhlkxzhSzjmJ+DNzPo6fHXVlvyhXVTzl6oDHAzmilaPMkBZ49zm/Wx/2q9G2FA4ZgAw6wXqfSaqDCvoM/M8Le9HmtDUcMNHg84PItkKXtNQg7EwMm9tQvnCkLYH9PcNY2pqZv2BZAgd6R/CmFbPREPGja20cR4fiyoJqrAnggRvOV3NBPvqG48r63dn1IeVWrQ35UBf249XOfowkbYujKRJwZdkBdursXWs7VPaWtK70FOnzlrXgjIWNLotCR14fmV6dC2mF6iNv+fzo+z19YQOe2dENn2Eoq0sK3ZmL3J23bh0Ui3EUY9X8ekc43Ne6sSaAz13hLW5h33f69SEiXHryHPVeCpnfJFdbVs6tU0IhJzSWIhwydqL/3kRTqsXxdgCnA3gZAIQQnUQ0cdMUJ4h8WVWmc2Es4Z6dbAuHvY0cMW05NICztRpGck6Daz1wZ4Qb9JkgZ6qK7MRlWYvzl7egezCBk+fW4+HNh1UA7ohmcdQEfPCZpDqIZEqoUbR3ImAu6p1Rod/jqgLsQKq8iYkIkYAPw4k00j57nz7TUBVgj59diycCR5zjCKhjHoxnxzjGSp2yOHK4qpxLo3dILSqrylmC1jn/ekpmwGcHPeVM9JDPPr7mmiC6B+0yFk01AVWptxxOX2iPxJ/fddQlHN1DcSRSFuY1hNV53nZowLXNWGYCz9ZrM4X8qAv5cdgZdNQEfPjERcuySqq8acVsmPcQ/rihE4ubI0owZFxpTn0IPrPwRDPpSixscTjioqcFO8Kg/z957vYcHVLXNFJCIDjqsThK+T8673n9QtRH/OpZGCvyGagPZ6yph//lAsxtCKkZ9XJeSinCMRko1cZJOGmwAgCIqKbI9lOSfK4q2femLeFZKlao76TrJJGyXKPSoM8umpar+KE9ATATG5C/AQDvO2sRHrzhfDUSkzeU1+LwGYaqWptIW/D7MlleQGmuKv24Z9eGlNjobqBwwMRIIo2RRErFM2Tu+vGzaz1ZVbYwDuWIcYyV5pogfAYpy0hvo7w2+jyF2pAPQZ+hZrHLcvcRLYUxlXZfd1kVWZYnn1MfRmONvc9yLY7jWqNY2BTB4876HRKZsjq3Iaw6qGPDyZyVBEaDPk9GxgvkvVUT9OHkefV4gzPnQNJYE8A5x9mDns9cfpL6vMGxHtrqih+7FJnGmvwWR2s0iIAzV0UiLQqXxeFMFOwdTmYsjhJST+U8lkxwvDyLozkaxN+fvXjcRvLSy6C7GU+YU4vakF8913LCYCkxDkk1YhuSUp/qu4joRwAaiOjDAD4I4CeVa1Z1yLI4tGwjwO5gvXX5pSvIFWTUAqhBv4n+kaQSDrkW+YDzuiUaRG3Qp0bGN77lBLx6sF+lq8o2yBvMHeMwVVqlTAWW7rVMjCN/RyTdCiOJzDEZBmF+Uxi7uoZcbqBIwLQtDicoD2RiA231ITWqa4jI4LidDDDeFkdjTQAP/csFWOx0tvocEYkUDnJKMUSDPlUkMOj46PWHTl9kSmdOXQh7u4cRDfrUiLeU4LgOEeGiE2fhty/tw0gijaFECl+7f4tytcxrCLs6+VCO4xkNui+/NuRzdVqFOt9/uWQ5Vi9uwpu1dOlr2udjcXPEZSHkQ8+qykd9xI9nPvVG1+Q2+fw0aO1srAngXe3z8cYTZikhKsV6kMfnDY5Xi9ocwiGRQvGqU06lHIvj+ZsuzlkaaCIo6akWQnybiN4EoB92nOMLQohHK9qyKuCdWSir28o4gyWEy1UFIBMcj2ZGSvqoVPr7XVVznc4h6Dfx1lVzcdFJs/Bbp2xzSzSA9521KPP/HctB5rXr7oVI0KeEIWVZ7rLqcgJggZz1eiUc7iqbZy9txq6uIVfqcdhvC0ciZakR3EM3nI/O3hiISLmM6sN+1IbsXP6UJbLmcYwHcuIVgJwjdDlqDfvtGcPRkE8FxwOmXVqllA7oylVzVZE+GQwu1+IAgHOXteAXz+3Bqwf78avn9+D/1ndig1MeZV5jGHUhnxOrEjmzxEaDdFUFTAMhv6kGCUD+Qo8AcOaiJpy5yL3OSEs0iMtOaSvpd+XvFJsk6Q1Sy/PrtVT+852rAGSWESjF4vCZhopJ+c3sRc0mmqhKBsg+J42RAEyDsLNrCAa5500Voz7iRz1KS5YYb0o6o0T0TSHEo0KIfxdC3CiEeJSIKlrivBpkuapMd4zD66oCMm6s2qBPjf5driq/neqpLymrL7ZkGHapD282lLcN3gWCADtALbdPpkTu4HgJrqphT3mRd545H4B7DYZIwMRIMoV1+47hZFVhN6oyR6T1EwmYrgKG4+2q8pKro5XHJd1ndfpsfb/dqZQiaFedNk+5bFbOrcfKuXWj6tjnaKL/rFMyfE/3EGodFxIRqZFncJxdVbWh7NFuJfP+S4lx5KI+7Me5y5qxenH24lgAVEHQUuMVUU8aezWRWVS5LA5Dc7s2R4MTWm9qLJR6B70JwKc9n12W47OpTZ55HHKSn2W5S44AUMFxIkJzTRCH+mOuAKpKx1VLypquGIfk5Hn1OKmtLmvEIQOJs2qDrvIlgG1xSLF6wam6KTOkMsJRPKvK66KTqZW62RwJ+PDyvmMYTqRxnlNkT+fKVXMxy8nScQlHnqVax4tcHa10d4Q98SPAFuIvvXWly33xn+88FbFkGl+4d3Pe3/nwBUvxYWeWebnIelB/3nZEuU8sAdfM5Na6EDr7YuNncThiJTvQOs1lWMmZxrL8RikTJXUMg/CbD51VZN8hlyurELVBH7oG4mUHxitBLvHWufD4VuzuHnJdo8lOwZYS0UcBfAzAUiLaqH1VC+DZSjasGuSbOS4H7ZYnHRfIWByA7a46PBBzVX6V/n69+GFGODKdRPviJjx4w/lZbTrnuGb89sNn4aS2OtSGfIhr5ZhlyQ8A+Miv1rnaXIqrKt9DSER46f+72BUcPGV+vVrF7Nzl2cKxekmTmlmrW1zRCtfRydXRSv+6sji0yqhBn6Hmjkje1W7X2pzXEM673OpYkBl3L+6y63O21gbRNRB3DTBac6QXj4WQUzhSdVqR0mIcY+XMRY340z+dh5Od8v3jyX+/94yiLjBJtEhnPZFEC8Q4AHti4S+e21NwsazJRjGJ+18ADwL4Btx1ogaEEIVXl5+C5CurLi2OdI4Yh6l1rs3RIGbVBl0TAMN+n4oNyH0qV1UJbgnDIJztZLpEgz61QA1gC4d33Q7vzPFSLI5ceH3QN1y8HLc9tRN1IX/RlFRdOKPByj64Yb+JmoCpMrwALcbhsTgCplEwU0YWChxv7ElrPtUxnHNcM+5d3+myOGbVyfTi8fPHz64Lqgwj3V1XSfchEVVENIBMteFSkG7Wt5/hXeJn4qkP+3Hlqrm44PjsAReQWYJ4vOc8VZKCd5AQog/24knvBgAimgUgBCBKRFEhxL7KN3HiyJdVJQPglpXL4sg86O9qn4+Dve4MnWjIh8F4StV+cgXHy5zxGdXmMHQPxtEQCag0U4ksF5EpOZK/o6wrYzQW8pv46+fflHOBKS8BJ1usezA+7llVXvymgYf/9QLc8uRO/PYl+3bMsjhC2aVVJpqW2iD6Yyk0RPw4qa3OFg5NgPWSG+PFp95yoitpQVKN2kYTjV07LoX3OOXmq4lhEP7r3afn/T7kN/HbD59V9hyhalLqzPErAXwXwFwARwAsgl3RdmXlmjbxeF1VcrRuuiyO/K6qK07NLvkt/ZZy6c98rqpSkCPF1Usace3rFuLcZS14YNNB9f1ZS5vUg5JvCVl328vrpPLNFs7F3IaQLRwT0EnNb4y4XGJy3oG0OMoRyErREg1iV9cQ5taH1Wxw3eKQ4uZdR3ssXKKl1NpzfuzJopMhYFxp7rj+bPSNJF2W1mRGehWmCqUOwf4DwFkAtgshlgC4GNMwxuF1VUn0mePeEXehEs1AJjDWPeCegAaMwuIIypG0Dxcc3wrTcKca/uIDq1X6qJrHUaDkiOSSk2YV3aZc2jxZPZVGX63Rm1EjhSPX4lQThYxhzG0I4ezjmvH20+fhnOMyrgsZh9EX6xpPiAh1YT9qAmbRe3Y6sKSlRiV5MONPqU91UghxlIgMIjKEEE9Ox3TcLF+Vg3SLW5ZA73AS0aBPBVGLuRZkVtRRx6UU9Bk5s6pKQXbCenBTd0WFPGU0vN/nYtfXL8+q7z8eqEqoE+QWyUx8zKwN4o1xpPKNDCYAWVmgrT6M+rAf3/tb95oMcmSsr7E93tSH/RiaAaLBVJ5Sn+peIooCeBrAb4joCIDKDI2qSL5uRbmqLIG+kSSaowElHEaRXld29nJx+qDP1CyO0bmq9BIK+YLfQRUcL9y+So0+LzyhFXuODk3YrF29jLxcWCoT46i+T1/Okm7LM/P8whNacdGJs3DTZbkXphoP6kI+sGww40GpT9RVAGIA/hXAe2Gv/f2VSjWqWlh5RqRSOOJOxdfls6PYe9Qu8FeyxTGoxzjkzPHRBccj/sxly5duW0p13EryxhNmqfWXJwKv++8fLzwOr3PKX0+GlExZUiNfWfZIwIefvf91FW1DyxSaYMZMbkotOTIEAERUB+CPFW1RFclncUirQi5Hqi8qU+xBVDEOWWRvDK4qaXHorippcXj35TdLc1VNF7wLV91wyXL13WQIjsuCid5S3RPJF69ciUS6enEeZvpQalbVR2BbGCMALAAEu58d3VTaSUqeEEeWcOh1qcoJjpsGwTRo1K4qtRqgy1Vl/763cwz77SyamZB6CQDBAllkk8HiuGB5K2597xlZ60hMJN51HBhmtJTaq9wIYKUQoruSjak23lpVEtkXyTkTTTWZUhzFXFUq6BnLlBhvjgZd60OXivz/etzA77iivH78cMDEnR85u+hSrNOFQskAkyEl02caJRcKZJjJTqnCsRNA7oWTZwDS4pCztptq/Fnf5SPoM1TlUznB572vX4g3nTS7bH+zCo5rMQ5ZyiSXO6aao9uJRi1clcOKq+bEP4aZjpQqHJ8B8BwRvQhATVUWQvxzRVpVJfK5qmQHr2IcmsVRrPMnItSG/OgZSij/dshvjsptkFnsJiMS8vf1dbFnIpkaXTMjpsMw1aRU4fgRgCcAbIId45iW5HNVZcU4akoPjgN2bKJnKFFwyc1SOHNRI37y9+2u0tOnL2zEf7/3DFx04sRlME1GvMFxL1+6cgUWtUzLhSsZZsIpVThSQohPVrQlk4C8WVWOOBwdsgPcerXVUoRDCtJYM2qIyFVGQnI5+85ViZV8bqn3n7tkIpvDMNOaUp2/TxLR9UTURkRN8t9of9T5/48S0WvO35zOeCJKE9F65999o/29UsnrqkCXsCoAABBZSURBVNIsjoaw31XY0Cxh2nXvsD0beH4jZ7VUimIWB8Mw40epFsd7nL+f0T4bSzruTQAeF0LcTEQ3Oe9zLQo1IoQ4LcfnFSGvq8rpi44OJlAf8UM3MkqxOGT9oWrm8E93WDgYZuIodQLgeNv5VwF4g/P6lwD+jEm8mqCKcQwnsLilxpVJVU5m1FhjHEx+MllVHBxnmEpTbAXAi4QQTxDRO3J9L4T4wyh/d7YQ4qCzj4POOh+5CBHRWth1sW4WQvxfgbZeD+B6AFi4cHQ1+IVmcbx79QL1WoqDEPZksnKF429OacP9mw5Oiolo0xVZ5DDAFgfDVJxiFseFsLOprszxnQCQVziI6DEAc3J89dmSWwcsFEJ0EtFSAE8Q0SYhxM5cGwohbgdwOwC0t7ePqgyqLFW142uXuYoH6kIR9psusShFOL5/7Wn45jtPHU2TmBJhVxXDTBzFVgD8ovPyK0KI3fp3RFTQfSWEuCTfd0R0mIjaHGujDfbiULn20en83UVEfwZwOuzJiBVBGhze5UV1cQj4DNfiTaUEx/2mwR1ahVHCwZP9GKbilPqU/T7HZ3eP4XfvA3Cd8/o6APd6NyCiRiIKOq9bAJwL4NUx/GZR5AqAXinQjQrvutVcbXRyoGIcLNAMU3GKxThOhL08bL0nzlEHe+3x0XIzgLuI6B8A7ANwjfN77QD+UQjxIQAnAfgREVmwBe5mIURFhcNSFof7c91VFfAZLiuDhWNyoGaOs8XBMBWnWIzjBABXAGiAO84xAODDo/1RIcRR2MvPej9fC+BDzuvnAJwy2t8YZcMAlOKqyryfCctwTgV8BiHkN1ATmBnVgBmmmhSLcdwL4F4iOlsI8fwEtalqCGRbG0C2cOjbFKuOy0wMRITffvgsLG7msiIMU2lKtevfTkR1ROQnoseJqJuI/q6iLasClhA5q916Yxwui6MSC3Yzo+L0hY1o1OqIMQxTGUoVjjcLIfphu606ABwP4N8r1qoqIUR2YBxwi0PQb4x6AiDDMMx0oFThkDPXLgfwWyFET4XaU1VKclWZbuFgVxXDMDONUiOJfySirbCXjv0YEbUCiFWuWdXBEiIrMA54LA4OjjMMM8MpyeIQQtwE4GwA7UKIJOzVAK+qZMOqQh5XlTc4bnBwnGGYGUxB4SCiT2lvLxFCpAFACDEEYFqt/gfkd1V553EYHBxnGGYGU8ziuFZ7/RnPd5eOc1uqjsiXVaWdpYBpcnCcYZgZTTHhoDyvc72f8lj5XFUFZo6zq4phmJlGMeEQeV7nej/lESJ71jiQI8ahnTUOjjMMM9MollW1ioj6YQ/Ew85rOO/HUqtqUiIgcsY4dDHxpuOWUh2XYRhmOlGs5Ig5UQ2ZDOSbAJhVq0oXDpOFg2GYmQWXEtUQeeZxmJ55HLp7ii0OhmFmGiwcGgLuulQSPaYh132Q23FWFcMwMw0WDo1SZo7L9R6kYLBwMAwz02Dh0Mgb48ghHFJg2FXFMMxMg4VDw545nmsCoDurCrAFg4jTcRmGmXmwcGjYwfHC2+iuKrY2GIaZibBwaORzVelkXFVsbTAMMzNh4dAQonjRQuWqMojLjTAMMyNh4dCwSnBVyRiIQeyqYhhmZsLCoSFQeuVGg4hdVQzDzEhKXQFwRpCvyGEuTAMQgoWDYZiZR1UsDiK6hog2E5FFRO0FtruUiLYR0Q4iuqnS7Solq0rCFgfDMDOVarmqXgHwDgBP59uAiEwAtwC4DMAKAO8mohWVbFS+FQBzwTEOhmFmKlVxVQkhtgBF3UKrAewQQuxytr0D9jrnr1awXSUvBculRhiGmalM5uD4PAD7tfcdzmcVI98KgLkwiMWDYZiZScUsDiJ6DMCcHF99Vghxbym7yPFZ3lUHieh6ANcDwMKFC0tqY66dlxocNwyCMe3WQGQYhilOxYRDCHHJGHfRAWCB9n4+gM4Cv3c7gNsBoL29fVRdernB8em36jrDMExxJrOrag2A5US0hIgCAK4FcF8lf7CUkiMSk4PjDMPMUKqVjvt2IuoAcDaA+4noYefzuUT0AAAIIVIAPgHgYQBbANwlhNhcyXbZa46X7qriGAfDMDORamVV3QPgnhyfdwK4XHv/AIAHJq5duVcAlCxqjqjXBhWva8UwDDMd4ZnjGpYQoDzOqi1fudS1hKxpUMnWCcMwzHSChUPDLjmS+7twwHS9J+LquAzDzEwmc3B8wiknHdckXjaWYZiZCVscGkKIkrOqQn6Tg+MMw8xIWDg0CrmqvHzxypUFA+kMwzDTFRYODYHSM6VOmFNb2cYwDMNMUjjGoVHKCoAMwzAzHRYOjXJmjjMMw8xUWDg0ysmqYhiGmamwcGiUU+SQYRhmpsLCocGuKoZhmOKwcGgIlL4CIMMwzEyFhUPDskqfx8EwDDNTYeHQEMhf5JBhGIaxYeHQEAIc5GAYhikCC4dGsfU4GIZhGBYOF+yqYhiGKQ4Lh0Y5RQ4ZhmFmKiwcGpbgdFyGYZhisHBo2CVHqt0KhmGYyQ0Lh4YQ1W4BwzDM5IeFQ0Owq4phGKYoLBwa7KpiGIYpTlWEg4iuIaLNRGQRUXuB7fYQ0SYiWk9EayvdLi5yyDAMU5xqLR37CoB3APhRCdu+UQjRXeH2AOAihwzDMKVQFeEQQmwBJt+iSVzkkGEYpjiTPcYhADxCROuI6PqJ+DF2VjEMwxSmYhYHET0GYE6Orz4rhLi3xN2cK4ToJKJZAB4loq1CiKfz/N71AK4HgIULF46qzXZW1aj+K8MwzIyhYsIhhLhkHPbR6fw9QkT3AFgNIKdwCCFuB3A7ALS3t49qRgaXHGEYhinOpHVVEVENEdXK1wDeDDuoXjG4yCHDMExxqpWO+3Yi6gBwNoD7iehh5/O5RPSAs9lsAM8Q0QYALwG4XwjxUCXbJQRgTFopZRiGmRxUK6vqHgD35Pi8E8DlzutdAFZNZLsswRYHwzBMMXh8rSEATqpiGIYpAguHjgBPAGQYhikCC4eG7apiGIZhCsHCocFFDhmGYYrDwqEh2FXFMAxTFBYODXZVMQzDFIeFQ0MIcFYVwzBMEVg4PLCrimEYpjAsHBrsqmIYhikOC4cGFzlkGIYpDguHBq8AyDAMUxwWDg2LLQ6GYZiisHBoCC5WxTAMUxQWDhe8AiDDMEwxWDg02FXFMAxTHBYODcHrcTAMwxSFhUNDAOyqYhiGKQILh4ZlCRD7qhiGYQrCwqEhqt0AhmGYKQALhw6XVWcYhikKC4eGJQRnVTEMwxSBhUODq6ozDMMUh4VDQwjA4LQqhmGYgrBwaLxl5WycOKe22s1gGIaZ1FRFOIjoW0S0lYg2EtE9RNSQZ7tLiWgbEe0gopsq3a7vX3s63nHG/Er/DMMwzJSmWhbHowBOFkKcCmA7gM94NyAiE8AtAC4DsALAu4loxYS2kmEYhsmiKsIhhHhECJFy3r4AINcwfzWAHUKIXUKIBIA7AFw1UW1kGIZhcjMZYhwfBPBgjs/nAdivve9wPssJEV1PRGuJaG1XV9c4N5FhGIaR+Cq1YyJ6DMCcHF99Vghxr7PNZwGkAPwm1y5yfJZ3crcQ4nYAtwNAe3s7TwJnGIapEBUTDiHEJYW+J6LrAFwB4GIhRK6OvgPAAu39fACd49dChmEYZjRUK6vqUgCfBvBWIcRwns3WAFhOREuIKADgWgD3TVQbGYZhmNxUK8bxQwC1AB4lovVEdBsAENFcInoAAJzg+ScAPAxgC4C7hBCbq9RehmEYxqFirqpCCCGW5fm8E8Dl2vsHADwwUe1iGIZhikO5wwtTGyLqArB3FP+1BUD3ODenWvCxTD6my3EAfCyTlbEcyyIhRGspG05L4RgtRLRWCNFe7XaMB3wsk4/pchwAH8tkZaKOZTLM42AYhmGmECwcDMMwTFmwcLi5vdoNGEf4WCYf0+U4AD6WycqEHAvHOBiGYZiyYIuDYRiGKQsWDkz8uh/jARHtIaJNzgTKtc5nTUT0KBG95vxtdD4nIvqBc3wbieiMKrf9Z0R0hIhe0T4ru+1EdJ2z/WtOCZvJcixfIqIDzrVZT0SXa999xjmWbUT0Fu3zqt+DRLSAiJ4koi1EtJmIbnA+n1LXpsBxTLnrQkQhInqJiDY4x/Jl5/MlRPSic37vdKprgIiCzvsdzveLix3jqBBCzOh/AEwAOwEsBRAAsAHAimq3q4R27wHQ4vnsPwHc5Ly+CcA3ndeXw65ATADOAvBildt+AYAzALwy2rYDaAKwy/nb6LxunCTH8iUAN+bYdoVzfwUBLHHuO3Oy3IMA2gCc4byuhb1Wzoqpdm0KHMeUuy7OuY06r/0AXnTO9V0ArnU+vw3AR53XHwNwm/P6WgB3FjrG0baLLY7pte7HVQB+6bz+JYC3aZ//j7B5AUADEbVVo4EAIIR4GkCP5+Ny2/4WAI8KIXqEEMdgLw52aeVb7ybPseTjKgB3CCHiQojdAHbAvv8mxT0ohDgohHjZeT0Au9TPPEyxa1PgOPIxaa+Lc24Hnbd+558AcBGAu53PvddEXqu7AVxMRIT8xzgqWDjKXPdjEiEAPEJE64joeuez2UKIg4D98ACY5Xw+FY6x3LZP9mP6hOO++Zl07WAKHYvj4jgd9gh3yl4bz3EAU/C6EJFJROsBHIEtwjsB9IrMYnh6u1Sbne/7ADRjnI+FhaPMdT8mEecKIc6AvbTux4noggLbTtVjBPK3fTIf060AjgNwGoCDAL7jfD4ljoWIogB+D+BfhBD9hTbN8dmkOZ4cxzElr4sQIi2EOA320hKrAZyUazPn74QcCwvHFF33Q9gFISGEOALgHtg31GHpgnL+HnE2nwrHWG7bJ+0xCSEOOw+7BeDHyLgEJv2xEJEfdmf7GyHEH5yPp9y1yXUcU/m6AIAQohfAn2HHOBqISBap1dul2ux8Xw/blTqux8LCMQXX/SCiGiKqla8BvBnAK7DbLTNYrgNwr/P6PgB/72TBnAWgT7oeJhHltv1hAG8mokbH5fBm57Oq44kfvR32tQHsY7nWyXxZAmA5gJcwSe5Bxxf+UwBbhBDf1b6aUtcm33FMxetCRK1E1OC8DgO4BHbM5kkA73Q2814Tea3eCeAJYUfH8x3j6JjIDIHJ+g92dsh22L7Dz1a7PSW0dynsDIkNADbLNsP2ZT4O4DXnb5PIZGbc4hzfJgDtVW7/b2G7CpKwR0L/MJq2w16vfofz7wOT6Fh+5bR1o/PAtmnbf9Y5lm0ALptM9yCA82C7LzYCWO/8u3yqXZsCxzHlrguAUwH81WnzKwC+4Hy+FHbHvwPA7wAEnc9DzvsdzvdLix3jaP7xzHGGYRimLNhVxTAMw5QFCwfDMAxTFiwcDMMwTFmwcDAMwzBlwcLBMAzDlIWv+CYMw4wFIpLprAAwB0AaQJfzfrWw6yAxzJSB03EZZgIhoi8BGBRCfLvabWGY0cKuKoZhGKYsWDgYhmGYsmDhYBiGYcqChYNhGIYpCxYOhmEYpixYOBiGYZiy4HRchmEYpizY4mAYhmHKgoWDYRiGKQsWDoZhGKYsWDgYhmGYsmDhYBiGYcqChYNhGIYpCxYOhmEYpixYOBiGYZiy+P8BTCAp66z/PywAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mean = numpy.zeros(300)\n", "T = numpy.zeros(300)\n", "Tval = 10\n", "i = 0\n", "while Tval < 3010: \n", " E = numpy.random.normal(0, 4, Tval + 500)\n", " X = .5 * E[500:Tval + 500] + E[499:Tval + 499] \n", " mean[i] = statistics.mean(X)\n", " T[i] = Tval\n", " Tval = Tval + 10\n", " i = i + 1\n", " \n", "import matplotlib.pyplot as plt\n", "plt.plot(T, mean)\n", "plt.ylabel('Estimate')\n", "plt.xlabel('T')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Generating Non-Stationary AR(1) Data:\n", "$$X_t = 1 X_{t - 1} + e_t$$\n", "$$e_t \\sim N(0, 4)$$\n", "Or, alternatively, \n", "$$X_t = \\sum_{j = 0}^\\infty e_{t-j} + lim_{j \\rightarrow \\infty} X_{t-j}$$\n", "$$ (\\sum X_t)/T \\not\\rightarrow E[X_t] $$" ] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEKCAYAAAA8QgPpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXe8JFWZ//95qqq7b/fNk/MMYQhDhkFBQVRAwQTmuLCGxYBf3fW7X4V1v6vr6n7dNaGrqyLLgj/T4qoLCAoMIiMShzCEyTnP3JmbQ6eq8/uj6lSdqq7qru7pcO+d5/163dftrq6uOlXd/TzniYeEEGAYhmGYatFaPQCGYRhmasIKhGEYhqkJViAMwzBMTbACYRiGYWqCFQjDMAxTE6xAGIZhmJpgBcIwDMPUBCsQhmEYpiZYgTAMwzA1YbR6AI1k1qxZYtmyZa0eBsMwzJTi6aefPiyEmF1pv2mtQJYtW4Y1a9a0ehgMwzBTCiLaGWc/dmExDMMwNcEKhGEYhqkJViAMwzBMTbACYRiGYWqCFQjDMAxTE6xAGIZhmJpgBcIwDMPUBCsQhmGaTtG0cMdTu2FavKT2VKZlCoSIFhPRQ0S0noheIqJPO9tnENEDRLTZ+d/rbCci+g4RbSGi54no3FaNnWGYo+OpHQP47K+ex7O7Blo9FOYoaKUFUgTwv4UQpwK4AMD1RLQCwA0AHhRCLAfwoPMcAK4EsNz5uw7A95s/ZIZh6kHBtAAA+aLV4pEwR0PLFIgQYr8Q4hnn8QiA9QAWArgKwO3ObrcDuNp5fBWAHwubxwH0ENH8Jg+bYZg6YAnbdVVgF9aUZlLEQIhoGYBzADwBYK4QYj9gKxkAc5zdFgLYrbxtj7ONYZgphqM/YFpsgUxlWq5AiKgDwK8A/LUQYrjcriHbSqYvRHQdEa0hojV9fX31GibDMHVEWiBFky2QqUxLFQgRJWArj58KIX7tbD4oXVPO/0PO9j0AFitvXwRgX/CYQoibhRArhRArZ8+u2I2YYZgWID1XRXZhTWlamYVFAP4DwHohxDeVl+4CcK3z+FoAdyrbr3GysS4AMCRdXQzDTC1cC4QVyJSmlRbIKwH8BYDXEtFzzt8bAHwVwOVEtBnA5c5zALgXwDYAWwD8CMAnWjBmhmHqgHAUSJwYyKNbDmM8X2z0kJgaaNmCUkKIRxAe1wCAS0P2FwCub+igGIZpCk4WLwoVYiBD4wW8/z+ewP976xl4z8uWNGFkTDW0PIjOMMyxh+VaIOUVSLZoQghgLG82Y1hThrFcEbli6+8JKxCGYZpO3BiIfJ3Tff184D+ewNd+v7HVw2AFwjBM85F1IEWzvGIwTQ62h3FoOIdDI7lWD4MVCMMwzSeuC6vgWB5cL+LHtMSkaETJCoRhmKYTtw5ECkm2QPwUWYEwDHOs4lWil3dhScuDYyB+LCEmhVJlBcIwTNMRMYPorgXCLiwftgur9UqVFQjDME3HdWFVUAxFGQOZBLPtyYRlCUwGncoKhGGYphM3jde04gXbW8Ha3YN49w8fi1WPseHAsGt11YMiWyAMwxyrSH1QSQgW3SB664VlkBf2DuGJ7f3oH8uX3e+PGw/hipv+hF+u2VO3c5tCTAq3HisQhmGajpyNV2plMpljICIkFVkIgRf3Dvn22354DACwbn+51Sqqw7KEa8W1ElYgDMM0HSuma6o4idN45YhU4+iWP23Hm/7tETy1o9/dppHd8q+ebjiTs7AYhjlWiV8HYjn/Wy8sg7hKULEEpOI4rFSJa5qtQOplMViWgBCT456wAmEYpulUWwdSqLBfK/DiOJ4gH3PazrenvEbnjv6omwIxY1bxNwNWIAzDNB0RInzDmMxZWFIhqIphNGsrkLaE7m7THRdWvfIAJtM9YQXCMEzTqbYbb6P9/fsGJ/DSvqHKOyqEKcHRnK1AVKXixkDqZYFMorgQKxCGYZqO6SqQCt14rXiurqPl6/dvxPU/faaq94Q1hHQViLLNjYHUSeDLe1ev4x0NrEAYhmk6Xjv3Ct14zeZUoh8czqKvyvbobhaWYlmM5eyiQnW8uiNl62WBWGyBMAxzLBM3jbdZ/v4jo3mM5U3ki/EtnXIWiBnmwqqXBcIxEIZhjmWk7Cu0MAaybt8wntxup90ecarJByfKV5WrSB0Rll1lmqoFQr79jxZWIAzDTHt+/cweLLvhHhweLXUNebP31sVAblq1CV+46yVYlnDbkQyNF2K/37OiSl9rqAUSMwGhGbACYRimIfzXU7sBAJsOjpS8FreVSbGBs+1s0cJ4vojhbME9/uBEFQokkIWlNktUxysVSN3qQCbROvGsQBiGaQiyFiJXKBV0YUV4QQ4OZ2MH0S/75sP49qrNVY2vaFrIFkwcHvXcVgMVGiOqCPjrQMbyXldevwKBb7+jhV1YDMNMe1KGLV6yhdJ255XqQAbG8rj4Xx7C7188YO9XwVLZcmgU31q1qarxFUwLE3nT1033aCyQIeW9qnAn1wKpaniRsAJhGGbaIy2QbMh6Gd6CUuFumIHxPPKmhf1DE/Z+DXDXFEyBbNHCESVGU00MxO3G6/xX3xvs0BvcVo5b/rTN7eAbRtwizGbACoRhGADAXWv34c7n9tbteG0JaYGUCv9KQjXvKJYJxy3UiNl2wbSQL1roUxTIwHh8F5YVKOiLskDC6kWimMib+PI963H32n2R+0idOxnauRuVd2EY5ljgJ4/vBARw1dkL63K8tLRAanBhyXqMCee9lYLttSDdYnsHbSunO52onwtL1GaByNUNy61yOJmW+WUFwjAMgPrP8qULazxfvQsr5ygQqTjKja3Wlh4yQL93YALd6QRmdSTxyzW7Masjhc9cflLk+w6NZPHQhkMlzRRVoV/0ubBQ8Rok8rrDEg8k0psnhH3tslVKK2ipC4uIbiWiQ0T0orJtBhE9QESbnf+9znYiou8Q0RYiep6Izm3dyBlm+lG06rtIkRRsY051tkpcC0QdWxS1jrngSOL9Q1nMaE+iO51AwRT4zoPls7k+fNsafO5XL6Bv2HZ9SR2YU8asKjX5MI7HSSqOXJmKeNW6UR8/uP4gVm/qq3ySOtLqGMhtAK4IbLsBwINCiOUAHnSeA8CVAJY7f9cB+H6TxsgwxwSWJeraoE/OuMMUSKVZeVCBlKt5qNVyKhTt9/WN5NDZZmDjAa9eRZSR9rsHxu0HzsRfCnF1zD4LBP5gezniuLDU61Uff/j2Nbjm1icrnqOetFSBCCFWA+gPbL4KwO3O49sBXK1s/7GweRxADxHNb85IGWb6U7REXd1YMsYwmouOgUgX1XDWH3sIzsDLpfEWFOVSzfilC+vwaA7tSQOnLex2X8uXqXyXgX0ZHZdKN1/BAokT9HZdWOUsEOXYrY6DtNoCCWOuEGI/ADj/5zjbFwLYrey3x9nmg4iuI6I1RLSmr6+55hzDTGUsS9Q1s0cGe0dzpYFptZXJ2t2DOPOL9+N3L+x3Xw8K8HKCUu07JVcEjINUION5Ex1tBn70FyvxoVceZx8nROlJpHAPVsmrY/bHQPzZWuWIEwOJskBawWRUIFGERYpK7p4Q4mYhxEohxMrZs2c3YVgMMz0oWlZ9LRDXhVUqjKWsLVoCL+0bBgD8caM34St1YcWzQOSKgHFQM7s6Uga6MwmcMr/TGXPl47gFfSEurDBFHOfWStdVOQtIPTYrkFIOSteU8/+Qs30PgMXKfosARCdLMwxTFaYl6rZmBeBZBiOhMRAniG4KJJ2KdXXd86ACKcSMgYzGEPwStTixPWVnjHU4a5nHOY58f5gLS3W5hbV9j8JzYZVL41VdWK3thzUZFchdAK51Hl8L4E5l+zVONtYFAIakq4thmKPHFPUNokuhXy4Ly7Q8BZJTFEhQgMqU1TBUYT0S0wIRQgQskAQAoN1RIOMxXGElFohpucWT/joQ+3+sGEihsgvLF19pcT/FVqfx/hzAYwBOJqI9RPRhAF8FcDkRbQZwufMcAO4FsA3AFgA/AvCJFgyZYerOkdEcvnLPuoYv21oJ06yzBVImC8tbD8RCUre90+oMPmxhp6g4SLEGCyR4rA7HAmlP6s5xoi2A4DFUCySpa9A18mWNVRdEl1lY4d+Fv/3lWrfLsT2G1n5nWlpIKIR4b8RLl4bsKwBc39gRMUzz+dPmw/jRn7bjbecuwqnzu8ruO54v4i//8yn801Wn4+R5nXUdhykEyKpfUZqbhRViFUhhqsrUci4sINoFZNYQAykElLW0POT/qmIgliwktJA0NGSLlm+NkLiV6JsOjrgWVJQL648bD/ksp1bHQLgSnWFajGz1EWc51Z1HxvHk9n48u2ug/grEEqDSvJSakbPjsXwRQgi3Ky3gVxyyV5ZPgYRYYwXLQhp66XZTtUDitSIJtkaRsY+qYiCySt45lGuBkN8C8VxY0ccaGi/gdd9a7T6PskCKlnDbuwCtVyCTMQbCMMcUUiCUy/2XDDv9lsLagxwt9Q6iSwFrCfiEnr3NO0+YAg21QCJqQVQhGjcGErRAOgIWyHgVFojrwjJtC8TQyG+BoLIFkjP99ycqBmKawndvWIEwzDFO1m1fUVkpSAEZFMj1oFjnSnRfbCIg2NXTTIQokDBlGhUDUZVB7BhIQBlJxZFxYiBjMRS0VLbyf8FxYWlHEQORRH0XgveACwkZ5hhHzsDLZd5IZMV2OR/91+/biOd2D1Y9jrqn8SrCLSiQ/RaIfd15RaiHubCiAsZmGUUVRYkF0mYrkJRjQUQpIrXFSVghoWuBhGVhlevnFVBoUdZo0OIIs0CaaZWwAmGYFpONUTwmkRZIlAuraFr47kNbcPX3/lz1OMw6tzJRhXSwpbsqiL2W7YoFEqJMo9qZqLPww8raHnHHBnguLCJCe8qIdGHlfO6j0jqQpC4tkNI6kHK3NnjfoxRIsB4m7PMKXlsjYQXCMC0mV5ULy7ZAJiIUiKqEthwarWocZoOaKQKlCkSVg6ExkBAhGKXcpGJJGlrZlfxUgkF06cICbGUSlcarKm4viK4oEDcGolggcvxlrLugK8q0RElat2WJko6+YS6sOLG0esEKhGFajFQG8VxY9sw4queTeox7Xyits/3ne9fjN8/uCX2vKerrwipYwl0XvfogeqkAj3Jhye2nzOvE1r6xsp103bFFWCCAXZUe5SJUtweD6DnTQtLQoRGF9sIqN66wbsO/eGo3dh0Zd5+HKYuwuApbIAxThm19ozjji/dhd/945Z2nANkKxWMqlSwQ9RhDIavr3fXcPvxhQ2mTUTm7rWddmmlZ6HRiC0HlqMpCV4FUqAOJLCR0LIGT5nZiNFfEoZHKbqySOpCklx6cSRqRCjoshdZngegaDJ18llycBaXCVlz8+/95Ee+5+bGS86mEufXipIPXC1YgzJTj50/uwki2iHtCZthTETeIHsOFFbRAsgUTy264B7f9eXvJMcKOlzet0CVmiwFhGOS2P2/Hp37+bMXx+Y5pCndmXz4G4tSBRLiwZPlIpRjIyXPtupitfZVdd6oySid0GLonCjtSRiwLxAuiO2MumkgZdh1ImAVSToFEvTacVc8Xz63HFgjDlEH6p1W3w1QmG6P/kUTWgUgLRD7/4t3r7GMoQjgbcrx8MVyBqL2pwlwtT+8axBPbj1Qcn0rREm5sIY4LK6soPHUW3Wbo7tjCkNtlYeVvn9+PgbF82bFJZZU0NF/8A5AurHBlrlp+rgtL6YUl03jV6/PSeKPHE2VdSQtOPZ9KmMJnBcIwZZCzwOmjQKpxYfmzsILvUZVQmKLIm1Z4hlPIAkj+45pVZ2iZlmqBVHZhFZQiOfW6Uk6DwkoxkAU9aQDAz57Yha/fv7Hs2ArOADpThtsHS9KeMiLTeNV0ZHleN41XurA08llL8lG5OpCoPmiqAglTMmGxEw6iM0wZ5I/b0OvXt6mVVOfC8leiq+/ZOzjhex4U2kLYAjqsCLHSIkXZohXqpy9HwfRiICVZWKoLSxHKsgtumAUS6cKSWVi6hm+9+ywA9jK1gL1O+FM7goueehbI+y9Yio9cfLzvtUxSj+zGq24P9sKSWVga+S2QOC6saAskUXI+3/tC7km1n9PRwAqEmXJIBVJumdOphFtIV5UFImMg3nvW7Oj3zT6DCkkKljDLJKxuwT9Gs+puwaoFUtaFVVQVSGlGlrRAItN4nVm4rhPees4irJjf5R7/w7evwTt/8Fjke648fR4+cMFS32sJXYv8bvnSeIMuLJnGqwdjIKXXHES9Nk2ZF3VVsEDCjslBdIYpg6w2jlN4Vy1fuPNF/PixHXU/bjlqycIKs0C2Hx5zn7cn9RJFIe9XNsTSqWSB5Aqm6/aJS8EUSCcNENnvV/HXgXhPpGL0ubDkglMV2rknHMmbSmgV76Wsek+EWLEJXYtcwGo8JAYi/xecxbE0qr6QUFUOXWnP6vBZICFKLUypNDMGMj2cyMwxxViIm6Ne3P7YTgDANRcuq/uxo4gbA7ED4PYsN1e0l59V4xmHRnLu8+50osSFJe9XWHDdp0BCLRCrBgvEQkIntBk6soFrs4SARrZQVau+ZfBanRy0JXT3eOHnscerSwViVFYg8loSeukcOhjDUFGD82oQXQhhB9H1MoWEZVuZeOPtaktgcNyeKMjFtoD4WVhsgTBMGaQF0syZViOZCLEmwpDWx9yuFAB7tq4Kyr6RnPu8K50osUDk/QpP41Wa/4XGQExYonw/p9JjCugaoS2hldStCGGnzwL+JW/VyYF0f0kLJEqoS9ecodn7JQ29ogKR98KIUiAR2WiHRrK+6wNsIS4VntdMsbxLMIjfAvHm9VYFy1BuU8faCMs8ClYgzJRDCpzpokDk7LxSGq+sCZjX1QbAdqdIZTC7M+UoEPt5dzpR4qryLJCQNF7l1OEuLKdWo4pKw6IpkNA1tCVK3WmWsN1bgL+2YjznxUCkApEWSFSgWVomMqkiZWglLrMg5VxYUqmE3Qe1SNFzYXn3Vlog/iB69Die3N6PbCDDrTPlua38659Hu7DU93MaL8NEYFpeqmczs00ahWWFp66GscPp83SSUzA3njfd9yzqTfsskO50wqeQNh0cwYFhe/ZcMEubJqoWSKgLqyjrToo4OJwteT0M07FA0gk9NIguW6dbwisWzLn3wlQ65HoKRAiBL/92HdYq3Ybl90B1YQVn4UG3juvC0kIsEEepqAJ7y6ER3Lx6Kw4N59DtxChUF5arQAx7SduwQkLAb1EcHM7iXT98DJd+42Hf+FQLRP0sylkg6kvswmImNT97Yhfuf+lA2X02HRwJTZ88WoaV9hyN/KFEtQqpN+WypoJsOjgCADhrcQ8A6cKy37O4N4O+kZw70w+6sP7y1ifxDaU2olxabZiRIfe/adUmvPuHpVlNYRQsCwmNkErooXUgGaV9SIdjjeSKdraXJbw6H5mF9amfP4tHtx7BLY9sx4MbDrnvlUJUxjNShl5izUn3nzs2qUCMUhEolYo6k//VM3vxz/duwLbDo5jf3ebbP+jC0rXwViaAXyHIzsF7Bydw19p97vauiMB5eB2IPxMsOO5GwwqEqZpvP7gJ3394a+TrQgi87lurQ9Mnj5aBcS+IWc0PRQiBj9y+Bqs3lfaBkqg/+iNj8dqCHy2qIK+kEDcdHMWczpQrwGwXlv2exTPSyJuWW//Q1ZZwXWNCCPSN5nzul6ACUYVT0AIRQrjnOTicxeHR8lXegNdbS9c0tCW0EuUohEBaUSCyGjxXtFxhHIyBAMAX73rJ3k8Zv7QmZPqrTDJQZ/7S/ff0zn48ub1fiZuEubAcC0QR3jJ4ni1YmNsVUCCqBaKXtjJR5f6b/+0R7DxiW5LDE57rTrXq1Cws3+dSJoiuWidsgTCTllzRxMHhHNbvH47Mynlml+deGA7M/CxLHNUMqVYFki1YWLX+ID5421OR+6huj8e2HsGegcY3a1RdO7mihdFcEV+/b2NoIdvmQyM4aW6nO3O3XVj2+xf1ZgAAewYmQGRXMOeLFixLIOcUAaoCK5gVpQrLYKBctZLUc5ajoMQl0gm9xKKzlCA6YLcPAWzhJwWgp0C8/ZbP7bDHryoQS8DQyF1zPWVoyBdN3+cpLZC3f/8xvOuHj3kWSFgQXZdpw977+5Xsq3kBBaK6IaUF4s/C8h5vODCCHziTL/W3oX4PfBaIMgb1M5KKL1iLAvgX5mo0rECYihwayeKmVZuwZkc/9g5MALAF8raItRd+uWa3+3iz43aR/PV/PYfln/9dzWPpH/N+dNXEQKQwKbsutSIo/89/P4/Xfv3hGkYYn219o7ju/1vjO/89z+/Ddx/agm/cv8m3r2UJbD44iuVzO5BJeut2S1fNol67jceegQmkDM0NPOeKllt8qLr/SgV6tK9ddQdNFMzQGEoQ+bqhkR1EL4YE0RUF0pHyXFhSGMsq9raEJ6akwPcv7CR8XQlkHYj6/VCVJ2ArQo28uImKFM7qNaoTl7khLqxcOQUSuFXtzuenrt+edT6PN545H5etmIPPXH4S2hIa1K+4ekyZ3iuVvWqcsAXCTCre/6MncNOqzbjlT9uxx1EgAPDi3qHQ/V/cN4TjZ7UDsN0ukqJpub7esEygOBxRVpyTSqFoWrhp1SYMjZe2L3f3jfGjCu7T6HTIb9y/CS/uHQZgC8tcwXQDtMG1PPYOTmCiYOKkuZ2uABrNFZEtmtA1ct1auwfGkTJ0V+hmC6Y7+1avpxoXlir8J5RK8Se2HcE/37s+9Nrk8QzddmGFxUA0jVxB6LqwChZ2O5af3JZUrAQpdNXxF0zhpvACQEq3FYjq5gpawuv3D4daH4Ayu1ekd1kLRJTGQPxL2vrvZ8a5LlWhjztj/fJVp+O0Bd341KXLcfqCbr8FEqJAwiwQjoEwZXl21wCu+t6fSwKDR4PMbnl214Bv+6HhLDY7K9uN5oruj5sIrvALHmd73xguWj4LmaSOjQc8C+RPmw+7j+MuPSopmhYe2nAIR5wf8sz2pCvwH9/Wj5tWbcbf/c8Lke+P86MKKoxkhICpF7M7U+7j7nTCFnrONe0fyvrW85Dxi3ndbZjbnYKuEXYeGUeuYKHN0DC7wxZqg+MFnwWSLZq+ma4k6IZS3VYlLqyC34UF2AL83Tc/jptXb/PtPzRRwMOb+lzhKy2Q0joQu5CwzRGEbQkdGtnfiw/dtgazO1N49cmzAQCzOlJ427kLAXiuqODSsn4LxHPxSYYDa6Os3T0Y+flKxTKcLWDfoD1hGlAmJ/O6U7791czAlB7mwvIj1x7xubCcsarXoQUKGtVjJtxUY8eyZgXCxOWBdQexdvcg1uwcqLxzTA4O53DLI9tx53P7fNtfcKyMzjYD/WN57BmYQEInnLagy80KUukbyWEsb+KE2R1YPrfTt899SuZWnEV/VP64sQ8fvO0pPLyxDx0pA51thvtDkSmoz5a5H3EskEJgn+McK6pRqMFhqUBUt8uft3gKVwqb7nQCKUPHspkZbDo4glzRQiqhoyttuDPnVEJTLBArtLNs0CKIa4GEdQEeV2b6//uO53DtrU+6gld3FEiJwhICGpGr6AyNkDJ0bDs8hqGJAr7w5hVYNtO+/wmd8M13nY0F3W2uK8pngTgxEIm8r6riVO+f/bwY2YxTbv/Bw9tw9ff+DNMSGBzPY0F3G2a2JzGnMxhER9kYSLCQ0LNAiuhIGU6SgRMzUiypYEW7qhjkNcpNvhgIu7CYcmxwZvXlBKaKEAKPbj1cdmYiBb3MENl4YARD4wW8sHcIRMDFy2dhYDyP3f3jWNCTxtKZ7dg7OFFyHBkXOX52O5bNzPhcXs/tHsTM9iQAr1uqyot7h/A3//Ucfv/iAbzzB4/6Zq1S4azdM4gZ7Um7X5FzPXKmvm8o62QMmXjzvz2C1379j65FFccdFdyn1pnc3sGJstleEnWxoK62hM//D9iBfHdf5xplgPXkeZ2OArEXMSIi9Dr3NmXobuBZdWGpBF1Y5XphqftOuL2qTLd2Y1S5jl3OKpHy83VbmYS5sIjcFN2EriGV0NzPsiNluG4aXabnJvRwC8QUvliGFK6q4gxzfUW7sOztB4eyGBwvYHiiAEsAH7n4eDx642t9MRnAnxiSdBaUKhcDkUMdzhbQ1Wb4LCH1OoKusHIWiC8GwhYIU471+23X0bNKMZVkJFvAl3+7zmcef+uBTXjfj57AXQHrQsVVIP3jyBZMXP29P+PfH96CF/YM4YTZHVjcm0H/WB67ByawqDeNRT1p7B2Y8Lkv7lq7z02zPG5WOzJJw52xTuRNbD40iktPnQMgXIHctXYffvPsXnzsJ0/jqR0DeHGfF2ORQcxc0cLMDluB7Oofxzlfuh9PbvfqTfYMTGDPwARe2DuEbYfH3NeqiYF84IIlOGVep5sZ86+/34Abfx3tHgvygVuewDW3Pol80cKvn9mD139rdWgLEFWwt6cMFEzhCrizFnXj0a2HUTQt3Lx6q5ee6xSZLZ/TiZ39467LCgBmZKQC0XwxkOEQF1awsM83Yw7cKl+zQ6VvV5iglsHwwQn789I1DemkFlpISOS1ak/ohKSuuT2g0gkdczpTeMd5i3Dh8TMB2C7FMEVgZ2F5oizpjsu7v8MTBTcjTaJRhAXiCPHxQhFFy0K/892b2ZFEytDdbC+J6sKS3Xj9Foj/+PK7MDxRQFc64et3pVpSQQtEtRI18lqu2OdgCyQWRHQFEW0koi1EdEOrx9NsBsfz2D+URVLX8NyuwRLB9Nvn9+OWR7bjX3+/AYDdofU7f9gCANh0aARCCPzPs3sxOJ7HP9z5IrYcshWHjFXs6Z/A83uGMFEwsevIOF7aN4wzFnajtz2JXNHC1kOjWNybwaJeu+7gwQ2HsG6frdA+9fNnXetoQXcamaTuzljX7R+CaQm89pQ5IPIUyL0v7MclX3sIRdMqyYh5SQnSq03sZrankDA0vLh3GAPjBfz0iV3uay/uHcLQhLev9F2rs7JvPrAJN63a5AuMAt6s9vIV8/Dy42a4yu/f/7gVP39yl09RlUNe29a+Ufx5yxFsPDjiWiX3PO8Fx1XBvqDHdovIth6XnDQbW/vG8MC6g/jnezfgbif5QLVAhABe2jfsuoF62+3XUobmCuZc0fJeJecgAAAgAElEQVRZCJKgRVCumaIqrOVLuYLlzpxVBSID3zKhQVogZiB9WwQsECNggaST9jKzX3/nWThxjp26m0pobhsbdfzFYAzEuXa/C6voKrEV87sAwK3MDyKPNZ63+38dcepeeh0FHVQ8phDud6XN0Eu68YpAFEQK/ZFsEV1tCdea0MiOe0j0MjEQQ/MH67mVSQyISAfwPQBXAlgB4L1EtKK1o2ou6xzr441nzsdIroi1e/xWiCxI+sP6Q7j1ke24w0mp7c0ksPngKF7cO4y//q/n8A93voQfP7YT9zxvxyU2OYHyvGnht8/bwmr74TEcGM7i+Fnt7ux2NFe0LRCn7uBjP3ka/3i3bXXImVRStxvKZZI6xgsmhBB4fo+tDM5Z0ouZ7UnXJfXFu17CziPj2HBgBMMTBcxsT2LtF16Hme1J91oBuLNAAJjVkUQywn89kiu6s1jAUzxqfOM7D27GTas240d/2ua7b4POORI6oS3ptd+QsZAflCmeVFk2y743Gw4MY/th+75++Z51uObWJ3H9z55xrcPhiQIuXj4LL/3j611/vxSQ5x83AwDw1A7bBbe1b8wXHD/JqYfYOzjhWSCqCyuhurDCFEhtLixJrmi651AVlMwQG3QUgYyBBI9juUF0zwJJGbp7b9QUX4kaM1JjKrIOJLifqtiGJwqukL9o+aySY6tIgS5dqDLhQ97fYOavZQm38HRmR9K2HMr0wpL3dzhbQFfacM9nBNqqBGMpPgtEI/t1U+AHD291f+dAc1v8TLV27i8DsEUIsQ0AiOgXAK4CsK5VAyqYFnQijBdMJHXNFaKrN/WhK53AvK42fPV367FvKItrLlyKN525AID95cwWTNdvDdhupK/dtxFb++x2Cde/5kQMTxRwxenz3X3kbP+Trz0R97ywH3c+tw/nLOl1X5cZU/uGsvjSb+3bcu6SHiyekcGaHQN4dKsdmJVpojv7x5wagxGcNLcDmw6O4ldP7wEAbHTcWotnZHzrRi+ekXHrDkxLYO/gBIazBeSLFq69cCk++MrjANizSCHs2eLDm/qwsCeNuV1tmNWRcmfpy2a249BIDk/t6MdwtoiudALd6QRWLOjCS/s8BaIqhZkdSezq9//YZKpotmC6+2aSuqt4wvzCh50xZAsmXv7PD7qCJ2VoSCd05J2W6dIykp1YV607iMGJAl532lzkixZmdfizcuZ3p/Hi3mFs2D+CHUfsmMD96w66r0/kTXS1JTCcLWBhbxrtKcOdiUthLzO09g7a7x/NFX1ZW3L5Vnu8jgUiXViBIHqcGEixnAIJcYnkiuUtEPkZGJqGNifraKJguutbmFYwiK4hZWiusG0LUSCqq8dngZiWP43XuXZVsRUszz243LFoonBdWI4Ckd/V3vZoC+TwaB4pQ0NHyrC78apL2opwC2Q4W8DJczvdho5BC9zQtEAMxLtmXYPbc+t/nt3r+34304U11RTIQgC7led7ALy83icZzhbw8Z88DY3s6lad7B+sJQTaEjp2HRmHphH2D01gw/4RzOxIon8sj7ldbXjzWQtweCSHXz69BzPak2hP6Tg0nMO87jb87S/X4uXHzcTh0Rw+/pOnMTBecJfg7E4n8G9/2IKndwzgouWz8MeNfXjfj54AAPzu0xfjVMfsfm73IBb2pHHC7A5cduoc3L12Hz7/xlPdWcyWg6N41Umz8cYz5uHJ7QP41TN7cMXp81AwBe58bp8ryOSXeNeRcewdnMB43sTlK+Zi08FRd91n+d1d1JuG+ptZ1JvGwl5PgB0YymLnYVvQXXD8TCxzZuwZRwisPzCMhzf14dOXLgfgdI51ZnXSXbBm5wBGs0V3BbYVC7pw6yPbMTieR08m6XM3zWhPlQRAF3Snse3wGMbzpjsDO25WOwbG8vjzlsMl7irAc2/J2JBbDKbrbrX3RMELbEvhfvPqbdjVP45/+u06DE0UsOOrb/QdV1boP77tiHteVYZIQTY8UXRdUlIYj2QLSOqau11NVOhWWlxkkoaz9KrpCkw5QzY0tZAwngUSLCRcvakPY7kirjxjfmhn22zBjIiB2OeV6daGRm6qbjbvd2ERkXuMhDL5AsIViFqRXq6QMKnrvnFpZAtfqRAySQNnLuqOFLRGwAKRCkRa4cHQiWUJHB7JYVZHCkRUYoEEYyCuBTJhT5hcCyRgVZezQHTnPJYQsITfPdjMIPpUUyBhfgvfx0NE1wG4DgCWLFlS00mEZft47Q/H9rE+s2sQCZ0wmitiQU8aOhFmdaTw4YuOw84j45jVmcTj2/px8+ptaDM0XHbqXDy08RCGJgq446MXoDeTxGXffBhXfvtPODyaQ28mgfakjg/dZlciJw0NEMA1Fy7F379pBZ7eOYBV6w/iP/+8HT9+bAf+39vOBGArkLMWdwMA3nneYtz7wgF85Z71eHHvEL73/nOx/fAYXn3KbLz7/CV445kLMLcrhXetXIw1jivk6Z0D6EgZ7o9rZ/+4G0C/5KQ5+K+nbOvj6rMX4JZHtgOwLQ615fai3gwySQMz2m3BXrQEnnGynRbPyLj7yYrp2x/dAZ0I73uZ/Xn0ZJL40+bDuPybD7sf3pod/VjQk3b7AL3m5Dm4efU2vPYbD+Nnf/Vy170E2C6soAKZ29WGbYfHMJE3MS6KIAKWzszgjxv78P5bniiZ3Z00twNDE3k8t3sQ331oi++1pGOBALYQkYJGzmj3DIxH+s8BL0C9do+/0PKcJT14dtegO3uWLgzAq10YyRaRNDT3PuxVstjU5U3t+5DCrv7xEgskVzR9biNVwCecAG+5NN6iZeGaW58EAOz46hujLRCpQBQLR16HdB0aOrnbVLeT68JKqC4s7zNV+2S5x1ZdWCWFhP5KdMBT+OmEjqIp3PYwmaSOO69/ZUkwXCItAimIByfy0MgbU/C7ZDq9xmY5FqJW0o3Xf/yiKWBZAiMyC8uQLqxSBeLrkqwqEMeFVbTsrgCFYmuC6FNNgewBsFh5vgiAL7VICHEzgJsBYOXKlTU5A7szCfz3x19R9fuEEL4v5a+f2YOkoeG8pbY/+8MXHYcntvfjY5ccj6vPWQidCM/sGoBpCXz8p8/AtARef/o8AMB5S3tx3tJeDI7ncceaPVi3fwQHhiZwcDiHay5cCgB49cmzsWJ+F257dAcAe2acNy0sn2O3++5IGfjsFacAAM5c3I10QocpBK5/zYn4l99vQCapo28k57bGPmV+J1Z/9tVIGTpWb+rDLY9sR9LQMLsj5f545XMAOGVeJ9bvtwPZj2+zU06XzPQUiPzBbdg/giUzMpjjVPCevbgHd6/dh82HRt0f48HhHEwLeLnj+7/g+Jn47f+6CB+67Sl8+LY1ODSSRWebgZFsETPbU0ga/h/bjPak48YyMeFUdM9sT7mzTvXHl07oOGF2BzYfGsVHbl+DIEm1EK9gurO7kWwRBdMqqzzs9/h/wGcs7MYLe4fwsuNmOArEdl/mi5ZraXi1CwUkdEJ70i6sUwvY1CZ7gK1Id/WPl1gg43nTm/UXLF9GXlLXIPSwGEh4+jBQuhwtADd9GIBrsQJehpFngWiudaVaDbIORB7DcGIgkjbDP0EA/C6sXNHCgaEsHt16GHsHJ3zKNVgHkk7aQXyp2NPJ0kwqlaCCyBUs38JTJS4sCzg8msdCJxHCCHTjDdaBmJaFsXwRloDPAtHDYiBqN16zVIGYpj3JLUx2FxYRXQRguRDiP4loNoAOIcT2xg0tlKcALCei4wDsBfAeAO9r8hgiCX4p33buIt/zz7+xNN5/6alzAQDvOX8x/rixD+cq8QwAuOGKU5FJGli7exAHh21T+uzFve75PnflKfg/v1yLgmnhJ4/by7GetqCr5DxzOtuw9guvc83ew6M5zGhP4mv3bcQD6w9hfnebr4nbPKc1xqKeNDSN0NWWgEbecwC4+ZqV2HhgGG///mN4fNsR9GQSvmNIN9DBkSyWKpbJB1+xDMfPascHb3sKpiWwfI4tzA+P5tz+R/Z1dOOrbzvTbYD4tnPmY9/QBE5f2FVigXRnEkgnbJfO0EQBPemEL76ksqg3jd72JPY5rrsbrjwFj209goed2o2ETl6/KcUCyZsWdh4ZD3VJ6Brh6Z39WLd/BNmC7Q58cns/hiYKuOL0edhyaBTnOZ+tGtjucvs9+S0QIkJXOuGL/aj3FoAbe5ECU17vWK4YGUSXxw72plK9Hr9+Zq/vuKFB9ILl+nKe2N6PgrkJn750ueu66XeCyrrmZVqprhW7lQmQ0rwYSNJ1Z1HoSoGqBZI3LVzx7dXu/XmZM/GQ1wh4abzppI6iJVyXVFiAXiX43coWLXe9dSDchTUwnseZC23PQLAbb5CiJdzPpLPNcC2eoAUSdIWpkyCNyI2RqGnE9vEnmQIhoi8AWAngZAD/CSAB4CcAXtm4oZUihCgS0ScB3AdAB3CrEOKlZo6hUfzjW05DPiSVtTuTwP99k614fveCnaJ75qJu9/VLTpqNJ/7uUvzj3etw26M7cPHyWW68JIj8YWkg/N83rcDzewbxtfs2Yv3+YVxy0mzfvrK30iJH8GsaoTeTdJ8DtoVz4mzb2hkYL/jGBXgurMHxAs5Y6Ak/TSMcP9ur8j5tQZcb/A/Oss9XBMOp8zvxL++wXXnBH3lPOuEuXDQ4UUB3JokZGeWc5PmiF/Wm0ZNOuNbJ3K6Ub22KpGHXLgBODMS03HiD2ppFMpotojuTwNu/b7evP25WO9IJHas/+xps6xvFigVduPqchTgwlHWPKa0Ceb1SqA1nC+61dbYZPgXSHbRAHJeJVD7SRx+0QEZzRXedcFnoVprG6z2XFpacOYetoZ4rWm5m2+pNfVi9qQ8ffdUJrpCTsZ+EThDCUSBFu938l367DqPZIjQi91pVF1ZY/APwx0AAf2KFPwvLHwNJJ2wLxIuBlFcgQUGeK5i+32XQAilaFo6M5TGr077/0pKwLAFNo5Igutp8MWXoNcVADJ2gafax1F5cQHOzsOKm8b4VwFsAjAGAEGIfgM5GDaocQoh7hRAnCSFOEEJ8pRVjaASGrrkCN4orz5iPX338FSU/MCLCW85egM6Ugc85Lqs4LJvV7mbryLRQSXc6gY6U4bMc3rFyEd5y1gLffurqaZc71pRE/aF2Bvz3C3vS7g/1tAWe4gn6+TtSXosO1aJIBlwcPZkE0k7q7dB4vsQCUSeEi2dk3HgBAMztbPPd+5Suu/d4PF9E3rRc99CGA3ZmmNpQbyTnz3KayJtoS2joTidwzpJepAwdC3vSvsyoYGW5fG14ouheW9DiUO814Fkg3v1JuGM2nKVV7V5YBcxx1lFPGhpSidI1w8PirgVTIFc0I9N4g/UGA+N5V+CpqwTK68kXLTyzawB3r92HiYLpy8JK6Jor+KMshOBnrqJaLG5w33VhGX4LpIICCbVAyriw+sfs65afh9z17C/dj939pRarHbewnH3JdfGFx0DCs7CkBVK07HiKX4FMvjqQvLDVqAAAImpskyCmas5d0ovnv/g6nL6wu/LODl1tCdzx0Qvx/pcvwTvOW+x7jYhw+4fOx/967YnuthuvPBXvOG9RyX6Sd5/vP4ZPgaT8wtDQNTfgvkJxuQUtEMAOhgN28F0SbITXnXYUSN62QHoyCZ+SkFx26hy8bsU89CjWydzuthILRCqU0WwRQsBtwbJ+v22BfONdZ+Gt59gN/oK9poazhdBZtJoZJYsIpVKQr0046eBAiAIJPJ/d4VkcgBdEn9luC7JMUsd4roiRbNHt3yTX7A6u5RK2WBEAjOXMEncXYLuwgoJqcLxQ4u9Xs6vyRf977CC6V0go94sS8KlyCiSsF5ZrgWgwLcuLgVRwYQUtgVILxL+/jFN5CkQ2YyzipX3DEBDozSTwr28/07WG3G7FGikxkBALJGJFQl0j26q2BEwh3EB90tB8sZJGEzcGcgcR/RBADxH9FYAPAbilccNiaqFcYDCKMxf14MxFPaGvyeB/JT7+6hMwliu6QXKJKgiCs2cAWDYzg+2Hx7BkRsbNDAsKSQC48ISZ2No35gusJgI/8lkdKWQShq1Axu0YyIyQGMh333cu2hI6Vil1GXO72pBJBVxYjpCRldGqBTK3K4VXnjgLpiXwm2f3YiRb9LkpxvNmWQWSLZiuBSLrItoM//nVeybdb0EX1gxHUcg6j/aUgW+/52w3HtDZlsCgU0AnEx+SyvriKlE++zFnzZHguubZolniKhkcz5cIL10jGM73MhdQIHYab2kWVpsRz4UF2D3XtvWN+YRvMmiBJOwYiJrGW46SIHrAAon6nQUtEAA4MDQBS9jHfNf5i/HV329A0bLc+6Rr5C6rGywkLBcD0YlcF5eqx1OG1lQLJJYCEUJ8nYguBzAMOw7yD0KIBxo6MmbKEOU2U3+onSGKYenMdgB9mN2ZwpzOFEZzxRJXFwD8/RtXYOXSGb5AqfxBn7W4B598zYm45KTZuP2xnRgcz2M4a8dAwoLocnYvLZCOlIGOlIFMwhPWukYhCsQWDnsGJnDe0l7nmjwrJbhgUZgCUVODZRBVKia1QV/SjYHYY1zUm8Gu/vES60wqGDVIftXZC93HnW0G9jtxF9eFpZOdml5igUQokHwRY/kiutKGf/XEglVSbzA4UShpg5LQyXX55E3Tl25qu7A84SmD7W0RFkiYC+uE2R3Y1udf2CzYysTNwnLSeIPNEIMkAoI8W8ECkczq8MdAAGD/cNaxDrxiQdPyFuQydPKC6CUxENtFdcVNq/Hely1B0RJI6ISC0zxSI3KsD++epgzdt5Jio4kbRP8XIcTnADwQso1hQikXAwGAD1ywFMfNakdbQsecrhS2HR4LdWG1JXRcfc5C3zapQDIJHZevsGMv6YSGjc4PtsfpAvD+ly/B3Wv32e27NXIzyKQ7TApWuaRq0I0iLYWZHZ4yWuhUgctrGskVS1J7w4SUGwMpWnhh7xAWdLe5s1ZVaCYCLqyLls/Cwp50SaKDbFwY1q5djm//kF1H4lkgtlBSLY5vr9qMnzyx0/fe7nQCQxMFjOWKGBgrYHZnys0CBEqtCcCOgQR7s+la0IXlbwgohb2axpuOEPBhLqwTZnfgARz0rfeR0AlEXj1Om2HXgUwUTKQT5VN45VhU8kV/r62oJozSQlTffmAoi7Sz1gngNUiU919X0pxLYyC2Yt/aN4odR8aQ0DUYmgaCcNN4pQtLvUeqkm40cWMgl4dsu7KeA2GmH3abcftxmAVy4pwOXPuKZQDg+ujDXFhhSKGkKqlM0nCrhnsyCega4StvPcN10akz2F7HApnrnFcqDPljls+DLizAWz5WXtNLe4dw19q9vvGFuWHalBbrz+4cwLlLe0P3D7qwejMJXP+aE33tZADg1PldWLm0F19482kl55LjOxC0QAwZA/GEzP3rDpR0R5b3ZzRnYmA87yogSa5olqyfMjheKHGFGUqQuCQGonkWSFLXXAURFaNIhSgW2WhRzchSK9yTutcddzxvVszAAkKC6AXTJ9yjFIj8fHTl/fuHsk7Vvf1cLf4DKsVANKcJpV1pXnQKJhM6eb2whH954VRCmzxpvET0cQCfAHA8ET2vvNQJ4M+NHBgz9SEiZBI6xvJmSXZVkDlOSmpYrCQMafarcZa2hO5mvKgCXwopVTDI2aKsd5FNAJOBTKDhgKsJgNvGRVoAP1ztNWVUxxJEc4TpziPj2DeUxUeUmh/b1WPHOoJZWGHKV56jXMFrh9MiHvB6ayV1DZblrxUIZmQBdsbbjiPjGHeaUwaTM3JFvwtLIzsGYgmBRb1pdx0YQ/eysHJFyxdkJ6US3YiRxhtMnNA1chMsBifyJftmCxYSOtmCWNhZWJUysICQNN6ihR7FLRWmP4i8yYyu7HBgyK6BkkrHs0C8LKyoZorqOApFAZ0s6DpBB0Encrv+qlZfytCRL1a3WNvRUMkC+RmANwO4y/kv/84TQnygwWNjpgFpRzBHCUHJSXM70ZEyQjOnwnBdWD4LxHssM5EATyD51l3QNZy7pMeNZ0jBIoWYTD8NW+BIdiLOJPVIf3iUn70tobkNLVULhJSUVi8GYrjnrgXVbTizPQUib70K1VIIrhYIeBldo7kiBsbz6M0kfIkLdoW+wNmLe3DFafMwt6sNA+MFmJZAe9LAO51svfaU4d7TvBl0YSmV6FocC8T/OXanE65lpFogwX2l0JYurEoELYFswfS5sIKvA0BH0nBdY6rgPzCUhSW8HkxhFojbyiSkDkRSsCy367Ch2dekO0XBQRdWuSLGelNWgQghhoQQO4QQ7xVC7AQwATuVt4OIams0xRxTSKEeFgNReft5i/DI514TOfsM4ikQ77iqcFBjFkHBLPn1J16JD1xgt4WRFogqJNMJXWnxrrlWlIyBEFGJWylsLCptCd2NJZwwu73kNUB1YUkLpDYF0tGmJjEYSCd0JA29xIUV1vpCJhkcGcvbM/BM0jdDHsvZSueyU+fgB39xHnoySQw6CkTXCP/6jjPx+I2Xoqst4XNhqecK9sKS1x0VRE+5lpl9Xd3phFtMGbSi1CaNupO2HNeFJRsiSnJFq6ILS/0eqGt65E0LR8ZyinLRYJpKGq+uuUsThKXxSgqmcO/trI4kejJJd+VDVV80OwsrVgyEiN5MRJsBbAfwMIAdAH7XwHEx0wT5gw2moAbRNfLVeVQiGKsIPg5zYZUrRHNjIGpDv4TuurCSTqtuwIuBAF5/oo9ecjx+8mGvMXSUIlQFZtCySCuvAd49C0ssiIMaT+pssyv17ZiAFrBAQlxYzmexZ2Dcfa4KUZk6LO9XTzqBQaeQUNfsLtbSPahptkAOBt41IldRybXe1fsQJOja607bDUkBYEG3P4Vc3tugBRJ3gqJaA8EVD8OsznYlDTzoAts3OFEaAzFLYyBhrUwkBScBwdA0/PhDL8NnXncSNM1ePdLnwkroKJiipPq9UcQNon8ZwAUANgkhjgNwKTgGwsQgHdMCqZaE00wxowgEKXg6UoZPUKiLFkURzMICbOUnXVhJXUNnWwKzOpK+Y8tMn/MDacZhAV/AU2Y9mWRJNlAqoOjsAPkKvPKE8gsgRdEZsEDeeOZ8XLR8ZkkhYS6kVUm70wFAxjJ6MwmfUJVFelL49bYnnDRe/wxckjS0kiA6EeG0Bd34n+tfifOW9ioxkPJZWJ1pT4EQEX75sQvxm+v9XZXkeilJXYPuuOzySgfhSgRTedVrD8vi6lCUtVTOslvBSLboKhA7oG8pWVhKHUhIjMc7pgXTslsdzelqQ0fK8ILoAReWOoZGE1eBFIQQRwBoRKQJIR4CcHYDx8VME6QFUqsfP4pEGQskWEAYFgMpGadTB6K6udoSurs0a9KwXQdydcIgS2dmkDS8lMyoma5Ucr2ZkHRlQ7raZGBZwwdfeVxsoRdE3nNZpPelq07HW89Z5LNAhBChMZCUoSGT1N128j2ZpCvg0gndLdKTn0N3OulYIFbJTFoeLyyNF7C7MxN5LqzIGIhzf6QLS1ov5y+bgbmBIlbpZlQtkKJllWRYRVEuHqGOXdKhWCAypViOTy6eJY8TlYUV1spEknfcXkFXWtESvnbxrgJpUjV63F/1IBF1AFgN4KdEdAhAePI5wyikEwbanfWt60lYDEQqKzX+AXiKpZzwyIRYIB0pw800Suo6vvzWM0rqHCSyLUt7Skd+3IqupnYVSKm7zs0WM6ItpWroVLK41FmzoXnrTBQDPnRJUrdddnJBq972hNuRtqPNW09GKszejN05uGgKXxaSezxHgehl6imkgohupuiPDZVzi8pMuVzRcgvyZBpsHILf1+B3RyPyZZS1K99D2ShTVSDyrGFZWLFiIIHVMeXrwYJQeQ/zpoU04rnrjoa4v+qrYAfQ/wbA7wFshZ2NxTBl6UjpNfvwyxFWByIFz8yABaLWBESRCVEyql87YRAW9qR9C2ap45DnlgotKl20rawCcTLB6qRspQUSdB+qQfSw+AdgWy3tKcPXZ0sK1U5lu1R27Sm7YWG2YEILGX7S0OwsrEAQXUW68Cr1wupMGejNJLAk8FmoSAukbyTnKo1gS5JylLMG7LH7n6sWtrw30hIuKhaItBo8C0SLFQMpOm4vnwIhKrE05D0MKpZGEbeVyRgAEFEXgLsbOiJmWvHRS07wtdeoF6FB9ETtLqw2QweRv9pZzayJUj5//NtXu3ESwBMkkWm87tod4RX3QHlLqRqi0oBVF1bYYlGAXQ+jXn+PEgNRs7vUduyALaQ7QuJdSd22QKwysYS4vbBShoZVn7mkbGq4THQYzRVdoZstmGXjYCpBF1bwfUEjS/0efuySEzAwlsfV5yzEvS8csJWFEgPJFpQYiK64sEpiIN7zvJOF5auI16gk40rew2a1dI/byuSjAL4E2wqxYN8OAeD4xg2NmQ6cOr8rcn2So0Fm4qgzec+F5a+ajkrjVdGc/ldBF5YkSqgv6Em7AVvAs1qihGA5CyQdQ9FVg1eIGGaB2IKnvAUirSrdXrdC89en2PvJuhmvWDC4sh5gK6Rc0YRQnB7BWf3MjhQ0Qkk8wzuG/d5UQi/5jIMsVDLl5Ex+omDGdqUGg+jBawpaIJlABuDX3nmWu9Jn0bKUGIiGomX6YyCRS9p6jwtFaYEo94+opB+ZVLLNSuWNGwP5WwCnCSEON3IwDBOX0xd24Y6PXojzlyntQCJcWHHSeAHb/ZTUwy2Qcq3EVdpdC6RSEL1UgQSzsI4WaQkEZ+qGXtmFlTQ09KSTvvHI2a9apJkMWiAFE2GT/KSzoJVXUlfqwlrYk8afPvfakpRcSRxXpGSWMsa6WCAVguhhgX+pNIIxEMvXC6tcDMS7TpmFlQjEQKIskGZlYcVVIFsBjDdyIAxTDUTkS5sF7OA5kb9OA1DrK8oLnqvPXuBr2eFzYcVVIE4MJErheGm8jXdhuTGQoAtL8yrRgxlYRIAQtiY3oGEAABsqSURBVJC+4cpTsHRmBstm2plncvY+S5n9B9tw2BZIdBaWOnMPK8hb2JMu2eYeowoFK1OJE7pXFJgtWCXtQqII7lcpBhJW/CgPUQzLwnIEf+w6ENNOAvB1BdaiYyCTzQK5EcCjRPQEALfRihDiUw0ZFcPUwPzuNH7/6Vdh+Rz/6opxYiAA8Pdv8q9Zr6ZmxhXqsnVHWC2EOpbQILpRXxeWXNckmMRgx0AcF1agBkQnQlEIJJwFvz6rtOqXs1+5dCvgWR6GEgOJUiBqPAKofv2adELHpy9djitPnxdr/4f+9tXIJHXcr6z9EvdzDFoqQddXcOiZMhaI2kzRzsIK1IFExkCUNF4nC0utL9IJU8aF9UMAfwDwAuwYCMNMSk6eV7rSclgzxTjUYoEs7GmL9OEDShpvyFolamfaenHTe87GKYF7kvBZIP6fs+Z0dAy7Xqkk1M680n8vBW4+QoHIILq/jqG6ayEi/M3lJ8XeX9bsqOeM68IKXkOlrKywzDHVSqHIOpDoLCw1HVpmYWXUivipEkQHUBRCfKahI2GYBqFm71SDGkSPWz/wideciGucFvVhlCskrHcQHQBef1rpbF3XNAhh++aDfbCk0ApTttKtI/tPAZ6yk/76vGmFuqZkHYg6y45qi15vVGEfjG1EEbQGgu8rcWGFWCDqIUrrQNQ14yNiIMo5ZS8s9XuoE5UoCq+QcHJZIA8R0XWwU3hVF1Z/Q0bFMHVE7T9VDV6Ldy22u6UtoZfttzS3K4W2hOa2Vw8bZz0tkDCkMCxaVkkMRAqxuBaI3E8N7oYpWxlET+jRdSCNwvCNrUYXVmCwwe9DWBCdfBaI/V/X7DXL41SiB3thldSBOMpIRVq4QddWo4irQN7n/L9R2cZpvMyUIKxRYhykC6ueAv0tZy3AK06YFVrD4FWiN1iBOEKoaIpSF5Yjn8KuWQpftdI/zH8f2gvLcWEVDH8vrGYQtl56JUrW5iipRPfvH7bOuh6SMOBaIKYFIvtexYmBFJwsLJ8LMCLWBEyyViZOA0WGmZK0uX766gSzrHeop0vJ0DW3S22QVNMsEE/IVGOBJJz1u9UEgLD1vCNbmZil3XibgaoM4roiK1kgcuxzOlM4NJLDCXNKe6RpYRaI09hR7WsVuSKh8v6C0wsrWIkexIuBTAILhIheK4T4AxG9Lex1IcSvGzMshqkfmaQBjfz9iuLQrjQjbAZetlhjz+daIJZVmoWl+YPivvfpdsBXddGFrecdFmewV8oLb6bYaPwxkHjKuTSIHm6B/MUFS3HdJce7cTbfPmoMxGeB+PtaeYtqBTO/vOf20rZ+CyQ8280ex9o9Q5jX3YbTFnSX7FNPKv2iLoGdfRXW90oAYAXCTHrSSR23f+hlOGNhdT+msBbvjUTGFuKuylgrXgyk1IUl5WuYtZZQ3C3utpDlWMsF0fMtsEBUQRu/kLB8EF0qBE2jUOUBBCwQZZtngfgt43KtTAC7jkWvcJ9lmu93HtyMhzf14c5Am/t6U1aBCCG+4Dz8khBiu/oaEbFbi5kyXLx8dtXv6WhADKQcFxw/A6s+cwmOn91ReeejQLbpsBVIwIVF/lmxiuG4sHzHMkotltA0XseFpWZ9NUl/BNJ447YyqeDCcg5T7hr8RZPecUxL+CyQdFKHrpGv7ijsnNm8GbBASs+pfm7B4zWCuDb9rwCcG9j23wDOq+9wGGbykE7Ya57XqzK8EkSEE+c0VnkAnoAvmp4L64m/uxQJXcNbvvsIgPBrftfKxThrcY9vmxcDqTAzdgTbWK5Ydr9G4HNhVdnOPZ3QMVEwIyvRy11DmAvLi4F47qjudAL//bELS3rGBY+dLZr+dvhlXFhA/dfgCaNSDOQUAKcB6A7EQboARFdLMcw0gIjQnjSqrh+Z7Eh3TMHJwtI1cosfywXRVy6bgZXL/O1jpDVjVBDS0opTezTFzKg9amqyQHRyajQ0TBTM0PVA7P/Rx4hrgQDAOUt6g28vcZsVTOGzhssF0QF/IWyjqHQ3TwbwJgA9sOMg8u9cAH9V60mJ6J1E9BIRWUS0MvDajUS0hYg2EtHrle1XONu2ENENtZ6bYaqhPWU0zQJpFtL3bjouLFXolCskDEPOgo0YLqyS907iQkK7xYjXRyt4TXLo5a5B98VA5HE0JyBeeXGrsPuoWhVRLWPC9m0UlWIgdwK4k4guFEI8VsfzvgjgbbBbpLgQ0QoA74Ft9SwAsIqIZO+C7wG4HMAeAE8R0V1CiHV1HBPDlNCe0psWRG8WUoh+7b4NWLX+kK8qXm1CWI7FM9LY3T/hHTPQYiNImBXXtDRevXoLZH53GvO7067LLXg/pHIoV8sSlsYrlUbB9K/OGDrukPuYUeIa4a5C7/XJYIFI3kpEXUSUIKIHiegwEX2g1pMKIdYLITaGvHQVgF8IIXJO0H4LgJc5f1uEENuEEHkAv3D2ZZiGsnhGBvPK9LaaikhhuGr9IQB+oaMTIalXrry/8/qLcPcnLyo5JhBdiR6keRaId+64WVjXvep43PXJVyoWSLgLq9zR/DEQORb7QS5GZ+Cw+6M2baxk6bXcAlF4nRDis0T0VtgWwDsBPATgJ3Uez0IAjyvP9zjbAGB3YPvL63xuhinhe+87t2mCrlkEhaEqdDStNNMqjBntSd/Kj7424yH3K6y9y+RuZWLXvEgrIWpFwrgxEOnC8pbXLQ3Ml4w75HPIlHFh2Qkf3rb2iKWB60lcBSJt3DcA+LkQor/SDIWIVgEI67v8ecc1Fvq2kG0C4ZZSaK2+07PrOgBYsmRJ2TEyTCWa4QZoNsEUVbmGN2CnhtbislNdQ2GCMbjIF9CaViZxYyDu/lp4lbgbRC+jBPxrd/i35QKdiSu9X6IWwwYVtUbke08zvrtxz3A3EW2AvaTtJ4hoNoBsuTcIIS6rYTx7ACxWni8CsM95HLU9eN6bAdwMACtXrmxOQxiGmUIEBdPh0Zz3GpUWC8ahUoV0WPPIVlgg1db0yGspqUR360CiL0J9qdQCCW97rxJmLakxkOClaBr5xtMMF1asuymEuAHAhQBWCiEKsFcnbEQM4i4A7yGilFOouBzAkwCeArCciI4joiTsQPtdDTg/w0x7yrXz0JzU1aM5ZphgnBMSR2qNBVLdtRmaX/BLqk3j9Xph2efPFc3KFojzJnW99fIWiP/9LQ+iE9FnlaeXCSFMABBCjAGoeTVCInorEe2BrZTuIaL7nOO+BOAOAOsA/B7A9UIIUwhRBPBJAPcBWA/gDmdfhmGqpFyMQwbRqyW4TkWQMH+8JZrUMbaGZooS1wKJaGVCZcLovjReClgghcoWiIy/qIpAVSblmi8CQEdb6y2Q9yiPbwy8dkWtJxVC/EYIsUgIkRJCzBVCvF557StCiBOEECcLIX6nbL9XCHGS89pXaj03wxzrlBNctVsg5V1YYWtjBNeyaBR6DWm8kmDHXIm8xHI6QJXncj9/DKT8WOS5VeXbXjaIHlAgkyAGQhGPw54zDDMFUIXhBy5Ygr9UVlCsNQaS0Mq7sOzz2ivopRM6xvOmryq9kdSypK1EXkstrUyICBoBlvCEpbQS4mRhydfVtUZUZVLiwgocr+UuLPgznYKfNgeoGWYKogqu85fNwIlzvDXTZ3WmMCck4F0JTSNvVh4hGHucLsNy9T6rWRZIDa1MJEZEe3t5yEphnKCikZZarmhVzAiTyka1JCql8ap0VLl8QS1UOsNZRDQMW4Gmncdwnk+v6iqGOUZQrYXgUqxffdsZNccmDGfVwbAYCAD0pBPoG8m5NSGtsECqTeP1LBC/4qEYFgjgKFNLhBYSVrJA7Kwqf+aV+nkF73PweO2t7sYrhGj8CBiGaSpqTCAdCG4fjdvD0Ah5RAeq5TonXgykOavm+bvxVmmByG7DEbP9yhaI3C+YxmvGGouhkS/zyl9bEh7Y98be+BY806vJD8MwFVELCTN1rFaWwjHKhdXt9NwadfpLNc8C8cRctRlmbhpvsBeWFs8CcXtmue+zz2+JeBlhn7viFLxj5aLwYwcuJcryaySsQBjmGEOdmYa1GKkVb23v8Nc/+/qTcdqCLrzKWdyrFTGQ6l1YpSsuAkovrJgxkGAzRQAVmykCwEcuPh5nLeoJfS2ovCq5xBoBKxCGOcZQBU0wBnI0SOEcjBdIls/txD2futjt/jsVYiBRhYRxYyDBtu+1LG4VlTkW1WK+mbACYZhjjESZGMjR4PaNqiDJ5KJUJ8/rLLtfvZDBaMCfQBAHaSUEFY8WUAyR79fKWCCxFUj4mKOC6D2ZRNXpyrUy/TrFMQxTlsZbIOX3e/NZC7ByWS/md6frdu5KGBrBEuWbH0a9z/5f/YqE6n50VBZI+A0NXos81xN/dymaVOTPCoRhjjXUWXg9YyBRa2eE0UzlAdiCW6tBqEa1MomdhSUtEOe5v2I/njUUZakELRC5m7q+S6NhBcIwxxjqzLWe671XCqK3kmrTd733lY+BVGoIGXR16UfRlytIpVYmzYAVCMMcw9SzI27c1NZWoFH17itAycKK7IUVM433KGIgUQSvpxVZWKxAGIapC1LI1jrbbySGrtUkYI+mnTtQmq2lKpxqLJCb3n02Vizo8m0LurCa1R5fhRUIwzB1QRYoTkL9AV2jkpUY474PKKdAYmZhOc/bEkrTySoypa4+Z2HJtuB9boXrcBJ+1AzDTEXcLKxJ6MIyNKqptYeXGBBRc1FlK5OZ7V6jyqOOgSj3OaFTS1yHrEAYhqkLUeuHTwZ0jaouIgRsK8EILBUrjwfEaKYYiIF0pY2qstXKji3QZZiD6AzDTFm8OpDJp0AMrbaVFi86cRaGJwol22PXgQRcWESEnkwCh0fzR22BqEH0pKE1bY15FVYgDMPUhelogVy8fDYudnp3qQRblEQRlq3VnbYVyNHeJ9WFlawxSeBoYRcWwzB1QbbPmIxpvIam1TU7rNZmioC3sFY960ASusZZWAzDNIeb3n12TSsPlsNN421SH6Zq0DWqa38oNzheIYoebGUC2BaIHNPRjcFfENqK5AVWIAxzDBKWFnq0uMHhyWiB6LWt9R5F3BhIsJkiYK/MCNTPAiGyr68V6dOsQBiGqQuVFpRqJUtntvvWFj9a3ALBSsvSkv8/AHSl69POXupDnQi6xllYDMNMYbxK9MmnQP7tvefU9XhhiiF0PzcLy9uxx1kPJSy7q7oxeDEnQ2tNHQgrEIZh6sJkDqLXm7DYRrn9tBAX1uBRKhC3FkUD/vIVy45qPftaYQXCMExdkAJtMqbx1hsviB5zPzWI7lggg+P1sUB0Irz9vPB10xsNp/EyDFMXEpPYhVVv4i5pG5bGe76zIuObzpx/VGPQJ0HMiS0QhmHqwmQOotebuM0UXQWi2CqLejPY8dU3Nm0MjYQtEIZh6oIMok/GNN56IzOgKl1qWBpvvZBpu610GbZEgRDR14hoAxE9T0S/IaIe5bUbiWgLEW0kotcr269wtm0hohtaMW6GYaJhC6QUr+VJ/cegH8MWyAMAThdCnAlgE4AbAYCIVgB4D4DTAFwB4N+JSCciHcD3AFwJYAWA9zr7MgwzSTgmYyAVJKhngdT/nnhJC3U/dGxaEgMRQtyvPH0cwDucx1cB+IUQIgdgOxFtAfAy57UtQohtAEBEv3D2XdekITMMU4E3nDEPGqEl6aTNpvpWJg0YwyRYQngyxEA+BOB3zuOFAHYrr+1xtkVtL4GIriOiNUS0pq+vrwHDZRgmjKUz2/HRS05o9TCaQux27jEVTS1MBhdWw6YKRLQKwLyQlz4vhLjT2efzAIoAfirfFrK/QLiiC+0DIIS4GcDNALBy5cqj6xXAMAwTQlh9R/h+8RRNLUyGupuGKRAhxGXlXieiawG8CcClQggp6PcAWKzstgjAPudx1HaGYZimQrEtkAa6sKj1CqRVWVhXAPgcgLcIIcaVl+4C8B4iShHRcQCWA3gSwFMAlhPRcUSUhB1ov6vZ42YYhgHiZ2HFXfq2FhqZIhyXVkW7vgsgBeABR5M/LoT4mBDiJSK6A3ZwvAjgeiGECQBE9EkA9wHQAdwqhHipNUNnGOZYx3Nhld+P3BhI48bQyrqbVmVhnVjmta8A+ErI9nsB3NvIcTEMw8QhbgZUI9N4iQgaHYMuLIZhmKkMxbRAGhkDAWzlcayn8TIMw0wp4qbQhvXCqicatWYlQvf8rTs1wzDM1CR+M0X//3qja9TSGAgrEIZhmCqZNC4sopb2HmMFwjAMUyVxFYPWwCC6PD5bIAzDMFMIzzUVz4XFQXSGYRgGAHDS3E6cNLcDnW3lKyEaWUgoj9vKIPr0b5vJMAxTZ15x4izc/zeXVNzPy8JqDLrGdSAMwzDTkkYvO6sTu7AYhmGmJY2OgWgcA2EYhpmeNLKViTw+u7AYhmGmIdToGAi7sBiGYaYnzXBhHXNrojMMwxwLNDqN923nLsS8rraGHDsOrEAYhmEaBDW4lcknXh25MkZTYBcWwzBMg9CpsUH0VsMKhGEYpkG4MZDWDqNhsAJhGIZpEHFXLpyqsAJhGIZpEI1u595qWIEwDMM0CJli28Jav4bCCoRhGKZBNHpJ21bDCoRhGKZBNDqNt9WwAmEYhmkQuluJPj01CCsQhmGYBuEuadvicTQKViAMwzANwl0PZJpK2ml6WQzDMK2Hg+gMwzBMTcg03mkaAmmNAiGifyKi54noOSK6n4gWONuJiL5DRFuc189V3nMtEW12/q5txbgZhmGqgbgXVkP4mhDiTCHE2QB+C+AfnO1XAlju/F0H4PsAQEQzAHwBwMsBvAzAF4iot+mjZhiGqQLdXRO9xQNpEC1RIEKIYeVpOwDhPL4KwI+FzeMAeohoPoDXA3hACNEvhBgA8ACAK5o6aIZhmCqRwfPpGgNp2XogRPQVANcAGALwGmfzQgC7ld32ONuitjMMw0xaNLZAaoOIVhHRiyF/VwGAEOLzQojFAH4K4JPybSGHEmW2h533OiJaQ0Rr+vr66nEpDMMwNTHdmyk2zAIRQlwWc9efAbgHdoxjD4DFymuLAOxztr86sP2PEee9GcDNALBy5cpQJcMwDNMMNA6i1x8iWq48fQuADc7juwBc42RjXQBgSAixH8B9AF5HRL1O8Px1zjaGYZhJi5vG29phNIxWxUC+SkQnA7AA7ATwMWf7vQDeAGALgHEAHwQAIUQ/Ef0TgKec/b4khOhv7pAZhmGqY7qn8bZEgQgh3h6xXQC4PuK1WwHc2shxMQzD1BNO42UYhmFqQuNKdIZhGKYWprsLixUIwzBMg9CJ27kzDMMwNeAVEk5PFcIKhGEYpkFwDIRhGIapCZ0tEIZhGKYWzlrcg4++6nics6Sn1UNpCC1rpsgwDDPdaUvouPENp7Z6GA2DLRCGYRimJliBMAzDMDXBCoRhGIapCVYgDMMwTE2wAmEYhmFqghUIwzAMUxOsQBiGYZiaYAXCMAzD1ATZazhNT4ioD/aKh7UwC8DhOg6nlUyXa5ku1wHwtUxW+FpslgohZlfaaVorkKOBiNYIIVa2ehz1YLpcy3S5DoCvZbLC11Id7MJiGIZhaoIVCMMwDFMTrECiubnVA6gj0+Vapst1AHwtkxW+lirgGAjDMAxTE2yBMAzDMDXBCiQAEV1BRBuJaAsR3dDq8cSBiHYQ0QtE9BwRrXG2zSCiB4hos/O/19lORPQd5/qeJ6JzWzz2W4noEBG9qGyreuxEdK2z/2YiunYSXcsXiWiv89k8R0RvUF670bmWjUT0emV7S7+DRLSYiB4iovVE9BIRfdrZPuU+lzLXMhU/l7b/v737CY2riuI4/j2kSRRbTOrf0Aom0k0XUoOUgtKFSrTdRKGLrBQVBP8sXLgoFKQuFXUnFkShithqVexGtPgHV6aiNjVSrNNWsDQ0Qm21G//U4+KeVzPDzCRzSX3vwe8Dj7nvvrc4J/dOTua+m8TMDprZTOTyTPSPmtl0fI33mtlA9A/GeSOu37hYjj1zdx1xAH3AMWAMGABmgPVlx7WEuH8Crm7pew7YHu3twLPR3gp8CBiwCZguOfbNwDgwmxs7sBo4Hq/D0R6uSC47gafa3Ls+5tcgMBrzrq8KcxAYAcajvQo4GvHWbly65FLHcTFgZbT7gen4er8NTEX/LuDRaD8G7Ir2FLC3W445MekTSLONQMPdj7v7n8AeYLLkmHJNArujvRu4d0H/6558CQyZ2UgZAQK4+xfAmZbuXmO/Gzjg7mfc/VfgAHDPpY++WYdcOpkE9rj7H+5+AmiQ5l/pc9Dd59z9m2j/DhwB1lDDcemSSydVHhd39/Nx2h+HA3cA+6K/dVyK8doH3GlmRucce6YC0mwN8POC85N0n2xV4cDHZva1mT0Sfde5+xykNxFwbfTXIcdeY696Tk/E0s5rxbIPNckllj1uIf20W+txackFajguZtZnZoeAeVJBPgacdfe/28R1Mea4fg64imXMRQWkmbXpq8M2tdvcfRzYAjxuZpu73FvXHKFz7FXO6WXgJmADMAe8EP2Vz8XMVgLvAk+6+2/dbm3TV/Vcajku7n7B3TcAa0mfGtr9w/UirkueiwpIs5PADQvO1wKnSoplydz9VLzOA++TJtbpYmkqXufj9jrk2Gvslc3J3U/Hm/4f4BX+WyqodC5m1k/6hvumu78X3bUcl3a51HVcCu5+Fvic9AxkyMxWtInrYsxx/UrSEuuy5aIC0uwrYF3sahggPXjaX3JMXZnZFWa2qmgDE8AsKe5i18sDwAfR3g/cHztnNgHnimWJCuk19o+ACTMbjqWIiegrXcvzpftIYwMpl6nYKTMKrAMOUoE5GOvkrwJH3P3FBZdqNy6dcqnpuFxjZkPRvhy4i/RM5zNgW9zWOi7FeG0DPvX0FL1Tjr37P3cR1OEg7Sg5Slpb3FF2PEuId4y0o2IG+L6ImbTW+QnwY7yujn4DXor8vgNuLTn+t0hLCH+RfjJ6OCd24CHSw8AG8GCFcnkjYj0cb9yRBffviFx+ALZUZQ4Ct5OWNA4Dh+LYWsdx6ZJLHcflZuDbiHkWeDr6x0gFoAG8AwxG/2Vx3ojrY4vl2Ouh30QXEZEsWsISEZEsKiAiIpJFBURERLKogIiISBYVEBERybJi8VtEZLmYWbEVFuB64ALwS5xv9PR3lkRqQdt4RUpiZjuB8+7+fNmxiOTQEpaIiGRRARERkSwqICIikkUFREREsqiAiIhIFhUQERHJom28IiKSRZ9AREQkiwqIiIhkUQEREZEsKiAiIpJFBURERLKogIiISBYVEBERyaICIiIiWf4FCp/U0C0Fkv8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mean = numpy.zeros(300)\n", "T = numpy.zeros(300)\n", "Tval = 10\n", "i = 0\n", "while Tval < 3010: \n", " E = numpy.random.normal(0, 4, Tval)\n", " X[0] = 0\n", " for j in list(range(Tval)):\n", " X[j] = X[j - 1] + E[j]\n", " mean[i] = statistics.mean(X)\n", " T[i] = Tval\n", " Tval = Tval + 10\n", " i = i + 1\n", " \n", "import matplotlib.pyplot as plt\n", "plt.plot(T, mean)\n", "plt.ylabel('Estimate')\n", "plt.xlabel('T')\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "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.7.1" } }, "nbformat": 4, "nbformat_minor": 2 }