{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Venus\n", "\n", "OK. Venus is at its maximum angular separation from the Sun when Earth-Venus-Sun [form a right triangle](https://en.wikipedia.org/wiki/Elongation)—which is when the Earth-Venus line-of-sight is tangent to the orbit of Venus... \n", "\n", " " ] }, { "cell_type": "code", "execution_count": 32, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import matplotlib.pyplot as plt\n", "plt.style.use('seaborn-whitegrid')\n", "\n", "import numpy as np\n", "import pandas as pd" ] }, { "cell_type": "code", "execution_count": 33, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/anaconda3/lib/python3.6/site-packages/ipykernel_launcher.py:21: FutureWarning: from_items is deprecated. Please use DataFrame.from_dict(dict(items), ...) instead. DataFrame.from_dict(OrderedDict(items)) may be used to preserve the key order.\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 33, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD0CAYAAACLpN0/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvDW2N/gAAIABJREFUeJzs3XdcVFf+//EXvYPSi6JIERBRwQIau0Zj710TU81uys/9riUxJiZ2syabuDEmuhpjjBWMGo1rTEw0FlQQaYM0Ox2kDcwwzMzvj7thY6zUGeA8H488ogxz7+eCvr2ce87nGGi1Wi2CIAhCk2eo6wIEQRCE+iECXRAEoZkQgS4IgtBMiEAXBEFoJkSgC4IgNBMi0AVBEJoJY12dODo6WlenFgRBaNJCQ0Mf+HGdBTo8vKgnIZPJCAgIqMdq9J+45uavpV0viGuuqUfdDIshF0EQhGZCBLogCEIzIQJdEAShmRCBLgiC0EyIQBcEQWgmRKALgiA0EyLQBUEQmgmdzkPXR0qlkkOHDpGdnY2joyPTp09/7HuKioo4ffo0o0ePboQKBUFoqpIyS/j3/kO4mSkbZO69uEP/k7y8PPbt21ej91y9epWff/65gSoSBKGpU1apWX/8KuP/9St/L3yPyaqDDXKeOt2hX7lyhX/84x/s2LHjno9v27aN/fv3Y29vD8D7779Phw4danTsiOjb7L1066Gvl5eXY3mqqEbHnNK9LRND2zzyczZt2kRaWhpxcXE89dRTHDt2jKKiIt58800GDRrEDz/8wFdffYWhoSGhoaH8/e9/Z9OmTSQnJ7Nnzx66devGmjVr0Gg0lJSU8M477xASElKjOgVBaD6ib9xlUUQcabllvOdzDbfb+dzu+EyDnKvWgb5582YOHTqEhYXFfa8lJiaydu1agoKC6lScLsybN4+UlBT69u1LdnY2K1euJCoqii1bthASEsKGDRuIiIjAwsKCBQsWcObMGebNm8fu3buZOnUqR48eZdGiRXTs2JHDhw8TGRkpAl0QWiC5sop/HL/KV2ev425nwVdzezDg3Odg50mp+1MNcs5aB7qnpycbNmxg4cKF972WmJjIl19+SV5eHgMGDOCVV16p8fEnhrZ55N10Y/R/6NSpEwCOjo4oFApu3rxJYWEhL7/8MgByuZxbt27h5eVV/R5nZ2c2btyIubk5crkca2vrBq1REAT9czo1j7ci47l9t4I54e1YONwf67vJcP00DP0ADBvm8WWtjzps2DBu3779wNdGjhzJjBkzsLa25rXXXuPkyZMMHDjwvs+TyWS1PT0KhaJO73+YvLw8ysvLycvLo6qqCplMxu3bt5HL5ZSXl9O6dWsWLlyIsbExP/30EzY2Nty6dYuioiJkMhlLly5l/vz5tG3bll27dpGbm1tvdTbUNeuzlnbNLe16oXldc6lSzZZLhRxPK8XD1oQPh7sR5GLErYxUXC+uw87IjFSrXg12zfX+z4RWq+XZZ5/FxsYGgP79+5OUlPTAQK/LHXZD3aF36NABIyMjbGxscHV1JSAgAFNTU6ysrAgLC2PevHmsXLkStVqNh4cHL7zwAiUlJfz73/8mKiqKKVOmsH79ehwcHHB1dUWprL+n2aIrXfPX0q4Xms81H0vIZumRBArllbw6wJs3B/tibmIkvVheCBHHocs0OnYNa7Bui/Ue6GVlZYwaNYqjR49iaWlJVFQUEydOrO/TNBgzMzMOHrz3CbS3t3f1g9+xY8cyduzYe163sLDghx9+qP793LlzG75QQRD0Ql6pkmWHEjkSn0Wgmy3bnutBkIfdvZ8U8zVUVUCvmg8/10S9Bfrhw4cpLy9n6tSpzJ8/nzlz5mBqakp4eDj9+/evr9MIgiDoBa1WS2TMHT74PokKlZoFwzrycr8OmBj9aTa4ugouboH2fcGlU4PWVKdAb9OmDXv37gW4Z1HNuHHjGDduXN0qEwRB0FO375az5EACv6bkEdquNWsnBuPj/JAJEFePQvEtGL6mwesSK0UFQRCekEaj5ZuoG6z9IRkt8P6YTswOa4ehocHD33R+I7TyhAaae/5HItAFQRCeQHpeGYsj4rh4/S59fR1ZNb4zbe0tH/2mOzFw8xwMWwWGRg1eowh0QRCER1CpNWw+ncE/T6RiYWLEPyZ3YWKIBwYGj7gr/935jWBqA91mN3yhiEAXBEF4qIQ7xSyKiCMxs4Rnglx5f2wnnG3Mn+zNJZmQeAB6vgzmtg1b6H+J5lx/olQq2bdvHxs2bGDXrl1P9J6ioiIOHz5cp/Pm5eWxbNmyOh1DEIT6oVCp+fA/yYz97Aw5JUo+nxnC57NCnzzMAS58CVpNg09V/CMR6H+iq26LTk5OItAFQQ9cul7IiE9P89nJdCZ08+Cnv/Xnmc5uNTtIpRwubQP/kdC6fYPU+SD6O+QSuwsuf/PQlz3L5XDeqmbH7DYLuj66v3lDdVscPHgwXbp04ebNm/j6+rJy5Uo+++wzLl++THl5OStXruStt95i7969jB49mu7du5OSkoKXlxcODg5cunQJlUrFzp07USgULFmyhLt37wLwzjvv0LFjx5p9LQRBuEeZsooPjyXz9fkbuNtZ8PXzPenn51S7g13ZBYoiCPtr/Rb5GOIO/U/mzZuHj48Pf/3rX3FxcWH79u28/fbb7Nq1i6KiIjZs2MBXX33Frl27yMnJqe62GBYWxtSpU0lLS2PRokV89dVXzJ07l8jISABycnJ488032b9/P+Xl5Zw4cQKQWg3s3r0bMzOz6hrkcjmjRo1i586dXLp0iZCQEHbu3ElVVRVpaWls2rSJsLAwduzYwfLly8WdvSDU0a8peQz7+BRfn7/Bs+HtOT6/X+3DXKOB85+DezfwDKvfQh9Df+/Qu05/5N30zSbWbdHNzY127doB0K1bN65duwZwz3sfdG5bW1u8vb0BsLKyQqlUkpKSwvnz56vbDZSUlNT3pQtCi1BUXsny72VExNzG28mKfa+E0729fd0OmvYjFKTBhC3wJDNh6pH+BrqOGBoaotFoAO6bltSmTRvc3NzYunUrJiYmREZGEhAQQFlZWfV7Vq5cyT/+8Q+8vb359NNPuXPnDiDdoefl5eHk5ERMTAxjx44lKSkJQ8MH/5D0qClRHTp0YMyYMYwePZqCgoIaj/kLggA/xGex9GAid8sreW2gD68N8vlfM626OPcZ2LhDp8ZfLS8C/U8cHBxQqVQoFIr7XrO3t+e5555j9uzZ1d0Wn3nmGUpKSkhJSeGrr75izJgx/OUvf6nutvj7OLepqSnLly8nKyuLLl26MGjQIJKSkmpV47x581iyZAl79+6lrKyM1157rU7XLAgtSW6JgncPJnIsMZsgD1u2P9+DTu52j3/jk8iMhWu/wuD3wMikfo5ZAwZarVbb6GdFagEZGhpa6/c3tZabffr04cyZM3U6RlO75vrQ0q65pV0vNN41a7Va9kffZvn3SSiqNMwf4sdLfb0w/nMzrbrY9xyk/QTzE8D84f9I1LV97sOyU9yhC4LQ7N0qLOftA/GcTs2nZ3t71kzsTAenet5NrCAdkg5C79cfGeYNSQR6I6nr3bkgCDWn1mj5+tx1PvzPVQyA5WM7MbPXY5pp1dbZDdLWcmF/qf9jPyER6IIgNEtpuaUsiogn+sZd+vs5sWpCZzxa3b+pfb0ozYHYb6HrDLBxbZhzPAER6IIgNCsqtYYvfk3n05/SsDQz4qMpXRjf7QmbadVW1OegroTebzTcOZ6ACHRBEJqNhDvFLNgfhyyrhJHBbiwb3QknG7PHv7EuFMVw8d8QOBYcvBv2XI8hAl0QhCZPoVLzzxOpbD6dgYOVKV/MDmVYp0Ya+ri0DZQl8NT/a5zzPYII9D9RKpUcOnSI7OxsHB0dmT790b1fQOq2ePr06Xu24XuU+fPnM23aNJRKJVlZWUydOpWPP/6Y06dPs2jRIiIiIrh+/TqrV6+uXiUqCMKDRWUUsDgynmv5cqZ2b8vbIwOws2ikOeCqCqnnuVd/aam/jolA/5Pfuy327dv3id/ze7fFJw303/Xr16/610ePHuXAgQNYW1szf/58zp49W6NjCUJLU6pQse7YVXacv0Fbewt2vtiLPj6OjVtE9HYoy4GJ/27c8z6E3gb6ofRDHEg98NDXy8vLsbzxmO2f/mS873jGeI955Oc0VLfFnTt3sm/fPpycnCgoKAAgMjKSjIwMzM3Nyc7O5pVXXsHT05OSkhJeffVVPv/88xpdnyC0FCev5rIkMp6sEgXP9/Hi78P8sDRt5DhTKeC3j6HdU+D15DeADUlvA11X5s2bR0pKCn379iU7O5uVK1cSFRXFli1bCAkJYcOGDURERGBhYcGCBQuquy3u3r2bqVOncvToURYtWkTHjh05fPgwkZGR+Pr68vXXX3P48GEMDAyYMGHCPed87bXXiIyMZOvWrZiZmXHq1CkR5oLwAHfllSz/PonIy3fwdbYm4tXehHi21k0xMduhLBsmbtHN+R9AbwN9jPeYR95NN8Zy4frqtpiRkYGPjw+mpqYABAcHN2jdgtDcaLVajsRn8d7BRIorVLwx2Je/DvTGzLjhN15+oOq78z56c3cOdeyHfuXKFWbPvn/z059//pmJEycydepU9u7dW5dTNLon7ba4Y8cOZs2aRZcuXe55z8qVK3njjTdYu3Ytfn5+aLVa2rZtS1paGgqFArVajUwma/TrEoSmKqdEwcs7onnt28t4tLbg8OtP8behfroLc4CYr6E0CwYs1l0ND1DrO/TNmzdz6NAhLCzuXXmlUqlYvXo1+/fvx8LCgunTpzNw4ECcnGrZLL6RNUS3RXt7e958802mTZuGvb39fV8zQRDup9Vq2XvpFiuOyKis0vD2CH+e71PPzbRqQ6WA3z4Cz97QXn/uzqEOge7p6cmGDRtYuHDhPR9PT0/H09MTOzupOU1oaCiXLl3imWeeqVuljcTMzIyDBw/e8zFvb2927NgBwNixYxk7duw9r1tYWFRvNgEwd+7c+447YsQIRowY8dDz/nFPUtH3RWjpbhaUszgyjrPpBfTysmftxGDaO9Zwy8mGcmmrdHc+/otG38DicWod6MOGDeP27dv3fbysrAwbG5vq31tZWVFWVvbAY9Rl6EGhULS4oQtxzc1fS7teuPea1Roth5JL2H65EEMDeD3MkeF+NlTk3USWp+NCAUOVHO9f1qJw6cEtpTPU8nvVUN/nen8oam1tjVwur/69XC6/J+D/qC4PNUXf6JahpV1zS7te+N81p+SU8tb+OGJvFTHI35mV44Nws9Oz4cmfV4KyCOsx6wjw0E1+RUdHP/S1eg90b29vbty4QVFREZaWlly6dIkXXnihvk8jCEIzoVJr+eREKv86mYq1mTGfTOvKmC7uDdtMqzbKcqXt5QLHgUeIrqt5oHoL9MOHD1NeXs7UqVNZvHgxL7zwAlqtlokTJ+Li4lJfpxEEoRm5cquIN4/c4frdSsZ0cee90YE4WDdwM63aOvUhVClg0FJdV/JQdQr0Nm3aVE9L/OOy90GDBjFo0KC6VSYIQrNVUanm4xMpbDmdQWsLI7bM6c6QQD2+8Su8JjXhCpkDjj66ruah9HZhkSAIzdO59ALeiozjekE503t6MtHbkO76HOYAP30g7UbUf5GuK3kkEeiCIDSKEoWKNT8k823UTdo5WPLtS73o7e2o/7N6bpyFxEjovxhs3XRdzSOJQBcEocH9nJzD25EJ5JYqeKmvF38b2hELUx2u9HxSGg0cWwy2HtDnTV1X81gi0AVBaDAFZUo++D6Jg7GZdHSxYdPsULq2baXrsp5c7E7IuiK1xzWtWXdXXRCBLghCvdNqtRy6ksn7h5MoVaiYP8SPVwd4Y2qs42X7NaEogZ/eh7ZhEDRR19U8ERHogiDUq6ziCt45kMBPybl0aduKdROD6ej64MWFeu3UhyDPh5n79G6J/8OIQBcEoUY0Wg2V6kqUaiUqjQqlWolSrUSrgaMJOXz5SwYqDbzxtDfTe3phY2qARqvB0KAJ3Z3nJsP5z6HrTL3YWu5JiUAXhBZMpVZRoCggrzyPvIo88ivyya/I567iLiWVJZRUllBaWSr9Win9ulJT+eiDeoIJsO2W9N/vLIwtsDKxwtrEGgcLB5wsnHC0cERTquGa2TXa2bWjnU07rE2tG/SaH0ujge//H5hZw9D3dVtLDYlAF4RmSqvVUlJZwp2yO2SWZXKn7E71rzPlmeSX53NXefe+9xlggI2pDbamttia2WJraouzpbP0e1NbzI3NMTUyxczIDBNDU6LSizmWkI+xgSFju7nT06sVWrRotBpUGhUVVRXIVfLq/0orS8mvyCepIIm8ijwqqir49va31ee3N7fHy86LAPsAAh0C6eTYifa27RvvDj/2G7h5DsZ+BlaNvEdpHYlAF4QmrkhRxPWS61wrvsb1kuvcKLnBrdJbZJZlUqa6t9OptYk1HtYeeFh50NWpq3SXbOmIk4VT9R2zvYU9JoYmjz1vcnYJi/bHceW2NUMCQlgxLghXO/Ma1x+TEIONhw03Sm5U/5denM6+lH0o1UoAbExsCHUJpadbT3q69sS3tW/DBHxZHhxfKu1E1HVm/R+/gYlAF4QmQKPVcLv0NmlFadXBfb34OtdLrlOkLKr+PGNDYzxtPPG08aS7S3fcrd1pY90Gd2t33K3dsTW1rXPTK2WVms9OprPxZBp2FiZsmN6NUcFutT6uhZEFvq198W3te8/HqzRVZBRnkFSQRGxuLBezL/LL7V8AcLZwZpDnIIa0G0KoSyjGhvUUZceXQKUcRv2zyTwI/SMR6IKgZ/Ir8km9myr9VyT9P6M4g4qqiurPcbRwpL1te4a0G0J72/Z42XnR3rY97tbu9RduD3D55l0WRcSRklPG+G4eLB0ViL2VaYOcy9jQGL/Wfvi19mOczzgAsuXZnM86z6+3fuW7tO/YfXU3rcxaMcJrBBN8J9DRvmPtT5h2AuL2SMv7nfzq6Soalwh0QdARlVpFenE6sgIZKXdTSL2biixfRsmFkurPsTe3x7e1LxN9J+LX2g+fVj60t2uPjWnjTgMsr6xi/fEUtp65hqutOVuf684g/8bvv+Jq5co4n3GM8xlHRVUFZ++c5dj1Y+xL2ce3yd/SyaETUztOZVSHUZgYPX7YqFpFERx8HRw7wlN/a7gLaGAi0AWhESjVSlIKU5AVykgqSEJWKCP1bioqjQqQZoD4tPIhtFUoPbx6SEMQrXxxsHDQceVwNi2fxZHx3CwsZ1aYJ4uG+2NjXoOwbCAWxhYMbjeYwe0GU6Qo4si1I+xP2c+7Z9/lX5f/xezA2Uzym/Rks2aOvQVlOTBtJ5jU/DmAvhCBLgj1rFxVztW7V6XgLpAhK5SRXpSOWqsGwM7MjgD7AGYFziLQPhB/e388bT0xNDDUqx2LiitUrD4qY/fFW7R3sGT3y2GEddD9PzAP0sq8FTMDZjLDfwbnMs+xNWEr66PXszl+My91fonpAdMxM3pIn/Xko3DlW+i3UG83rnhSItAFoQ6qNFWkF6UTlx9HfF488fnxpBelo0ULSEMmgQ6BDGg7gED7QAIcAnCzqv0DxMbyY1IO73wXT16pklf6d2D+ED/MTfS/mZaBgQG9PXrT26M3CfkJ/Cv2X6yPXs+3yd/yerfXGdVh1L1fe3kBHH4TXDtDvwW6K7yeiEAXhCek1WrJlmcTny8Fd1xeHLJCWfXDylZmrQhyDGJou6F0cuhEgEMAThZOeh/ef5RfpmTZoUS+j8vC39WGzXO6E9ymCTXT+oMgxyA2DdnE+azzfHTpI97+7W2+S/uOpWFLaW/XHrRaaQFRxV2Y8x0YN8zD3cYkAl0QHqKssoyEgoTqO+/4/HjyK/IBMDE0IcA+gAm+E+js2Jlgx2Da2LRpUuH9R1qtlu9i7/D+4STKlWr+b6gf8wZ4Y2LUhJbrP0SYWxi7R+0mIjWCjy99zMRDE3mlyys8rzTEWHYIhn4ALp10XWa9EIEuCEiBdrv0NpfzLhObG8vl3Mv3DJ20t21PuFs4QY5BBDsF07F1x5rNotBjmUUVLDkQz8mreXTzlJpp+bo0wWZaj2BoYMhkv8kMaDOANRfWsOHyBk4pVaz17o9H+Ou6Lq/eiEAXWqRKdSWyQll1eMfmxlKgKACk1ZRdnLrwdPunCXYMJsgxCDszOx1XXP80Gi07L9xk7Q/JqDVa3h0VyLO922Nk2DR/yngSTpZOrO/9AUeTT7LcTMskw1zeu3mc4e2H67q0eiECXWgR7iruEpsbS2xeLLG5sSTkJ1Q3mWpj3Ybe7r3p6tyVrs5d8bbzxshQ/x8A1sW1fDmLIuK4cK2Qp3wcWT2hM23t9X8Dh3rxwyJGZGcQPGULi69/x4JfFyArkPFGtzea/PddBLrQ7Gi1Wq6XXOdy7uXqu+/rJdcBafVhoH0g0/yn0c25G12cuuBk6aTbghtRlVrDlt+u8fGPKZgaG7JuYjCTuzfdsf8ai/kaLu+Avn+nTeAktnUcy+oLq9masJXUu6ms7be20Rdt1ScR6EKTp9FqSL2byqWcS0TnRBOdE02hohCQ5nx3derKWJ+xdHPuRieHTpgbN92FI3WRlFnCwogrJNwp4elAF5aPC8LFtgV9LW5Hw5H/gw4DYeDbAJgYmfBu+Lv42/uzOmo1zx57lk1DNuFs6azjYmtHBLrQ5Kg0KmQFsurwjsmNobSyFAB3K3f6uPchxCWEEJcQvGy9Ws7d50Moq9T86+c0Pv8lnVaWJmycGcIzQa4t6+tSlgt7ZoGNK0zaCn8aWpnScQptbdry5sk3mfPDHL4Y+gXtbNvpqNjaq3WgazQali1bxtWrVzE1NWXFihW0a/e/L8CKFSuIiYnBysoKgI0bN2Jj03R/lBF0R6lWEp8XXx3gsXmx1XO/29u25+l2TxPqEkqoSyju1u46rla/RN+Qmmml5ZYxIcSDpSMDad1AzbT0lloF+56T5pu/cBws7R/4aeHu4WwdtpW/nPgLc36Yw5dDv6xbsy8dqHWgnzhxgsrKSvbs2UNsbCxr1qzh888/r349MTGRLVu2YG//4C+eIDxMuaqc2LxYLmVf4rfrv5F2KQ2VRoUBBvi29mW8z3hCXUIJcQnB0aJpbUDQWOTKKv5x/Cpfnb2Ou50FX83twYCOTXMYoU60Wji6AG6cgQmbwS34kZ8e5BjE9me28+LxF3np+EtsHbYVn9Y+jVRs3dU60KOjo+nbty8AXbt2JSEhofo1jUbDjRs3ePfdd8nPz2fSpElMmjSp7tUKzVJpZSkxOTHVY+BJBUmotWqMDIzwsvRiZsBMQl1C6ebcrVlOH6xvp1PzeCsyntt3K5gT3o6Fw/2xNmuho6tnN0D0NnhqPgRPeaK3eNl5sXXYVuYem8uLx19k6/CtdLDr0MCF1o9af5fLysqwtv5fFzMjIyOqqqowNjamvLycWbNmMXfuXNRqNXPmzCEoKAh/f/97jiGTyWpduEKhqNP7m6Lmcs0KtYLk0mQSShNILEkkQ56BFi0mBib4WPsw1m0sgTaB+Fn7gQrMzc2hDGnrNDJ1XX6Dqsv3uFSpZsulQo6nleJha8KHw90IcjHiVkZqPVdZvxrqz7XNrZ9pc3YpJW0Hc8d9MtTwHG/7vM0y2TKeO/IcywOX42xWfz/hNNQ11zrQra2tkcvl1b/XaDQYG0uHs7CwYM6cOVhYWAAQFhZGcnLyfYFel65y+tSVrrE01WuuqKqo3nHmQvYFEvMTqdJWYWxoTLBjMC97v0xP1550ce5yX0e8pnrNtVXb6z2WkM3SIwkUyiv5ywBv3hjs2ySaaUEDfY9vX4ILH0CbntjO2YmtiUWNDxFAAO282vHssWf56PpHfP3M1/X2E2Jdrjk6Ovqhr9U60ENCQjh58iQjRowgNjYWP7//7fBx/fp15s+fz4EDB9BoNMTExDB+/PjankpoYpRqJXF5cVzIvsCFrAvE5cdRpanCyMCITo6deC7oOXq49qCbczcsjGv+F034n9xSBcsOJXI0PptAN1u2PdeDII8WPixVkA67pkkzWqbvglqE+e98W/vyycBPeOXHV3jj5zf48ukvH96GVw/UOtCHDh3KmTNnmDZtGlqtllWrVrFt2zY8PT0ZPHgwo0ePZsqUKZiYmDB27Fh8fX0ff1ChSVKpVcTlSwF+MfsiV3KvUKmpxNDAkAD7AGYHzKaHaw9CXEKwMrHSdbnNglarJTLmDh98n0SFSs2CYR15uV+HZtFMq06K78DXY0GrgZn7waruD817uPZg5VMrWXhqIe/89g7r+q3T2ymftQ50Q0NDPvjgg3s+5u3tXf3rl156iZdeeqn2lQl6S6VRkZifWD2EEpsbi0KtwAADOtp3ZKr/VHq69iTUJbRJr7rTV7fvlvP2gQROpeQR2q41aycG4+P8BLvyNHfyfNgxDhTF8OxhcKy/m8hnvJ7hTtkdPon5BH97f17o/EK9Hbs+tdBH30JNqDVqku8mE5UVxYXsC1zOuUx5VTkAPq18mOA7gZ6uPenu2l3MQmlAGo2Wb6JusPaHZLTA+2M6MTusHYbNuJnWE1OUwDcToegmzIoE9671fooXgl4guTC5OtT7ePSp93PUlQh04T5arZYbJTeIyorifNZ5LmRfoKRS2rjYy86L0d6j6eHagx6uPbA3F+sMGkN6XhmLI+K4eP0ufX0dWTW+BTXTepxKuTRmnpMA03ZB+4YJWgMDAz7o/QEZxRksPLWQ3aN209ambYOcq7ZEoAsA5JXncT7rPFFZUURlR5EtzwakXdYHeQ6il1svern2alGNrPSBSq3hy1MZfPJTKhYmRvxjchcmhnjo7Rhuo1OWwbdT4OY5aeGQ39MNejpLE0s+GfAJ045MY/7J+ewcuVOvHpKKQG+hSitLuZR9iajsKM5nnie9OB2Qmln1dO3JS51fopdbLzxtPEV46EjCnWIWRcSRmFnCiM6uLBvTCWebFtRM63GUpbBzMty6IIV558ZZvNjWti2r+67mrz/9lY8ufcRbvd5qlPM+CRHoLUSlupIreVc4l3mOqOwoEvMTUWvVmBuZE+ISwhifMYS5heFv74+hQQufKaFjCpWaT39K5YtTGbS2NGWZ4uELAAAgAElEQVTTrBCGB7npuiz9oiiBnZOk+eYTt0DQhEY9fb82/ZgVMItvZN8Q7h7OgLYDGvX8DyMCvZn644PMqKwoYnJiUKgV1XPBX+j8AmFuYXRx6oKpUQtr1qTHEnMV/PXoaTLy5EwObcM7IwOxs2weW93Vm/JCaZgl87LUObHTOJ2UMT90PtE50Sw9s5T9o/fjYuWikzr+SAR6M6HVarlZevOeB5nFymJAmoky0W8ivVx70d21u5hKqIfKlFV8eCyZr89l4t7Kgq+f70k/P/G84j4lmbBjAhSmw+SvIGC0zkoxNTJlXb91TPl+Cm//9jabn96s859uRaA3YfkV+dUPMs9nnb/nQeaANgMIcw8TDzKbgF9T8ng7Mp7M4grGBNiyalo4Vi21mdaj5KfBjvFQUSgtGurQX9cV0d6uPYt7Lua9s++xK3kXMwNm6rQe8aemCSlXl/PLrV+qAzytKA0AW1Nbern14sWgFwlzDxMPMpuIovJKPvg+iciYO3g7WbF/XjiW5TkizB8k87I0zxwDeO57cO+m64qqjfcZz483fuSTmE/o69EXT1tPndUi/uTosSpNFQn5CZzLOse5zHNcyb2CBg1mRmaEOIcw2ns0vdx64d/av8lvbtvSHI3P4t2DCRSVq3htoA+vDfLB3MQImSxH16Xpn9QTsO9ZsLCH2QfAUb/6kxsYGPBe+HtMODiBpWeWsm34Np0NvYhA1yNarZZbpbc4l3mOs5lnuZh9kVJVKQYYEOgQyBi3MYzuPPqBXQmFpiG3RMG7BxM5lphNkIct25/vSSd3sbr2oS5shh8WgksnmLEXbPVzRypXK1cW9lzI0jNL+Vb2LbMCZ+mkDhHoOlasLCYqK6r6LvxO2R1A2hvz6fZPE+4eTi/XXrQybyW13HRrOa1kmxOtVsu+6Nus+D4JRZWGRcP9eamvF8YtvZnWw6ir4D9vw4UvwO8ZaWqimX73qxnrPbZ66KVfm346GXoRgd7IVGoVsXmxnMs8x/ms8yQWJKLRarA2saaHaw+e6/Qc4e7hYhy8GblVWM7bB+I5nZpPz/b2rJnYmQ5O+h1OOqUshf3PQ+pxCPsrPL38vk2d9dHvQy9jvxvL8vPL+XLol43+d1gEegPTarVkFGdwLvMc57LOcTH7IhVVFRgZGNHZsTOvBL9Cb/feBDkGYWwovh3NiVqj5etz11l37CqGBrB8bCdm9hLNtB6pIB12z4T8FBj1MXR/XtcV1YizpTNvhLzBqqhVHLl2hFEdRjXq+UWCNICCigLOZ52vDvHc8lwA2tm2Y4z3GHq796aHaw8xH7wZS8stZeH+OGJuFjGgoxMrx3fGo5XYzOORrv4AkS+DoTHMigDvgbquqFam+E3hcPphPrz4IX09+jZqB1IR6PVAUaUgJjeG85nnOZd1juTCZEDqi9LLtRfh7uGEu4fjYe2h40qFhqZSa/ji13Q+/SkNSzMjPp7ahXFdRTOtR9Jq4OQq+HUtuHWBqd9AK91N/asrI0Mj3g1/l2nfT+Pj6I9Z1ntZo51bBHotaLQaUu+mcjbzLOcyzxGTG4NSrcTY0Jhuzt14o9sb9Hbvjb+9mE7YksTfLmbB/iskZ5cyMtiN98d0wtFazEZ6pIq7tD39d8g6C11nwsj1ddoyTl/42/szO3A2XyV+xRjvMYS4hDTKeUWgP6EceU71TJTzWecpVBQC0rL6yX6T6e3em1CXUCxNRI/qlkahUvPPE6lsPp2Bg5UpX8wOZVgnV12Xpf9uXYT9z2NVmgkjP5LGy5vRTzKvdnmV/1z/Dx+c+4B9o/dhYtTwPXlEoD9EuaqcSzmXpHHwzHPV7WUdzB2kIRS3cMLcwvSiIY+gO1EZBSyOjOdavpyp3dvy9sgA7CxEM61H0mjg7Cfw8wqwdef6oE149Zis66rqnaWJJUt6LeG1n1/jG9k3zA2a2+DnFIH+X2qNGlmhrHoYJTYvlipNFWZGZoS6hDLOZxzh7uH4tfYT46ECpQoVa48l8835m7S1t2Dni73o41P3DYmbvbJc6cFnxkkIHAejP0FxPUvXVTWY/m37M6DNADZd2cQIrxENfgPYogP9Ttmd6gCPyoqq3mYtwD6A2YGzCXcLJ8QlRKzKFO5xMjmXJQfiySpR8MJTXvzf035Ymrbov0pPJu0nODAPlCUw+hMIefa/QyzNN9ABFvZcyLjvxrE+ej3r+q1r0HO1qD+FpZWlXMi+UD2McrP0JgAuli4M8hxEuFs4vdx64WDhoONKBX1UKK9k+fdJHLh8B19nayJe7U2IZ2tdl6X/KuXw43twcTM4+cOcg+ASqOuqGk1bm7Y83/l5Nl3ZxGS/yfRw7dFg52rWga7SqKTmVv/tjZKQn4Baq8bC2IKerj2ZETCDcPdwvGy9xDCK8FBarZYj8Vm8dzCR4goVbwz25a8DvTEzFjOYHuvWBTjwChRmQNhfYPC7zWIWS009H/Q8h9IOsSpqFXtH722w89Q60DUaDcuWLePq1auYmpqyYsUK2rVrV/363r172b17N8bGxrz66qsMHNjwiwR+363+99koF7IvIFfJMTQwJMghiBc6v0C4WzhdnLo0yhNnoenLKVHwzncJ/JiUQ3AbO755sRcBbra6Lkv/VSnhlzVw5p9g6wHPHgavfrquSmcsjC1Y2HMh/+/k/2NP8h5CDUIb5Dy1DvQTJ05QWVnJnj17iI2NZc2aNXz++ecA5OXlsWPHDiIiIlAqlcyYMYM+ffpgalr/W50VKYo4n32e85nnOZt5liy5NB7nYe3BCK8RhLuH09O1Z6Ou1hKaPq1Wy56Lt1h5VEZllYa3R/jzfB/RTOuJZMdLY+U5CdBtFgxbDebiH8FBbQfRx70Pn8V+xkedPmqQc9Q60KOjo+nbty8AXbt2JSEhofq1uLg4unXrhqmpKaampnh6epKcnExwcHDdKwaSCpLYdWsXqempJBUkoUWLjYmNtMlD5xcJdwunrW3bejmX0PLcLChncWQcZ9ML6OVlz9qJwbR3tNJ1WfpPVQG/roOzn4JFa5i+Gzo+o+uq9IaBgQGLey5m/KHxRGRGEN4lvN7PUetALysrw9r6fx3jjIyMqKqqwtjYmLKyMmxs/tenxMrKirKysvuOIZPJanxejVbD3Ji5VKor8bPxY7LHZILtgvG28sbIwAg0UHanDNmdmh9b3ykUilp9zZqyxrxmtUbLIVkx2y/fxdAQXg9zZLifDRV5N5HlNUoJTfZ7bJkbg+vF1ZiV3aKo/Uhyu76BWmMHT3AtTfWaa+vV9q+iVCkb5JprHejW1tbI5fLq32s0GoyNjR/4mlwuvyfgfxcQULve3ofbHebOtTuEBDXOclp9IZPJav01a6oa65pTckp5a38csbeKGOTvzMrxQbjZNf7Duyb3Pa4ogh/fhZjt0KodzP6OVt4DaVWDQzS5a66jAALqdM3R0dEPfa3WA4IhISGcOnUKgNjYWPz8/KpfCw4OJjo6GqVSSWlpKenp6fe8XlcuVi5YGLW8J+VC/aus0vDJiVRGfnqam4XlfDKtK/9+trtOwrxJ0Woh6SB81hMu74Deb8BfzjfZDonNRa3v0IcOHcqZM2eYNm0aWq2WVatWsW3bNjw9PRk8eDCzZ89mxowZaLVa5s+fj5mZWJwj6Jcrt4pYFBFHcnYpY7q4897oQBxEM63HK0iHowsg/Sdw7Qwz9ujVps0tWa0D3dDQkA8++OCej3l7e1f/esqUKUyZMqX2lQlCA6moVPPxiRS2nM7A2cacLXO6MyRQ9OR5rEo5nF4PZzeAkZk0e6XnSyCmAOuNZr2wSBD+7Fx6AYsj47hRUM70np68NcIfW3MRSI+k1YLssLTHZ/EtCJ4KQz8AG9FRUt+IQBdahBKFitVHk9l14SbtHCz59qVe9PYWzbQe64/DK86d4Lmj0L6PrqsSHkIEutDs/STLYcmBBHJLFbzcrwPzh/hhYSqW7T9SRRGc+hCivpCW6g9fAz1eAiMRGfpMfHeEZqugTMn7h5M4dCWTji42bJodSte2NZlQ1wKpVXBpG/yyGiruSis9By0FG/GMoSkQgS40O1qtlkNXMnn/cBKlChXzh/jx6gBvTI3Fsv2H0moh9TgcfwfyU6B9Xxi2CtzqZ3W30DhEoAvNSlZxBe8cSOCn5Fy6tG3FuonBdHS9f1Gb8AfZCXB8CWT8Ag4+0pJ9v+HNaju4lkIEutAsaDRadl28yeqjyVRpNLwzMoC5fbwwMhSh9FDFt6WhldhvwcwWhq+V9vU0rv8mekLjEIEuNHnX8+UsjozjfEYhvb0dWD2hM+0cRDOthyovlOaTX9gMaKHXPOi3ACztdV2ZUEci0IUmq0qtYeuZa6w/noKpkSFrJnRmao+2YrOSh1GWwfnPpW6IlWXQZToMWAytPHVdmVBPRKALTVJydgmL9sdx5XYxQwJcWDEuCFc7c12XpZ+qKqXmWb+uA3kudBwJg5eCc8tpiNVSiEAXmhRllZrPTqaz8WQadhYmbJjejVHBbuKu/EE0akiIgJMr4e51aNcHpu2Etj11XZnQQESgC03G5Zt3WRQRR0pOGeO7ebB0VCD2VuIB3n00Gkg6AL+shfyr4NIZZu4HnyFi5kozJwJd0HvllVWsP57C1jPXcLU1Z+tz3RnkLxa63EejAdlBKcjzZODkD5O2QeA4MBRz8FsCEeiCXrucVcHLh09xq7CCWWGeLBruj41opnUvjQaSD0tBnpsIjn4w8d/QaTwYihYHLYkIdEEvFVeoWH1Uxu6LWXg5WrHn5TB6dXDQdVn6RauF5O+lIM+JlxYFTdgMQRNFkLdQItAFvXM8MZt3vksgv0zJpCA7VkwNx9xEBFQ1jQZSfpAWBWXHg30HGP8FBE0SzbNaOPHdF/RGfpmSZYcS+T4uC39XG7Y82x2T0iwR5r9TV0HSd9KioNwkaO0F4z6HzlNEkAuACHRBD2i1Wr6LvcP7h5MoV6r5v6F+zBvgjYmRITJZlq7L072qSriyC377GO5eA8eO4o5ceCDxp0HQqTtFFSw5EM8vV/Po5ik10/J1Ec20AKgslxYEnd0AJXfArStM2QH+o8SsFeGBRKALOqHRaNl54SZrjsrQaOG90YHMCW8vmmkBKIqlPivnN0J5gbQgaMyn4D1YzCMXHkkEutDoMvLKWBwRz4XrhTzl48jqCZ1pa2+p67J0T56PU/wm+O4AKIvBZyj0/T9oF67ryoQmQgS60Giq1Bq2/HaNj39MwczYkHWTgpkc2kYs2y/MgHOfweWdOFQpIHAMPPU3cO+q68qEJkYEutAokjJLWBhxhYQ7JQzr5MLysUE427bwZlq3o+HsJyA7DIbGEDyFDNeRePcaoevKhCZKBLrQoBQqNf/6OY1Nv6bTytKEjTNDeCbIteXelWs00lZvZz+FG2fAzA76vCn1JLdxpVIm03WFQhNWq0BXKBQsWLCAgoICrKysWLt2Lfb29zbHnzdvHkVFRZiYmGBmZsaWLVvqpWCh6Yi+UcjC/XGk58mZEOLB0pGBtG6pzbSqlBC3V5qxkn8VbNtIe3aGzAEzMatHqB+1CvRdu3bh5+fH66+/zpEjR9i4cSPvvPPOPZ9z8+ZNjhw50nLvxFowubKKD/9zle3nruNuZ8FXc3swoKOzrsvSjYoiuLQVor6Asmyp8+GEzVKfFSPRk0aoX7UK9OjoaF588UUA+vXrx8aNG+95PT8/n5KSEubNm0dJSQkvv/wyAwcOrHu1gt47nZrHW5Hx3L5bwbPh7Vgw3B9rsxY4snf3hhTiMdul3YE6DITxn0v/Fzc5QgN57N+0ffv2sX379ns+5uDggI2N9GOilZUVpaWl97yuUql4/vnnmTNnDsXFxUyfPp3g4GAcHO5triSrw3ihQqGo0/ubIn2+5lKlms2XCvgxrYw2tiZ8ONydIBcjbmWk1um4+nzN99Fqsci/gn3KHmzu/AoYUOI5lIKOM1G29oVKIDn5kYdoUtdbT8Q115/HBvrkyZOZPHnyPR977bXXkMvlAMjlcmxtbe953dHRkWnTpmFsbIyDgwMBAQFcu3btvkAPCKj9FlgymaxO72+K9PWajyVksfRIIoXySv4ywJs3BvvWW/8Vfb3me1RVSj1Wzm+EzMtg3kp60NnjRezs2mBXg0M1ieutZ+KaayY6Ovqhr9XqZ+GQkBB+/fVXgoODOXXqFKGhofe8fvbsWXbu3MmXX36JXC4nNTWVDh061OZUgh7LLVXw3sFEfkjIJtDNlm3P9SDIoybx1cTJCyB6K1zYIo2PO/rByI+gyzQwtdJ1dUILVKtAnz59OosWLWL69OmYmJiwfv16ANatW8fw4cPp378/v/32G1OmTMHQ0JC//e1v982CEZourVZLRMwdln+fRIVKzYJhHXm5XwdMjFpIf5GcJIj6XJq1UqWQluSP/Qy8B4keK4JO1SrQLSws+PTTT+/7+MKFC6t/vWTJktpXJeit23fLeftAAqdS8ujerjVrJgbj42yt67IankYDaT9KwyoZv4CxBXSZLs0fd/bXdXWCAIiFRcIT0mi07Dh/g7XHpId674/pxOywdhg292ZayjKpdW3UJihIAxt3GPwehD4HluKnTkG/iEAXHis9r4xF++O4dOMu/fycWDU+iDatm3kzrYJ0uLgFLu+UGmW5h0j7dAaOFfPHBb0lAl14KJVaw5enMvjkp1QsTIz4x+QuTAzxaL6LxTRqSP0RLnwJ6T+BoYkU4D1fhrY9xfxxQe+JQBceKOFOMQv3x5GUVcKIzq4sG9MJZ5tm2kyrvBAufyPdkRfdABs3GLgEQp4FGxddVycIT0wEunAPhUrNJz+l8uWpDFpbmrJpVgjDg9x0XVbDyIqT7sbj90mzVdr1gaHvSzsCiWEVoQkSgS5Uu3i9kEX748jIlzM5tA3vjAzEzrKZBVtVJcgOSTsC3ToPJpbSvPEeL4FrkK6rE4Q6EYEuUKasYt2xZL4+d4M2rS3Y8UJP+vo66bqs+lWSBdFfQfQ2KMuB1l5St8OuM8Cita6rE4R6IQK9hfs1JY+3I+PJLK7gud7tWTCsI1bNpZmWVgs3z0vDKrJD0kNP36HSQ07vwWIRkNDsNJO/uUJN3ZVXsvxIEpExd/B2smL/vHBC2zWTedWKEojbA5e2QW4imNtJC4C6Pw8O3rquThAajAj0Fkar1fJDQjbvHkygqFzFawN9eG2QT70109KprCtS7/G4faCSg2swjP4EOk8WvVWEFkEEeguSW6Jg6cEE/pOYQ5CHLduf70kn9ybeTEtVAQmRUpDfuSQtyQ+aCD2elxYDibnjQgsiAr0F0Gq17Iu+zYrvk1BWaVj8jD8vPuWFcVNuppWfKoV47E5QFEudDoevkWasiIecQgslAr2Zu1VYzluR8fyWlk/P9vasmdiZDk5NtJlWVSVcPQIX/w3XT0srOQNGQ48XpDnk4m5caOFEoDdTao2Wr89dZ92xqxgawPJxQczs6dk0m2kV3ZSmHMbsAHku2HnC4Heh22ywbqF7lQrCA4hAb4ZSc0pZFBFHzM0iBnR0YuX4zni0stB1WTWjUUPaCdqc+gSyzkp3377DpLtx70Fg2Awe4gpCPROB3oyo1Bo2/ZLOhp/TsDQz4uOpXRjXtYk10yrNgctfQ/R2KL6FhbkD9Pu71FelVVtdVycIek0EejMRf7uYBfuvkJxdyshgN94f0wlHazNdl/VkNBrIOCkNq1w9Cpoq8OoPw1aSqu1AQKfOuq5QEJoEEehNnEKl5uMTKWw+lYGjtRlfzA5lWCdXXZf1ZEoypX7jl7+WxsktHaQFQKFzwdFH+pwWthu8INSFCPQmLCqjgMWR8VzLlzOtR1veGhGAnYWeN9P6ved4zHZIOQZajXQ3PmSZ1OXQuIn8VCEIekgEehNUqlCx9lgy35y/SVt7C3a+2Is+Po66LuvRim7B5R1S3/GSO2DlDH3ehJA5YN9B19UJQrMgAr2JOZmcy9sH4skuUfDCU17839N+WJrq6bdRrZLuwqO3Q9oJ6WM+g6UFQB2fET3HBaGe6WkSCH9WKK9k3elcTmZk4OtsTcSrvQnx1NMVkYXXIOZraRVnWY60sXK/BRAyG1p56ro6QWi2RKDrOa1Wy/dxWSw7lEhReSVvDPblrwO9MTPWs3nYVZWQ/L00Np7xCxgYSvPGQ58Fn6FgJP6oCUJDq9Pfsh9//JFjx46xfv36+17bu3cvu3fvxtjYmFdffZWBAwfW5VQtUk6JgiUHEjghyyG4jR3LBzkxoo+frsu6V36qFOKx30J5gbSKc+AS6DoT7Dx0XZ0gtCi1DvQVK1bw22+/ERAQcN9reXl57Nixg4iICJRKJTNmzKBPnz6YmprWqdiWQqvVsufiLVYelVFZpWHJiADm9mlPaspVXZcmUSmkDSOiv4IbZ8DQWBoTD30OOgwUqzgFQUdqHeghISEMGTKEPXv23PdaXFwc3bp1w9TUFFNTUzw9PUlOTiY4OLhOxbYENwrkvBUZz9n0Anp52bN2YjDtHfWkl3dOkjQ2fmUXKIqkbdyGLIMuM8DGRdfVCUKL99hA37dvH9u3b7/nY6tWrWLEiBFERUU98D1lZWXY2NhU/97KyoqysrL7Pk9Wh0UjCoWiTu/XN2qNloOyYr6+fBdDQ3g93JHhvjZU5N1Elid9ji6u2VAlx/bmj7TKOIRFYRJaQ2NKPAZQ5D2OcucQaaz8diFQ2CDnb27f58dpadcL4prr02MDffLkyUyePLlGB7W2tkYul1f/Xi6X3xPwv3vQcM2TkslkdXq/PrmaXcriiDiu3CpikL8zK8cH4WZ3fzOtRrtmrRZunJXmjCd9B6pycAqAYaswCJ6KnZUjjbUtRnP6Pj+Jlna9IK65pqKjox/6WoNMPQgODuaf//wnSqWSyspK0tPT8fPTs4d5eqCySsPGX9L47GQaNuYmfDKtK2O6uOuumVZptvRw8/I3UJgOpjYQPAW6zQEPsfuPIOi7eg30bdu24enpyeDBg5k9ezYzZsxAq9Uyf/58zMzEku4/unKriIX747iaU8rYru68OyoQB10001KrIPW41Gs89Tho1eDZW+pwGDhW7MUpCE1InQK9V69e9OrVq/r3c+fOrf71lClTmDJlSl0O3yxVVKr56Mer/Pu3azjbmLNlTneGBOrggWJ+qrQUP3aXtGmEtQv0fl3aNOL3xliCIDQpYrVHIzqbns9bkfHcKChnRi9PFj/jj615Iy5/r5RD4ndSkN88BwZG4DdMCnHfoWIpviA0cSLQG0GJQsXqo8nsunCTdg6WfPtSL3p7N1IzLa0Wbl+SWtQmREJlGTj4/He64XSwaSKtdgVBeCwR6A3sRFIOS76LJ69Uycv9OjB/iB8Wpo2w8EaeD1d2Sw8482RgYgmdxkO3WeAZLh5wCkIzJAK9gRSUKXn/cBKHrmTS0cWGL2Z3p2vbVg17Uo0a0n+WFv9c/QE0KvDoDqM/gU4TwNy2Yc8vCIJOiUCvZ1qtlkNXMll2KJEyZRXzh/jx6gBvTI0NG+6khdekzoax30q9xi0doOfL0t24S2DDnVcQBL0iAr0eZRVX8M6BBH5KzqVr21asmxSMn8v9C6rqhaoCZN9LY+PXTgEGUq/xYaug4wgwFn1zBKGlEYFeDzQaLbsu3mT10WSqNBreGRnA3D5eGBk2wDh11hVpznj8XlAUS/3FBy6BrjPArk39n08QhCZDBHodXcuXszgijqhrhfT2dmDNhGA8HSzr9yTyfIjbi1fUVihKBSMzCBgtbRjRvh8YNuBwjiAITYYI9FqqUmvYeuYa64+nYGpkyJoJnZnao239LdtXV0Haj9IslZT/gEaFtrU/jPgHBE0ES/v6OY8gCM2GCPRakGWVsCgijrjbxQwJcGHFuCBc7czr5+C5yRD7DVzZI63gtHT87wPOmVwvNGxxTYwEQXhyItBrQFml5rOT6Ww8mYadhQn/mtGNkZ3d6n5XXlEECRHSTJU70dKGEb7DpHFxv2H/W8FZ2LJajAqCUDMi0J9QzM27LNofR2puGeO7efDuqEBaW9VhJolGA9d+gcs7pb04qxTgHAhPr5Q6HFo711vtgiC0DCLQH6O8sor1x1PYeuYarrbmbHuuBwP96xC2hRnSfPHYXVByG8ztpPniXWeCezexglMQhFoTgf4IZ9LyWRwZx63CCmaFebJouD82tWmmpSyDpIPSkMqNM4ABeA+Cpz+AjiPBpJ7G3wVBaNFEoD9AcYWKVUdk7Ll0Cy9HK/a8HEavDg41O4hWK3U0vLwTEg+ASg72HWDQUqkplp1HwxQvCEKLJQL9T44nZvPOdwnklyl5pb/UTMvcpAbNtIpvS5sox34rDa+YWkPQeOg6CzzDxJCKIAgNRgT6f+WVKll2OJEjcVn4u9qw5dnuBLd5wmZaKoX0YDN2J6SfBLTQ7inotwACxoCZdYPWLgiCACLQ0Wq1fBd7h/cPJ1GuVPP3p/14pb83JkaPWX2p1UJmjDSkkrBfWoZv11YK8a7TpeEVQRCERtSiA/1OUQVLDsTzy9U8QjylZlo+zo9pplWWC3F7pCDPk4GxubQMv+tM8OovluELgqAzLTLQNRotO6NusOaHZDRaeG90IHPC2z+8mZZaJS2/j90pbaSsqZL6jI/6WOozbtHAfc4FQRCeQIsL9Iy8MhZHxHPheiFP+TiyekJn2to/pJlWTqJ0Jx63B8rzpY2Uw/4i3Y07+zdu4YIgCI/RYgK9Sq1h8+lrfHwiBXNjQ9ZNCmZyaJv7l+2XF0rL8C9/A1mxYGgCHYdLs1R8hoBRi/mSCYLQxLSIdErKLGFhxBUS7pQwrJMLy8cG4Wz7h8U8GrU0OyX2G0g+AupKcOkMw9dA5ylgVcM56IIgCDrQrANdoVLzr5/T2PRrOq0sTfl8ZgjPdHb73yfkp0nj4ld2Q2kmWNhD9+elIRW3YN0VLgiCUAt1CvQff/yRY8eOsX79+sGmjysAAAd3SURBVPteW7FiBTExMVhZWQGwceNGbGwaaDu2B4i+UcjC/XGk58mZGNKGpaMCaGVpCspSaeXm5Z1w6zwYGILPUHhmDfgNB2OzRqtREAShPtU60FesWMFvv/320P7ciYmJbNmyBXv7xt2IQa6s4sP/XGX7ueu421mw/fme9PdxkHqoxO6UeqqoysHBF4Ysg+BpYOv2mKMKgiDov1oHekhICEOGDGHPnj33vabRaLhx4wbvvvsu+fn5TJo0iUmTJtWp0CdxKiWPtyLjuVNUwbPh7VgYboVV0lY4uhOKboCZrdSatussaNNdLMMXBKFZeWyg79u3j+3bt9/zsVWrVjFixAiioqIe+J7y8nJmzZrF3LlzUavVzJkzh6CgIPz9753qJ5PVfsMGhUJR/f5SpZrNFwv4Mb0MbxsNe7rJCMr8CMuN0YAWuUt3inrNpbTNALTG5lAGJCfX+ty68sdrbila2jW3tOsFcc316bGBPnnyZCZPnlyjg1pYWDBnzhwsLCwACAsLIzk5+b5Ar8t2ajKZjICAAI4lZLH0+wTaVSRy0DOG4OKf/3979xcSVRqGAfyZPyrorJksXbShNKmQuWEaSrtmQdAfSCKpyGAKA9E0TK3QJoaMZgOrmzIW6qIE6yayLiODiGnW9GIqYURqg2pxMqkU1hkCz8x592Jxtkh3G51tnM/nd+U534ff+5yRV/k45wjD4DiQlgmsPwbkVyAlLQMpM15p7pjMPJ/Mt8zzLS/AzJHyeDzTjv0vd7m8evUKjY2NuH37NnRdx+PHj7F9+/aorjH6MYjmjm4s/L0LXUluZCQMAWPJQO62v+9SyfyZj+ET0bwS1YZ+9epVZGRkYMOGDSgrK8OuXbuQkJCAbdu2ITs7O2rr9HTfQPC3X/EL+mFO0KH/UAysagZWbAeSvt2dNEREc8msGnpxcTGKi4vDx5WVleGvq6qqUFVVNZtvP6WgpiGvpwETxiT8uaoW6T9Vwvh9VtTXISKKN3H3YJE5IQHGpgGM/DGEFXk/xrocIqI5Iy43mS2pC2HkO1WIiD4Tlw2diIi+xIZORKQINnQiIkWwoRMRKYINnYhIEWzoRESKYEMnIlKEQUQkFgv/2wtmiIhoeoWFhVOej1lDJyKi6OKWCxGRItjQiYgUEVcvRNF1Ha2trXj27BkSExPhdDqRmZkZ67KiRtM02O12+Hw+TExM4MCBA8jKykJLSwsMBgOys7Nx4sQJGI1GXLx4EQ8ePIDZbIbdbsfKlStjXf6MffjwAeXl5bhy5QrMZrPyeQHg0qVLuH//PjRNQ0VFBYqKipTNrWkaWlpa4PP5YDQacerUKaU/5/7+fpw7dw6dnZ14/fr1V+ecbm5EJI7cvXtXmpubRUTkyZMnUlNTE+OKouvmzZvidDpFRGR0dFTWrVsn1dXV0tvbKyIiDodDuru7xev1is1mE13XxefzSXl5eSzLnpWJiQmpra2VjRs3yosXL5TPKyLS29sr1dXVEgqFxO/3y4ULF5TOfe/ePamvrxcREbfbLQcPHlQ27+XLl2Xr1q2yc+dOEZGIck41N1JxteXi8Xiwdu1aAEB+fj68Xm+MK4quzZs349ChQ+Fjk8mEgYEBFBUVAQBKS0vR09MDj8eDkpISGAwGLF68GKFQCKOjo7Eqe1ba2tqwe/duLFq0CACUzwsAbrcbOTk5qKurQ01NDdavX6907qVLlyIUCkHXdfj9fpjNZmXzZmRkoL29PXwcSc6p5kYqrhq63++HxWIJH5tMJgSDwRhWFF0pKSmwWCzw+/2or69HQ0MDRAQGgyE8Pj4+/sV1mDwfb27duoX09PTwL2kASuedNDY2Bq/Xi/Pnz+PkyZM4cuSI0rmTk5Ph8/mwZcsWOBwO2Gw2ZfNu2rQJZvM/O9mR5JxqbqTiag/dYrEgEAiEj3Vd/+ziqWB4eBh1dXXYs2cPysrKcPbs2fBYIBBAamrqF9chEAjgu+/i71/vdXV1wWAw4NGjRxgcHERzc/Nnf5GplndSWloarFYrEhMTYbVakZSUhLdv34bHVcvd0dGBkpISHD58GMPDw9i3bx80TQuPq5b3U5/ugf9XzqnmRrze7Mr9tgoKCuByuQAAT58+RU5OTowriq73799j//79OHr0KHbs2AEAyM3NRV9fHwDA5XJh9erVKCgogNvthq7rePPmDXRdR3p6eixLn5Hr16/j2rVr6OzsxPLly9HW1obS0lJl804qLCzEw4cPISIYGRnBx48fsWbNGmVzp6amhhvzggULEAwGlf65/lQkOaeaG6m4erBo8i6X58+fQ0Rw+vRpLFu2LNZlRY3T6cSdO3dgtVrD544fPw6n0wlN02C1WuF0OmEymdDe3g6XywVd13Hs2LEZffhzic1mQ2trK4xGIxwOh/J5z5w5g76+PogIGhsbsWTJEmVzBwIB2O12vHv3DpqmYe/evcjLy1M279DQEJqamnDjxg28fPnyq3NONzcScdXQiYhoenG15UJERNNjQyciUgQbOhGRItjQiYgUwYZORKQINnQiIkWwoRMRKYINnYhIEX8B02Kb2UrU0QUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "theta = []\n", "thetaprime = []\n", "thetadiff = []\n", "\n", "b = 0.723\n", "pi = 3.141592653\n", "\n", "for i in range(-500,500):\n", " thetai = pi*i/1000\n", " y = np.sin(thetai)\n", " x = np.cos(thetai)\n", " yprime = y - b\n", " xprime = x\n", " thetaprimei = np.arctan(yprime/xprime)\n", " thetadiffi = thetai - thetaprimei\n", " theta = theta + [thetai]\n", " thetaprime = thetaprime + [thetaprimei]\n", " thetadiff = thetadiff + [thetadiffi]\n", "\n", "data = [('theta', theta),('thetaprime', thetaprime), ('thetadiff', thetadiff)]\n", "venus = pd.DataFrame.from_items(data)\n", "venus.plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Venus\n", "\n", "\n", "\n", "## Catch Our Breath—Further Notes:\n", "\n", "* Questions\n", "* Comments\n", "* Readings \n", "\n", "
\n", "\n", "----\n", "\n", "* Weblog Support \n", "* nbViewer \n", "\n", "\n", " \n", "\n", "----" ] } ], "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.6.6" } }, "nbformat": 4, "nbformat_minor": 2 }