{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Exoplanets In My Solar System" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### This notebook is based on the November 2019 transit of Mercury." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Planetary transits in the solar system are perfect analogies of explanetary transits. If sufficiently intelligent aliens were monitoring the brightness of our star on Monday 11th November 2019, they may have detected something strange." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "![caption](files/transit_white.png)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The above image is a simplified diagram of how planetary transits can be detected by measuring the dip in brightness of the host star, known as the \"Transit Method\". Below is an image I captured of the transit on the 11th. Note the tiny black dot. Thats Mercury!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "![caption](files/my_mercury.jpg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The change in brightness $\\Delta F$ is related to the size of the planet $R_P$ through a very simple formula." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$\\Delta F / F = (R_P/R_*)^2$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets now import some `python` libraries and get `Juypter` to do the calculation for us." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'The radius of Mercury is 2530000.0, and the radius of the Sun is 695000000.0'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import astropy.units as u\n", "F = 1370 # Stellar Flux\n", "R_p = 2.53E6 # Radius of Mercury in m\n", "R_s = 6.95E8 # Radius of the Sun in m\n", "\n", "F\"The radius of Mercury is {R_p}, and the radius of the Sun is {R_s}\"" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'The drop in flux due to the transit of Mercury is 1.3251695046840224e-05'" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "delta_f_f = ((R_p/R_s)**2)\n", "F\"The drop in flux due to the transit of Mercury is {delta_f_f}\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Using the [BATMAN](https://www.cfa.harvard.edu/~lkreidberg/batman/index.html) python code, it is possible to simulate what the Mercury transit lightcurve would look like observer from a large distance." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import batman\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "params = batman.TransitParams()\n", "params.t0 = 0. #time of inferior conjunction\n", "params.per = 88. #orbital period\n", "params.rp = R_p/R_s #planet radius (in units of stellar radii)\n", "params.a = 5.7E7/R_s #semi-major axis (in units of stellar radii)\n", "params.inc = 87. #orbital inclination (in degrees)\n", "params.ecc = 0. #eccentricity\n", "params.w = 90. #longitude of periastron (in degrees)\n", "params.u = [0.1, 0.3] #limb darkening coefficients [u1, u2]\n", "params.limb_dark = \"quadratic\" #limb darkening model\n", "\n", "t = np.linspace(-30, 30, 100)\n", "\n", "m = batman.TransitModel(params, t) #initializes model\n", "flux = m.light_curve(params) #calculates light curve\n", "\n", "plt.plot(t, flux)\n", "plt.xlabel(\"Time from central transit\")\n", "plt.ylabel(\"Relative flux\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The purpose of code is to make life easier. In the sections above, I have defined the variables `R_s` and `R_p` but going forward I want to be able to convert planet radii in meters to solar radii easily, so I will create a small subroutine to do this." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "def cvnt_sol(radius):\n", " rad_sol = radius/R_s\n", " return rad_sol " ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'In this example I have converted the radius of Venus into solar radii. The result is 0.008705035971223022'" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "example = cvnt_sol(6.05E6)\n", "F\"In this example I have converted the radius of Venus into solar radii. The result is {example}\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets now use `PANDAS` to do some data manipulation." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "import pandas as pd # Import the PANDAS library" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "file = \"solar_system_data_1.csv\" # Assign the data file to a variable." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "df = pd.read_csv(file) # Read in the data file into a dataframe." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMass(10^24kg)Radius(km)Diameter(km)Density (kg/m3)Gravity (m/s2)Escape Velocity (km/s)Rotation Period (hours)Length of Day (hours)Distance from Sun (106 km)...Orbital Period (days)Orbital Velocity (km/s)Orbital Inclination (degrees)Orbital EccentricityObliquity to Orbit (degrees)Mean Temperature (C)Surface Pressure (bars)Number of MoonsRing System?Global Magnetic Field?
0MERCURY0.33002439.5487954273.74.31407.64222.657.900...8847.47.00.2050.03416700NoYes
1VENUS4.87006052.012,10452438.910.4-5832.52802.0108.200...224.735.03.40.007177.400464920NoNo
2EARTH5.97006378.012,75655149.811.223.924.0149.600...365.229.80.00.01723.4001511NoYes
3MOON0.07301737.5347533401.62.4655.7708.70.384...27.31.05.10.0556.700-2000NoNo
4MARS0.64203396.0679239333.75.024.624.7227.900...68724.11.90.09425.200-650.012NoNo
5JUPITER1898.000071492.0142,984132623.159.59.99.9778.600...433113.11.30.0493.100-110Unknown79YesYes
6SATURN568.000060268.0120,5366879.035.510.710.71433.500...10,7479.72.50.05726.700-140Unknown82YesYes
7URANUS86.800025559.051,11812718.721.3-17.217.22872.500...30,5896.80.80.04697.800-195Unknown27YesYes
8NEPTUNE102.000024764.049,528163811.023.516.116.14495.100...59,8005.41.80.01128.300-200Unknown14YesYes
9PLUTO0.01461185.0237020950.71.3-153.3153.35906.400...90,5604.717.20.244122.500-2250.000015NoUnknown
\n", "

10 rows × 22 columns

\n", "
" ], "text/plain": [ " Planet Mass(10^24kg) Radius(km) Diameter(km) Density (kg/m3) \\\n", "0  MERCURY  0.3300 2439.5 4879 5427 \n", "1  VENUS  4.8700 6052.0 12,104 5243 \n", "2  EARTH  5.9700 6378.0 12,756 5514 \n", "3  MOON  0.0730 1737.5 3475 3340 \n", "4  MARS  0.6420 3396.0 6792 3933 \n", "5  JUPITER  1898.0000 71492.0 142,984 1326 \n", "6  SATURN  568.0000 60268.0 120,536 687 \n", "7  URANUS  86.8000 25559.0 51,118 1271 \n", "8  NEPTUNE  102.0000 24764.0 49,528 1638 \n", "9  PLUTO  0.0146 1185.0 2370 2095 \n", "\n", " Gravity (m/s2) Escape Velocity (km/s) Rotation Period (hours) \\\n", "0 3.7 4.3 1407.6 \n", "1 8.9 10.4 -5832.5 \n", "2 9.8 11.2 23.9 \n", "3 1.6 2.4 655.7 \n", "4 3.7 5.0 24.6 \n", "5 23.1 59.5 9.9 \n", "6 9.0 35.5 10.7 \n", "7 8.7 21.3 -17.2 \n", "8 11.0 23.5 16.1 \n", "9 0.7 1.3 -153.3 \n", "\n", " Length of Day (hours) Distance from Sun (106 km) ... \\\n", "0 4222.6 57.900 ... \n", "1 2802.0 108.200 ... \n", "2 24.0 149.600 ... \n", "3 708.7 0.384 ... \n", "4 24.7 227.900 ... \n", "5 9.9 778.600 ... \n", "6 10.7 1433.500 ... \n", "7 17.2 2872.500 ... \n", "8 16.1 4495.100 ... \n", "9 153.3 5906.400 ... \n", "\n", " Orbital Period (days) Orbital Velocity (km/s) \\\n", "0 88 47.4 \n", "1 224.7 35.0 \n", "2 365.2 29.8 \n", "3 27.3 1.0 \n", "4 687 24.1 \n", "5 4331 13.1 \n", "6 10,747 9.7 \n", "7 30,589 6.8 \n", "8 59,800 5.4 \n", "9 90,560 4.7 \n", "\n", " Orbital Inclination (degrees) Orbital Eccentricity \\\n", "0 7.0 0.205 \n", "1 3.4 0.007 \n", "2 0.0 0.017 \n", "3 5.1 0.055 \n", "4 1.9 0.094 \n", "5 1.3 0.049 \n", "6 2.5 0.057 \n", "7 0.8 0.046 \n", "8 1.8 0.011 \n", "9 17.2 0.244 \n", "\n", " Obliquity to Orbit (degrees) Mean Temperature (C) \\\n", "0 0.034 167 \n", "1 177.400 464 \n", "2 23.400 15 \n", "3 6.700 -20 \n", "4 25.200 -65 \n", "5 3.100 -110 \n", "6 26.700 -140 \n", "7 97.800 -195 \n", "8 28.300 -200 \n", "9 122.500 -225 \n", "\n", " Surface Pressure (bars) Number of Moons Ring System? \\\n", "0 0 0 No \n", "1 92 0 No \n", "2 1 1 No \n", "3 0 0 No \n", "4 0.01 2 No \n", "5 Unknown 79 Yes \n", "6 Unknown 82 Yes \n", "7 Unknown 27 Yes \n", "8 Unknown 14 Yes \n", "9 0.00001 5 No \n", "\n", " Global Magnetic Field? \n", "0 Yes \n", "1 No \n", "2 Yes \n", "3 No \n", "4 No \n", "5 Yes \n", "6 Yes \n", "7 Yes \n", "8 Yes \n", "9 Unknown \n", "\n", "[10 rows x 22 columns]" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df # Display the dataframe." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "# Convert the Mass from units of 10^24kg to kg.\n", "def to_24(mass):\n", " mass_24 = mass * 1E24\n", " return mass_24" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMass(10^24kg)Radius(km)Diameter(km)Density (kg/m3)Gravity (m/s2)Escape Velocity (km/s)Rotation Period (hours)Length of Day (hours)Distance from Sun (106 km)...Orbital Velocity (km/s)Orbital Inclination (degrees)Orbital EccentricityObliquity to Orbit (degrees)Mean Temperature (C)Surface Pressure (bars)Number of MoonsRing System?Global Magnetic Field?Mass
0MERCURY0.33002439.5487954273.74.31407.64222.657.900...47.47.00.2050.03416700NoYes3.300000e+23
1VENUS4.87006052.012,10452438.910.4-5832.52802.0108.200...35.03.40.007177.400464920NoNo4.870000e+24
2EARTH5.97006378.012,75655149.811.223.924.0149.600...29.80.00.01723.4001511NoYes5.970000e+24
3MOON0.07301737.5347533401.62.4655.7708.70.384...1.05.10.0556.700-2000NoNo7.300000e+22
4MARS0.64203396.0679239333.75.024.624.7227.900...24.11.90.09425.200-650.012NoNo6.420000e+23
5JUPITER1898.000071492.0142,984132623.159.59.99.9778.600...13.11.30.0493.100-110Unknown79YesYes1.898000e+27
6SATURN568.000060268.0120,5366879.035.510.710.71433.500...9.72.50.05726.700-140Unknown82YesYes5.680000e+26
7URANUS86.800025559.051,11812718.721.3-17.217.22872.500...6.80.80.04697.800-195Unknown27YesYes8.680000e+25
8NEPTUNE102.000024764.049,528163811.023.516.116.14495.100...5.41.80.01128.300-200Unknown14YesYes1.020000e+26
9PLUTO0.01461185.0237020950.71.3-153.3153.35906.400...4.717.20.244122.500-2250.000015NoUnknown1.460000e+22
\n", "

10 rows × 23 columns

\n", "
" ], "text/plain": [ " Planet Mass(10^24kg) Radius(km) Diameter(km) Density (kg/m3) \\\n", "0  MERCURY  0.3300 2439.5 4879 5427 \n", "1  VENUS  4.8700 6052.0 12,104 5243 \n", "2  EARTH  5.9700 6378.0 12,756 5514 \n", "3  MOON  0.0730 1737.5 3475 3340 \n", "4  MARS  0.6420 3396.0 6792 3933 \n", "5  JUPITER  1898.0000 71492.0 142,984 1326 \n", "6  SATURN  568.0000 60268.0 120,536 687 \n", "7  URANUS  86.8000 25559.0 51,118 1271 \n", "8  NEPTUNE  102.0000 24764.0 49,528 1638 \n", "9  PLUTO  0.0146 1185.0 2370 2095 \n", "\n", " Gravity (m/s2) Escape Velocity (km/s) Rotation Period (hours) \\\n", "0 3.7 4.3 1407.6 \n", "1 8.9 10.4 -5832.5 \n", "2 9.8 11.2 23.9 \n", "3 1.6 2.4 655.7 \n", "4 3.7 5.0 24.6 \n", "5 23.1 59.5 9.9 \n", "6 9.0 35.5 10.7 \n", "7 8.7 21.3 -17.2 \n", "8 11.0 23.5 16.1 \n", "9 0.7 1.3 -153.3 \n", "\n", " Length of Day (hours) Distance from Sun (106 km) ... \\\n", "0 4222.6 57.900 ... \n", "1 2802.0 108.200 ... \n", "2 24.0 149.600 ... \n", "3 708.7 0.384 ... \n", "4 24.7 227.900 ... \n", "5 9.9 778.600 ... \n", "6 10.7 1433.500 ... \n", "7 17.2 2872.500 ... \n", "8 16.1 4495.100 ... \n", "9 153.3 5906.400 ... \n", "\n", " Orbital Velocity (km/s) Orbital Inclination (degrees) \\\n", "0 47.4 7.0 \n", "1 35.0 3.4 \n", "2 29.8 0.0 \n", "3 1.0 5.1 \n", "4 24.1 1.9 \n", "5 13.1 1.3 \n", "6 9.7 2.5 \n", "7 6.8 0.8 \n", "8 5.4 1.8 \n", "9 4.7 17.2 \n", "\n", " Orbital Eccentricity Obliquity to Orbit (degrees) Mean Temperature (C) \\\n", "0 0.205 0.034 167 \n", "1 0.007 177.400 464 \n", "2 0.017 23.400 15 \n", "3 0.055 6.700 -20 \n", "4 0.094 25.200 -65 \n", "5 0.049 3.100 -110 \n", "6 0.057 26.700 -140 \n", "7 0.046 97.800 -195 \n", "8 0.011 28.300 -200 \n", "9 0.244 122.500 -225 \n", "\n", " Surface Pressure (bars) Number of Moons Ring System? \\\n", "0 0 0 No \n", "1 92 0 No \n", "2 1 1 No \n", "3 0 0 No \n", "4 0.01 2 No \n", "5 Unknown 79 Yes \n", "6 Unknown 82 Yes \n", "7 Unknown 27 Yes \n", "8 Unknown 14 Yes \n", "9 0.00001 5 No \n", "\n", " Global Magnetic Field? Mass \n", "0 Yes 3.300000e+23 \n", "1 No 4.870000e+24 \n", "2 Yes 5.970000e+24 \n", "3 No 7.300000e+22 \n", "4 No 6.420000e+23 \n", "5 Yes 1.898000e+27 \n", "6 Yes 5.680000e+26 \n", "7 Yes 8.680000e+25 \n", "8 Yes 1.020000e+26 \n", "9 Unknown 1.460000e+22 \n", "\n", "[10 rows x 23 columns]" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df[\"Mass\"] = df[\"Mass(10^24kg)\"].apply(to_24)\n", "df" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "# Define a function to calculate the Planet - Star Radius Ratio\n", "def RR(x):\n", " RR = (x / R_s)\n", " return RR" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMass(10^24kg)Radius(km)Diameter(km)Density (kg/m3)Gravity (m/s2)Escape Velocity (km/s)Rotation Period (hours)Length of Day (hours)Distance from Sun (106 km)...Orbital Inclination (degrees)Orbital EccentricityObliquity to Orbit (degrees)Mean Temperature (C)Surface Pressure (bars)Number of MoonsRing System?Global Magnetic Field?MassRadius Ratio
0MERCURY0.33002439.5487954273.74.31407.64222.657.900...7.00.2050.03416700NoYes3.300000e+230.000004
1VENUS4.87006052.012,10452438.910.4-5832.52802.0108.200...3.40.007177.400464920NoNo4.870000e+240.000009
2EARTH5.97006378.012,75655149.811.223.924.0149.600...0.00.01723.4001511NoYes5.970000e+240.000009
3MOON0.07301737.5347533401.62.4655.7708.70.384...5.10.0556.700-2000NoNo7.300000e+220.000003
4MARS0.64203396.0679239333.75.024.624.7227.900...1.90.09425.200-650.012NoNo6.420000e+230.000005
5JUPITER1898.000071492.0142,984132623.159.59.99.9778.600...1.30.0493.100-110Unknown79YesYes1.898000e+270.000103
6SATURN568.000060268.0120,5366879.035.510.710.71433.500...2.50.05726.700-140Unknown82YesYes5.680000e+260.000087
7URANUS86.800025559.051,11812718.721.3-17.217.22872.500...0.80.04697.800-195Unknown27YesYes8.680000e+250.000037
8NEPTUNE102.000024764.049,528163811.023.516.116.14495.100...1.80.01128.300-200Unknown14YesYes1.020000e+260.000036
9PLUTO0.01461185.0237020950.71.3-153.3153.35906.400...17.20.244122.500-2250.000015NoUnknown1.460000e+220.000002
\n", "

10 rows × 24 columns

\n", "
" ], "text/plain": [ " Planet Mass(10^24kg) Radius(km) Diameter(km) Density (kg/m3) \\\n", "0  MERCURY  0.3300 2439.5 4879 5427 \n", "1  VENUS  4.8700 6052.0 12,104 5243 \n", "2  EARTH  5.9700 6378.0 12,756 5514 \n", "3  MOON  0.0730 1737.5 3475 3340 \n", "4  MARS  0.6420 3396.0 6792 3933 \n", "5  JUPITER  1898.0000 71492.0 142,984 1326 \n", "6  SATURN  568.0000 60268.0 120,536 687 \n", "7  URANUS  86.8000 25559.0 51,118 1271 \n", "8  NEPTUNE  102.0000 24764.0 49,528 1638 \n", "9  PLUTO  0.0146 1185.0 2370 2095 \n", "\n", " Gravity (m/s2) Escape Velocity (km/s) Rotation Period (hours) \\\n", "0 3.7 4.3 1407.6 \n", "1 8.9 10.4 -5832.5 \n", "2 9.8 11.2 23.9 \n", "3 1.6 2.4 655.7 \n", "4 3.7 5.0 24.6 \n", "5 23.1 59.5 9.9 \n", "6 9.0 35.5 10.7 \n", "7 8.7 21.3 -17.2 \n", "8 11.0 23.5 16.1 \n", "9 0.7 1.3 -153.3 \n", "\n", " Length of Day (hours) Distance from Sun (106 km) ... \\\n", "0 4222.6 57.900 ... \n", "1 2802.0 108.200 ... \n", "2 24.0 149.600 ... \n", "3 708.7 0.384 ... \n", "4 24.7 227.900 ... \n", "5 9.9 778.600 ... \n", "6 10.7 1433.500 ... \n", "7 17.2 2872.500 ... \n", "8 16.1 4495.100 ... \n", "9 153.3 5906.400 ... \n", "\n", " Orbital Inclination (degrees) Orbital Eccentricity \\\n", "0 7.0 0.205 \n", "1 3.4 0.007 \n", "2 0.0 0.017 \n", "3 5.1 0.055 \n", "4 1.9 0.094 \n", "5 1.3 0.049 \n", "6 2.5 0.057 \n", "7 0.8 0.046 \n", "8 1.8 0.011 \n", "9 17.2 0.244 \n", "\n", " Obliquity to Orbit (degrees) Mean Temperature (C) Surface Pressure (bars) \\\n", "0 0.034 167 0 \n", "1 177.400 464 92 \n", "2 23.400 15 1 \n", "3 6.700 -20 0 \n", "4 25.200 -65 0.01 \n", "5 3.100 -110 Unknown \n", "6 26.700 -140 Unknown \n", "7 97.800 -195 Unknown \n", "8 28.300 -200 Unknown \n", "9 122.500 -225 0.00001 \n", "\n", " Number of Moons Ring System? Global Magnetic Field? Mass \\\n", "0 0 No Yes 3.300000e+23 \n", "1 0 No No 4.870000e+24 \n", "2 1 No Yes 5.970000e+24 \n", "3 0 No No 7.300000e+22 \n", "4 2 No No 6.420000e+23 \n", "5 79 Yes Yes 1.898000e+27 \n", "6 82 Yes Yes 5.680000e+26 \n", "7 27 Yes Yes 8.680000e+25 \n", "8 14 Yes Yes 1.020000e+26 \n", "9 5 No Unknown 1.460000e+22 \n", "\n", " Radius Ratio \n", "0 0.000004 \n", "1 0.000009 \n", "2 0.000009 \n", "3 0.000003 \n", "4 0.000005 \n", "5 0.000103 \n", "6 0.000087 \n", "7 0.000037 \n", "8 0.000036 \n", "9 0.000002 \n", "\n", "[10 rows x 24 columns]" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df[\"Radius Ratio\"] = df[\"Radius(km)\"].apply(RR)\n", "df" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "def cvnt_m(km):\n", " m = km * 1000\n", " return m\n" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMass(10^24kg)Radius(km)Diameter(km)Density (kg/m3)Gravity (m/s2)Escape Velocity (km/s)Rotation Period (hours)Length of Day (hours)Distance from Sun (106 km)...Orbital EccentricityObliquity to Orbit (degrees)Mean Temperature (C)Surface Pressure (bars)Number of MoonsRing System?Global Magnetic Field?MassRadius RatioRadius (m)
0MERCURY0.33002439.5487954273.74.31407.64222.657.900...0.2050.03416700NoYes3.300000e+230.0000042439500.0
1VENUS4.87006052.012,10452438.910.4-5832.52802.0108.200...0.007177.400464920NoNo4.870000e+240.0000096052000.0
2EARTH5.97006378.012,75655149.811.223.924.0149.600...0.01723.4001511NoYes5.970000e+240.0000096378000.0
3MOON0.07301737.5347533401.62.4655.7708.70.384...0.0556.700-2000NoNo7.300000e+220.0000031737500.0
4MARS0.64203396.0679239333.75.024.624.7227.900...0.09425.200-650.012NoNo6.420000e+230.0000053396000.0
5JUPITER1898.000071492.0142,984132623.159.59.99.9778.600...0.0493.100-110Unknown79YesYes1.898000e+270.00010371492000.0
6SATURN568.000060268.0120,5366879.035.510.710.71433.500...0.05726.700-140Unknown82YesYes5.680000e+260.00008760268000.0
7URANUS86.800025559.051,11812718.721.3-17.217.22872.500...0.04697.800-195Unknown27YesYes8.680000e+250.00003725559000.0
8NEPTUNE102.000024764.049,528163811.023.516.116.14495.100...0.01128.300-200Unknown14YesYes1.020000e+260.00003624764000.0
9PLUTO0.01461185.0237020950.71.3-153.3153.35906.400...0.244122.500-2250.000015NoUnknown1.460000e+220.0000021185000.0
\n", "

10 rows × 25 columns

\n", "
" ], "text/plain": [ " Planet Mass(10^24kg) Radius(km) Diameter(km) Density (kg/m3) \\\n", "0  MERCURY  0.3300 2439.5 4879 5427 \n", "1  VENUS  4.8700 6052.0 12,104 5243 \n", "2  EARTH  5.9700 6378.0 12,756 5514 \n", "3  MOON  0.0730 1737.5 3475 3340 \n", "4  MARS  0.6420 3396.0 6792 3933 \n", "5  JUPITER  1898.0000 71492.0 142,984 1326 \n", "6  SATURN  568.0000 60268.0 120,536 687 \n", "7  URANUS  86.8000 25559.0 51,118 1271 \n", "8  NEPTUNE  102.0000 24764.0 49,528 1638 \n", "9  PLUTO  0.0146 1185.0 2370 2095 \n", "\n", " Gravity (m/s2) Escape Velocity (km/s) Rotation Period (hours) \\\n", "0 3.7 4.3 1407.6 \n", "1 8.9 10.4 -5832.5 \n", "2 9.8 11.2 23.9 \n", "3 1.6 2.4 655.7 \n", "4 3.7 5.0 24.6 \n", "5 23.1 59.5 9.9 \n", "6 9.0 35.5 10.7 \n", "7 8.7 21.3 -17.2 \n", "8 11.0 23.5 16.1 \n", "9 0.7 1.3 -153.3 \n", "\n", " Length of Day (hours) Distance from Sun (106 km) ... \\\n", "0 4222.6 57.900 ... \n", "1 2802.0 108.200 ... \n", "2 24.0 149.600 ... \n", "3 708.7 0.384 ... \n", "4 24.7 227.900 ... \n", "5 9.9 778.600 ... \n", "6 10.7 1433.500 ... \n", "7 17.2 2872.500 ... \n", "8 16.1 4495.100 ... \n", "9 153.3 5906.400 ... \n", "\n", " Orbital Eccentricity Obliquity to Orbit (degrees) Mean Temperature (C) \\\n", "0 0.205 0.034 167 \n", "1 0.007 177.400 464 \n", "2 0.017 23.400 15 \n", "3 0.055 6.700 -20 \n", "4 0.094 25.200 -65 \n", "5 0.049 3.100 -110 \n", "6 0.057 26.700 -140 \n", "7 0.046 97.800 -195 \n", "8 0.011 28.300 -200 \n", "9 0.244 122.500 -225 \n", "\n", " Surface Pressure (bars) Number of Moons Ring System? \\\n", "0 0 0 No \n", "1 92 0 No \n", "2 1 1 No \n", "3 0 0 No \n", "4 0.01 2 No \n", "5 Unknown 79 Yes \n", "6 Unknown 82 Yes \n", "7 Unknown 27 Yes \n", "8 Unknown 14 Yes \n", "9 0.00001 5 No \n", "\n", " Global Magnetic Field? Mass Radius Ratio Radius (m) \n", "0 Yes 3.300000e+23 0.000004 2439500.0 \n", "1 No 4.870000e+24 0.000009 6052000.0 \n", "2 Yes 5.970000e+24 0.000009 6378000.0 \n", "3 No 7.300000e+22 0.000003 1737500.0 \n", "4 No 6.420000e+23 0.000005 3396000.0 \n", "5 Yes 1.898000e+27 0.000103 71492000.0 \n", "6 Yes 5.680000e+26 0.000087 60268000.0 \n", "7 Yes 8.680000e+25 0.000037 25559000.0 \n", "8 Yes 1.020000e+26 0.000036 24764000.0 \n", "9 Unknown 1.460000e+22 0.000002 1185000.0 \n", "\n", "[10 rows x 25 columns]" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df[\"Radius (m)\"] = df[\"Radius(km)\"].apply(cvnt_m)\n", "df" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We'll take our data, and extract only the columns we're interested in, storing them in a new dataframe which we will call `df2`" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMassRadius (m)
0MERCURY3.300000e+232439500.0
1VENUS4.870000e+246052000.0
2EARTH5.970000e+246378000.0
3MOON7.300000e+221737500.0
4MARS6.420000e+233396000.0
5JUPITER1.898000e+2771492000.0
6SATURN5.680000e+2660268000.0
7URANUS8.680000e+2525559000.0
8NEPTUNE1.020000e+2624764000.0
9PLUTO1.460000e+221185000.0
\n", "
" ], "text/plain": [ " Planet Mass Radius (m)\n", "0  MERCURY  3.300000e+23 2439500.0\n", "1  VENUS  4.870000e+24 6052000.0\n", "2  EARTH  5.970000e+24 6378000.0\n", "3  MOON  7.300000e+22 1737500.0\n", "4  MARS  6.420000e+23 3396000.0\n", "5  JUPITER  1.898000e+27 71492000.0\n", "6  SATURN  5.680000e+26 60268000.0\n", "7  URANUS  8.680000e+25 25559000.0\n", "8  NEPTUNE  1.020000e+26 24764000.0\n", "9  PLUTO  1.460000e+22 1185000.0" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cols = ['Planet','Mass','Radius (m)']\n", "df2 = df[cols]\n", "df2" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets now create another subroutine which will take the `Radius (m)` column from the `df2` dataframe, and calculate the delta flux" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "def flux(Radius):\n", " flux = (Radius)**2 / (R_s **2)\n", " return flux" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\mwood\\AppData\\Local\\Programs\\Python\\Python36-32\\lib\\site-packages\\ipykernel_launcher.py:1: SettingWithCopyWarning: \n", "A value is trying to be set on a copy of a slice from a DataFrame.\n", "Try using .loc[row_indexer,col_indexer] = value instead\n", "\n", "See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy\n", " \"\"\"Entry point for launching an IPython kernel.\n" ] } ], "source": [ "df2[\"Delta Flux\"] = df2[\"Radius (m)\"].apply(flux)\n", "#dft = df2.drop(columns=\"Flux\")" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PlanetMassRadius (m)Delta Flux
0MERCURY3.300000e+232439500.00.000012
1VENUS4.870000e+246052000.00.000076
2EARTH5.970000e+246378000.00.000084
3MOON7.300000e+221737500.00.000006
4MARS6.420000e+233396000.00.000024
5JUPITER1.898000e+2771492000.00.010581
6SATURN5.680000e+2660268000.00.007520
7URANUS8.680000e+2525559000.00.001352
8NEPTUNE1.020000e+2624764000.00.001270
9PLUTO1.460000e+221185000.00.000003
\n", "
" ], "text/plain": [ " Planet Mass Radius (m) Delta Flux\n", "0  MERCURY  3.300000e+23 2439500.0 0.000012\n", "1  VENUS  4.870000e+24 6052000.0 0.000076\n", "2  EARTH  5.970000e+24 6378000.0 0.000084\n", "3  MOON  7.300000e+22 1737500.0 0.000006\n", "4  MARS  6.420000e+23 3396000.0 0.000024\n", "5  JUPITER  1.898000e+27 71492000.0 0.010581\n", "6  SATURN  5.680000e+26 60268000.0 0.007520\n", "7  URANUS  8.680000e+25 25559000.0 0.001352\n", "8  NEPTUNE  1.020000e+26 24764000.0 0.001270\n", "9  PLUTO  1.460000e+22 1185000.0 0.000003" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cols = ['Planet','Mass','Radius (m)','Delta Flux']\n", "data = df2[cols]\n", "data" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAFDCAYAAADyPJMJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xu8HfO9//HXW4LE/dIoEsS1btXQUHVpkdal1UZbKrROtMrRclo/SnFK0eqhF3WK9lSPlmoVpUidtHUUvR0VSUUJVUEQ17hFqIjE5/fH97uT2Strrb129qzZayfv5+OxH3vu85m1ZtZn5jvf+Y4iAjMzszIs198BmJnZ0sNJxczMSuOkYmZmpXFSMTOz0jipmJlZaZxUzMysNE4qfSDpDEk/LWlZb5X0B0lzJH27jGWWRdInJN3Uh/lHSgpJg8uMq7D8kLRZi9NK0o8lvShpUjviKYOk2yR9pqRlXSrpa2Usq49xbCjpFUmDSlpeacdfH2Joed9bVixVSUXSOEl3SHpV0rO5+3OS1N+xteAo4DlgtYg4oXZkf/4wRMTPImLvQiz9diCV8GO7G/B+YERE7FRSWEbP301EPBYRq0TEgiVY9h6SZvYtwmWLpMMl/anq9S41SUXSCcB/At8E1gXeChwN7Aqs0GCeUs6YSrIRcF8s4dOo7boKWAptBMyIiFd7O6M/YyvLUr0vRcSA/wNWB14FPtbDdJcC3wcm5unfB3wQuAt4GXgcOKMw/UggSFcRTwJPAScUxp8BXA38BJgDTANGN1n/LsCdwOz8f5dCXG8A84BXgPfVzHdUzfhf5eEzgC8BfwNeBwYDJwMP5XjuAz5SWM7hwJ+AbwEvAo8A+9WMfzjP+wjwieJ8ufsP+TN5NcdycJ3tHJTX8Vxe3jF5nsGF7+uS/Hk+AXwNGFRY15+BC/Ln9HdgTB53NrAAmJvXfWEeHqQTiAfzdl0EqE5cR+R5F+T5z8zDjwSmAy8AE4D1C/NEjv9B4JEG3+vOwP8BLwF3A3sUxn0KuD9/pg8D/1oz71hgKmn/ewjYNw+/Dfhq/izmADcBb2mw/j2AmcCp+TOf0fXdFfavr+XuNYEbgVn5s7qRdNVGK+tttK2NvpuaOEfW7ActbSOwMvAa8GZe9ivA+vRw/OVprs3b+gjw+SbH5gdIx8sc0j75xcK4nvaPzXJ3K78lRwCPAX+oE8Nb8vfxUl7XH0kn/icC19ZMewFwfqPjFtiK7vv6S3naFUnH5mPAM8B/AUNr9qOTgGdJx+cB+bP5R47p1B5/j8v6Ye/PP2BfYH7XztpkuktJP1S75i9rSP4g3577t8sf9AE1O8LP84799ryDvi+PPyN/cR8g/ZD+B/CXButei3QQH0b68T8k969de+A3if1rNcNmkH6QNijsGAeRDqblgINJP/7rFXa+N0gHySDgs6Rkqbx9LwNvy9OuB2xTmO9P9Q6kBrEeTUoGG+TtvpXuPybXAz/I61wHmET+sc3rmg/8P2D5vA2zgbUKP0SfqVlfkA7GNYAN83e0b4PYardlL9IP8Q6kA+4CCgd8Xvb/5u0YWmd5w4Hn8z6wHKlo7XlgWOGHZtP8Gb8X+CewQx63U9629+d5hwNbFrbzIWALYGjuP6fBNu2RP7Pz8ja8N3/vXd/lpSxKKmsDHwNWAlYFfgFcX1hWw/W2sK2LfTc1cY5k8aTSm22cWTPsDBocfzm+KcDppJKKTUg/vPs0WP5TwO65e83Cd9TK/rFZIcaefkt+Qtrv6+1L/0H6kV8+/+1O2m/Wy9/nGnm6waQf/XfSi+M2DzuflBjXyt//r4D/qNmPTs/rP5J0LF2Rp90mf96bNP2dXdIf8k76Az4JPF0zrOts6jXgPYWD6yc9LOt84Ds1O8KWhfHfAC4p7NQ3F8ZtDbzWYLmHAZNqht0OHF574DeYf7HxpKTy6R62ZyowtrCTTS+MWylv37p553yJ9IMztGYZ3XZOek4qtwBHF/r3zvMMJhVLvl5cBynB3lpY15MUrjRISeew3H0b9ZPKboX+q4GTG8RWuy2XAN8o9K9CSrwjC8veq8m2fgm4vGbYb4HxDaa/HvhC7v5B175WZ7rbgC8X+j8H/KbBtHuQfgxWrvkMTutp3wJGAS+2st6etrXed1Mz7UgWTyq92cZ6SaXu8Qe8C3isZvpTgB83WP5jwL+S7mkWh7eyf9Q9Fqj/W9LwBxk4C7ih3vKAXwNH5u79SUXl0LvjVqTktGlh2LvJV+D5M36NRaUGq+aY31WYfgo5UTb6W1ruqTwPvKVYThkRu0TEGnlccTsfL84o6V2SbpU0S9Js0ln2W2qWX5znUdKVQJenC93/BIY0KC9dP89b9Cjp7K8varfnXyRNlfSSpJeAbem+PQvjjYh/5s5VIt1jOJi0/U9J+h9JWy5hTOuz+GfWZSPSWdBThRh/QLpi6fJE5D24MH/xM6+n9ntYpRexLowvIl4h7TPF7+Xx2pkKNgIO6tqWvD27kc4YkbSfpL9IeiGP+wCLvo8NSGfqZWzTi9H9PlHdz0zSSpJ+IOlRSS+TijPXqLm/2Gi9Tbd1CS3p99Zo/q7jbyNg/ZpYTyWd1NTzMdJ386ik30t6dx7eyv4BLNFvSa1vkorZbpL0sKSTC+MuI508k/9fnuPpzXE7jHQiOaXwmfwmD+/yfCyqSPFa/v9MYfxr9PAdLS1J5XbS2e/YFqaNmv4rSJeDG0TE6qTLz9raYhsUujcknUn31pOkHb1oQ1L5bStq415suKSNgB8Cx5KK1dYA7mXx7am/oIjfRsT7ST8Sf8/LWhJPsfhn1uVx0nf1lohYI/+tFhHbFKYZXlNjr/iZN/ocllS370XSyqQiouL30mydj5PO3tco/K0cEedIWpFUpv8t4K35+5jIou/jcVLRWBnWzLF3abSfngC8jXT2uRrwnjy8lX2k4bbm8WV/N0W9XfbjpDPwYqyrRsQH6i484s6IGEs6ubmedKUHre0fXVr5LWm4HRExJyJOiIhNgA8Bx0sak0dfD2wnaVvSlcrPCvM1Om5r1/UcKSlsU/hMVo+I3ibyppaKpBIRLwFnAt+TdKCkVSQtJ2kU6fKwmVWBFyJirqSdgEPrTHNaPsPbhnTj9aolCHMisIWkQyUNlnQw6XL9xhbnf4ZULtzMyqQdaRaApE+RrlR6lJ+T+XA+aF4n3dxrVPWzp1iuBj4vaYSkNUmVBwCIiKdIN2S/LWm1/D1tKum9hfnXyfMvL+kg0k3HiS2uu7euAD4laVROAl8H7oiIGS3O/1PgQ5L2kTRI0pBc/XUEqSx/RdL3MV/SfqSiwC6X5HWPyZ/D8D5cHQKcKWkFSbuTfnh+UWeaVUk/LC9JWgv4Si+W32xbofzvpugZYG1Jq7c4/STgZUlfkjQ0x7utpB1rJ8yf2SckrR4Rb5DuUXTt+73ZP1r5LWlI0v6SNssnVF0xLACIiLnANTmeSRHxWJ6n2XH7DDBC0gp5GW+SEs53JK2T5x8uaZ/exNmTpSKpAETEN4DjWVRz4RlSscqXSPdXGvkccJakOaQbVFfXmeb3pMvS3wHfiohePwgYEc+TDvQTSJfPJwH7R8RzLS7iEmDrfNl6fYN13Ad8m3Tl9gzppuGfW1z+cjm2J0m1PN5L+mzqOQO4LMfy8Trjf0gqa78b+Cvwy5rx/0L6wb2PVFnhGroXodwBbE46szobODB/fpCqjR+o9PDid1vctoYi4nfAaaQriqdIVw7jejH/46Qr5FNJyeNxUm2d5SJiDvB50j71IulHZkJh3kmkk5TvkG7Y/57Fr2Zb9XRex5Oks9ijI+LvdaY7n3RT/DngL6Tij5Y029Y8SanfTc26/06qMPNw3u+aFofmIpwPke4ZPULa3v8m1Tys5zBgRi4SPJpc1NTL/aOV35JmNgduJiWG24HvRcRthfGXkY7pywvDmh23t5BqxD0tqet35kuk37K/5G29mXTlWhp1L7q2IkkjSTvk8hExv3+jWTZIOpx0s3e3/o5loJC0B/DTiBjR07Q2cEnakFS8tW5EvNzf8TSy1FypmJktrSQtRyqJubKTEwq0OalI2lfSA5Km19Rk6Bq/oqSr8vg78pVB17hT8vAHimV+PS3TzGxpku+XvEx6Lqg398D6RduKv5SqKP6D9EHMJD1Bfkgu9++a5nPAdhFxtKRxpKe/D5a0Nan8dCdSlb6bSQ9I0dMyzcys/7TzSmUn0oN2D0fEPOBKFq/yO5Z08wnSzdoxuebDWNJl3usR8QjpxtJOLS7TzMz6STsbNRtO9wd9ZpKecq07TUTMV3pgaO08/C8183Y9bNTTMgGQdBSpzSxWXnnld265ZV9qalbgybv6voz1t+/zIqY9P63Py9hm7W16nmgAePbROX1exjobrdqn+Z95eHqfY3jrJn1vUHreE6/0eRkrDO/b4xBPPrkkj4d1t/76PT1D27M5c+7p8zJWXfXtfV5Gu02ZMuW5iBjW85TdtTOp1HuYqrasrdE0jYbXu7KqW34XERcDFwOMHj06Jk+e3DjSTnBGq9Xvmy2j79v49sv6vrNPHt/hn/UA8u2D9+/zMk64qtVHoRqbefIf+7yMEefs3qf5zzjjjD7HUMYyfndL359XHbNX5x8jkmpbAGlJO5PKTLo/VT2CxZ/w7ZpmplLTCquT6lo3m7enZQ5MZ8zu7wgAuGd838/CzGzZ1c57KncCm0vaOD/ROY7Cg1/ZBGB87j4QuCW3+TQBGJdrh21MeihoUovLNDOzftK2K5V8j+RY0pPVg4AfRcQ0SWcBkyNiAukp8csldb2rYFyed5qkq0lPXM8Hjulq5KzeMtu1DWZm1jttfftYRExkUZtNXcNOL3TPJb3/o968Z5Oa6OhxmUvijTfeYObMmcydO7evi7JsyJAhjBgxguWXX76/QzGzfrL0vtKyBzNnzmTVVVdl5MiRaEC8wr6zRQTPP/88M2fOZOONN+7vcMysnyyzzbTMnTuXtdde2wmlJJJYe+21feVntoxbZpMK4IRSMn+eZrZMJxUzMyvXMntPpdbIk/+n1OXNOOeDPU4zaNAg3v72t/PGG28wePBgxo8fz3HHHcdyyzXO9TNmzGD//ffn3nvvZerUqTz55JN84AN1X2bXcP6tttqKt71t0SsUJk2axBVXXMHkyZO58MILW16WmVktJ5V+NHToUKZOnQrAs88+y6GHHsrs2bM588wzW5p/6tSpTJ48uVdJBWDTTTdduF4zszK5+KtDrLPOOlx88cVceOGFRAQLFizgxBNPZMcdd2S77bbjBz/4Qbfp582bx+mnn85VV13FqFGjuOqqq5g0aRK77LIL22+/PbvssgsPPPDAEsVy+OGHc8011yzsX2WV1GbTddddx/ve9z4igqeeeootttiCp59+esk32syWOr5S6SCbbLIJb775Js8++yw33HADq6++OnfeeSevv/46u+66K3vvvffCm+ErrLACZ511Vrciq5dffpk//OEPDB48mJtvvplTTz2Va6+9drH1PPTQQ4waNQqAXXfdlYsuuqil+D7ykY9w7bXXctFFF/Gb3/yGM888k3XXXbekrTezpYGTSofper/NTTfdxN/+9reFVwyzZ8/mwQcfZIsttmg47+zZsxk/fjwPPvggknjjjTfqTteX4q8LLriAbbfdlp133plDDjlkiZZhZksvJ5UO8vDDDzNo0CDWWWcdIoILLriAffbZp9s0M2bMaDj/aaedxp577sl1113HjBkz2GOPPZYojsGDB/Pmm28CKcnNmzdv4bgnnniC5ZZbjmeeeYY333yzaaUCM1v2+BehQ8yaNYujjz6aY489Fknss88+fP/73194tfGPf/yDV199tds8q666KnPmLHrnx+zZsxk+PL125tJLL13iWEaOHMmUKVMAuOGGGxbGMH/+fD71qU9xxRVXsNVWW3Heeect8TrMbOnkK5WslSrAZXvttdcYNWrUwirFhx12GMcffzwAn/nMZ5gxYwY77LADEcGwYcO4/vrru82/5557cs455zBq1ChOOeUUTjrpJMaPH895553HXnvttcRxHXnkkYwdO5addtqJMWPGsPLKKwPw9a9/nd13353dd9+dUaNGseOOO/LBD36Qrbbaask/BDNbqrTtHfWdpN5Luu6//37/GLaBP9dy+SVdiyxdL+l6qM/LaDdJUyJidG/nc/GXmZmVxknFzMxKs0wnlWWh6K9K/jzNbJlNKkOGDOH555/3D2FJut6nMmTIkP4Oxcz60TJb+2vEiBHMnDmTWbNm9XcoS42uNz+a2bJrmU0qyy+/vN9QaGZWsmU2qZiZLYmBUB24Py2z91TMzKx8vlIxswGhjAcXrf18pWJmZqVxUjEzs9I4qZiZWWmcVMzMrDROKmZmVhonFTMzK42TipmZlcZJxczMStPWpCJpX0kPSJou6eQ641eUdFUef4ekkYVxp+ThD0jaJw/bQNKtku6XNE3SF9oZv5mZ9U7bkoqkQcBFwH7A1sAhkraumewI4MWI2Az4DnBunndrYBywDbAv8L28vPnACRGxFbAzcEydZZqZWT9p55XKTsD0iHg4IuYBVwJja6YZC1yWu68BxkhSHn5lRLweEY8A04GdIuKpiPgrQETMAe4HhrdxG8zMrBfamVSGA48X+meyeAJYOE1EzAdmA2u3Mm8uKtseuKPeyiUdJWmypMl+Z4qZWTXamVRUZ1jtaxYbTdN0XkmrANcCx0XEy/VWHhEXR8ToiBg9bNiwFkM2M7O+aGdSmQlsUOgfATzZaBpJg4HVgReazStpeVJC+VlE/LItkZuZ2RJpZ1K5E9hc0saSViDdeJ9QM80EYHzuPhC4JdJL4ycA43LtsI2BzYFJ+X7LJcD9EXFeG2M3M7Ml0Lb3qUTEfEnHAr8FBgE/iohpks4CJkfEBFKCuFzSdNIVyrg87zRJVwP3kWp8HRMRCyTtBhwG3CNpal7VqRExsV3bYWZmrWvrS7ryj/3EmmGnF7rnAgc1mPds4OyaYX+i/v0WMzPrAH6i3szMSuOkYmZmpXFSMTOz0jipmJlZaZxUzMysNE4qZmZWGicVMzMrjZOKmZmVxknFzMxK46RiZmalcVIxM7PSOKmYmVlpnFTMzKw0TipmZlYaJxUzMyuNk4qZmZXGScXMzErjpGJmZqVxUjEzs9I4qZiZWWmcVMzMrDROKmZmVhonFTMzK42TipmZlcZJxczMSuOkYmZmpXFSMTOz0jipmJlZaZxUzMysNE4qZmZWmrYmFUn7SnpA0nRJJ9cZv6Kkq/L4OySNLIw7JQ9/QNI+NfMNknSXpBvbGb+ZmfVO25KKpEHARcB+wNbAIZK2rpnsCODFiNgM+A5wbp53a2AcsA2wL/C9vLwuXwDub1fsZma2ZFpKKpLWlbSTpF26/lqYbSdgekQ8HBHzgCuBsTXTjAUuy93XAGMkKQ+/MiJej4hHgOl5eUgaAXwQ+O9WYjczs+oM7mkCSV8HPgn8HViQBwfwgR5mHQ48XuifCbyr0TQRMV/SbGDtPPwvNfMOz93nAycBq/YUu5mZVavHpAJ8DNgiIub2ctmqMyxanKbucEn7A89GxBRJezRduXQUcBTAhhtu2HO0ZmbWZ60Ufz3S4nS1ZgIbFPpHAE82mkbSYGB14IUm8+4KfFjSDFJx2l6Sflpv5RFxcUSMjojRw4YNW4LwzcystxpeqUj6DumqYQ5wl6Sbgde7xkfE8T0s+05gc0kbA0+QbrwfWjPNBGA8cDtwIHBLRISkCcAVks4D1gc2ByZFxO3AKTm+PYAvRsQnW9xWMzNrs2bFX/fm/9OA39SMqy3GWky+R3Is8FtgEPCjiJgm6SxgckRMAC4BLpc0nXSFMi7PO03S1cB9wHzgmIhYUHdFZmbWMRomlYi4BEDSsRFxYXFcThY9ioiJwMSaYacXuucCBzWY92zg7CbLvg24rZU4zMysGq3cK/l0nWFHlB2ImZkNfM3uqRxMKo7aWNIvC6NWA15qd2BmZjbwNLunMgl4nlTz6qLC8DnAXe0MyszMBqZm91QeIVUnvlnSMGBH0g36RyLijYriMzOzAaTHeyqSDgf+SqoO/ElgsqTxbY7LzMwGoFaeqD8F2CEiZgFIegvwZxa12WVmZga0VvvrCbrfmJ9NeuLdzMysm1auVB4Dbpd0PemeygHAnZI+DxAR321jfGZmNoC0klQez38r5v6up+vdoJaZmXXTY1KJiNMgvaUxIl7vaXozM1t2tVL7aydJ9wAP5v53SLqg7ZGZmdmA08qN+u8C+5MehCQi7gb2bGdQZmY2MLWSVJaLiEdrhrnFYDMzW0xLN+ol7UR68+Ig4N+Af7Q3LDMzG4hauVL5LHA8sCHwDLBzHmZmZtZNK7W/niW/PMvMzKyZplcqknaXdLWku/PflZJ2qyo4MzMbWBomFUn7AZcD/0t6UdcRwC3ATyTtU014ZmY2kDQr/joJ+EhEFN+dMlnSJOB80rvnzczMFmpW/LV+TUIBICKmAuu2LyQzMxuomiWVV5qMe7XsQMzMbOBrVvy1ac276bsI2KRN8ZiZ2QDWLKl8rMm4C8sOxMzMBr5m76j/XZWBmJnZwNfKE/VmZmYtcVIxM7PStPI+lY+2MszMzKyVK5Uv1xn272UHYmZmA1/DG/W5KZZ9geGSziuMWg14s92BmZnZwNOsSvGzwL3AXGBaYfgc4OR2BmVmZgNTsyrFdwF3SfppRLxeYUxmZjZANWul+Oe58y+S/lr718rCJe0r6QFJ0yUtdnUjaUVJV+Xxd0gaWRh3Sh7+QLFVZElrSLpG0t8l3S/p3S1vrZmZtVWz4q8T8/8Dl2TB+dXDFwHvB2YCd0qaEBH3FSY7AngxIjaTNA44FzhY0takF4NtA6wP3Cxpi4hYAPwn8JuIOFDSCsBKSxKfmZmVr+GVSkTMzJ1PAg9HxEO5/23Aoy0seydgekQ8HBHzgCuBsTXTjAUuy93XAGMkKQ+/MiJej4hHgOnATpJWA94DXJJjnBcRL7UQi5mZVaCVKsV/BIZKWg/4Pen99D9qYb7hwOOF/pl5WN1pImI+MBtYu8m8mwCzgB9LukvSf0taud7KJR0labKkybNmzWohXDMz66tWkspyEfFPUgOTF0bEh4DtWphPdYZFi9M0Gj4Y2AH4fkRsT2qCv25NtIi4OCJGR8ToYcOGtRCumZn1VUtJRdKOwKHAjXnYoBbmmwlsUOgfQSpKqzuNpMHA6sALTeadCcyMiDvy8GtIScbMzDpAK0nleOBM4H8i4l5Jm5CKxHpyJ7C5pI3zDfVxwISaaSYA43P3gcAtERF5+LhcO2xjYHNgUkQ8DTwu6W15njHAfZiZWUdoVvsLgIi4Bbil0P8w8LkW5psv6VjSu+wHAT+KiGmSzgImR8QE0g33yyVNJ12hjMvzTpN0NSlhzAeOyTW/AP4N+FlOVA8Dn2p5a83MrK16TCqSNiNdrYwsTh8Re/c0b0RMBCbWDDu90D0XOKjBvGcDZ9cZPhUY3dO6zcysej0mFdJ9i0uAnwILepjWzMyWYa0klTcj4oK2R2JmZgNeKzfqb8jPfAyTtFrXX9sjMzOzAaeVK5XP5P+nFYYFsGH54ZiZ2UDWSu2vDXqaxszMDFq7UkHSlsDWwJCuYRFxRbuCMjOzgamVKsVfBvYGtiQ9c7IP8CfAScXMzLpp5Ub9wcCewFMRcRjwDlq8wjEzs2VLK0nltfw0+3xJqwJPk1oLNjMz66aVK467JK1Bau5+MvAy0NKbH83MbNnSNKnkF2adkV+EdZGk3wKrRYSTipmZLaZp8VduMfjGQv90JxQzM2uklXsqkyT5nSVmZtajhsVfkgbnV/zuBhwp6SHSmxZFuohxojEzs26a3VOZRHqr4gEVxWJmZgNcs6QigIh4qKJYzMxsgGuWVIZJOr7RyIg4rw3xmJnZANYsqQwCViFfsZiZmfWkWVJ5KiLOqiwSMzMb8JpVKfYVipmZ9UqzpDKmsijMzGyp0DCpRMQLVQZiZmYDXytP1JuZmbXE70Uxsx6NOGf3/g7BBghfqZiZWWmcVMzMrDROKmZmVhonFTMzK42TipmZlcZJxczMStPWpCJpX0kPSJou6eQ641eUdFUef4ekkYVxp+ThD0japzD8/0maJuleST+XNKSd22BmZq1rW1KRNAi4CNgP2Bo4RNLWNZMdAbwYEZsB3wHOzfNuDYwDtgH2Bb4naZCk4cDngdERsS2pJeVx7doGMzPrnXZeqewETI+IhyNiHnAlMLZmmrHAZbn7GmCMJOXhV0bE6xHxCDA9Lw/SA5tDJQ0GVgKebOM2mJlZL7QzqQwHHi/0z8zD6k4TEfOB2cDajeaNiCeAbwGPAU8BsyPipnorl3SUpMmSJs+aNauEzTEzs560M6nUazo/Wpym7nBJa5KuYjYG1gdWlvTJeiuPiIsjYnREjB42bFgvwjYzsyXVzqQyE9ig0D+CxYuqFk6Ti7NWB15oMu/7gEciYlZEvAH8EtilLdGbmVmvtTOp3AlsLmljSSuQbqhPqJlmAjA+dx8I3BIRkYePy7XDNgY2ByaRir12lrRSvvcyBri/jdtgZma90LZWiiNivqRjgd+Samn9KCKmSToLmBwRE4BLgMslTSddoYzL806TdDVwHzAfOCYiFgB3SLoG+Gsefhdwcbu2wczMeqetTd9HxERgYs2w0wvdc4GDGsx7NnB2neFfAb5SbqRmZlYGP1FvZmalcVIxM7PSOKmYmVlpnFTMzKw0TipmZlYaJxUzMyuNk4qZmZXGScXMzErjpGJmZqVxUjEzs9I4qZiZWWmcVMzMrDROKmZmVhonFTMzK42TipmZlcZJxczMSuOkYmZmpXFSMTOz0jipmJlZaZxUzMysNE4qZmZWGicVMzMrjZOKmZmVxknFzMxK46RiZmalcVIxM7PSOKmYmVlpnFTMzKw0TipmZlYaJxUzMytNW5OKpH0lPSBpuqST64xfUdJVefwdkkbm4WtLulXSK5IurJnnnZLuyfN8V5LauQ1mZta6tiUVSYOAi4D9gK2BQyRtXTPZEcCLEbEZ8B3g3Dx8LnAa8MU6i/4+cBSwef7bt/zozcxsSbTzSmUnYHpEPBwR84ArgbE104wFLsvd1wBLxz7oAAAUa0lEQVRjJCkiXo2IP5GSy0KS1gNWi4jbIyKAnwAHtHEbzMysF9qZVIYDjxf6Z+ZhdaeJiPnAbGDtHpY5s4dlAiDpKEmTJU2eNWtWL0M3M7Ml0c6kUu9eRyzBNEs0fURcHBGjI2L0sGHDmizSzMzK0s6kMhPYoNA/Aniy0TSSBgOrAy/0sMwRPSzTzMz6STuTyp3A5pI2lrQCMA6YUDPNBGB87j4QuCXfK6krIp4C5kjaOdf6+hfghvJDNzOzJTG4XQuOiPmSjgV+CwwCfhQR0ySdBUyOiAnAJcDlkqaTrlDGdc0vaQawGrCCpAOAvSPiPuCzwKXAUODX+c/MzDpA25IKQERMBCbWDDu90D0XOKjBvCMbDJ8MbFtelGZmVhY/UW9mZqVxUjEzs9I4qZiZWWmcVMzMrDROKmZmVpq21v4ys7454aob+zsEs17xlYqZmZXGScXMzErjpGJmZqVxUjEzs9I4qZiZWWmcVMzMrDROKmZmVhonFTMzK42TipmZlcZJxczMSqMmb+9dakiaBTzah0W8BXiupHD6ohPi6IQYoDPi6IQYoDPi6IQYoDPi6IQYoO9xbBQRw3o70zKRVPpK0uSIGO04OiOGTomjE2LolDg6IYZOiaMTYujPOFz8ZWZmpXFSMTOz0jiptObi/g4g64Q4OiEG6Iw4OiEG6Iw4OiEG6Iw4OiEG6Kc4fE/FzMxK4ysVMzMrjZOKmZmVxknFzMxK46RiZmalcVIxM7PSOKmYmVlpnFQKJE2UNLID4viQpI0K/adLulvSBEkbVxTDr/L66v5VEUOnkrS8pO0lrdMP6x4k6S2F/hUkHSXp/qpj6U8dcozsVejeuGbcR6uIoWadQyRtK2kbSUOqXv/COPycyiKSPg58DbgM+EZEvNFPcfwN2Dki/ilpf+A84BBge+CgiNinghje29UJ/BD4THF8RPy+3THkOG4FGu2kERFjKojhv4ALImKapNWB24EFwFrAFyPi5+2OIccxDvgB8CrwIHAGcDlwJ/DViPhrBTHMYdH3ofw/gMHAChExuN0x5Dg64Rj5a0TsUNtdr7/NcQwGvg58mtRw7nLACODHwL9X/TtWyQ4wUETE1ZL+BzgdmCzpcuDNwvjzqgsl/pm7PwpcEhFTgCmSPldRAAuThqRXqkoidXyxzrCdgZOAZyuKYfeIODp3fwr4R0QcIGld4NdAJUkF+DLwzoiYLmkHUnIbFxHXVbR+ImLVYr+kVYHPAf8KVBYHHXCMsCip1nbX62+nbwKrAhtHxBwASasB38p/X6gwFieVOt4gnQmuSPqi3mw+eVtI0irAP4ExwPcK4/rjsrbfLmfzDwWw8OrpNNJ3c3RE/LqiMOYVut8P/CLH9rRU5W8H8yJiel73XyU9UmVCKZK0BnAc8C/AFcCOEfF8tSH0+zESDbrr9bfT/sAWUSh2ioiXJX0W+DtOKv1H0r6ky+gJwA6FM6GqnQ9MBV4G7o+IyTm+7YGnqghA0lqF3kGS1qRw9hURL1QRR45lH1IymQucHRG3VrXu7KVcxPIEsCtwRI5rMDC0wjjWkXR8oX+VYn8VV9L5fs4JwMHAj4DtI2J2u9dbR78fI8Am+f6iCt3k/kru62RRTCiFgQskVX5C6HsqBZL+SDoDntYBsQwH1gHujog387D1gOUj4rEK1v8I6Wyr3ql4RMQm7Y4hx3EnMIx0iX97nUCquI+wBfBdYF3g/Ii4NA/fB9g7Ik5odwx5fV9pNj4izqwghleBWaTy+jl1YqiqiLgTjpH3Nhtf4X3H64FfRsRPaoZ/Evh4RHy4ijgWrtdJpWf5Uv+YiDi7ovXV3uAL4LmIeLyK9ecYNoqIvrwts6w4bqP5jfq9GoyrhKSVI+LV/oyhSpLOoEnRThWJLcfR78dIp8jJ9ZfAa8AU0mexI+kq+iMR8USl8TipLCJpA9LN0OHA9aSy4q8ChwE/j4hKyiZzjadaawErAIdExNQKYqis9spAkA/c9YC/RcS8XJ34OODwiFi/ohiujoiP5+5zI+JLhXE3RcTeVcTRSJUJtkOOkX6vmZjjGBwR83MV521IpQvTIuJ3Vay/lu+pdPcT4PekrL8v8BdgGrBdRDxdVRARsWe94ZJGk4ph3lNBGJXegW4m/4AfQzpgArgPuCgiKqn9Jek44N+B6cCKkv6TdO/tJ8A7q4gh27zQ/X7gS4X+Xr9LfEk1S7BAJQm2Q46RTqiZCDCJdA/4FuCWCtdbl5NKd2tFxBm5+7eSniHVanm9H2NaKCIm5xovVRgu6btNYvl8FUFI2pV0xXgp6UdcwA7AJEmfiIg/VxDGUcDbIuIFSRuSkst7IuIvFay7qFmxQiVFDh2UYOuq8hjpkJqJ0EEngOCkspiaWk5PAytJWhmqrfFUj6S3Ul1Vxa7y2f72beCAiLirMOwGSdeRHgR8VwUxzO367iPiMUn/6IeEAmlf3J70cNvQ3K38V1UttE5JsHVVfIx0Qs1EgGE1tQK7qbLyBDip1Fqd9ENazPxdtYsCqKrG0wUsfmCsBexCdXXOn4+IyypaVzOr1SQUACJian7wrgojaq7a1in2V3XVRqoq2/UD8XShu6u/Ch2RYDvhGKlXM7FYgaCKmonZIGAVOuSKxUmlICJG9ncM2eSa/gCeB46v6j4C3R/460+StGZEvFgzcC2qa7vuxJr+/rqCO6UDrgg6JcF2wjHyKvAKcGD+q42nqpqJT0XEWRWtq0eu/VXgaoo9k/Q2UntXR1a0vqOAI0k3RbvO/N4JnAv8KCJ+UEUcDWIbAnwoIn5R0fr6vUaepPHNxnfI1e0yRdJdEbF9f8fRxUmloBOqKRbi6O9GFLcjtRu0Pql69QWkpjDeBXw7Ir7T7hgKsexPqlGzTR40DfhmRPyqqhgKsQwC9iY1XrgP8MeIqD1Lbde6O+rHoz91yDFS2xJxAM8BU7va4KqCpGGkBk6LcbxU7yn7SuJxUulZrqZ4XkRUUU0RSfVq0SysqhgRO1YQwx3A90llxfvmdV8BnBYRc9u9/k4j6T3AocAHSVU4dwU2qbIpH0kvAX9oNL6KJ6cl/YrF27x6Drg1In7a7vUX4uiEY+THdQavBWwHHJGr+LZdg9YvVgHuBj4TETOqiGNhPE4qremvooeaqopfr6qqoqSpETGq0P84MDIiFjSZrV2x7AecTPfnVM6NiIkVrX8m8BgpyV4fEXNyY45Vtu+EpAepeQVBURXNgjRommQt4JPAgxFxcrtjqNVfx0iTeDYCro6IKmomNovjo8BREbFvlev1jfoWVF1NMa+zv6sqDilUWYV0Q3I7KTXLW1XNFklHkppVP4lFN2dHA+dIGhERF1cQxrXAAaRGFBdIuoH+abm5P19BADROXLkxxSmk5F+JDjhG6oqIRyUt3wFx/FLSl6ter69UCnqqplhVGX6HNKLY7ACtrM0tSfcBu9U+IyRpbeBPEbFVRXEI2JN0L+UDwGqk1oonRsQrFcXwy4io/I2Craq9um3zuvr9GGkkV2a5NCLe3V8x5DhWIR0jlXwnC9frpLJInZotXdUU76ywmuJAaERx+ajobXKS7m+UOJqNa3NMywP7AeNIrRS/pYdZylrvx2jemOMvK4hhrTqD1yS9V2WziPhEu2PIcdzGos+i9n5CJcdInftLkE5C1wM+GRGLJbs2xVHvwcc1gQ+TmjOq4mp+IRd/dTcRGBYR9xUHKr3zOSJiVhVBRMQeVaynNwpn6ocCHwLeWtGqX5b0joi4uyaed1Cn6fUq5IQ6AZggqcr3qezfZFyQ2qxrt65WcIuvEn4OuA34bAXrTyttcoxUWPT0rZr+rpPQByOiyue8ah8CDtLDsJ+MiHsqjAPwlUo3kq4Evl9bbpzLbsdHxKEVxXFSRHwjdx9UfA5C0tcj4tQq4sjrexcpkXyEdBZ2DDCh9mHENq5/N+BnpPd3FJv1Hk86aP5UQQx/azY+IrZrdwydQtK7qzoD743ak56IaPtJjzqgZeieSPpWRNRr+LJ963RSWUTStIjYpsG4eyNi24riWFjTrLbWWVW10CSdDXycVOvp56T3j0+uusZTjmVd0nvQFzbrTbqsr6RpEklTScnsCuBXpHbRFoqK3juTE/zFwKbAPcCnI+L+KtZdiKHfH8As6s+Tnk77LOqR9FhEbFjlOl381V2zy+Yqa3OoQXe9/nY5CniAVI32xoiYq354NSmkd8EDp/fHuvP6R0naknST/gpSleYrgJsiYn6FoVxEalngD6Ty8vNJD2BWqSPal6pz0nMW6aSnyif616jzAORCVdzjakHl35eTSncPSvpA7fMP+TmJhyuMo/bhskbj2mldFj05fn6uDTZU+YVAFcXQrOhJpBuylRQ9RcTfga8AX5F0MKmp93NJtY+qslxE/G/u/oWkUypcd5eNtehd7Iup4gHMrBNOelYn3eeq+8ptqrnH1ajyBCxqwbpSTird/T/gRkkfZ1GjgaOBd9P8JmnZ3iHpZdIOMTR3k/uHVBFAfsjx18CvcxtX+wMrAU9I+l1V95eAN2lS9FQVpRdTjSMVs7xI2leuqziM2jPjbv0VnRnPIr2OoL91wknPoxHx6YrW1Uxt5YmiyhuG9T2VGpJWJJXRdt0/mQZcsSw2TVKPpNVI772urJihUPT0Ifqh6EnS70k1bK4GrgG6PTNT+wxNG+Oo1yxIIYz2/8B14n2EwknPIcBuQCUnPW6LrT4nlRYoNSI4LiJ+VtH6biT9cF5fZdtSNTGcHxHH5e4vRMR/FsZdGhGH91NcB5PuLZwbEZUUPUmaQfdnIhaOIv2YV/KenU7Q7AHMKp9fakTpHTsfi4hLK1jXdqR3LG0G3BMRv233OhvEsQ5wao7jb8A5EfFy87naGI+TyiL5LPwYYDhwA3Bz7j+R1PLo2IriGEsqankf6Z3TPyc9uV3ZpWwn1EArrK+26Olq4LqqnmTvFHUecut6RuRPEfFIP4TUL1V5C+seBKwZEc/l/hVIVc2Pr+KhWEnfI9VI/D9gDPCriPhqu9dbJ47fkIrA/kC6Ylu1v076wEmlm9ym04ukZh/GkJ5KXYHUREslzd7XxDOUVMtnHOm+zkTg54Wbte1c98JL+9rL/CqTSqcUPXUCSV+pM3gtUg2wMyLiygpj6e/nl8aRXif9KvAgcAZwOXAn8NWKmjK6F3hHRCyQtBLpNQj1Wk9udxy1jb/2axGlk0qBpHsi4u25exDpLHDDqPDdCI3kS+3LgO0iYlAF67sb2IP0dsVbcnfXjcBbI+Id7Y4hxzEDFz01lWv/3LwsPb+Uf9APiIjpSi/Xu51URF1Z5Yn+voIvrLfrWF14fBb7qz7xcu2v7haWB+ezj0f6M6EotY78cdKVynrAL4BPVbT61UmX1F07avHMr7IzkeicVzx3rIh4IRdDVaETqvICzIuI6ZAaj8zHatW18bYsVHkXsGnur7S6O4sfq7DoeA3SfZ/KOKl011WVF7pX5+3aSVarIgil5t4PAd5Gqut+UkT8uYp1d/GP+cAhaS9SsW0VOqEqL8A6NfeYVin2R8R5FcRQeWOm9XTaseqkUlBFsVKLdgHOIRVpvNnfwVhnkHQP9VvFfZLUSnDbNXh+aWWqf37ph3RvSLG2vwpD80OxSFoxIl7vGiFpZ6CS5ns6je+pFEjaK/IrQCVtXKxRI+mjHdLsgi2jlN4oWBTA8xHxan/EU5Sr8n604mZS+lUn1ZDsJL5S6e5bQNeOcG2hG+DLVNTsglkD9e7vrZgf2H29iuQi6UPA37oa0ZR0OvAx0ln5F9q9/kIc360Z1FW9+taooOXqrjAadNfrX2Y4qXTnncQ6WbPmOAbne/Unt/kh3bOBnQEk7U96N/0hwPbAf1FdA5dT6gxbC/impKsi4vwKYuiENvo6jpNKd95JrGP1VG1X0jDg96T3z7QxjIWtPHwUuCQipgBTJH2ujeutDaJuMZuk/yI9jFhFUhmRr5hU6Cb3D69g/R3JSaW7TXILrCp0k/srf4+IWW9ExCxJX2rzaqT07vN/kh4Q/l5hXCWNnTYTEa9VV7uaEwvdk2vG1fYvM5xUuis2w1L7qtDafrOOExG/avMqzgemAi8D90fEZABJ2wNPtXndTUkaDBwGzKxifctSpYTecO0vM+uV3BbbOsDdXVXeJa0HLB8Rj1UUwxwWL5J+jVT8d1xEPFlBDD+uE0OXiIgj2h1DJ3JSKZDfRW5mLZL0sTqDNwSOAwZFxIiKQ+oILv7qriNeCGVmnS8iru3qlrQJqfn595AeXL6kv+Lqb75SqdHfL4Qys4FD0lbAv5OqVH8T+Omy/lvhpNJEf7wQyswGBkm/IL1u/FukVzMsKI5fll7LUOSkUsMvhDKzVjR4LUNXfeZl9rUMTioFfiGUmVnfOKkU+IVQZmZ946RiZmalWa6/AzAzs6WHk4qZmZXGScXMzErjpGJmZqVxUjEzs9I4qZiZWWmcVMz6SNICSVMl3SvpF5JWysNLb4VB0gGSti57uWZlcVIx67vXImJURGwLzAOObuO6DgCcVKxjOamYleuPwGbFAZJWkfQ7SX+VdI+ksXn4SEn3S/qhpGmSbpI0NI/bVNJvJE2R9EdJW0raBfgw8M18ZbRp5Vtn1gMnFbOS5NfZ7gfcUzNqLvCRiNgB2BP4tha9SH1z4KKI2AZ4Ceh68dPFwL9FxDuBLwLfi4j/AyYAJ+Yro4fau0VmveeXdJn13VBJU3P3H1n8BU0Cvi7pPaQXwQ0H3prHPRIRXfNOAUZKWgXYBfjFotzDiu0K3qxMTipmffdaRIxqMv4TwDDgnRHxRm64dEge93phugXAUFIJwks9LNOsI7n4y6z9VgeezQllT2CjZhNHxMvAI5IOAlDyjjx6Dun1DGYdyUnFrP1+BoyWNJl01fL3Fub5BHCEpLuBacDYPPxK4ERJd/lGvXUiN31vZmal8ZWKmZmVxknFzMxK46RiZmalcVIxM7PSOKmYmVlpnFTMzKw0TipmZlaa/w9aJLacs7tqAQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data.plot(kind='bar',x='Planet',y='Delta Flux')\n", "plt.title(\"Graph of transit depth for each planet in the solar system\")\n", "plt.xlabel(\"Planet\")\n", "plt.ylabel(\"Transit Depth\")\n", "plt.gca().invert_yaxis()\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc4AAAEWCAYAAADvi3fyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3XmcHFW5//HPd5KQhD2GqJBBgySiiBhhBBUXFMHAVYKClyhqUATlynW7XhaVn4g7915xQ70IyCYCgkjUq8gi4MKSCXtAJWxmCEtIAiQIIck8vz/Oaag03T1dk+npmcz3/Xr1a7qqTlU9VXP6PLWXIgIzMzNrTke7AzAzMxtOnDjNzMxKcOI0MzMrwYnTzMysBCdOMzOzEpw4zczMShiUxClpd0k96zD+jyQdO5Ax5ekeJ+nsgZ7ucCEpJE1tdxwDQdKVkj6Svx8k6fftjqmRYrxDiaQ3SvrbOk6jZfVK0r2S3jbA0/ycpFMaDD9Y0p8aDP+tpNkDGVM7jfR2sRlNJ85cYZ+UtELSg5JOl7TxQAdUq5JGxMci4ssDPa914crVHElTckM6erDmGRE/jYi9Bmt+A62ddSsi/hgR2xViGfBE1azcxnyl1fOJiK9FRGWjq3R9jYi9I+KM1kVoRe2skxVl9zjfGREbA9OBVwPHDHxINtINZpIdjpT4NIs9Y6j8ZiSNancMgyIimvoA9wJvK3SfAPym0D0W+G/gH8BDwI+A8XnY7kBPoezRwF3AcuB24F25/8uBp4A1wArg0dz/dOArhfEPBRYAS4E5wFaFYQF8DLgTWAacBKjOMh0HXACcl2O5AXhVYfhWwIXAYuAe4BO5/wzgaWBVjvNm4C3ArYVxLwOuL3T/Cdiv0XTzsI7C+lkCnA88Lw+bkpdvdl7PjwCfb/A/Oz3/Hy7Ny3cV8OKqdTU1f/8X4EbgcWAhcFyhXMP59hHzP/K4K/LndQ3+D2fn+X8E2AW4BngUeAD4PrBBYZw9gb8Cj+VhVwEfycMOBv5UFfvowrhXFspOzeM+lpfrvAbr8+fAg7ns1cArqtb1ScBv8rq+Dti2mXir5vGculWI+avAn4Enc9wfAu7I87sb+GhhOrsDPcB/AA/ndfihwvB9SL+95cD9wGerf6vAWUBvnt8K4Mg66+U/8/QXAR9m7XrVZ7sAfC6v+3uBg/Kww/I6eDrP+1eFduizwC15XZ4HjKsT133Azvn7+3Nc2+fujwC/LNS/s+vVV3J9ysuxjPSb3btOfWpYtkaML8/jPwrMB/atNd3qel34/X6c1NbdU2Pa40i/qSV5+nOBFxTaoDmkNnQBcGjV7/HsEvX+h8D/AU9QyBFVcd9Nqmv3AAflerEUeGWh3PNJdW0SsAXw6xz3UuCPpHamZp0EXgv8JZe/Gdi9aj1+JQ9fAfwKmAj8lNTezAWm1Psf1fy/NV2wkDiBTuBW4DuF4d/O/4jnAZvk4L5e/WPM3e/J/7gO4MC8wresVTkK/5yv5O9vJf3Idsor/3vA1VWV6dfA5sCLSMlpRp1lOo704zwAGEP6Qd6Tv3cA84D/B2wAvCT/899ep3KNy//MLYDRuaItyutifB42sYnpfgq4Nq/jscD/Aj+rSgI/ztN8FbASeHmd5TudVFnflKf1HZ77w5ta+B+9Mse3I6mR26+Z+TYZ8+haMVb9H/bL8x8P7Ez6MYzO07gD+FQuvwWpwlf+b58GVtO/xPkz4PN5vuOANzSI88P5/zmWVN9vqlrXS0kJfzTpR3luM/HWWR9nV/W7ktSovyJPfwxpY2dbQMCbgX8COxX+n6uB43PZffLwCXn4A8Ab8/cJVeMVf6v3UqMxLAyfkevKDsBGwDmsXa/6ahdWA9/K6/TNpLZgu+rffVU815Paj+flevGxOrGdCfxH/n4yacPu8MKwT1ev7zr15WBS/TwUGAUcTvptq0Z9ali2Kr4xpKT1OVJb8FbS73W76unWahtznJfm9TC+xvQ/mtf3hjmWnYFN87CrgB+Q6vx0Uju5R522ra96/xiwG/k3VBXDRqS6X1mmLcmJN8//m4Wyn+TZDaSvkzayxuTPGwvr+17W3ombTNo42CfHsGfunlRYjwtIv5XNSBuMfwfeRvotnQn8pF4dr1m3mi6Ygl2R/7EBXA5snoeJVOGLW9ivI28FUfVjrDHtm4CZtSpH9Q8IOBU4oTBsY1JFnVKoTG8oDD8fOLpBA3VtobuD3KAAuwL/qCp/TGUFV1eu3O+PwLtJDf7v87xnkPZGb8ll+pruHeQKXKhoq3g2gQTQWRh+PTCrzvKdTm68C+tqDbB1YV1NrTPut4ETqxqTmvNtMua+EufV9YbnMp8CLsrfP1j1fxNpz6U/ifNMUqPa2Wj+NeLZPE93s8K6PqUwfB/gr83EW2d91Eqcx/cR0y+BTxZ+c09WLffDwGvz93+QGtZNq6axO+US52nANwrdL63UK5prF1YDG1X9Xo+t/t1XxfP+QvcJwI/qxHYIMKdQRz/Csxsz9/HsxsIz67tOfTkYWFDo3jCXeWGN+tSwbFV8byRtYHcU+v2MfLSH5hLnWxv8bz5M2svasar/1qR2YJNCv68Dp9erf33U+zMbxLARaS9wf6qSO6ktXFhZfqAb+Nf8/XjgYmq0Tzw3cR4FnFVV5hJgdmE9Fo+Q/Q/w20L3OylsDDTzKXueZL+I2IRU4V9G2pKGtGu9ITBP0qOSHgV+l/s/h6QPSrqpUHaHwrT6shWp0gMQEStIWxeTC2UeLHz/Jylh1LOwMK1eUoO2FfBiYKtKjDnOzwEvaDCtq0jr5k35+5Wkreg3526amO6LgYsKw+4gVfLifPu7fCtIe0VbVReStKukP0haLOkx0uHu6v9Jvfk2E3NfFhY7JL1U0q/zhWiPA18rxLNV1XJF9fglHElq4K+XNF/Sh2sVkjRK0jck3ZXjuTcPKq6jeutnoOKtXkd7S7pW0tK83vepimdJRKyuE9P+ufx9kq6S9Lp+xANVy0bht0lz7cKyiHiiavzn1M8qzdb/q4A3SnohaY/rPGA3SVNIex439TGfmvOMiH/mr/Xm22zZrYCFud2puI+127K+NKpHZ5ESyLmSFkk6QdKYPN+lEbG8r/k2We/rxpD/tweS2pMHJP1G0svysOtIG1Zvzv2mko5OAPwXaS/x95LulnR0g+V8MfCeqjb1DaQN+IqHCt+frNFd6kLXfl1gEBFXkbY0/jv3eiTP/BURsXn+bBbpQqK1SHox6ZDfEcDEiNgcuI3UeEHammlkEWlFVaa3EekQ6P39WRbS1ldlWh2kw42LSJXhnsLybB4Rm0TEPg3irE6cV/HcxNnXdBeSzokUh4+LiIFYvo1Jh3UW1Sh3DqnSbh0Rm5EOk6hGuVoaxdzX/7OiutwPSecEp0XEpqSNi0o8D1Qtl4rdVSqN8oaFfi98ZqYRD0bEoRGxFWkP7Ad1bqV4HzCTdHhnM9KeCTS3jsrEC/XX2TP9JY0lnSf/b9J5q81J55ma+p9FxNyImEk6r/RL0p5emVgq1lo20umRimbahQn5N1wcv1I/m607NUXEAlJi/QTpiMZyUlI7jLTn1ltrtHWZZ0mLgK2rLvR6Ec+2ZU9Qp94W1I03IlZFxJciYnvg9cA7SEc/FgHPk7RJnfkWNVPvG66ziLgkIvYkJbK/ktr/ijNI558/AFwQEU/lcZZHxH9ExEtIe4SfkbRHnfktJO1xFtufjSLiG43iWhfrcmXet4E9JU3PFfDHwImSng8gabKkt9cYbyPSgi/O5T5E2uOseAjolLRBnfmeA3xI0vTceHwNuC4i7u3ncuws6d35qrRPkc7dXUs6FPm4pKMkjc9bXjtIek0hzilVlf4vwHak81zXR8R8UpLflXRSnSam+yPgq3kDA0mTJM3s57IB7CPpDXl9fpm0rmptIW5C2gp9StIupB9MsxrFvJh0Mv8lJePehHRuZEXeGj28MOw3wCsK/7dPULtRISIWkxqE9+d1/WHSuQ5yrO+R1Jk7l5Hq5po68awkHd3YkFTvmtV0vFmtulVtA9I5p8XAakl7A03dgiNpA6V7XTeLiFWk9VxrmSuxNPrfnQ8cLGl7SRsCX6wMKNEufCnH9EZS4/7zJufdjKtIG+mVDdcrq7qr9be+9kdlj+tISWMk7U5KEufm4TcB75a0Yd6YO6TMxCW9RdIrla50fZx0+mRN/v3/Bfi6pHGSdszT/mmNyaxLvUfSCyTtmzeOVpJO9xXr2lnAu0jJ88zCeO+QNDVvZFbqZ2W86npxNvBOSW/Pv/FxSs8O6KRF+p04c4N0JlB5MMFRpF3ra/Mu/WWkJFI93u2kY8zXkFbAK0lXClZcQbq67EFJj9QY//I8zwtJW7vbArP6uxyk4+gHkhrNDwDvzltqa0iVeDrpgqFHgFNIW13w7I97iaQbcmxPkK7MnR8RT+fh1wD3RcTDuUxf0/0Oac/v95KWk5L4ruuwfOeQGrOlpIsDDqpT7t+A4/M8/x/190BqqRtzPlT1VeDP+TDKa5uc5mdJyXs5qfE9rzIgIh4hXWD2DdIPehpr16Fqh5Ku/FxCurjmL4VhrwGuk7QiL8MnI+KeGtM4k3Q4637SxQXXNrkc/Yn3OXWrxjSXkxLw+aS6+z6ePczVjA8A9+bf6sdIDVctXwe+kP93n60Rx29JG9FXkH7/V1QV6atdeDDHv4jUcH8sIv6ah50KbJ/n/csSy1Z0Fanxv7pOd/Xy9Le+lpbbiH2BvUntwA+ADxaW/0TSVcUPkfbMaiW2Rl5Iulr9cdLpk6tISQbgvaS9x0XARcAXI+LSGtPod73POkhXdi8itUFvJrU1AERED6nNDNI1IhXTSHVlBakN/UFEXJmHrVUn84bATNJRqcWkPdD/pJ/5TemUTb12MpXJJ0dtPSTpdNKFHl9odyxm1fIe1tkR0bI9Axv6JJ0GLBpO7dSQuGnWzMxGHqULtd5NeqDOsOGnj5iZ2aCT9GXShaH/Vef0yJDlQ7VmZmYleI/TzMysBJ/jBLbYYouYMmVKu8MwMxtW5s2b90hE1HzQzfrMiROYMmUK3d3d7Q7DzGxYkXRf36XWPz5Ua2ZmVoITp5mZWQlOnGZmZiU4cZqZmZXgxGlmZlaCE6eZmVkJTpxmZmYlOHGamZmV4MRpZmZWghOnmZlZCU6cZmZmJbQ1cUqaIelvkhZIOrrG8LGSzsvDr8svPUXSnpLmSbo1/31rYZydc/8Fkr4rSYO3RGZmtr5rW+KUNAo4Cdgb2B54r6Ttq4odAiyLiKnAicA3c/9HgHdGxCuB2cBZhXF+CBwGTMufGS1bCDMzG3Hauce5C7AgIu6OiKeBc4GZVWVmAmfk7xcAe0hSRNwYEYty//nAuLx3uiWwaURcE+kN3WcC+7V+UczMbKRoZ+KcDCwsdPfkfjXLRMRq4DFgYlWZ/YEbI2JlLt/TxzQBkHSYpG5J3YsXL+73QpiZ2cjSzsRZ69xjlCkj6RWkw7cfLTHN1DPi5IjoioiuSZNG3HtYzcysn9qZOHuArQvdncCiemUkjQY2A5bm7k7gIuCDEXFXoXxnH9M0MzPrt3YmzrnANEnbSNoAmAXMqSozh3TxD8ABwBUREZI2B34DHBMRf64UjogHgOWSXpuvpv0gcHGrF8TMzEaOtiXOfM7yCOAS4A7g/IiYL+l4SfvmYqcCEyUtAD4DVG5ZOQKYChwr6ab8eX4edjhwCrAAuAv47eAskZmZjQRKF5+ObF1dXdHd3d3uMMzMhhVJ8yKiq91xDDY/OcjMzKwEJ04zM7MSnDjNzMxKcOI0MzMrwYnTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxKcOM3MzEpw4jQzMyvBidPMzKwEJ04zM7MSnDjNzMxKcOI0MzMrwYnTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxLamjglzZD0N0kLJB1dY/hYSefl4ddJmpL7T5T0B0krJH2/apwr8zRvyp/nD87SmJnZSDC6XTOWNAo4CdgT6AHmSpoTEbcXih0CLIuIqZJmAd8EDgSeAo4FdsifagdFRHdLF8DMzEakdu5x7gIsiIi7I+Jp4FxgZlWZmcAZ+fsFwB6SFBFPRMSfSAnUzMxs0LQzcU4GFha6e3K/mmUiYjXwGDCxiWn/JB+mPVaSahWQdJikbkndixcvLh+9mZmNSO1MnLUSWvSjTLWDIuKVwBvz5wO1CkXEyRHRFRFdkyZN6jNYMzMzaG/i7AG2LnR3AovqlZE0GtgMWNpoohFxf/67HDiHdEjYzMxsQLQzcc4FpknaRtIGwCxgTlWZOcDs/P0A4IqIqLvHKWm0pC3y9zHAO4DbBjxyMzMbsdp2VW1ErJZ0BHAJMAo4LSLmSzoe6I6IOcCpwFmSFpD2NGdVxpd0L7ApsIGk/YC9gPuAS3LSHAVcBvx4EBfLzMzWc2qwAzdidHV1RXe3714xMytD0ryI6Gp3HIPNTw4yMzMrwYnTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxKcOM3MzEpw4jQzMyvBidPMzKwEJ04zM7MSnDjNzMxKcOI0MzMrwYnTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxKcOM3MzEpw4jQzG6GWrFjJzQsfZcmKle0OZVgZ3e4AzMxs8F180/0cdeEtjOnoYFVvLyfsvyP7Tp/c7rCGBe9xmpmNMEtWrOSoC2/hqVW9LF+5mqdW9XLkhbd4z7NJbU2ckmZI+pukBZKOrjF8rKTz8vDrJE3J/SdK+oOkFZK+XzXOzpJuzeN8V5IGZ2nMzIaHnmVPMqZj7eZ/TEcHPcuebFNEw0tTiVPSCyXtIun1lc+6zljSKOAkYG9ge+C9kravKnYIsCwipgInAt/M/Z8CjgU+W2PSPwQOA6blz4x1jdXMbH3SOWE8q3p71+q3qreXzgnj2xTR8NJn4pT0NeB64CukZHUs8IUBmPcuwIKIuDsingbOBWZWlZkJnJG/XwDsIUkR8URE/ImUQIuxbglsGhHXREQAZwL7DUCsZmbrjYkbj+WE/Xdk3JgONhk7mnFjOjhh/x2ZuPHYdoc2LDRzcdD+wEsj4qk+S5YzGVhY6O4Bdq1XJiJWS3oMmAg80mCaPVXTrHm2W9JhpD1TXvSiF5WN3cxsWNt3+mR2m7oFPcuepHPCeCfNEppJnPfQmnOhtc49Rj/K9Kt8RJwMnAzQ1dXVaJpmZuuliRuPdcLsh7qJU9KJpKSzHLhR0mXAM5dcRcRn1nHePcDWhe5OYFGdMj2SRgObAUv7mGZnH9M0MzPrt0Z7nLflv/OB31UNG4g9tLnANEnbAPcDs4D3VZWZA8wGrgEOAK7I5y5riogHJC2X9FrgOuCDwPcGIFYzMzOgQeKMiFMBJB0REdW3fByxrjPO5yyPAC4BRgGnRcR8SccD3RExBzgVOEvSAtKe5qxCDPcCmwIbSNoP2CsibgcOB04HxgO/zR8zM7MBoQY7cKmAdENE7FTV78aIeHVLIxtEXV1d0d3d3e4wzMyGFUnzIqKr3XEMtkbnOA8k7eFtI+kXhUGbAo+2OjAzM7OhqNE5zuuBJaQLbE4q9F8O3NjKoMzMzIaqRuc47yHdinKZpEnAa0gXBd0TEasGKT4zM7MhpZknBx0M3EC64vX9QLek2S2Oy8zMbEhq5gEIxwA7RcRiAElbAH/m2UfhmZmZjRjNPBHofta+GOgx1n6snZmZ2YjRzB7nP4BrJP2SdI5zP2CupE8ARMR3WxifmZnZkNJM4lyYP5UHGlaeIjSpJRGZmQ0DS1as9APSR6g+E2dEHAvppdIR4deDm9mId/FN93PUhbcwpqODVb29nLD/juw7veaLmGw91MxVtbtIuhW4M3e/SpKf/2pmI9KSFSs56sJbeGpVL8tXruapVb0ceeEtLFnh/YqRopmLg74LvIP0MAQi4mbgLa0MysxsqOpZ9iRjOtZuOsd0dNCz7Mk2RWSDrZnE2RER91X1W9OKYMzMhrrOCeNZ1du7Vr9Vvb10ThjfpohssDWTOBdK2gUISaMkfQr4e4vjMjMbkiZuPJYT9t+RcWM62GTsaMaN6eCE/Xf0BUIjSDNX1R5OOlz7IuAh4LLcz8xsRNp3+mR2m7qFr6odoZq5qvZhCu/BNDOztOfphDkyNTxUK+mNks6XdHP+nCvpDYMVnJmZ2VBTN3FK2hs4C7gU+DBwCHAFcKaktw9OeGZmZkNLo0O1RwLviojiuze7JV0PfBu4pKWRmZmZDUGNDtVuVZU0AYiIm4AXti4kMzOzoatR4lzRYNgTAx2ImZnZcNDoUO22kn5Ro7+Al7QoHjMzsyGtUeLcv8Gw7w90IGZmZsNB3cQZEZcPZiBmZmbDQTOP3GsZSTMk/U3SAklH1xg+VtJ5efh1kqYUhh2T+/+teHuMpHsl3SrpJkndg7MkZmY2UjTzyL2WkDQKOAnYE+gB5kqaExG3F4odAiyLiKmSZgHfBA6UtD3paUavALYCLpP00oioPHz+LRHxyKAtjJmZjRjNvI/z3c3064ddgAURcXdEPA2cC8ysKjMTOCN/vwDYQ5Jy/3MjYmVE3AMsyNMzMzNrqWYO1X6hRr/PD8C8JwMLC909uV/NMhGxGngMmNjHuAH8XtI8SYfVm7mkwyR1S+pevHjxOi2ImZmNHHUP1ebzhjOAyZK+VRi0KdBbe6xSVKNfNFmm0bi7RcQiSc8HLpX014i4+jmFI04GTgbo6uqqnq+ZmVlNjc5xPgzcBjwFzC/0Xw4850KefugBti50dwKL6pTpkTQa2AxY2mjciKj8fVjSRaRDuM9JnGZmZv3R6HaUG4EbJZ0dEStbMO+5wDRJ2wD3ky72eV9VmTnAbOAa4ADgiogISXOAc/Ke8FbANOB6SRsBHRGxPH/fCzi+BbGbmdkI1ehQ7c8i4r3AtZKecygzInZalxlHxGpJR5AeFj8KOC0i5ks6HuiOiDnAqcBZkhaQ9jRn5XHnSzofuB1YDXw8ItZIegFwUbp+iNHAORHxu3WJ08zMrEgRtU/vSeqMiB5J29YaHhF3tTSyQdTV1RXd3b7l08ysDEnzIqKr3XEMtrpX1UZET/66CLi7kCi3A+5rdWBmZmZDUTO3o/wRGC9pS+Aq4HDgtJZGZWZmNkQ1kzg7IuKfpIe+fz8i3gns2NqwzMzMhqamEqek15CueP117jeqdSGZmZkNXc0kzs8AXwJ+ExG3SXoJ6fCtmZnZiNPnQ94j4grgikL33cC/tTIoMzOzoarPxClpKmmvc0qxfETs1bqwzMzMhqZmXit2AelBBGcDa/ooa2Zmtl5rJnH2RsT3Wh6JmZnZMNDMxUEX51dwTZK0aeXT8sjMzMyGoGb2OD+S/x5b6BfAiwY+HDMzs6Gtmatqt+6rjJmZ2UjRzB4nkl4GbA+Mq/SLiHNaFZSZmdlQ1cztKF8gvdfyZaRXgL0d+BPgxGlmZiNOMxcHHQi8BXggIj4AvIom91TNhoIlK1Zy88JHWbKiFe9jN7ORppkE+GR+SfRqSZsADwIvaXFcZv22ZMVKepY9SeeE8fxpwSMcdeEtjOnoYFVvLyfsvyP7Tp/c7hDNbBhrJnHeKGlz0qvEuoHHgRtaGpVZP1180/3PJMqn1/SypreX1b3wFL0AHHnhLew2dQsmbjy2zZGa2XDVMHFKEnBcRDwKnCTpEmDTiHDitCFnyYqVHHXhLTy1qveZRFltTEcHPcuedOI0s35reI4zIoJnXyVGRCxw0rShqmfZk4zpaHzaflVvL50Txg9SRGa2Pmrm4qDrJe3U8kjM1lHnhPGs6l17T3PMKDF2tNhk7GjGjenghP139N6mma2TuodqJY2OiNXAG4BDJd0FPAGItDPqZGpDysSNx3LC/jtyZNXFQLtN3eKZi4WcNM1sXTU6x3k9sBOw3yDFYrbO9p0+uWaidMI0s4HSKHEKICLuGqRYjGdvpdhog1E88fQa7yX1w8SNx3qdmVnLNEqckyR9pt7AiPjWus5c0gzgO8Ao4JSI+EbV8LHAmcDOwBLgwIi4Nw87BjiE9I7QT0TEJc1McyAV7xcEan6fuPHYtco1atArt1JEb7ByTTBuTDoF7XsPzcyGjkaJcxSwMXnPc6BJGgWcBOwJ9ABzJc2JiNsLxQ4BlkXEVEmzgG8CB0raHpgFvALYCrhM0kvzOH1Nc0AU7xd8avUaIoLxY0av9X1Vby//unMn58/r6fMG/OKtFBWV77730Mxs6GiUOB+IiONbOO9dgAURcTeApHOBmUAxyc0EjsvfLwC+n+8tnQmcGxErgXskLcjTo4lprrN69wsuX7n6Od/PvPYfQN834Fdupah1/6HvPTQzGzoa3Y7Skj3NgsnAwkJ3T+5Xs0y+wvcxYGKDcZuZJgD55dzdkroXL15cKvBm7hesp5IEq9W6laLC9x6amQ0djVr/PVo871qJOZosU7b/c3tGnBwRXRHRNWnSpIaBVmuU5PpSLwlWbqUYN6aDsaPSYowb0+F7D83Mhpi6h2ojYmmL590DFF+S3QksqlOmR9JoYDNgaR/j9jXNdVZ9v2DDc5xdnZzfvfY5znpJsHgrha+qNTMbmtr5erC5wDRJ2wD3ky72eV9VmTnAbOAa4ADgiogISXOAcyR9i3Rx0DTSfadqYpoDovp+Qah/Ve0n93hp0zfg+1YKM7OhrW2JMyJWSzqC9HLsUcBpETFf0vFAd0TMAU4FzsoX/ywlJUJyufNJF/2sBj4eEWsAak2zVctQneQafXcyNDNbPyg9x31k6+rqiu7u7naHYWY2rEiaFxFd7Y5jsPXv0lAzM7MRyonTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxKcOM3MzEpw4jQzMyvBidPMzKwEJ04zM7MSnDjNzMxKcOI0MzMrwYnTzMysBCdOMzOzEpw4zczMSnDiNDMzK8GJ08zMrAQnTjMzsxKcOM3MzEpw4jQzMyvBidPMzKyEtiROSc+TdKmkO/PfCXXKzc5l7pQ0u9B/Z0m3Slog6buSlPsfJ+l+STflzz6DtUxmZjYytGuP82jg8oiYBlyeu9ci6XnAF4FdgV0xK/0mAAAL20lEQVSALxYS7A+Bw4Bp+TOjMOqJETE9f/6vhctgZmYjULsS50zgjPz9DGC/GmXeDlwaEUsjYhlwKTBD0pbAphFxTUQEcGad8c3MzAZcuxLnCyLiAYD89/k1ykwGFha6e3K/yfl7df+KIyTdIum0eoeAASQdJqlbUvfixYv7uxxmZjbCtCxxSrpM0m01PjObnUSNftGgP6RDuNsC04EHgP+pN/GIODkiuiKia9KkSU2GZGZmI93oVk04It5Wb5ikhyRtGREP5EOvD9co1gPsXujuBK7M/Tur+i/K83yoMI8fA7/ub/xmZma1tOtQ7RygcpXsbODiGmUuAfaSNCEfct0LuCQf2l0u6bX5atoPVsbPSbjiXcBtrVoAMzMbmVq2x9mHbwDnSzoE+AfwHgBJXcDHIuIjEbFU0peBuXmc4yNiaf5+OHA6MB74bf4AnCBpOunQ7b3ARwdhWczMbARRujB1ZOvq6oru7u52h2FmNqxImhcRXe2OY7D5yUFmZmYlOHGamZmV4MRpZmZWghOnmZlZCU6cZmZmJThxmpmZleDEaWZmVoITp5mZWQlOnGZmZiU4cZqZmZXgxGlmZlaCE6eZmVkJTpxmZmYlOHGamZmV4MRpZmZWghOnmZlZCU6cZmZmJThxmpmZleDEaWZmVoITp5mZWQlOnGZmZiU4cZqZmZXgxGlmZlZCWxKnpOdJulTSnfnvhDrlZucyd0qaXej/VUkLJa2oKj9W0nmSFki6TtKU1i6JmZmNNO3a4zwauDwipgGX5+61SHoe8EVgV2AX4IuFBPur3K/aIcCyiJgKnAh8swWxm5nZCNauxDkTOCN/PwPYr0aZtwOXRsTSiFgGXArMAIiIayPigT6mewGwhyQNaORmZjaitStxvqCS+PLf59coMxlYWOjuyf0aeWaciFgNPAZMrFVQ0mGSuiV1L168uGT4ZmY2Uo1u1YQlXQa8sMagzzc7iRr9YqDGiYiTgZMBurq6+pqumZkZ0MLEGRFvqzdM0kOStoyIByRtCTxco1gPsHuhuxO4so/Z9gBbAz2SRgObAUvLxG1mZtZIuw7VzgEqV8nOBi6uUeYSYC9JE/JFQXvlfs1O9wDgiojw3qSZmQ2YdiXObwB7SroT2DN3I6lL0ikAEbEU+DIwN3+Oz/2QdIKkHmBDST2SjsvTPRWYKGkB8BlqXK1rZma2LuQdsnSOs7u7u91hmJkNK5LmRURXu+MYbH5ykJmZWQlOnGZmZiU4cZqZmZXgxGlmZlaCE6eZmVkJTpxmZmYl+HYUQNJi4L4mim4BPNLicAbKcIoVhle8jrU1hlOsMLzibVWsL46ISS2Y7pDmxFmCpO7hcs/ScIoVhle8jrU1hlOsMLziHU6xDgc+VGtmZlaCE6eZmVkJTpzlnNzuAEoYTrHC8IrXsbbGcIoVhle8wynWIc/nOM3MzErwHqeZmVkJTpxmZmYlOHE2SdIMSX+TtEDSkH3Pp6TTJD0s6bZ2x9IXSVtL+oOkOyTNl/TJdsfUiKRxkq6XdHOO90vtjqkvkkZJulHSr9sdSyOS7pV0q6SbJA3pd/xJ2lzSBZL+muvu69odUz2StsvrtPJ5XNKn2h3XcOdznE2QNAr4O+ml2z2kF2u/NyJub2tgNUh6E7ACODMidmh3PI1I2hLYMiJukLQJMA/YbyiuVwBJAjaKiBWSxgB/Aj4ZEde2ObS6JH0G6AI2jYh3tDueeiTdC3RFxJB/oICkM4A/RsQpkjYANoyIR9sdV19yO3Y/sGtENPPAF6vDe5zN2QVYEBF3R8TTwLnAzDbHVFNEXA0sbXcczYiIByLihvx9OXAHMLm9UdUXyYrcOSZ/huyWp6RO4F+AU9ody/pC0qbAm4BTASLi6eGQNLM9gLucNNedE2dzJgMLC909DOEGfjiSNAV4NXBdeyNpLB/6vAl4GLg0IoZyvN8GjgR62x1IEwL4vaR5kg5rdzANvARYDPwkHwI/RdJG7Q6qSbOAn7U7iPWBE2dzVKPfkN3TGG4kbQxcCHwqIh5vdzyNRMSaiJgOdAK7SBqSh8MlvQN4OCLmtTuWJu0WETsBewMfz6cchqLRwE7ADyPi1cATwJC95qEiH1LeF/h5u2NZHzhxNqcH2LrQ3QksalMs65V8rvBC4KcR8Yt2x9OsfHjuSmBGm0OpZzdg33zu8FzgrZLObm9I9UXEovz3YeAi0umRoagH6CkcabiAlEiHur2BGyLioXYHsj5w4mzOXGCapG3yltssYE6bYxr28sU2pwJ3RMS32h1PXyRNkrR5/j4eeBvw1/ZGVVtEHBMRnRExhVRfr4iI97c5rJokbZQvDiMf9twLGJJXhUfEg8BCSdvlXnsAQ/JitirvxYdpB8zodgcwHETEaklHAJcAo4DTImJ+m8OqSdLPgN2BLST1AF+MiFPbG1VduwEfAG7N5w0BPhcR/9fGmBrZEjgjX53YAZwfEUP6No9h4gXARWk7itHAORHxu/aG1NC/Az/NG9F3Ax9qczwNSdqQdEfAR9sdy/rCt6OYmZmV4EO1ZmZmJThxmpmZleDEaWZmVoITp5mZWQlOnGZmI0yZl0FIOrHwkPi/SxoujxhsGSdOG3EkrcmNwG2Sfp4v10fSir7G7ce89pO0fZ1hx0m6P8dyu6T39mP6K/LfrSRdsK7x5ml9StIHS5TfQNLVknx72/BxOk0+vCMiPh0R0/MTs74HDJsHlbSKE6eNRE/mhmAH4GngYy2c135AzcSZnZgbpJnA/+YnKZUWEYsi4oD+jFuUk9+HgXNKzPtp4HLgwHWdvw2OWi+DkLStpN/l5wX/UdLLaozqByngxGn2R2BqsYekjSVdLumG/I7Imbn/lPz+xR/n93H+Pj9BqGajI+n1pOeD/lfeq9y2XhARcSfwT2BCnt6hkuYqvfvzwsJe8TaSrsnDvlyIeUrlsJukgyV9vzDs15J2zw+oPz3vad8q6dM1Qnkr6dFsq/O4V+ZDdVfnZX+NpF9IulPSVwrj/RI4qOm1bkPRycC/R8TOwGeBHxQHSnoxsA1wRRtiG1KcOG3EyntXewO3Vg16CnhXfuj4W4D/yY8HBJgGnBQRrwAeBfbP/Z/T6ETEX0iPZvzPvId7V4NYdgLuzM9qBfhFRLwmIl5Fet3aIbn/d0gPGH8N8GDJRZ4OTI6IHSLilcBPapTZjfRe1KKnI+JNwI+Ai4GPAzsAB0uamMvcBrymZDw2ROQXLbwe+Hl+itf/kp6UVTQLuCAi1gx2fEONz0nYSDS+8Ii/P5LfrVgg4Gv5DR29pFfIvSAPuyciKuPOA6ZUNTqVaYxtMpZPSzqU9Lqq4jmnHfIe3ebAxqTHPUJKbJVkfRbwzSbnA+nxcC+R9D3gN8Dva5TZkpSoiyrPZb4VmB8RDwBIupv08oMlEbFG0tOSNsnvVrXhpQN4NJ82qGcWaaNpxPMep41ElXOc0yPi3/M5uqKDgEnAzrkheQgYl4etLJRbQ9r4fKbRKXxe3mQsJ0bEdqTzg2dKqszndOCIvGf4pcL8oe9X2q1m7d/2OICIWAa8ivRWl49T+wXXT1bNC55d5l7WXv5e1t74HkvaW7dhJr/O7x5J74H0AgZJr6oMzw+1nwBc06YQhxQnTrPn2oz0LstVkt4CvLhR4T4aneXAJn3NML9SrRuYnXttAjyQLxYqnjv8M2nLH+qfU7wXmC6pQ9LW5Fd0SdoC6IiIC4Fjqf06rDuoOufbjHzIdnFErCo7rg0+pZdBXANsJ6lH0iGk+nSIpJuB+aQL1ireC5wbfrg54EO1ZrX8FPiVpG7gJpp7ddhBwA8lfQEYQ3oH5s35748lfQI4oNF5TuB44BxJPyYltuuA+0iHSCvJ95O5zCdJ7zGt5c/APXm824Abcv/JwE8kVTaYj6kx7m9Jh4DLegswVN9qY1Uiot6tTzVvUYmI41oXzfDjt6OY2VokXQQcma/0bXacXwDHRMTfWheZ2dDgQ7VmVu1onntFZV1K76X8pZOmjRTe4zQzMyvBe5xmZmYlOHGamZmV4MRpZmZWghOnmZlZCU6cZmZmJfx/Wm3Lzx4lHekAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data.plot(kind = 'scatter',y='Delta Flux', x='Radius (m)')\n", "#plt.gca().invert_yaxis()\n", "plt.title(\"Relation between planet radius and transit depth within our solar system.\")\n", "plt.xlabel(\"Planet Radius (m)\")\n", "plt.ylabel(\"Transit Depth\")\n", "\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5,1,'Mass-Period Relation (Solar System)')" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEWCAYAAABliCz2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3XucVXW9//HXe2AYEFARposQkmmWGqJNV7uY1i+zxIrqaGXasTydOl0t7fbo3rGo0/VX50RmmZqmknlJTTtqV7VAASW1FC0QS0QQRmEYmM/5Y31H12zmsobZa/YM6/18POYx6/797O/e85nv+q61v0sRgZmZ7fyaGh2AmZkNDyd8M7OKcMI3M6sIJ3wzs4pwwjczqwgnfDOzinDCt7qS9HFJZ+zgvj+S9IV6x5SOfb2kd+zgvjMltUsaU++40vHPk/TaOhyntPobDSS1SLpD0hMaHctI5YTfIJLulbRF0rSa5UskhaRZwxhLSHokJbX7JH1tR5NbRPxnROxQYh0gxhMlbUsxbpC0VNJr6l1OKuteSS/vno+Iv0fEpIjYVkJZs4GDgEvS/DhJ/yVpVXqt90j6er3LLRjbiyT9QdLDkh6S9HtJzxniMXvUbT1FRAdwJnBaGcffGTjhN9Y9wHHdM5KeBUxoUCwHRcQk4AjgzcA7B3sASWPrHlVPN6QYdwe+C5wvafeSyyzbvwHnxuPfgPwY0AY8F5gMvAy4pcwAlGmqWbYrcDnwbWAPYDrwWaCjzFjq4CfACZJaGh3ISOSE31hnA2/LzZ8A/Di/gaRXS7oltWpXSvpMbt14SedIWitpvaQ/SXpiWneipBWSNqZW4luKBBQRdwC/BQ5Mx9lT0kJJa9Jx3pcr/zOSLkoxbABOTMvOyW0zV9LyFN/1kp6ZW3ewpJtTjD8FxheMsSvV3URg39zxnp9apOvTGcBhve0v6WmSrk319qCkc7v/cUg6G5gJXJZa2KdKmpXOgsbm6uTS1Oq9S9I7c8f+jKQLJP04va7lktr6eTmvAn6dm38OcHFErI7MvRHx2GdC0jNTPa5Px57bx2ucIuny9L6tS9Mzcuuvl/RFSb8HHgX2rjnE01NdnxcR2yJiU0RcHRHLUtfJQ6mB0n28J0jaJKlV0rRU3vq03W8lNfVWtwO9bynOL6T17ZIukzQ1vWcb0md+Vvf2EbEKWAc8v586r66I8E8DfoB7gZcDdwLPBMYAK4G9gABmpe0OA55F9s95NvBP4LVp3b8BlwG7pP2fDexKlgg3APul7Z4MHNBPLAHsk6b3B/4BnJTKXAx8ChhHlhRWAK9M234G6ARem7adkJadk9Y/HXgEeAXQDJwK3JWONQ74G/DBtO4N6Vhf6CPGE4HfpekxwHuALcAT0rLpwFrgqBTLK9J8a1p/PfCONL1PWt8CtAK/Ab5R+97k5melOhqb5n9NdoYxHpgDrAGOyNXJ5hTHGOB04MY+XtPEdNzW3LJPAn8H3p3ed+XWNaf6+3iqv8OBjbn3+Ufd9QdMBealz8Zk4ELg57ljXZ/KOQAYCzTXxLZrqr+zyP4pTalZ/13gy7n59wOXpenTgf9J8TYDL+5+Hb3UbZH37S7gacBuwJ+Bv5D97YwlayD9sCa2S4H3NfpvfCT+jLgWvqQzJT0g6bYC235dWZ/3Ekl/kbR+OGKss+5W/iuAO4D78isj4vqIuDUiuiJiGXAe8NK0upPsD3ufyFphiyNiQ1rXBRwoaUJE3B8RyweI42ZJ68j+gZwB/JCstdkaEZ+LiC0RsQL4PnBsbr8bIuLnKb5NNcf8F+AXEXFNRHQCXyX7p/BCshZYM1mi7YyIi4A/DRDj89N7vDkd660R8UBa91bgioi4IsVyDbCILJH0EBF3pZg6ImIN8DUer9N+SXoK8CLgtIjYHBFLUn0dn9vsdymObWTv70F9HK67O2pjbtnpwJeBt6T475N0QvfrByYBX0rvx7Vk3S7HUSMi1kbEwoh4NCI2Al/s5TX+KCKWR8TW9P7k99+QXmeQvedr0lnNE9MmZwFvznUFHZ9eK2SfyycDe6X39reRMnEvirxvP4yIuyPiYeBK4O6I+FVEbCX7R3ZwzTE38njdWs6IS/hkrZQji2wYER+MiDkRMYesr/FnZQZWkrPJ+sxPpKY7B0DS8yRdl07NHwbeBUzL7ftLsr7s1ZLmS2qOiEfIku27gPsl/ULSM9LxlqdT43ZJL84VdUhETImIp0XEJyPrNtkL2DOdaq9PyfbjwBNz+63s57XtSdaKBx7rillJ1qrbE7ivJhH8jf7dGBG7A1PIWnH5+PcC3lgT64vIEk8PqfvhfGUXqDcA5/B4nQ5kT+ChlETzcU/Pzf8jN/0oMF69X9/obqBM7l6Q/nF/JyIOJUtaXwTOTF1hewIrUz32VXb3a9xF0vck/S29xt8Au6vnxfj+3jsi4vaIODEiZpB18e0JfCOtu4ns7O2l6bO1D9l7AvAVslb51cq6FT/aTzFF3rd/5qY39TI/qeaYk3m8bi1nxCX8iPgN8FB+WepzvUrS4tQf+Ixedj2OrPU7qkTE38gu3h5F7/+wfkL2h/SUiNiN7FRZad/OiPhsROxP1mp+DemaQET8MiJeQfaHcwdZK42IOCCyO04mRcRvBwhvJXBPROye+5kcEfnWV3/Dra4m+4MGsouDwFPIzmLuB6anZd1mDhAP6TW0k3V5HC+pu3W3Eji7JtaJEfGlXg5xeop7dkTsStbKzMcx0GvaQ9Lk3LKZ1JyZFXwdjwB3k/rLe1m/KSK+Q9YnvX8q+ynqeYG1r7JPAfYDnpde40vS8qKvszaWO8gaYwfmFp9FVnfHAxdFxOa07caIOCUi9gaOBj4k6Yg+yhzM+1bUM4GlQ9h/pzXiEn4fFgDvjYhnAx8m6z98jKS9gKcC1zYgtno4CTg8JYBak8lalJslPZfsbAAASS+T9KzUattAdiq9TdITlV0snUh2V0U7sCO3FP4R2CDpNEkTJI2RdKCK35p3AfBqSUdIaiZLQh3AH4AbgK3A+ySNlfR6sjtTComItWRdKZ9Ki84Bjpb0yhTneEmH5S9U5kwmq5P1kqYDH6lZ/0+2v4jZXe7KFP/pqYzZZO/fuUVjr3EFua4WSR9IcU9I9XJCivcWoLtVfaqk5nRx82jg/D5e46b0GvcAPj2YoCQ9Q9Ip3fWXurKOA27MbXY28DqypJ+/sPwaSfukf+YbyD573Z+/2rodzPtWJO7pZHcV3TjQtlU04hO+pElkrdcLJS0Bvsf2p+nHkrUw6n6f9HBI/ZOL+lj9buBzkjaSJbcLcuueBFxE9kd1O9nFxHPI3tdTyFqED5EllHfvQFzbyBLKHLKzkAfJkuxuBfe/kywZfDvtezRwdOp/3gK8nqwrax1ZF9Rgu+S+ARwlaXZKxMeQdTmtIWs5foTeP+OfBQ4BHgZ+0Uu5pwOfTF0MH+5l/+PILuSuBi4GPp36nnfEAuAtuTOdTcB/kXULPUh2cXpeRKxIdTaX7CLqg2QNn7el1netb5BdL3mQLPldNci4NgLPA26S9Eg6xm1knyvgsTtibiZrtefPFvcFfkX2T/UG4LsRcX1a16NuB/m+FfFm4KzI7sm3Gt1XzkeUdJvV5RFxoLL7ge+MiO36YnPb3wK8JyL+MEwhmtWNpJ8AF0TEzxsdy2BJOhNYHRGfHAGxtJB15bwkdzHfcsr+osyQRcQGZfd/vzEiLkwtodkRsRRA0n5kF/FuaGigZjsoIt488FYjT2qYvZ7t75JpiNSq7+36niUjrktH0nlkyXs/ZV8vP4nsFrWTJC0FlpOdAnY7Dji/n9u+zKzOJH2erIvnKxFxT6PjsWJGZJeOmZnV34hr4ZuZWTlGVB/+tGnTYtasWY0Ow8xs1Fi8ePGDEdFaZNsRlfBnzZrFokV93Z1oZma1JA30DfXHuEvHzKwinPDNzCrCCd/MrCKc8M3MKsIJ38ysIpzwzcwaaG17B0tXrmdte/njvY2o2zLNzKrkkiX3cdrCZTQ3NdHZ1cX8ebOZO2e759nUjVv4ZmYNsLa9g9MWLmNzZxcbO7ayubOLUxcuK7Wl74RvZtYAq9ZtormpZwpubmpi1braR0PXT2kJX9J+evwB40skbZD0gbLKMzMbTWZMmUBnV1ePZZ1dXcyYMqG0MktL+BFxZ+4B488me5jzxWWVZ2Y2mkyd1ML8ebMZ39zE5JaxjG9uYv682Uyd1FJamcN10fYI4O70wG4zMwPmzpnOoftMY9W6TcyYMqHUZA/Dl/CPBc7rbYWkk4GTAWbOnDlM4ZiZjQxTJ7WUnui7lX7RVtI4sgcvX9jb+ohYEBFtEdHW2lpohE8zM9sBw3GXzquAmyPin8NQlpmZ9WE4Ev5x9NGdY2Zmw6fUhC9pF+AVwM/KLMfMzAZW6kXbiHgUmFpmGWZmVoy/aWtmVhFO+GZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFO+GZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFO+GZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFlP8R8d0kXSbpD0u2SXlBmeVa+te0dLF25nrXtHY0OxcwGqdSHmAPfBK6KiDdIGgfsUnJ5VqJLltzHaQuX0dzURGdXF/PnzWbunOmNDsvMCiqthS9pV+AlwA8AImJLRKwvqzwr19r2Dk5buIzNnV1s7NjK5s4uTl24zC19s1GkzC6dvYE1wA8l3SLpDEkTSyzPSrRq3Saam3p+XJqbmli1blODIjKzwSoz4Y8FDgH+OyIOBh4BPlq7kaSTJS2StGjNmjUlhmNDMWPKBDq7unos6+zqYsaUCQ2KyMwGq8yEvwpYFRE3pfmLyP4B9BARCyKiLSLaWltbSwzHhmLqpBbmz5vN+OYmJreMZXxzE/PnzWbqpJZGh2ZmBZV20TYi/iFppaT9IuJO4Ajgz2WVZ+WbO2c6h+4zjVXrNjFjygQne7NRpuy7dN4LnJvu0FkBvL3k8qxkUye1ONGbjVKlJvyIWAK0lVmGmZkV42/amplVhBO+mVlFOOGbmVWEE76ZWUU44ZuZVYQTvplZRTjhm5lVhBO+mVlFOOGbmVWEE76ZWUU44ZuZVYQTvplZRTjhm5lVhBO+mVlFOOGbmVWEE76ZWUU44ZuZVYQTvplZRTjhm5lVhBO+mVlFlPoQc0n3AhuBbcDWiPADzc3MGqTUhJ+8LCIeHIZyzMysH+7SMTOriLITfgBXS1os6eTeNpB0sqRFkhatWbOm5HDMzKqr7IR/aEQcArwKeI+kl9RuEBELIqItItpaW1tLDsfMrLpKTfgRsTr9fgC4GHhumeWZmVnfSkv4kiZKmtw9Dfw/4LayyjMzs/6VeZfOE4GLJXWX85OIuKrE8szMrB+lJfyIWAEcVNbxzcxscHxbpplZRTjhm5lVhBO+mVlFOOGbmVWEE76ZWUU44ZuZVYQTvplZRRS6D19SG/BiYE9gE9k3Zn8VEQ+VGJuZmdVRvy18SSdKuhn4GDABuBN4AHgRcI2ksyTNLD9MMzMbqoFa+BPJRrzc1NtKSXOAfYG/1zswMzOrr34TfkR8Z4D1S+objpmZlaVoH/5TgfcCs/L7RMTccsIyM7N6Kzp42s+BHwCXAV3lhWNmZmUpmvA3R8S3So3EzMxKVTThf1PSp4GrgY7uhRFxcylRmZlZ3RVN+M8CjgcO5/EunUjzZmY2ChRN+K8D9o6ILWUGY2Zm5Sk6tMJSYPcyAzEzs3IVbeE/EbhD0p/o2Yfv2zLNzEaJogn/0ztagKQxwCLgvoh4zY4ex8zMhqbfhC9Jkfn1QNv0c5j3A7cDu+5gjGZmVgcD9eFfJ+m9tQOkSRon6XBJZwEn9LWzpBnAq4Ezhh6qmZkNxUBdOkcC/wqcl4ZXWA+MB8aQ3ZP/9QHG0/kGcCowua8NJJ0MnAwwc6YH3jQzK8tAg6dtBr4LfFdSMzAN2BQR6wc6sKTXAA9ExGJJh/VTxgJgAUBbW1t/XUNmZjYERS/aEhGdwP2DOPahwFxJR5GdFewq6ZyIeOsgYzQzszoo7RGHEfGxiJgREbOAY4FrnezNzBrHz7Q1M6uIQglf0kRJTWn66ZLmpj79QiLiet+Db2bWWEVb+L8BxkuaDvwv8HbgR2UFZWZm9Vc04SsiHgVeD3w7Il4H7F9eWGZmVm+FE76kFwBvAX6RlhW+w8fMzBqvaML/APAx4OKIWC5pb+C68sIyM7N6K9RKT2Pp/Do3vwJ4X1lBmZlZ/RVK+JKuI3vCVQ8R4SdemZmNEkX74T+cmx4PzAO21j8cMzMrS9EuncU1i34vqc8hk83MbOQp2qWzR262CXg28KRSIjIzs1IU7dJZTNaHL7KunHuAk8oKyszM6q9ol85Tyw7EzMzKNdAjDg+PiGslvb639RHxs3LCMjOzehuohf9S4Frg6F7WBeCEb2Y2Sgz0xKtPp99vH55wzMysLAN16Xyov/UR8bX6hmNmZmUZqEun++Hj+wHPAS5N80eTDZlsZmajxEBdOp8FkHQ1cEhEbEzznwEuLD06MzOrm6KjZc4EtuTmtwCz6h6NmZmVpugXr84G/ijpYrK7c14H/Li0qMzMrO6KfvHqi5KuAl6UFr09Im7pbx9J48n6+VtSORd13/VjZmbDr/BTqyJisaSVZKNlImlmRPy9n106gMMjoj098Px3kq6MiBuHFrKZme2IQn34kuZK+ivZGDq/Tr+v7G+fyLSn2eb0s92Y+mZmNjyKXrT9PPB84C9pXJ2XA78faCdJYyQtAR4AromIm3rZ5mRJiyQtWrNmzSBCNzOzwSia8DsjYi3QJKkpIq4D5gy0U0Rsi4g5wAzguZIO7GWbBRHRFhFtra2tgwrezMyKK9qHv17SJLKLsOdKeoBBPPEqItZLuh44Erht0FGamdmQFW3hHwM8CnwQuAq4m94HVHuMpFZJu6fpCWTdQHfseKhmZjYURW/LfCRNdgFnSRoDHAuc289uT85t2wRcEBGXDyVYMzPbcQMNnrYr8B5gOtk4Otek+Y8AS+gn4UfEMuDgukVqZmZDMlAL/2xgHXAD8A6yRD8OOCYilpQcm5mZ1dFACX/viHgWgKQzgAeBmd2DqJmZ2egx0EXbzu6JiNgG3ONkb2Y2Og3Uwj9I0oY0LWBCmhfZl2l3LTW6kq1t72DVuk3MmDKBqZNaGh2OmVmpBhoPf8xwBTLcLllyH6ctXEZzUxOdXV3MnzebuXOmNzosM7PSFL0Pf6eytr2D0xYuY3NnFxs7trK5s4tTFy5jbXtHo0MzMytNJRP+qnWbaG7q+dKbEMtXb+hjDzOz0a+SCX/GlAl0dnX1WPZo5zbe+eNFXLrkvgZFZWZWrkom/KmTWpg/bzYtY9VjecdWd+2Y2c6rkgkfYO6c6Xz/bW3sMq7ndenmpiZWrdvUoKjMzMpT2YQPcMCeu9EVPZ/J0tnVxYwpExoUkZlZeSqd8Lu7dsY3NzG5ZSzjm5uYP2+278k3s51S4Wfa7qzmzpnOoftM8xewzGynV/mED1lL34nezHZ2le7SMTOrEid8M7OKcMI3M6sIJ3wzs4pwwjczq4jSEr6kp0i6TtLtkpZLen9ZZZmZ2cDKvC1zK3BKRNwsaTKwWNI1EfHnEss0M7M+lNbCj4j7I+LmNL0RuB3wE0bMzBpkWPrwJc0CDgZu6mXdyZIWSVq0Zs2a4QjHzKySSk/4kiYBC4EPRMR2TxiJiAUR0RYRba2trWWHY2ZWWaUmfEnNZMn+3Ij4WZllmZlZ/8q8S0fAD4DbI+JrZZVjZmbFlNnCPxQ4Hjhc0pL0c1SJ5ZmZWT9Kuy0zIn4HaMANzcxsWPibtmZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFO+GZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFO+GZmFeGEb2ZWEU74ZmYV4YRvZlYRTvhmZhXhhG9mVhFO+GZmFVFawpd0pqQHJN1WVhlmZlZcmS38HwFHlnh8MzMbhNISfkT8BniorOObmdngNLwPX9LJkhZJWrRmzZpGh2NmttNqeMKPiAUR0RYRba2trY0Ox8xsp9XwhG9mZsPDCd/MrCLKvC3zPOAGYD9JqySdVFZZZmY2sLFlHTgijivr2GZmNnju0jEzq4idJuGvbe9g6cr1rG3vaHQoZmYjUmldOsPpkiX3cepFyxjTJLZ1BV95w2zmzpne6LDMzEaUUd/CX9vewYcvXErH1i4e3bKNjq1dnHLhUrf0zcxqjPqEv3z1w3Ruix7LOrcFy1c/3KCIzMxGplGf8EGDXG5mVk2jPuEfsOeujK15FWObsuVmZva4UZ/wp05q4WtvmkPLWLFL8xhaxoqvvWkOUye1NDo0M7MRZae4S2funOkcus80Vq3bxIwpE5zszcx6sVMkfMha+k70ZmZ9G/VdOuAvXZmZFTHqW/iXLLmP0xYuo7mpic6uLubP85euzMx6M6pb+GvbOzht4TI2d3axsWMrmzu7OHXhMrf0zcx6MaoT/qp1mwa13MysykZ1wp84bgybO7t6LNvc2cXEcWMaFJGZ2cg1qhP+I1u2MabmC7VjlC03M7OeRnXCnzhuDDXD6LAtcAvfzKwXozrhP7JlW6/DKriFb2a2vVGd8CeOG8PWnl34bO1yC9/MrDelJnxJR0q6U9Jdkj5a7+M/smUb45t7voTxzU1u4ZuZ9aK0hC9pDPAd4FXA/sBxkvavZxkzpkwY1HIzsyors4X/XOCuiFgREVuA84Fj6lnA1EktzJ83m/HNTUxuGcv45ibmz5vtMXXMzHpR5tAK04GVuflVwPNqN5J0MnAywMyZMwddiEfKNDMrpsyE39sjp2K7BRELgAUAbW1t260vwiNlmpkNrMwunVXAU3LzM4DVJZZnZmb9KDPh/wnYV9JTJY0DjgUuLbE8MzPrR2ldOhGxVdJ/AL8ExgBnRsTyssozM7P+lToefkRcAVxRZhlmZlbMqP6mrZmZFaeIHboxphSS1gB/28HdpwEP1jGcenN8Q+P4hsbxDc1Ijm+viGgtsuGISvhDIWlRRLQ1Oo6+OL6hcXxD4/iGZqTHV5S7dMzMKsIJ38ysInamhL+g0QEMwPENjeMbGsc3NCM9vkJ2mj58MzPr387Uwjczs3444ZuZVcSIT/gDPTVLUoukn6b1N0malVv3sbT8TkmvbFB8H5L0Z0nLJP2vpL1y67ZJWpJ+ShtnqECMJ0pak4vlHbl1J0j6a/o5oUHxfT0X218krc+tK7UOJZ0p6QFJt/WxXpK+lWJfJumQ3LrhqLuB4ntLimuZpD9IOii37l5Jt6a6W9Sg+A6T9HDuPfxUbl2pT8wrGN9HcrHdlj5ve6R1pddf3UXEiP0hG4PnbmBvYBywFNi/Zpt3A/+Tpo8Ffpqm90/btwBPTccZ04D4Xgbskqb/vTu+NN8+QurwROD/97LvHsCK9HtKmp4y3PHVbP9esnGZhqUOgZcAhwC39bH+KOBKsuHAnw/cNFx1VzC+F3aXS/b0uZty6+4FpjW4/g4DLh/q56Ks+Gq2PRq4djjrr94/I72FX+SpWccAZ6Xpi4AjJCktPz8iOiLiHuCudLxhjS8irouIR9PsjWTDRA+noTx57JXANRHxUESsA64BjmxwfMcB59U5hj5FxG+Ah/rZ5Bjgx5G5Edhd0pMZnrobML6I+EMqHxrw+StQf30p/Yl5MOj4hvWzV4aRnvB7e2rW9L62iYitwMPA1IL7Dkd8eSeRtQa7jZe0SNKNkl5b59i6FY1xXjrtv0hS93MMRlQdpu6wpwLX5hYPRx32p6/4h6PuBqv28xfA1ZIWK3vyXKO8QNJSSVdKOiAtG1H1J2kXsn/YC3OLR0r9FVbqaJl1UOSpWX1tU+iJW0NUuAxJbwXagJfmFs+MiNWS9gaulXRrRNzdgBgvA86LiA5J7yI7Yzq84L7DEV+3Y4GLImJbbtlw1GF/Gvn5K0zSy8gS/otyiw9NdfcE4BpJd6QW73C6mWwsmHZJRwE/B/ZlhNUfWXfO7yMifzYwEupvUEZ6C7/IU7Me20bSWGA3slO04XjiVqEyJL0c+AQwNyI6updHxOr0ewVwPXBwneMrFGNErM3F9X3g2UX3HY74co6l5pR6mOqwP33FP2Ke+CZpNnAGcExErO1enqu7B4CLqX+X54AiYkNEtKfpK4BmSdMYQfWX9PfZa1j9DVqjLyL090N2BrKC7DS++8LNATXbvIeeF20vSNMH0POi7Qrqf9G2SHwHk1182rdm+RSgJU1PA/5KORelisT45Nz064Ab0/QewD0p1ilpeo/hji9ttx/ZRTI1oA5n0fdFx1fT86LtH4er7grGN5Ps+tULa5ZPBCbnpv8AHNmA+J7U/Z6SJcy/p7os9LkoO760vrsRObER9VfX19roAAq8GUcBf0lJ8xNp2efIWssA44EL04f6j8DeuX0/kfa7E3hVg+L7FfBPYEn6uTQtfyFwa/og3wqc1MA6PB1YnmK5DnhGbt9/TXV7F/D2RsSX5j8DfKlmv9LrkKxVdz/QSdbqPAl4F/CutF7Ad1LstwJtw1x3A8V3BrAu9/lblJbvneptaXrvP9Gg+P4j99m7kdw/pt4+F8MdX9rmRLIbQPL7DUv91fvHQyuYmVXESO/DNzOzOnHCNzOrCCd8M7OKcMI3M6sIJ3wzsxINNEBbzba9DrYo6WW5QdyWSNq8I98sd8K3SpIUks7OzY9NI4Ze3si4bKf0I4qPo3QL2a29s8nGBpsPj43JNSci5pB9C/5R4OrBBuKEb1X1CHCgpAlp/hXAfQ2Mx3ZS0csAbZKeJumqNA7PbyU9I217XQw82OIbgCtz2xXmhG9VdiXZN2WhZiRESc9N48ffkn7vl5YfIOmP6bR6maR9JU2U9Is0ANhtkv6lAa/FRpcFwHsj4tnAh4Hv9rJN7WB33bYb5qGokT54mlmZzgc+lbpxZgNnAi9O6+4AXhIRW9NYSP8JzCP7FuY3I+JcSePIxm0/ClgdEa8GkLTbML8OG0UkTSL7lviF2UjuQDYETH6b3gZbJA29/SzglztSthO+VVZELFP2hLTjgCtqVu8GnCVpX7JRGpvT8huF91ptAAABSElEQVSAT0iaAfwsIv4q6Vbgq5K+TPYwj98Oywuw0aoJWJ/647eTG2zxpZEbbDF5E3BxRHTuaMFmVXYp8FW2P0X+PHBdRBxINjTueICI+AkwF9gE/FLS4RHxF7IRRm8FTs8/ps+sVkRsAO6R9EZ47DGZB6Xpg4HvkY0j9UAvuw/pISxu4VvVnQk8HBG3Sjost3w3Hr+Ie2L3wjTu/oqI+Faani3pDuChiDhHUnt+ezNJ55E9ynGapFXAp4G3AP8t6ZNkZ4/nkw3E9hVgEo939/w9Iuam48wiGzL61zsciwdPsyqS1B4Rk2qWHQZ8OCJeI+kFZA+CWUP2hK3jI2KWpI8BbyUbXfEfwJuB55D9oXal5f8eEaPjodZWKU74ZmYV4T58M7OKcMI3M6sIJ3wzs4pwwjczqwgnfDOzinDCNzOrCCd8M7OK+D88S5I6F4toLwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data.plot(kind = 'scatter',x = \"Mass\",y = \"Radius (m)\")\n", "plt.title(\"Mass-Period Relation (Solar System)\")\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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.4" } }, "nbformat": 4, "nbformat_minor": 2 }