{ "cells": [ { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ] }, "execution_count": 19, "metadata": { }, "output_type": "execute_result" }, { "name": "stderr", "output_type": "stream", "text": [ "/ext/sage/sage-8.2/local/lib/python2.7/site-packages/sage/repl/ipython_kernel/__main__.py:11: RankWarning: Polyfit may be poorly conditioned\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "D-9 Esitmation for 2010: 1.0185356829006386e+44\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAw4AAAJICAYAAADINLOIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xl0lPWh//H3QAQRmCAqslnAClQQFAQFF9SKC1QixVonVRGXa70ubbW099bqr/bUFtfb4lJsBUGvmqhQEDWgYNgEooCiKC6AWlxYpEoCyCJkfn/MlYpsScjMd5b365w5Kck8z/MZtefkw3eLxOPxOJIkSZK0B3VCB5AkSZKU/iwOkiRJkvbK4iBJkiRprywOkiRJkvbK4iBJkiRprywOkiRJkvbK4iBJkiRprywOkiRJkvYqa4tDPB6noqICz7eTJEmS9l3WFod169aRn5/PunXrQkeRJEmSMl7WFgdJkiRJtcfiIEmSJGmvLA6SJEmS9sriIEmSJGmvsr44xGIxCgoKKCoqCh1FkiRJyliReJbuV1pRUUF+fj7l5eVEo9HQcSRJkqSMlvUjDpIkSZL2ncVBkiRJ0l5ZHCRJkqQM9NlnqX2exUGSJEnKMF9+CS1awJgxqXumxUGSJEnKMAsWwLZtcMwxqXumxUGSJEnKMC+/DA0awFFHpe6ZFgdJkiQpw5SVQc+ekJeXumdaHCRJkqQMEo/D3LnQu3dqn2txkCRJkjLIxx/Dp59Cr16pfa7FQZIkScogc+cmvlocJEmSJO1WWRm0bQvNm6f2uRYHSZIkKYOUlaV+tAFyoDjEYjEKCgooKioKHUWSJEnaJ5s3J85wSPXCaIAUbuAURnFxMdFoNHQMSZIkaZ8tXAhbtjjiIEmSJGkPysqgfv3Unhj9NYuDJEmSlCHmzoVjj4V69VL/bIuDJEmSlCFCLYwGi4MkSZKUEVasgH/+M8zCaLA4SJIkSRmhrCzx1REHSZIkSbtVVgatWkHr1mGeb3GQJEmSMsDcueGmKYHFQZIkSUp7X30F8+eHm6YEFgdJkiQp7S1aBBs3OuIgSZIkaQ/mzoX99oNu3cJlsDhIkiRJaa6sLHFadIMG4TJYHCRJkqQ0V1YWdpoS5EBxiMViFBQUUFRUFDqKJEmSVG2ffQZLl4ZdGA2QF/bxyVdcXEw0Gg0dQ5IkSaqRl19OfA1dHLJ+xEGSJEnKZHPnwqGHQtu2YXNYHCRJkqQ0VlaWGG2IRMLmsDhIkiRJaWrbNnjllfALo8HiIEmSJKWtt96C9evDr28Ai4MkSZKUtubOhbp1oUeP0EksDpIkSVLamjMncfBbw4ahk1SzOIwYMYKuXbsSjUaJRqP07t2bSZMm7fb9Y8aMIRKJ7PTatGlTlZ+5u3tEIhFWr15dnfiSJElSRpk9G044IXSKhGqd49C6dWtuu+02jjjiCAAefvhhzj33XF577TU6d+68y2ui0SjvvvvuDt/bf//9q/zMCy64gLPPPnuH7w0ZMoRNmzbRrFmz6sSXJEmSMsaqVbBsGZx4YugkCdUacRgwYAD9+/enQ4cOdOjQgT/+8Y80atSIsrKy3V4TiURo3rz5Dq9visfj3HHHHRx++OE0aNCAo48+mrFjx27/eYMGDXa4tm7dupSWlnL55ZdX86NKkiRJmWPOnMTXjCwO37Rt2zaKi4vZsGEDvfewP9T69etp06YNrVu35pxzzuG1117b4ec33XQTo0ePZsSIEbz11ltcf/31XHTRRcyYMWOX93vkkUc44IAD+NGPflTT6JIkSVLamzMHDjsMWrcOnSQhEo/H49W5YNGiRfTu3ZtNmzbRqFEjHn/8cfr377/L95aVlbF06VK6dOlCRUUFw4cPp6SkhNdff5327duzYcMGDj74YEpLS3coH1dccQVffvkljz/++E737Ny5M6eccgp//etf95izoqKC/Px8ysvLiUaj1fmIkiRJUnAnnABt2kBRUegkCdUuDlu2bGH58uWsXbuWcePGMXLkSGbMmEGnTp32em1lZSXdu3enT58+3HPPPcybN4/jjjuOht9aJr5lyxa6devGyy+/vMP3586dywknnMD8+fM59thj9/gsi4MkSZIy1aZNkJ8Pd90F110XOk1CtRZHA9SrV2/74ugePXowb948hg8fzt/+9re9XlunTh169uzJkiVLgESRAHjuuedo1arVDu+tX7/+TtePHDmSY445Zq+l4ZtisRh5eTt+zMLCQgoLC6t8D0mSJCmVFiyALVvSZ30D1KA4fFs8Hmfz5s1Vfu/ChQvp0qULAJ06daJ+/fosX76cU045ZY/Xrl+/nieffJJhw4ZVK19xcbEjDpIkScoos2cnzm7o2jV0kn+rVnG48cYb6devH4cddhjr1q2juLiY6dOnM3nyZAAGDx5Mq1attv9y//vf/55evXrRvn17KioquOeee1i4cCH3338/AI0bN2bo0KFcf/31VFZWctJJJ1FRUcGcOXNo1KgRl1xyyfZnP/HEE2zdupULL7ywtj67JEmSlJbmzIHjj4e8ff5r/tpTrSirVq3i4osvZsWKFeTn59O1a1cmT57MGWecAcDy5cupU+ffGzWtXbuWK6+8kpUrV5Kfn0+3bt2YOXMmxx133Pb3/OEPf6BZs2YMGzaM999/nyZNmtC9e3duvPHGHZ49atQoBg0axIEHHrgvn1eSJElKa/F4ojhcdVXoJDuq9uLoTOHiaEmSJGWiJUugQweYNAm+dQ5yUDU+x0GSJElS7Zs9GyIR6NUrdJIdWRwkSZKkNDJnDnTuDE2ahE6yI4uDJEmSlEZmz06vbVi/ZnGQJEmS0sTnn8PixYlTo9ONxUGSJElKE2Vlia+OOEiSJEnardmzoVkzOPzw0El2ZnGQJEmS0sScOYnRhkgkdJKdWRwkSZKkNPDVV/Dyy+k5TQksDpIkSVJaeP112LgxPRdGQw4Uh1gsRkFBAUVFRaGjSJIkSbs1ezbUrw/du4dOsmuReDweDx0iGSoqKsjPz6e8vJxoNBo6jiRJkrRHF1wAn3wCL70UOsmuZf2IgyRJkpTu4vH0PfjtaxYHSZIkKbDlyxOjDem6vgEsDpIkSVJwX09PcsRBkiRJ0m7NmgVHHgkHHxw6ye5ZHCRJkqTAZs2Ck08OnWLPLA6SJElSQP/6FyxebHGQJEmStAdfr2+wOEiSJEnarVmz4LDDoE2b0En2zOIgSZIkBTRrFpx0UugUe2dxkCRJkgLZsAFefTX9pymBxUGSJEkKpqwMtm61OEiSJEnag1mz4MADoVOn0En2zuIgSZIkBfL1+oY6GfBbeQZE3DexWIyCggKKiopCR5EkSZK2++qrxFSlTJimBJAXOkCyFRcXE41GQ8eQJEmSdvDqq/Dll5lTHLJ+xEGSJElKR7NmQYMG0L176CRVY3GQJEmSApg1C3r1gnr1QiepGouDJEmSlGKVlfDSS5kzTQksDpIkSVLKvfMOfP65xUGSJEnSHsyaBXXrJqYqZQqLgyRJkpRis2YlFkU3ahQ6SdVZHCRJkqQUmzUrs6YpgcVBkiRJSqnlyxMvi4MkSZKk3Zo1K/H1xBPD5qgui4MkSZKUQrNmwfe+B4ccEjpJ9VgcJEmSpBTKxPUNYHGQJEmSUmbNGli82OKQlmKxGAUFBRQVFYWOIkmSpBw3c2bi6ymnhM1RE3mhAyRbcXEx0Wg0dAxJkiSJGTOgXTv4zndCJ6m+rB9xkCRJktLFjBmZOdoAFgdJkiQpJT7/HN54w+IgSZIkaQ9mzYJ4HE49NXSSmrE4SJIkSSkwY0ZibUPbtqGT1IzFQZIkSUqB6dMzd5oSWBwkSZKkpFu7FhYuzNxpSmBxkCRJkpLupZcS6xsccZAkSZK0W9OnQ+vWcPjhoZPUnMVBkiRJSrKvz2+IREInqTmLgyRJkpRE5eXw6quZPU0JLA6SJElSUs2eDZWVmb0wGiwOkiRJUlLNmAEtWsARR4ROsm+yvjjEYjEKCgooKioKHUWSJEk56OvzGzJ5fQNAXugAyVZcXEw0Gg0dQ5IkSTlo3TpYsAAuuyx0kn2X9SMOkiRJUihz5sC2bZm/MBosDpIkSVLSTJ8Ohx4KHTuGTrLvLA6SJElSkmTD+Q1fszhIkiRJSbBhA8yblx3TlMDiIEmSJCXFnDmwdWvmn9/wNYuDJEmSlAQzZsAhh8CRR4ZOUjuSXhxGjBhB165diUajRKNRevfuzaRJk6p8/ZAhQ4hEIju9OnfunMTUkiRJ0r6ZPh369MmO9Q2QguLQunVrbrvtNubPn8/8+fP5/ve/z7nnnstbb71VpeuHDx/OihUrtr8++ugjmjZtyvnnn5/k5JIkSVLNbNgAr7ySPdOUIAXFYcCAAfTv358OHTrQoUMH/vjHP9KoUSPKysoAKC8v58orr6RZs2ZEo1G+//3v8/rrr2+/Pj8/n+bNm29/zZ8/ny+++IJLL7002dElSZKkGpk9G776Ck4/PXSS2pPSNQ7btm2juLiYDRs20Lt3b+LxOD/4wQ9YuXIlJSUlLFiwgO7du3P66afz+eef7/Ieo0aNom/fvrRp0yaV0SVJkqQqKy2F5s3he98LnaT25KXiIYsWLaJ3795s2rSJRo0aMX78eDp16kRpaSmLFi1i9erV1K9fH4C77rqLCRMmMHbsWK688sod7rNixQomTZrE448/norYkiRJUo2UlsL3v5896xsgRcWhY8eOLFy4kLVr1zJu3DguueQSZsyYwYIFC1i/fj0HHXTQDu/fuHEjy5Yt2+k+Y8aMoUmTJgwcODAVsSVJkqRqW7sWFiyAn/40dJLalZLiUK9ePY444ggAevTowbx58xg+fDiHH344LVq0YPr06Ttd06RJkx3+HI/Heeihh7j44oupV69elZ8di8XIy9vxYxYWFlJYWFj9DyJJkiTtxcyZUFmZGHHIJikpDt8Wj8fZvHkz3bt3Z+XKleTl5dG2bds9XjNjxgyWLl3K5ZdfXq1nFRcXE41G9yGtJEmSVHWlpdC2LbRrFzpJ7Up6cbjxxhvp168fhx12GOvWraO4uJjp06czefJk+vbtS+/evRk4cCC33347HTt25NNPP6WkpISBAwfSo0eP7fcZNWoUxx9/PEcddVSyI0uSJEk19vX6hmyT9OKwatUqLr74YlasWEF+fj5du3Zl8uTJnHHGGQCUlJTw29/+lssuu4zPPvuM5s2b06dPHw499NDt9ygvL2fcuHEMHz482XElSZKkGlu9GhYtgv/6r9BJal8kHo/HQ4dIhoqKCvLz8ykvL3eqkiRJklLiySfhggvgk0+gZcvQaWpXSs9xkCRJkrJZaSl07Jh9pQEsDpIkSVKtydb1DWBxkCRJkmrFRx/BkiUWB0mSJEl7MG1a4uuppwaNkTQWB0mSJKkWlJbC0UfDwQeHTpIcFgdJkiRpH8XjiRGHbJ2mBBYHSZIkaZ+9/z4sX25xkCRJkrQHpaVQty706RM6SfJYHCRJkqR9VFoKPXpANp87nPXFIRaLUVBQQFFRUegokiRJykLxeHaf3/C1vNABkq24uJhoNlc/SZIkBbV4Maxenf3FIetHHCRJkqRkKi2FevXghBNCJ0kui4MkSZK0D158EXr1ggMOCJ0kuSwOkiRJUg1t3Zo4v+GMM0InST6LgyRJklRD8+ZBRQX07Rs6SfJZHCRJkqQamjIF8vMTW7FmO4uDJEmSVENTpyZ2U8rL+r1KLQ6SJElSjaxbB3Pn5sY0JbA4SJIkSTUyc2ZicXQuLIwGi4MkSZJUI1OmwHe+A0ccETpJalgcJEmSpBqYMiUx2hCJhE6SGhYHSZIkqZo+/RQWL86d9Q1gcZAkSZKqberUxNfTTw+bI5UsDpIkSVI1TZ0K3brBIYeETpI6WV8cYrEYBQUFFBUVhY4iSZKkLBCPJ4pDLk1TAsj6oyqKi4uJRqOhY0iSJClLvPUWrFiRO9uwfi3rRxwkSZKk2jR1KtSvDyedFDpJalkcJEmSpGqYMiVRGho0CJ0ktSwOkiRJUhVt2QIzZuTeNCWwOEiSJElVVlYGGzbk3sJosDhIkiRJVTZ1KjRtmtiKNddYHCRJkqQqmjIlcehbnRz8LToHP7IkSZJUfWvXwiuv5Ob6BrA4SJIkSVUybRpUVubm+gawOEiSJElV8vzz0L49tGsXOkkYFgdJkiRpL+JxmDwZzj47dJJwLA6SJEnSXrz3Hvzzn3DWWaGThGNxkCRJkvbi+eehXj049dTQScLJ+uIQi8UoKCigqKgodBRJkiRlqMmT4eSToWHD0EnCyQsdINmKi4uJRqOhY0iSJClDbdoE06fD738fOklYWT/iIEmSJO2Ll16CjRtze2E0WBwkSZKkPZo8GVq2hKOOCp0kLIuDJEmStAfPPw9nngmRSOgkYVkcJEmSpN345BN4802nKYHFQZIkSdqt559PjDT07Rs6SXgWB0mSJGk3nn8eevaEgw4KnSQ8i4MkSZK0C9u2wZQpTlP6msVBkiRJ2oV58+CLL+Css0InSQ8WB0mSJGkXnn8e8vPhuONCJ0kPFgdJkiRpFyZPhjPOgLy80EnSg8VBkiRJ+pYvvoBXXnGa0jdZHCRJkqRvmToVKistDt+U9cUhFotRUFBAUVFR6CiSJEnKEJMnQ6dOcNhhoZOkj6yfsVVcXEw0Gg0dQ5IkSRkiHk8sjP7xj0MnSS9ZP+IgSZIkVceiRfDJJ9CvX+gk6cXiIEmSJH1DSQk0bAh9+oROkl4sDpIkSdI3lJRA375Qv37oJOnF4iBJkiT9ny++gDlzoH//0EnSj8VBkiRJ+j8vvADbtlkcdsXiIEmSJP2fkhLo2hVatw6dJP1UqzgMGzaMnj170rhxY5o1a8bAgQN5991393jNmDFjiEQiO702bdpU5efu7h6RSITVq1dX5yNIkiRJu1RZCZMmOdqwO9U6x2HGjBlcc8019OzZk61bt/Lb3/6WM888k8WLF9OwYcPdXheNRncqGPvvv3+Vn3vBBRdw9tln7/C9IUOGsGnTJpo1a1adjyBJkiTt0vz58Nln8IMfhE6SnqpVHCZPnrzDn0ePHk2zZs1YsGABffawX1UkEqF58+a7/Xk8HufOO+/kgQceYMWKFXTo0IGbb76ZH/3oRwA0aNCABg0abH//Z599RmlpKaNGjapOfEmSJGm3SkqgSRPo1St0kvS0T2scysvLAWjatOke37d+/XratGlD69atOeecc3jttdd2+PlNN93E6NGjGTFiBG+99RbXX389F110ETNmzNjl/R555BEOOOCA7cVCkiRJ2lclJXDWWZBXrb9azx2ReDwer8mF8Xicc889ly+++IJZs2bt9n1lZWUsXbqULl26UFFRwfDhwykpKeH111+nffv2bNiwgYMPPpjS0lJ69+69/borrriCL7/8kscff3yne3bu3JlTTjmFv/71r7t9bkVFBfn5+ZSXlxONRmvyESVJkpQjVq2C5s3hkUfg4otDp0lPNe5T1157LW+88QYvvfTSHt/Xq1cven1jvOfEE0+ke/fu3Hvvvdxzzz0sXryYTZs2ccYZZ+xw3ZYtW+jWrdtO95s7dy6LFy/mkUceqWl0SZIkaQeTJ0Mkkhhx0K7VqDhcd911TJw4kZkzZ9K6mntV1alTh549e7JkyRIAKisrAXjuuedo1arVDu+tv4vj+kaOHMkxxxzDscceW6XnxWIx8r413lRYWEhhYWG1ckuSJCl7lZRAz57gvju7V63iEI/Hue666xg/fjzTp0+nXbt21X5gPB5n4cKFdOnSBYBOnTpRv359li9fzimnnLLHa9evX8+TTz7JsGHDqvy84uJipypJkiRpt7Zuheefh+uvD50kvVWrOFxzzTU8/vjjPP300zRu3JiVK1cCkJ+fv33Xo8GDB9OqVavtv9z//ve/p1evXrRv356KigruueceFi5cyP333w9A48aNGTp0KNdffz2VlZWcdNJJVFRUMGfOHBo1asQll1yy/flPPPEEW7du5cILL6yVDy9JkiTNmQPl5W7DujfVKg4jRowA4NRTT93h+6NHj2bIkCEALF++nDp1/r1Z09q1a7nyyitZuXIl+fn5dOvWjZkzZ3Lcccdtf88f/vAHmjVrxrBhw3j//fdp0qQJ3bt358Ybb9zhOaNGjWLQoEEceOCB1YktSZIk7VZJSWKKUvfuoZOktxrvqpTu3FVJkiRJVdG1a6I0jBkTOkl626dzHCRJkqRM9tFHsGiR05SqwuIgSZKknFVSAnXrwrdOBtAuWBwkSZKUs559Fk48EZo0CZ0k/VkcJEmSlJO+/BKmToWCgtBJMoPFQZIkSTnpxRdh0yYYMCB0ksxgcZAkSVJOmjgROnRIvLR3FgdJkiTlnMrKxPoGpylVncVBkiRJOWf+fFi50mlK1WFxkCRJUs555hk48EA44YTQSTJH1heHWCxGQUEBRUVFoaNIkiQpTTzzTOLQt7y80EkyRyQej8dDh0iGiooK8vPzKS8vJxqNho4jSZKkNPHPf0LbtvDEE/DjH4dOkzmyfsRBkiRJ+qZnn4X99oOzzgqdJLNYHCRJkpRTJk6EU06B/PzQSTKLxUGSJEk5Y906mD7d3ZRqwuIgSZKknPHCC7Bli8WhJiwOkiRJyhkTJ8JRR0G7dqGTZB6LgyRJknLCtm1QUuJoQ01ZHCRJkpQTyspgzRooKAidJDNZHCRJkpQTJk6EZs3guONCJ8lMFgdJkiTlhK9Pi67jb8A14j82SZIkZb2lS+Htt52mtC8sDpIkScp6EybA/vvDGWeETpK5LA6SJEnKehMmwJlnQsOGoZNkLouDJEmSstqqVTBnDvzwh6GTZLasLw6xWIyCggKKiopCR5EkSVIATz8NkQicc07oJJktEo/H46FDJENFRQX5+fmUl5cTjUZDx5EkSVIg/fvDxo0wbVroJJkt60ccJEmSlLsqKuDFF52mVBssDpIkScpakybBli1w7rmhk2Q+i4MkSZKy1vjx0L07tGkTOknmszhIkiQpK23eDCUlTlOqLRYHSZIkZaXSUli3DgYODJ0kO1gcJEmSlJXGj4cjjoDOnUMnyQ4WB0mSJGWdbdsS5zf88IeJMxy07ywOkiRJyjplZbB6tdOUapPFQZIkSVln/Hg49FDo1St0kuxhcZAkSVJWicdhwoTE2Q11/G231viPUpIkSVnlzTdh2TK3Ya1tFgdJkiRllQkTIBqF738/dJLskvXFIRaLUVBQQFFRUegokiRJSoF//AP694d69UInyS6ReDweDx0iGSoqKsjPz6e8vJxoNBo6jiRJklJg2bLE2Q1jx8J554VOk12yfsRBkiRJuWPcODjgAOjXL3SS7GNxkCRJUtYYOzYxTemAA0InyT4WB0mSJGWFf/4T5s2DH/0odJLsZHGQJElSVhg3DvbfPzHioNpncZAkSVJWGDsWzj4bGjcOnSQ7WRwkSZKU8T7+GObOdZpSMlkcJEmSlPH+8Y/EuQ3nnBM6SfayOEiSJCnjPfUUnHkm5OeHTpK9LA6SJEnKaJ9+CrNnO00p2SwOkiRJymjjx0PdulBQEDpJdrM4SJIkKaONHQt9+8KBB4ZOkt0sDpIkScpYq1bBzJlOU0qFrC8OsViMgoICioqKQkeRJElSLZswASIROPfc0EmyXyQej8dDh0iGiooK8vPzKS8vJxqNho4jSZKkJDjjjMTXKVPC5sgFWT/iIEmSpOy0Zg1Mm+Y0pVSxOEiSJCkjTZgA8TgMHBg6SW6wOEiSJCkjFRfDaafBoYeGTpIbLA6SJEnKOKtWJaYpXXBB6CS5w+IgSZKkjDN2LNSpA4MGhU6SOywOkiRJyjhPPJHYUemgg0InyR1JLw7Dhg2jZ8+eNG7cmGbNmjFw4EDefffdat1jyJAhRCKRnV6dO3dOUmpJkiSlq48/hlmzIBYLnSS3JL04zJgxg2uuuYaysjKmTJnC1q1bOfPMM9mwYUOV7zF8+HBWrFix/fXRRx/RtGlTzj///CQmlyRJUjp68kmoV89D31It5QfAffbZZzRr1owZM2bQp08fAMrLy/nVr37FhAkT2LRpEz169ODPf/4zRx999C7vMWHCBAYNGsQHH3xAmzZtdvkeD4CTJEnKTscfDy1bwvjxoZPklpSvcSgvLwegadOmAMTjcX7wgx+wcuVKSkpKWLBgAd27d+f000/n888/3+U9Ro0aRd++fXdbGiRJkpSd3n8fXnnF3ZRCSGlxiMfj3HDDDZx00kkcddRRAEybNo1Fixbx1FNP0aNHD9q3b89dd91FkyZNGDt27E73WLFiBZMmTeKKK65IZXRJkiSlgSefhAMOgAEDQifJPXmpfNi1117LG2+8wUsvvbT9ewsWLGD9+vUc9K0l8Rs3bmTZsmU73WPMmDE0adKEgR4RKEmSlHOKi+Gcc6Bhw9BJck/KisN1113HxIkTmTlzJq1bt97+/crKSlq0aMH06dN3uqZJkyY7/Dkej/PQQw9x8cUXU69evSo9NxaLkZe348csLCyksLCw+h9CkiRJwbzzDrz+Ovzud6GT5KakF4d4PM51113H+PHjmT59Ou3atdvh5927d2flypXk5eXRtm3bPd5rxowZLF26lMsvv7zKzy8uLnZxtCRJUhZ44glo3Bj69QudJDclfVelq6++mscff5ynn36ajh07bv9+fn4+DRo0IB6P06dPH9atW8ftt99Ox44d+fTTTykpKWHgwIH06NFj+zUXX3wxS5YsoaysbK/PdVclSZKk7BGPQ6dO0LMnPPJI6DS5KemLo0eMGEF5eTmnnnoqLVq02P564oknAIhEIpSUlNCnTx8uu+wyOnToQCwW48MPP+TQQw/dfp/y8nLGjRtXrdEGSZIkZYdFixJTlTz0LZyUn+OQKo44SJIkZY8bb4QHHoCVKxOHvyn1Un6OgyRJklQd8XhifcOgQZaGkCwOkiRJSmtlZYmD3y68MHSS3GZxkCRJUlp77DFo1QpOOSV0ktxmcZAkSVLa+uqrxDSlwkKo42+uQfmPX5IkSWnrhRdgzRq46KLQSWRxkCRJUtp67DHo3Bm6dg2dRBYHSZIkpaV162DChMSi6EgkdBpZHCRJkpSWJkyAjRvhJz8JnURgcZAkSVKaeuwxOPn3RxMUAAAgAElEQVRkaNMmdBKBxUGSJElpaNUqmDLFsxvSSdYXh1gsRkFBAUVFRaGjSJIkqYqKi6FuXTj//NBJ9LVIPB6Phw6RDBUVFeTn51NeXk40Gg0dR5IkSdVw3HHQsmVinYPSQ17oAJIkSdI3LVkC8+bBk0+GTqJvyvqpSpIkScosjz0GjRvDOeeETqJvsjhIkiQpbcTj8OijcN550KBB6DT6JouDJEmS0sYrr8CyZe6mlI4sDpIkSUobDz8MrVrBaaeFTqJvszhIkiQpLWzaBEVFcPHFia1YlV4sDpIkSUoLzzwDa9fCJZeETqJdsThIkiQpLTz8MBx/PHzve6GTaFcsDpIkSQpu5UqYPNnRhnRmcZAkSVJwjz2WWNdwwQWhk2h3LA6SJEkKKh5PTFM691xo2jR0Gu2OxUGSJElBLVwIixY5TSndZX1xiMViFBQUUFRUFDqKJEmSduHhh+HQQ+Gss0In0Z5E4vF4PHSIZKioqCA/P5/y8nKi0WjoOJIkSdqFLVsSB74NHgx33x06jfYk60ccJEmSlL4mTYI1a5ymlAksDpIkSQrm4YehWzfo2jV0Eu2NxUGSJElB/Otf8OyzjjZkCouDJEmSgigqSmzF+pOfhE6iqrA4SJIkKYjRo6F/fzjkkNBJVBUWB0mSJKXcwoXw6qtw+eWhk6iqLA6SJElKuVGjoHnzxIiDMoPFQZIkSSm1cSM8+igMGQJ5eaHTqKosDpIkSUqp8eNh7Vq47LLQSVQdFgdJkiSl1MiR0KcPtG8fOomqw8EhSZIkpcyyZTBtGjzySOgkqi5HHCRJkpQyo0dDNArnnRc6iaor64tDLBajoKCAoqKi0FEkSZJy2tatieLwk5/AAQeETqPqisTj8XjoEMlQUVFBfn4+5eXlRKPR0HEkSZJy3nPPwTnnwPz5cOyxodOourJ+xEGSJEnpYeRIOPpo6N49dBLVhMVBkiRJSbdyJTz7LFxxBUQiodOoJiwOkiRJSrpHHoG6deHCC0MnUU1ZHCRJkpRU8TiMGgWDBsGBB4ZOo5qyOEiSJCmpXnoJ3nsPLr88dBLtC4uDJEmSkuqBB+CII+C000In0b6wOEiSJClp1qyBsWPhpz+FOv7mmdH81ydJkqSkGTMm8XXIkJApVBssDpIkSUqKykr429/g/PPh4INDp9G+ygsdQJIkSdmptBSWLoXRo0MnUW1wxEGSJElJ8cAD0LkznHhi6CSqDRYHSZIk1boVK2DChMSiaE+Kzg5ZXxxisRgFBQUUFRWFjiJJkpQzHnoI6tWDiy8OnUS1JRKPx+OhQyRDRUUF+fn5lJeXE41GQ8eRJEnKGdu2weGHQ9++iROjlR2yfsRBkiRJqTV5MixfDlddFTqJapPFQZIkSbXqgQege3fo0SN0EtUmi4MkSZJqzfLlUFKSGG1wUXR2sThIkiSp1jz4IDRsCIWFoZOotlkcJEmSVCu2bEkUh4svhkaNQqdRbbM4SJIkqVaMHQurVsG114ZOomSodnGYOXMmAwYMoGXLlkQiESZMmLDH90+fPp1IJLLT65133qnyM3d3j0gkwrx586r7ESRJkpQE994Lp58ORx4ZOomSIa+6F2zYsIGjjz6aSy+9lPPOO6/K17377rs7nKdwyCGHVPnaE044gRUrVuzwvZtvvpmpU6fSw+X6kiRJwc2fD2VlidOilZ2qXRz69etHv379qv2gZs2a0aRJk93+fPTo0dxxxx188MEHtG3blp/97GdcffXVANSrV4/mzZtvf+9XX33FxIkTufbaa4m4XF+SJCm4++6DNm3gnHNCJ1GypGyNQ7du3WjRogWnn34606ZN2+FnDz74IL/97W/54x//yNtvv82f/vQnbr75Zh5++OFd3mvixImsWbOGIUOGpCC5JEmS9uSzz6C4GK6+GurWDZ1GyVLtEYfqatGiBX//+9859thj2bx5M//7v//L6aefzvTp0+nTpw8Af/jDH7j77rsZNGgQAO3atWPx4sX87W9/45JLLtnpnqNGjeKss87isMMOS3Z8SZIk7cXIkYkzGy6/PHQSJVMkHo/Ha3xxJML48eMZOHBgta4bMGAAkUiEiRMn8tlnn9GsWTMaNGhAnTr/HgDZunUr+fn5rFq1aodrP/74Y9q0acOTTz65xzUWFRUV5OfnU15evsPaCkmSJNWerVvh8MPhjDNg1KjQaZRMSR9x2JVevXrx6KOPAlBZWQkkpisdf/zxO7yv7i7GukaPHs1BBx1EQUFBlZ4Vi8XIy9vxYxYWFlLoqSSSJEn7bOJE+Ogjt2DNBUGKw2uvvUaLFi0AOPTQQ2nVqhXvv/8+F1544R6vi8fjjB49msGDB7PffvtV6VnFxcWOOEiSJCXJfffBiSdCt26hkyjZql0c1q9fz9KlS7f/+YMPPmDhwoU0bdqU73znO/zmN7/hk08+4ZFHHgHgL3/5C23btqVz585s2bKFRx99lHHjxjFu3Ljt97jlllv42c9+RjQapV+/fmzevJn58+fzxRdfcMMNN2x/X2lpKR988AGXO4FOkiQpuLfegmnToKgodBKlQrWLw/z58znttNO2//nrX+wvueQSxowZw4oVK1i+fPn2n2/ZsoWhQ4fyySef0KBBAzp37sxzzz1H//79t7/niiuu4IADDuDOO+/k17/+NQ0bNqRLly784he/2OHZo0aN4oQTTuBITxWRJEkK7r77oEUL+L/9bZTl9mlxdDpzcbQkSVLyrF0LrVvDr34Fv/td6DRKhZSd4yBJkqTs8eCD8NVX8NOfhk6iVLE4SJIkqVq2boV774Wf/ASaNw+dRqlicZAkSVK1jBuX2IL1W8tRleVc4yBJkqQqi8ehVy9o3BimTg2dRqkU5BwHSZIkZaa5c+GVV+DZZ0MnUao5VUmSJElV9j//Ax07Qr9+oZMo1RxxkCRJUpV88AGMHw/33w91/OvnnOO/ckmSJFXJPfdAkyYweHDoJArB4iBJkqS9Ki+HUaPgqqvggANCp1EIFgdJkiTt1ahRsGkTXHNN6CQKJeuLQywWo6CggKKiotBRJEmSMtLWrYlpShdcAC1bhk6jULJ+cXRxcbHnOEiSJO2D8ePhn/+Ef/wjdBKF5AFwkiRJ2q14HI4/Hho2hGnTQqdRSFk/4iBJkqSamzED5s2DSZNCJ1FoWb/GQZIkSTV3++3QtSucdVboJArNEQdJkiTt0htvwOTJ8OijEImETqPQHHGQJEnSLt1xB7RpAz/+cegkSgcWB0mSJO3kww+huBhuuAH22y90GqUDi4MkSZJ28uc/Q34+XH556CRKFxYHSZIk7eBf/4KRI+HaaxPbsEpgcZAkSdK33H9/4vyG664LnUTpxOIgSZKk7b78Eu69NzFF6eCDQ6dROrE4SJIkabvRo+GLLxKLoqVvsjhIkiQJgK++grvuSmy/2q5d6DRKN1l/AFwsFiMvL4/CwkIKCwtDx5EkSUpbjz+e2IZ14sTQSZSOIvF4PB46RDJUVFSQn59PeXk50Wg0dBxJkqS0tm0bdOqUeI0fHzqN0lHWjzhIkiRp78aOhffeg8ceC51E6coRB0mSpBxXWQlHHw2tW8OkSaHTKF054iBJkpTjJk6EN9+EESNCJ1E6c8RBkiQph8Xj0LMnNGoE06eHTqN05oiDJElSDnv+eViwAKZMCZ1E6c4RB0mSpBwVj8NJJyV2VJo7FyKR0ImUzhxxkCRJylEzZsCcOfDMM5YG7Z0jDpIkSTmqb19YswZee83ioL1zxEGSJCkHzZkDL74ITz5paVDVOOIgSZKUg/r2hVWr4PXXoU6d0GmUCRxxkCRJyjEzZyZGG8aNszSo6hxxkCRJyjGnnQZr1ya2YbU4qKqyfsQhFouRl5dHYWEhhYWFoeNIkiQFNW1a4qC3p5+2NKh6HHGQJEnKEfE49OkDGzfCvHkuilb1ZP2IgyRJkhKmToWXXoLnnrM0qPoccZAkScoB8TiccELiq6dEqyYccZAkScoBkydDWVniq6VBNeGIgyRJUpaLx+G446BevcRUJYuDasIRB0mSpCz33HMwf35ijYOlQTXliIMkSVIWq6yEbt2gSZPENqwWB9WUIw6SJElZrKgI3ngDZs+2NGjfOOIgSZKUpbZsge99D7p0SRz4Ju0LRxwkSZKy1IMPwocfwsSJoZMoG3jQuCRJUhZavx7+8AcYPBiOOip0GmUDi4MkSVIWGj4cvvgCbrkldBJlC4uDJElSlvnXv+COO+A//xPatg2dRtki64tDLBajoKCAoqKi0FEkSZJS4rbbEtuw3nhj6CTKJu6qJEmSlEU+/hiOOAL++7+dpqTalfUjDpIkSbnkllugcWO44YbQSZRt3I5VkiQpS7z5JoweDXffDU64UG1zqpIkSVKW6N8f3nsPFi+GevVCp1G2ccRBkiQpC0yZApMmwdixlgYlhyMOkiRJGW7bNujWLTE9adYsiERCJ1I2csRBkiQpw40ZA4sWQVmZpUHJk5JdlWbOnMmAAQNo2bIlkUiECRMmVOv6W265hUgkstOrYcOGSUosSZKUGdavh5tvhlgMjj8+dBpls5QUhw0bNnD00Udz33331ej6oUOHsmLFih1enTp14vzzz6/lpJIkSZnlrrvg889h2LDQSZTtUlIc+vXrx6233sqgQYN2+fMtW7bw61//mlatWtGwYUOOP/54pk+fvv3njRo1onnz5ttfq1atYvHixVx++eWpiC9JkpSWPv0U7rwTfv5zaNs2dBplu7RY43DppZfy4YcfUlxcTMuWLRk/fjxnn302ixYton379ju9f+TIkXTo0IGTTz45QFpJkqT0cNNN0KAB3Hhj6CTKBcFPjl62bBlFRUU89dRTnHzyyXz3u99l6NChnHTSSYwePXqn92/evJnHHnvM0QZJkpTTFi5MLIq+5RbIzw+dRrkg+IjDq6++Sjwep0OHDjt8f/PmzRx00EE7vf8f//gH69atY/DgwamKKEmSlFbicbjuOjjySPjpT0OnUa4IXhwqKyupW7cuCxYsoG7dujv8rFGjRju9f+TIkZxzzjk0b968SvePxWLk5e34MQsLCyksLKx5aEmSpICeeAJeeilx6Nt++4VOo1wRvDh069aNbdu2sXr16r2uWfjggw+YNm0aEydOrPL9i4uLPQBOkiRljQ0bYOhQ+OEPoW/f0GmUS1JSHNavX8/SpUu3//mDDz5g4cKFNG3alA4dOnDhhRcyePBg7r77brp168aaNWsoLS2lS5cu9O/ff/t1Dz30EC1atKBfv36piC1JkpR2hg2DNWvg7rtDJ1GuSUlxmD9/Pqeddtr2P99www0AXHLJJYwZM4bRo0dz66238stf/pJPPvmEgw46iN69e+9QGiorKxkzZgxDhgzZaUqTJElSLnj//cS5Db/6FbRrFzqNck0kHo/HQ4dIhoqKCvLz8ykvL3eqkiRJygo//CHMnw/vvAMNG4ZOo1wTfI2DJEmS9m7KFJgwAYqKLA0KwxEHSZKkNPfVV3D00XDwwTBjBkQioRMpFzniIEmSlObuuw/efRcef9zSoHCCnxwtSZKk3fv4Y/h//w+uugqOOSZ0GuUyi4MkSVIau/76xJqGP/4xdBLlOqcqSZIkpanJk2HsWHjsMWjSJHQa5ToXR0uSJKWhjRvhqKMS5zVMmeLaBoXniIMkSVIaGjYssb6hpMTSoPTgGgdJkqQ08957cPvt8OtfQ8eOodNICU5VkiRJSiPxOJxxBnzwAbz5JjRoEDqRlJD1U5VisRh5eXkUFhZSWFgYOo4kSdIeFRfDiy8mpihZGpROHHGQJElKE59/Dp06wUknJXZTktKJaxwkSZLSxNChsGkT3HNP6CTSzrJ+qpIkSVImmDoVRo+Gv/8dWrYMnUbamVOVJEmSAtuwAbp0gTZtoLTU7VeVnhxxkCRJCuzmm2HFCnjhBUuD0pfFQZIkKaBXXoHhwxMHvh1xROg00u45VUmSJCmQLVvg2GOhXj14+WXI8690lcb8z1OSJCmQ22+Ht9+G+fMtDUp/bscqSZIUwFtvwa23wq9/DcccEzqNtHdOVZIkSUqxr76CXr1g40Z49VXYf//QiaS9c1BMkiQpxf70J3j9dSgrszQoczhVSZIkKYVefTUxRenGG6FHj9BppKpzqpIkSVKKbN6c2EVpv/0SuyjVqxc6kVR1WT9VKRaLkZeXR2FhIYWFhaHjSJKkHPa738F778GCBZYGZR5HHCRJklJg7lw46aTENKXf/CZ0Gqn6LA6SJElJ9uWXiS1XmzaFl17yzAZlJv+zlSRJSrJf/xo++gieecbSoMzlf7qSJElJ9OyzcP/9iVfHjqHTSDXnVCVJkqQkWbkSunRJHPY2cSJEIqETSTXnOQ6SJElJUFkJQ4ZA3bowapSlQZnPqUqSJElJcO+98PzzMGkSNGsWOo207xxxkCRJqmVvvJFYEP3zn8PZZ4dOI9UO1zhIkiTVoo0boWdPqFMHXnkF9t8/dCKpdjhVSZIkqRb96lewdCnMn29pUHaxOEiSJNWSsWP/vfXqUUeFTiPVLtc4SJIk1YJly+Dyy+H88+E//zN0Gqn2WRwkSZL20ebN8OMfwyGHwMiRbr2q7JT1U5VisRh5eXkUFhZSWFgYOo4kScpCQ4fCm29CWRm4J4uylbsqSZIk7YOnnkqMNvz1r05RUnZzqpIkSVINLV2aWNdwwQVw1VWh00jJZXGQJEmqgY0bE4Xh0EPh7393XYOyX9avcZAkSapt8ThcfTUsXgxz57quQbnB4iBJklRNI0bAmDHwv/8LxxwTOo2UGk5VkiRJqobZs+HnP4ef/Qwuuih0Gil13FVJkiSpij79FI49Ftq3hxdfhP32C51ISh1HHCRJkqpgy5bEqdB16sCTT1oalHtc4yBJklQFv/gFzJsHM2dC8+ah00ipZ3GQJEnaiwcfTCyI/tvfoFev0GmkMJyqJEmStAfTpye2Xr3qKviP/widRgrH4iBJkrQbS5fCeefBKafAPfd4yJtym8VBkiRpF9auhQED4OCD4amnXAwtZX1xiMViFBQUUFRUFDqKJEnKEFu3wgUXwMqV8MwzcOCBoRNJ4WX94uji4mLPcZAkSdXyy18mzml4/nno0CF0Gik9ZH1xkCRJqo4RIxLrGUaMgNNPD51GSh9ZP1VJkiSpqp55Bq69Fn72s8QuSpL+zeIgSZIEvPxyYl3DwIHwP/8TOo2UfiwOkiQp5y1dCuecA927w6OPQt26oRNJ6cfiIEmSctpnn8HZZ8NBB8HTT0ODBqETSenJxdGSJClnffllYqRh/XqYOzdRHiTtWtJHHG655RYikcgOr+bNm1frHvF4nLvuuosOHTpQv359DjvsMP70pz8lKbEkScoFW7dCLAZvvQXPPQft2oVOJKW3lIw4dO7cmalTp27/c91qThz8+c9/zgsvvMBdd91Fly5dKC8vZ82aNbUdU5Ik5YjKSrjsMpg0KbGT0rHHhk4kpb+UFIe8vLzdjjJs2bKFm266iccee4y1a9dy1FFHcfvtt3PqqacC8PbbbzNixAjefPNNOnbsmIq4kiQpi8Xj8ItfJBZBFxUl1jdI2ruULI5esmQJLVu2pF27dsRiMd5///3tP7v00kuZPXs2xcXFvPHGG5x//vmcffbZLFmyBIBnnnmGww8/nGeffZZ27drRtm1brrjiCj7//PNURJckSVnmllvg3nvhgQcS269KqppIPB6PJ/MBkyZN4ssvv6RDhw6sWrWKW2+9lXfeeYe33nqLtWvX0r59ez7++GNatmy5/Zq+ffty3HHH8ac//YmrrrqKMWPGcMwxx3DnnXeybds2rr/+eg488EBKS0t3+9yKigry8/MpLy8nGo0m8yNKkqQM8ec/ww03wG23wX/9V+g0UmZJ+lSlfv36bf/fXbp0oXfv3nz3u9/l4Ycf5rDDDiMej9OhQ4cdrtm8eTMH/d+2BpWVlWzevJlHHnlk+/tGjRrFsccey7vvvuv0JUmSVCUPPZQoDf/935YGqSZSvh1rw4YN6dKlC0uWLKFVq1bUrVuXBQsW7LRgulGjRgC0aNGCvLy8HcrFkUceCcDy5cv3WhxisRh5eTt+zMLCQgoLC2vj40iSpAzwxBPwH/8BV10Fbswo1UzKi8PmzZt5++23Ofnkk+nWrRvbtm1j9erVnHzyybt8/4knnsjWrVtZtmwZ3/3udwF47733AGjTps1en1dcXOxUJUmScthTT8GFFyZe990HkUjoRFJmSvoah6FDhzJgwAC+853vsHr1am699VZmzJjBokWLaNOmDRdddBGzZ8/m7rvvplu3bqxZs4bS0lK6dOlC//79qayspGfPnjRq1Ii//OUvVFZWcs011xCNRnnhhRd2+1zXOEiSpHHjEgugYzF4+GGo5o7wkr4h6bsqffzxxxQWFtKxY0cGDRpEvXr1KCsr2z5aMHr0aAYPHswvf/lLOnbsSEFBAS+//DKHHXZYImCdOjzzzDMcfPDB9OnThx/84AcceeSRFBcXJzu6JEnKYOPHJwrD+efDmDGWBmlfJX3EIRRHHCRJyl1PPw0/+hEMGgSPPQZ5KZ+cLWWflJzjIEmSlCoTJyZGGc49N3HIm6VBqh0WB0mSlDWeeALOOw8GDEicCr3ffqETSdnD4iBJkrLC6NHwk58k1jU88YSlQaptFgdJkpTx7rsPLrsscVbDww87PUlKBouDJEnKaLfdBtddlzgVesQIqONvN1JS+H8tSZKUkeJxuOkm+M1v4He/g7vu8nA3KZkcyJMkSRln61a4+mp48EG44w741a9CJ5Kyn8VBkiRllC+/TCyALilJHOx2ySWhE0m5weIgSZIyxr/+ldhq9fXX4ZlnoF+/0Imk3GFxkCRJGeHDD+Hss+Hzz2H6dOjZM3QiKbdk/eLoWCxGQUEBRUVFoaNIkqQaWrgQTjgBtmyB2bMtDVIIkXg8Hg8dIhkqKirIz8+nvLycaDQaOo4kSaqhp5+GCy+E730Pnn0WmjcPnUjKTVk/4iBJkjJTPA533gk//GFiitLM/9/evYdFVed/AH8PyEUURlm5iBaaueoTeYMSUjNvqKui1j4hGY/7WJvmUqZbednKfm4+5uaSXVZ37WZZgRe0fMpQXFAeA2vFQUWCvOUtkESYMShg4PP749sMjqBcnAszvF/Pc54znPnOmc93vp5z/JzvOd+TyaSByJGYOBAREVGbU10NPP448Pzz6jkNW7YAPj6OjoqofePN0URERNSmlJYCDz0EZGcDH30ExMc7OiIiApg4EBERURty7Ji6NEmvB/77X2DECEdHREQmvFSJiIiI2oTkZCAyEujcGfj2WyYNRG0NEwciIiJyqJoaYNEiIC5O9TZkZQG9ezs6KiK6Hi9VIiIiIoe5dAmIjVXPZnjzTSAhAdBoHB0VETWGiQMRERE5RFYW8PDDgNGo7me4/35HR0REN8NLlYiIiMiu6uqAV19ViUJoKJCTw6SByBkwcSAiIiK7uXQJmDQJWLZMPaNh3z6gRw9HR0VEzcFLlYiIiMgu0tOBWbNUj8Pu3cD48Y6OiIhagj0OREREZFM1NcCLLwLjxgF33QUcOcKkgcgZuXziMHPmTMTExCApKcnRoRAREbU7BQXAffcBq1YBf/+76mkIDnZ0VETUGhoREUcHYQsGgwFarRZ6vR5+fn6ODoeIiKhdqasD3noLWLIE6NUL2LQJiIhwdFREdCtcvseBiIiI7Ov8eSA6GnjmGeCJJ4DDh5k0ELkC3hxNREREViECbNwILFwI+PoCaWnqvgYicg3scSAiIqJbdvq0uuF5zhxg2jTg2DEmDUSuhokDERERtVptLZCYCISFASdPAqmpwIcfAl26ODoyIrI2Jg5ERETUKkePAlFRwLPPqnsZ8vKACRMcHRUR2QoTByIiImqRq1dVshAeDlRWAllZwNq1QOfOjo6MiGyJiQMRERE1iwjw6adAv37AunXA//2fGjEpMtLRkRGRPTBxICIioibl5QGjRwOzZqkHuhUUAMuWAZ6ejo6MiOyFiQMRERHdUGmpeh7D4MFAcTGwZw+wbRtw++2OjoyI7I3PcSAiIqIGqqqAt98GXnlFjZy0cqV6PgN7GIjaLyYOREREZCYCbN4MLF2qngD9xBPAyy8DgYGOjoyIHI2XKhEREREAYP9+daNzXBwwcKC6r2HdOiYNRKS4fOIwc+ZMxMTEICkpydGhEBERtUnZ2eopzw88oC5LSk8HPv8c6N/f0ZERUVuiERFxdBC2YDAYoNVqodfr4efn5+hwiIiI2pxDh4CXXgK++ko9+XnFCmD6dECjcXRkRNQWuXyPAxEREVnS6VSCcM89wJkzQHIycOQIMGMGkwYiujEmDkRERO2ACLBvHzBxIjB0qLp/4aOP1Dw2FnDj/wiIqAncTRAREbmwujp1v0JUlHqAW3ExkJSkHuAWHw+4uzs6QiJyFkwciIiIXNCvvwIffADcfbe6LMnTE9i1S12mNHMm0IEDshNRC3G3QURE5EIuXADWrwc2bFBPfZ4yRb0ePtzRkRGRs2PiQERE5OREgKws4M03gZQUwMcHmDMHSEgA7rzT0dERkatg4kBEROSk9Hrg00+Bd95RlyD17QusXQvMng34+jo6OiJyNUwciIiInIgI8PXXKlnYuhWorgYmTwZeeUWNmMTRkYjIVpg4EBEROYEff1S9C++9p0ZEuuMO4MUXVe9CSIijoyOi9oCJAxERURul16t7Fj75BMjIUCMjPfgg8K9/AQ88wN4FIrIvJg5ERERtSFWVGjb1k0+AL75QlyKNHg28+65KGrp0cXSERNReMXEgIiJysKtXga++ArZvV0nD1avAkCHAypXqmQs9ejg6QiIiQCMi4uggbMFgMECr1WLSpEno0KED4uLiEBcX5+iwiIiIAACXLwM7dwI7dgBpaaqnYcgQ1avw0EPAgAGOjpCIyJLLJw56vR5+fn6ODoX/KJYAABMISURBVIeIiNq5ujogN1f1LHz1FZCdrUZIGj5cJQvTpwO9ezs6SiKiG+OlSkRERDZy5QqwZ49KFHbvBi5dAjp3BsaNU093jokBgoMdHSURUfMwcbCC7Gzgyy+B++8HoqL40B0iovZKrwcOHAD27VOjIOl0qqfh7rvVsKmTJgH33adGRyIicjZMHKzg5En1IJ6VKwF3d2DoUGDUKJVIjBgBdO3q6AiJiMgWysuBrKz6ROHwYZUo9OihRkKaPx+YMIE3NxORa+A9DlYiAnz/PbB/P5CZqeYXLgAajTrTFBkJ3HsvMGyYuuHN3d3mIRERkRXV1gL5+aqX+eBBNX33nXqve3eVKIwerZ6v0KeP2v8TEbkSJg42IgKcPauSiMxM4NtvgePH1ZmoTp2AiIj6RCI8HAgN5UGGiKitqKsDTp1SNzPrdMA336j9+M8/q4euDRyoTghFRqpLj+68k/twInJ9TBzs6OefgZwcdfAxHYTOn1fvabXqQDRoUP08LAzw8XFszERErq6qSp3YMSUJubnAkSPqWQqA6k249976RCEiQt3gTETU3jhd4rBu3Tq89tprKCoqwl133YW1a9di5MiRDcrZO3FISkpq1XMiiorUNbFHj6oD1ZEj6pKnujp1VqtvX+Cuu4D+/dXUr5+atFobVKIda237UdvA9nNu9mo/gwEoKFCXF+Xnq/l33wGnT6t9rkaj9q+DB6vnKQwerKbAQJuH5vS4DTo3tp9zs2f7OVXisHnzZsTHx2PdunUYPnw4/vOf/+Ddd99Ffn4+br/9douy9k4cYmJisHPnTqusq7JSHdRMiUR+PlBYqO6ZMOneXR3g+vdXyUXv3vUTk4qWs2b7kf2x/ZybNdvvyhWVCJimU6fUvLAQuHixvlxoqLrfzDSFhane3k6drBJGu8Nt0Lmx/ZybPdvPqUZVSkxMxGOPPYbHH38cALB27Vrs3r0b69evx6pVqxwcnfX4+Kiu8IgIy+VXr6reiMJCddassBD4+mvgww+BX36pL9e1q2Ui0bs30LOnGtWjRw8gIIA3ZxOR86muVr20Fy6oJMA0P3euPlEoL68vr9Wqm5T79FFDoZqShP79mSAQEbWGm6MDaK7q6mrk5OQgOjraYnl0dDSysrIcFJV9JCUlAVDPhwgPBx55BFixAti8WV3iVFEBFBerkT6SkoDnnlPl9Hpgxw5gwQJg2jSViHTvDnh5AbfdBtx5ZxIefBBISABWrQLefx/YuVOt5+RJdQBubX+UKWZrs9V6bckZfwtnjNlWnPG3cKaYjUZ1j8GxY0B6OpCcDLz1FvDii8C8efX7ruBgwNsb6NVLDXMdG6vKfPEFUFam9nmLFwNbtgCHDqmeh/Jy4Nlnk7Blixou+9FHVTlrJA1sP/twxt/CGWO2FWf8LZwxZntymh6Hy5cvo7a2FkFBQRbLg4KCUFxc7KCo7KOpa9c0GiAoSE2RkQ3fr60FfvpJnZn78Uc1v3gR+PDDJFRWxiEzU/195UrDz3p4AN26qSkgQM27dFGTVls/Xf/3Rx8lITY2Dm5WTk2d8TpMW8Vsy9/CGWO2FWf8LewZc3W1OnlRWakGgNDr1VRefvPpyhW1X7pyRZ2gGDiwfp1eXmp/Exio9mtDh6onLPfoUd972rOn2tc0NZIR28/267UlZ/wtnDFmW3HG38IZY7Ynp0kcTDTXHSVEpMGyWyEiuGoaSqMFjEYjDAaD1eKw9rp9fNS9EH371i/LyTEiObl+vTU16iBeWgpcvqwm0+vSUjUVF6tLpAwG9Z8Dg0GdMWwkari7G+Dlpb7bxwfo2LF+3qmT5d8dO6qziZ6eKlnx9Kx/bfrbwwM4f17FbFrWoYO6iby1E6D+02IwGHH0qAEi9b0s186vf33tdKPldXVqfvmyEenpBvPfpvm1rxt7r6ky588b8cEHhkY/19zvAhovU1hoxMsvNx5zc2K/0Xfk5hoxZ47B4j2NRrXF9XOg8eWNLfvuOyOWLTPccF2NzU1TU+XOnjXinXcMN4znZvObfaakxIjduw1N1lmjUb+VaTL9djebioqM2LrVYP67tlZt3zU16j/61dXNe11drRKCykqVHBQUGNGnj8GcKFRWqnXfjOlkgp9f/evQUDVynOlkxMaNRqxebTCfpPD1bd7Qps3ZVdtq39zW9/n2XK8t1+1s67Xluhmz7ddry3W3xZh9fX1b/H9op7k5urq6Gj4+Pti6dStmzJhhXr5gwQLk5uZi//79FuVNN0dPmjQJHTpY5kdxcXE3zPpMnyMiIiIiclWtGUDIaXocPD09ER4ejrS0NIvEIS0tDdOmTbvh55KTk1v0o/j6+kKv199SrGRdIqpXw3Q29Nozokajer+2tvEzr9efnb22XG1t/VnNa89Am6bGljdV1jRv6Vlva53Fvtn6qXla0pvSnB6Yxj7b0s+4u9+850yjuXkZjcayF4//HoiIyNfXt8WfcZrEAQAWLVqE+Ph4REREICoqChs2bMC5c+cwb948q32HRqNpMw+MIyIiIiJqK5wqcYiNjUVpaSlWrFiBoqIihIWFYdeuXQgNDXV0aERERERELs1p7nFoKXs/AI6IiIiIyJU5zXMciIiIiIjIcZg4EBERERFRk5g4/CYzMxNTp05FSEgINBoNPvvsM4v3L126hD/96U8ICQmBj48PJk6ciBMnTliUOXXqFGbMmIGAgAD4+fnh4YcfxqVLlyzKlJWVIT4+HlqtFlqtFvHx8SgvL7d5/VydvdqvV69e0Gg0FtOSJUtsXj9Xt2rVKtxzzz3w9fVFYGAgpk+fjsLCQosyVVVVeOqpp9CtWzd06tQJMTExuHDhgkWZc+fOYerUqejUqRO6deuGp59+GtXV1RZl9u/fj/DwcHh7e+OOO+7Av//9b5vXz9XZq/327dvXYPvTaDQoKCiwSz1dlbXab8GCBQgPD4eXlxcGDx7c6HcdO3YMo0aNQseOHdGjRw+sWLECLnrFtN3Yq/1++OGHRre/1NRUm9bP1Vmj/Y4cOYK4uDjcdttt6NixIwYMGIA33nijwXdZ4/jHxOE3FRUVGDRoEN5+++0G74kIpk+fjtOnT+Pzzz+HTqdDaGgoxo0bh4qKCvPno6OjodFokJ6ejq+//hrV1dWYOnUq6kxPugLwyCOPIDc3F6mpqUhNTUVubi7i4+PtVk9XZa/2A2C+Od80vfDCC3apoyvbv38//vKXv+DgwYNIS0uD0WhEdHS0uX0A4JlnnsGOHTuQnJyMAwcO4Oeff8aUKVNQ+9sTyGprazF58mRUVFTgwIEDSE5ORkpKCv7617+a13HmzBn84Q9/wMiRI6HT6bBs2TI8/fTTSElJsXudXYm92s+ksLDQYhvse+2TLanFrNF+gNrXzpkzB7GxsY1+j8FgwPjx4xESEoL//e9/eOutt7BmzRokJibavI6uzF7tZ7J3716L7W/MmDE2q1t7YI32y8nJQUBAAD7++GMcP34cf/vb37B06VKL/xNZ7fgnLkqv1wsA0ev1Lf4sANmxY4f578LCQgEgeXl55mVGo1H8/f3lnXfeERGR3bt3i5ubm8X3XblyRQBIWlqaiIjk5+cLADl48KC5THZ2tgCQgoKCFsdJjbNV+4mIhIaGyuuvv26HWrRvJSUlAkD2798vIiLl5eXi4eEhycnJ5jIXL14UNzc3SU1NFRGRXbt2iZubm1y8eNFcJikpSby8vMzt+vzzz0v//v0tvmvu3LkSGRlp6yq1K7Zqv4yMDAEgZWVldqxN+9Oa9rvW8uXLZdCgQQ2Wr1u3TrRarfz666/mZatWrZKQkBCpq6uzQU3aJ1u135kzZwSA6HQ62wVPt9x+JvPnz5fRo0eb/7bW8Y89Ds1QVVUFAPD29jYvc3d3h6enJw4cOGAuo9Fo4OXlZS7j7e0NNzc3c5ns7GxotVoMGzbMXCYyMhJarRZZWVn2qEq7ZK32M1m9ejV+97vfYfDgwVi5cmWDS2Ho1pkewujv7w9AnU2pqalBdHS0uUxISAjCwsLM2052djbCwsIQEhJiLjNhwgRUVVUhJyfHXObadZjKHDp0CDU1NTatU3tiq/YzGTJkCLp3746xY8ciIyPD1tVpd1rTfs2RnZ2NUaNGWexnJ0yYgB9//BE//PCDdYInm7WfSUxMDAIDAzF8+HBs27bNOkGTmbXaT6/Xm9cBWO/4x8ShGfr374/Q0FAsXboUZWVlqK6uxquvvori4mIUFRUBUAlAp06dsHjxYlRWVqKiogLPPfcc6urqzGWKi4sRGBjYYP2BgYEoLi62a53aE2u1H6CuAU1OTkZGRgYSEhKwdu1azJ8/31FVc0kigkWLFmHEiBEICwsDoLYdT09PdO3a1aJsUFCQedspLi5GUFCQxftdu3aFp6fnTcsEBQXBaDTi8uXLtqpSu2LL9uvevTs2bNiAlJQUbN++Hf369cPYsWORmZlph5q1D61tv+a40fZneo9unS3br3PnzkhMTMS2bduwa9cujB07FrGxsfj444+tWof2zFrtl52djS1btmDu3LnmZdY6/jnVA+BawtfXF3q9vlWP076eh4cHUlJS8Nhjj8Hf3x/u7u4YN24cJk2aZC4TEBCArVu34sknn8Sbb74JNzc3xMXFYejQoXB3dzeX02g0DdYvIo0uJ+uwZvstXLjQ/HrgwIHo2rUr/vjHP5p7IejWJSQk4OjRow16ehpz/bbTnO3r+jLy242Z3Aatw5bt169fP/Tr18/8XlRUFM6fP481a9bg/vvvt0L0dCvt1xzc/mzLlu3XrVs3i2NgREQEysrK8I9//AOPPvpoq+IlS9Zov+PHj2PatGl46aWXMH78eIv3rLH9uWyPg0ajgZ+fn9V2RuHh4cjNzUV5eTmKioqQmpqK0tJS9O7d21wmOjoap06dQklJCS5fvoxNmzbh4sWL5jLBwcENRukBgJ9++qlBFkjWZY32a0xkZCQA4OTJkzavQ3vw1FNPYefOncjIyEDPnj3Ny4ODg1FdXY2ysjKL8iUlJeZtJzg4uMHZl7KyMtTU1Ny0TElJCTp06MDEzwps3X6NiYyMbDBCGrXOrbRfc9xo+wPAY6AV2Lr9GsPtz3qs0X75+fkYM2YM/vznPzcYuMVaxz+XTRxsRavVIiAgACdOnMChQ4cwbdq0BmW6deuGLl26ID09HSUlJYiJiQGgzo7p9Xp8++235rLffPMN9Ho97rvvPrvVoT27lfZrjE6nA6AuoaDWExEkJCRg+/btSE9Pb5CshYeHw8PDA2lpaeZlRUVFyMvLM287UVFRyMvLs7i0bM+ePfDy8kJ4eLi5zLXrMJWJiIiAh4eHrarn8uzVfo3R6XTc/m6RNdqvOaKiopCZmWlxX9iePXsQEhKCXr163XI92it7tV9juP3dOmu13/HjxzF69GjMnj0bK1eubPA9Vjv+tehWahd29epV0el0otPpBIAkJiaKTqeTs2fPiojIli1bJCMjQ06dOiWfffaZhIaGyoMPPmixjvfff1+ys7Pl5MmTsmnTJvH395dFixZZlJk4caIMHDhQsrOzJTs7W+6++26ZMmWK3erpquzRfllZWeb1nj59WjZv3iwhISESExNj17q6oieffFK0Wq3s27dPioqKzFNlZaW5zLx586Rnz56yd+9eOXz4sIwZM0YGDRokRqNRRNRIWWFhYTJ27Fg5fPiw7N27V3r27CkJCQnmdZw+fVp8fHxk4cKFkp+fL++99554eHjItm3b7F5nV2Kv9nv99ddlx44d8v3330teXp4sWbJEAEhKSord6+xKrNF+IiInTpwQnU4nc+fOld///vfmfXJVVZWIqNFhgoKCJC4uTo4dOybbt28XPz8/WbNmjd3r7Ers1X4bN26UTz75RPLz86WgoEBee+018fDwkMTERLvX2ZVYo/3y8vIkICBAZs2aZbGOkpIS8zqsdfxj4vAb0zB/10+zZ88WEZE33nhDevbsKR4eHnL77bfLCy+8YN6YTBYvXixBQUHi4eEhffv2lX/+858NhpgrLS2VWbNmia+vr/j6+sqsWbM4tKAV2KP9cnJyZNiwYaLVasXb21v69esny5cvl4qKCntW1SU11nYA5IMPPjCX+eWXXyQhIUH8/f2lY8eOMmXKFDl37pzFes6ePSuTJ0+Wjh07ir+/vyQkJFgM/Sgism/fPhkyZIh4enpKr169ZP369faookuzV/utXr1a+vTpI97e3tK1a1cZMWKEfPnll/aqpsuyVvuNGjWq0fWcOXPGXObo0aMycuRI8fLykuDgYHn55Zc5FOstslf7bdy4UQYMGCA+Pj7i6+sr4eHhsmnTJjvW1DVZo/2WL1/e6DpCQ0Mtvssaxz/Nb0ETERERERHdEO9xICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJjFxICIiIiKiJv0/dXZmPCdjL5MAAAAASUVORK5CYII=" }, "execution_count": 19, "metadata": { }, "output_type": "execute_result" } ], "source": [ "import pandas as pd\n", "\n", "df = pd.Series({'Year': [1994,1995,1996,1997,1998,1999,2000,2001,2002,2003], \"bbl/day*10^6\": [67.052,68.008,69.803,72.024,73.400,72063,74.669,74.487,74.065,76.777]})\n", "\n", "A = [1994,1995,1996,1997,1998,1999,2000,2001,2002,2003]\n", "b = [67.052,68.008,69.803,72.024,73.400,72063,74.669,74.487,74.065,76.777]\n", "show(df)\n", "\n", "import numpy as np\n", "\n", "c1, c2, c3, c4, c5, c6, c7, c8, c9, c10 = np.polyfit(A,b,9)\n", "t = var('t')\n", "\n", "f(t) = c10*t^9+c9*t^8+c8*t^7+c7*t^6+c6*t^5+c5*t^4+c4*t^3+c3*t^2+c2*t+c1\n", "print \"D-9 Esitmation for 2010: \", f(2010)\n", "plot(c1*t^9+c2*t^8+c3*t^7+c4*t^6+c5*t^5+c6*t^4+c7*t^3+c8*t^2+c9*t+c10, (t,1990,2020))\n" ] }, { "cell_type": "markdown", "metadata": { "collapsed": false }, "source": [ "It appears that the D-9 polynomial is a bad fit because of the Runge phenomenon. A linear regression would be more appropriate." ] }, { "cell_type": "code", "execution_count": 0, "metadata": { "collapsed": false }, "outputs": [ ], "source": [ ] } ], "metadata": { "kernelspec": { "display_name": "SageMath 8.2", "name": "sage-8.2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.14" } }, "nbformat": 4, "nbformat_minor": 0 }