Tle Spé - Loi de Kepler.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ce programme permet de montrer la validité de troisième loi de Kepler"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 0 : On importe la bibliothèque nécessaire"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np #Bibliothèque de calcul scientifique\n",
    "import matplotlib.pyplot as plt #Bibliothèque d'affichage de graphiques\n",
    "import scipy.stats as sc #Bibliothèque de calcul scientifique"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 1 : On affiche les données utiles pour tester le programme"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Vous trouverez ci-dessous les rayons et les periodes des orbites de Venus, la Terre, Mars et de Jupiter\n",
      "Venus : rayon orbite = 108 millions de km ; periode de revolution = 225 j \n",
      "Terre : rayon orbite = 150 millions de km ; periode de revolution = 365 j \n",
      "Mars : rayon orbite = 228 millions de km ; periode de revolution = 687 j \n",
      "Jupiter : rayon orbite = 778 millions de km ; periode de revolution = 4333 j \n"
     ]
    }
   ],
   "source": [
    "print(\"Vous trouverez ci-dessous les rayons et les periodes des orbites de Venus, la Terre, Mars et de Jupiter\")\n",
    "print(\"Venus : rayon orbite = 108 millions de km ; periode de revolution = 225 j \")\n",
    "print(\"Terre : rayon orbite = 150 millions de km ; periode de revolution = 365 j \")\n",
    "print(\"Mars : rayon orbite = 228 millions de km ; periode de revolution = 687 j \")\n",
    "print(\"Jupiter : rayon orbite = 778 millions de km ; periode de revolution = 4333 j \")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 2 : On saisie les demi-grand axes des orbites ainsi que des périodes des planètes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entrer le demi-grand axe de la 1ere planete en millions de km 108\n",
      "Le rayon est 108.0 millions de km.\n",
      "Entrer la periode de la 1ere planete en annees 0.6461\n",
      "La periode est 0.6461 annees.\n"
     ]
    }
   ],
   "source": [
    "a1 = float(input(\"Entrer le demi-grand axe de la 1ere planete en millions de km \"))\n",
    "print (\"Le rayon est\",a1,\"millions de km.\")\n",
    "T1 = float(input(\"Entrer la periode de la 1ere planete en annees \"))\n",
    "print (\"La periode est\",T1,\"annees.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entrer le demi-grand axe de la 2eme planete en millions de km 150\n",
      "Le rayon est 150.0 millions de km.\n",
      "Entrer la periode de la 2eme planete en annees 1\n",
      "La periode est 1.0 annees.\n"
     ]
    }
   ],
   "source": [
    "a2 = float(input(\"Entrer le demi-grand axe de la 2eme planete en millions de km \"))\n",
    "print (\"Le rayon est\",a2,\"millions de km.\")\n",
    "T2 = float(input(\"Entrer la periode de la 2eme planete en annees \"))\n",
    "print (\"La periode est\",T2,\"annees.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entrer le demi-grand axe de la 3eme planete en millions de km 228\n",
      "Le rayon est 228.0 millions de km.\n",
      "Entrer la periode de la 3eme planete en annees 1.88\n",
      "La periode est 1.88 annees.\n"
     ]
    }
   ],
   "source": [
    "a3 = float(input(\"Entrer le demi-grand axe de la 3eme planete en millions de km \"))\n",
    "print (\"Le rayon est\",a3,\"millions de km.\")\n",
    "T3 = float(input(\"Entrer la periode de la 3eme planete en annees \"))\n",
    "print (\"La periode est\",T3,\"annees.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entrer le demi-grand axe de la 4eme planete en millions de km 778\n",
      "Le rayon est 778.0 millions de km.\n",
      "Entrer la periode de la 4eme planete en annees 11.87\n",
      "La periode est 11.87 annees.\n"
     ]
    }
   ],
   "source": [
    "a4 = float(input(\"Entrer le demi-grand axe de la 4eme planete en millions de km \"))\n",
    "print (\"Le rayon est\",a4,\"millions de km.\")\n",
    "T4 = float(input(\"Entrer la periode de la 4eme planete en annees \"))\n",
    "print (\"La periode est\",T4,\"annees.\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 3 : On crée des listes pour le calcul et l'affichage"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "a = np.array([a1, a2, a3, a4])\n",
    "T = np.array([T1, T2, T3, T4])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 4 : On crée des listes $a^3$ et $T^2$ , on modélise la répartition des points par une droite et on calcule la masse du corps central"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "acube = (a*1e9)**3\n",
    "Tcarre = (T*24*3600)**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "droite = np.polyfit(acube,Tcarre,1) #On crée un array contenant les coefficients d'ajustement du modèle affine"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "G = 6.67e-11 #On définit la valeur de la constante de gravitation\n",
    "M = (4*(np.pi)**2)/(G*droite[0]) #On calcule la masse du Soleil grasse à la troisième loi de Newton"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Etape 5 : Affichage du graphique"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "La masse du Soleil est  2.650e+35 kg\n",
      "Coeff directeur :  2.233e-24\n",
      "Ordonnee à l'origine :  4.885e+07\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEfCAYAAABF6WFuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAnrUlEQVR4nO3de7xVc/7H8ddHyl0hjkuNGDHdFB3RuBUZuZUZBhn3IbdmkMsv18pcDGFcKoka5BIVKaUknYqK7ndyhJQQCUcl1ef3x3c3s89x6lzaZ619eT8fj/Ow19prrf3pq/bnfNf3uz5fc3dEREQ22SbuAEREJL0oMYiISDFKDCIiUowSg4iIFKPEICIixSgxiIhIMUoMIiJSjBKDSAlm1sLMJpvZBDN7wcyqp0FMeWY2yczGm9lbZrZP3DFJ9jI94CZSXOJLd5W7rzGze4Dp7j445piqAe7uG83sEqCOu/89zpgke6nHIFKCuy939zWJzXXAxlRc18z2NLMxZvatmfU3s3vM7PpyxrTB3TfFsQswP+m675lZo1TEKALqMUgOMrOLgV6JzW2B6sCmRLAeyHP3n8xsf2AgcJy7/5yCz30Q2N7drzGzPYFZwEFJSais85sBjwO1gN+5+6eJ/ecA57r7WVsbowgoMUiOM7O7gObu3r7E/l2B14Ar3P2DFH3WHOA6dx9nZjcDB7v7FZW4zjnACe5+VWJ7e+BzoKG7f5GKWCW36VaS5LpmwOzkHWa2LaGn0D0VScHMapjZd0ATYLiZzQVOAcaXOG4vMxtmZl+a2Q9mNjyRoDCzGkmHfges3rTh7muB6cDJWxurCCgxiDQj3NJJ1gE4ErjTzArM7NySJ5nZa2a2ajM/ryUf6+7rgJbAV+6+s7s3ISSJkklnV+BR4FfA/kBt4MpNcSZmSY0Drgd6lDh3IdC0In9wkc3ZNu4AROKS+G28HiUSg7sPAAZs6Vx3P72CH9eM4j2TWsAPJa5ZCBQmNn8yszHAbon33gOO28L1fwA0hVVSImt6DIlZHl+Z2bxyHHucmc0ws/VmdnbS/maJ+evzzWxOab8pSlZpSvhC/TiCz2pG8cTwLWF20X+Z2R/N7J3E3+NVQBdgUTmvvwuwauvDFMmixAA8BbQt57FLgEuA50vsXw1c5O6NEtd6yMxqpSg+ST/NgDleiRkYZva6mRVt5uf1Uk5pSvHEMAc4OOl6JwD3Em4T7Uu4jfQVv7zNtTkNSlxfpNKyJjG4+wRgZfI+M/u1mY0ys+lmNtHMfpM49hN3n0OJ+enuvsjdP0y8/pzwD3PPaP4EEoNmlP+Ltxh3PyUxXlDazymlnFIyMYwEji/x/meJY3YD+gN7AQvKiiUxK6k5MKYyfxaRkrImMWxGX+Av7t4cuAnoXd4TzawFUAP4qIpik/g1pZKJoSLMbG/Cl/37SbufAU41sx0S288RnqdYSZgm+yGwIDFwXZYzgILELzMiWy1rB5/NbGfgt8AgM9u0e7tynrsPYfDx4qSnTSXLuHt+RJ/zBSX+7rn712b2DGHW0UPu/hXh72uyv5XzI24C/rzVgYokZG1iIPSGVrl7s4qclJipMgK43d2nVEVgIgDufluKrnNkKq4jsknW3kpy9++Bj83sjwAWbHGed+IholeAZ+IumiYiEpesKYlhZi8ArQizOb4EugJvAY8R5ndXBwa6+91mdgQhAewGrAW+cPdGZnYB8B+SCpQBl7j7rKj+HCIiccuaxCAiIqmRtbeSRESkcpQYRESkmKyYlVS7dm2vV69ehc/78ccf2WmnnVIfUAZRG6gNQG0AudkG06dP/9rdf/EQb1Ykhnr16jFt2rQKn1dQUECrVq1SH1AGURuoDUBtALnZBmb2aWn7dStJRESKUWIQEZFilBhERKQYJQYRESlGiUFEJFN161Yll1ViEBHJVN27V8lllRhERDLQzz/DM1xIVVQ1UmIQEckk3box3ZqTX2M2F/MMY7dpA2Ypva2UFQ+4iYjkgjVroNuabty/TTfy8uCV5WfSxt9M+eeoxyAikgEmTICmTeG+++DSS2HBAjiTV6vks5QYRETS2PffwzXXwPHHw/r18Oab8OSTUKsW0LVrlXymEoOISJoaORIaN4Y+feCGG2DuXDjxxKQDqmi6qsYYRETSzNdfh0Tw7LPQsCFMmgRHHRXd56vHICKSJtzhxRdDMhg4EO66C2bMiDYpQMSJwcz6m9lXZjZvM++bmT1iZoVmNsfMDo8yPhGRuHz+OZx5Jpx3Huy/P0yfHp5f22676GOJusfwFNB2C++fAtRP/HQEHosgJhGR2LiHweSGDeGNN6BHD5g8GQ49NL6YIk0M7j4BWLmFQ9oDz3gwBahlZvtEE52ISLQWL4Y2beCKK6BZszC4fNNNsG3Mo7/pNvi8H/BZ0vbSxL7lJQ80s46EXgV5eXkUFBRU+MOKiooqdV42URuoDUBtANG2wYYN8PLLdejX7wCqVXM6d/6I005bztKlsHRpJCFsUbolhnJz975AX4D8/HyvzJJ8ubiUX0lqA7UBqA0gujaYNw8uvxzefRdOOy1MRa1T5xDgkCr/7PJKt1lJy4C6Sdt1EvtERDLaunVhMPnww+Gjj+D552H4cKhTJ+7IfindEsMw4KLE7KSjgO/c/Re3kUREMsnUqdC8eXge7eyzQzmLDh1C7bt0FOmtJDN7AWgF1DazpUBXoDqAu/cBRgKnAoXAauDSKOMTEUml1avDswj//jfssw8MGwZnnBF3VGWLNDG4e4cy3nfg2ojCERGpMgUFYSzho4/gyivh3nuhZs24oyqfdLuVJCKS0b77LiSC1q3D9ltvhQHmTEkKoMQgIpIyw4dDo0bhgbWbboI5c/6XIDKJEoOIyFZasQLOPx/atYPddgtPLvfoATvuGHdklaPEICJSSe5h2mmDBjB4cJiOOn06tGgRd2RbJ2MfcBMRidPSpXD11fDaayER9OsX1k7IBuoxiIhUwMaN8Pjjoejd2LHw4INhvYRsSQqgHoOISLkVFoaCdwUFcMIJ8MQTcOCBcUeVeuoxiIiUYf16uP9+aNIkLJzzxBNh7eVsTAqgHoOIyBbNmQN//jNMmxZmHfXuDfvtF3dUVUs9BhGRUvz0E3TtGmocffppWGpz6NDsTwqgHoOIyC9MmRJ6CQsWwAUXhFpHtWvHHVV01GMQEUn48Ufo3Bl++1v4/nsYMQIGDMitpADqMYiIADB9ei0uuww+/jg8n/Cvf8Guu8YdVTyUGEQkp61aFeoa9evXjPr1Yfx4OO64uKOKl24liUjOevXV8KDaf/4D5523hNmzlRRAiUFEctCXX8K558KZZ8Kee4b1l6+8cjE77BB3ZOlBiUFEcoY7PPts6CUMHQp//3t4PiE/P+7I0ovGGEQkJyxZAlddBa+/Di1bhqJ3DRrEHVV6Uo9BRLLaxo3haeVGjcLA8sMPw8SJSgpboh6DiGStRYvCussTJ0KbNtC3LxxwQNxRpT/1GEQk66xfD/feC4ceCnPnQv/+8MYbSgrlpR6DiGSV2bPhsstCFdTf/x569YJ99ok7qsyiHoOIZIW1a+GOO8IMo6VLYdAgGDJESaEy1GMQkYw3aVIoevf++3DxxWFVtd13jzuqzKUeg4hkrKIiuO46OOYYWL0aRo2Cp55SUthaSgwikpHeeCOss/zII3DttTBvHpx8ctxRZQclBhHJKN9+C5deGpLA9tuHqaiPPgq77BJ3ZNlDiUFEMsbLL4dyFgMGwK23wqxZ4TaSpFbkicHM2prZB2ZWaGZdSnn/V2Y2zsxmmtkcMzs16hhFJL188QWcfTacdRbsvTdMnQr//GfoMUjqRZoYzKwa0As4BWgIdDCzhiUOuwN4yd0PA84DekcZo4ikD3d4+unQS3jttZAM3nsPDjss7siyW9Q9hhZAobsvdvd1wECgfYljHNi0blJN4PMI4xORNPHJJ9C2LVxySUgMs2aF20fVq8ccWA6IOjHsB3yWtL00sS9ZN+ACM1sKjAT+Ek1oIpIONm4Mg8mNG8M774TXEybAb34Td2S5w9w9ug8zOxto6+6XJ7YvBI50905Jx3ROxPWAmbUE+gGN3X1jiWt1BDoC5OXlNR84cGCF4ykqKmLnnXeu9J8nG6gN1AaQPm2wZMmO9OhxCPPm1eSII1bSufMH7L33T5F8drq0QZRat2493d1/uRqFu0f2A7QERidt3wrcWuKY+UDdpO3FwF5bum7z5s29MsaNG1ep87KJ2kBt4B5/G6xb5/6Pf7jXqOG+227uTz/tvnFjtDHE3QZxAKZ5Kd+pUd9KmgrUN7MDzKwGYXB5WIljlgAnAphZA2B7YEWkUYpIZGbOhBYt4PbboV07WLgQLroIzOKOLHdFmhjcfT3QCRgNLCTMPppvZnebWbvEYTcCV5jZbOAF4JJEZhORLLJ2bRhMPuKIMB11yJBQ+C4vL+7IJPIieu4+kjConLzvrqTXC4Cjo45LRKLz9tuh6N2iReEp5gcegN12izsq2URPPotIZH74ATp1gmOPhXXrQr2j/v2VFNKNEoOIRGLUqDAFtXdv+Otfw8pqJ50Ud1RSGiUGEalS33wT1kg45RTYccdwG+nhhyHHZoZmFCUGEakS7jB4cHhq+fnnw+pqs2bBb38bd2RSFq3gJiIpt3x5WCPhlVegefMwltC0adxRSXmpxyAiKeMeBpMbNIDXX4d774UpU5QUMo16DCKSEh9/DB07wptvhllHTz4JBx8cd1RSGeoxiMhW2bAhDCY3bgzvvhtmHRUUKClkMvUYRKTSFiyAyy+HyZPDrKPHH4e6deOOSraWegwiUmE//wx//3tYMGfRInj2WRgxQkkhW6jHICIVMm1aKGcxZw6cey488gjstVfcUUkqqccgIuWyZg3ccgsceSSsWAFDh8LAgUoK2Ug9BhEp0/jxYSyhsDD8t0cPqFUr7qikqqjHICKb9f33cPXV0KpVWHLzzTfhiSeUFLKdEoOIlGrkSGjUCPr2hc6dw5jCiSfGHZVEQYlBRIr5+mu44AI47TTYdVeYNCmsl7DTTnFHJlFRYhARIJSzePHFUPTuxReha1eYMSMMNktu0eCziLBiRQ3OPBOGDYP8fBg7Fpo0iTsqiYsSg0gOcw81jW64oQUbNsD998N118G2+mbIafrfL5KjPvoIrrgCxo2DZs1+YNCg3TjooLijknSgxCCSYzYVvbvjDqhePdQ3Ouig2Rx0UKu4Q5M0ocFnkRwyb15YQe3GG8PU0/nzQ6nsbfRNIEn010EkB6xbB927w+GHw+LFYanNYcOgTp24I5N0pFtJIlnuvfdC0bt58+D88+Ghh2DPPeOOStKZegwiWWr1arjpJmjZEr79FoYPh+eeU1KQsqnHIJKFxo0Lxe4WL4YrrwxrL9esGXdUkinUYxDJIt99FxLBCSeAWUgQffooKUjFKDGIZInhw0M5iyefDLeQ5swJVVFFKkqJQSTDrVgRBpXbtYM99oApU8J6CTvuGHdkkqkiTwxm1tbMPjCzQjPrspljzjGzBWY238yejzpGkUzgHqadNmgAgweH6ajTpsERR8QdmWS6SAefzawa0As4CVgKTDWzYe6+IOmY+sCtwNHu/q2ZaeFAkRI++ywsoDNiRKh+2q9fWDtBJBWi7jG0AArdfbG7rwMGAu1LHHMF0MvdvwVw968ijlEkbW3cGEpYNGoUBpb//W945x0lBUmtqBPDfsBnSdtLE/uSHQwcbGbvmNkUM2sbWXQiaezDD8Nso6uughYtYO5cuP56qFYt7sgk21TqVpKZ7QSsdfcNKY4HQkz1gVZAHWCCmTVx91UlYugIdATIy8ujoKCgwh9UVFRUqfOyidog/dtgwwZj8OA69O9fj+rVnZtuKuTUU79gyRJYsiQ1n5HubRAFtUESdy/zh9CzOB8YAXxF+K3/K2AB0AM4qJzXaQmMTtq+Fbi1xDF9gEuTtscCR2zpus2bN/fKGDduXKXOyyZqg/Rug9mz3fPz3cG9fXv3Zcuq5nPSuQ2ikottAEzzUr5Ty3sraRzw68QX+d7uXtfd9wKOAaYA95rZBeW4zlSgvpkdYGY1gPOAYSWOGUroLWBmtQm3lhaXM06RrPDTT3DXXdC8OXz6aVhq85VXYN99445MckF5byW1cfefS+5095XAEGCImVUv6yLuvt7MOgGjgWpAf3efb2Z3EzLXsMR7vzOzBcAG4GZ3/6accYpkvClTQtG7BQvgwgvDAPMee8QdleSSMhODmZ0EnGNmPd19tpl1dPe+JY8rLXGUxt1HAiNL7Lsr6bUDnRM/Ijnjxx/D4jkPPxzKYY8cCaecEndUkovK02O4DLgauMPM9gCaVWlEIjlo7NiwzObHH8M118A998Cuu8YdleSq8owx/ODuq9z9JuB3gJ6rFEmRVatCFdQ2bWDbbWH8eOjVS0lB4lWexDBi0wt37wI8U3XhiOSOoUND0bunnoL/+z+YPRuOOy7uqETKcSvJ3V8tsf1o1YUjkv2+/BL+8hcYNAiaNg1VUZs3jzsqkf+p0JPPZjbBzHZNvL7KzK5PTDsVkTK4w4ABoZfw6qvwj3/A1KlKCpJ+KloSo6a7f29mzQk1jXYDnkh9WCLZZckSOO00uOgiOOQQmDULbrsNqpc5yVskehVNDD+b2bbARcC97t4VUPkukc3YuBF69w5F7iZMgEcegYkTQ6lskXRV0VpJjwCzge2BTWsp7JzSiESyxAcfhCmoEyfCSSdB375Qr17cUYmUrUI9Bnd/BjgSaOzua8zsIGBylUQmkqHWr4d//SsMLM+dC//5D4weraQgmaNcPQYzs8QTybh70ab97l4IXFryGJFcNWtWKGcxYwb84Q/Qsyfss0/cUYlUTLmL6JnZX8zsV8k7zayGmZ1gZk8DF6c+PJHMsHYt3H475OfDsmVhqc0hQ5QUJDOVd4yhLaE0xgtmdgCwijDOUA14A3jI3WdWSYQiaW7SpNBLeP99uPhiePBB2H33uKMSqbxyJQZ3Xwv0BnonqqjWBtZ4icVzRHJJUVGYctqzJ9StC6NGwcknxx2VyNar8ApuiSqqy6sgFpGM8cYb0LFjeD7h2mvhn/+EXXaJOyqR1ChzjMHMTjKzJ8ysWWK7Y5VHJZKmVq6ESy8NPYPttw9TUR99VElBskt5Bp8vA24GLjCzE1DZbclRQ4aEchYDBoRbSLNmwdFHxx2VSOqp7LZIGb74As4+O/zsuy9MmxbqHG2/fdyRiVSNcpXdNrMHQGW3Jbe4h5LYDRvCa6+FxXPefReaNYs7MpGqVWZiSJTdbp20rbLbkvU++QTatg3jCY0ahbUSunRR0TvJDRUtoieS1TZuDIPJjRuH5xN69gyrqh1ySNyRiUSnvNNVm5rZx8BcYF7Sfxe6+/qqCk4kSgsXhmU2J00KvYU+fWD//eOOSiR65e0xzAGOBnoC3xAGof8DfGNm86ooNpFIrF9v/POfYezg/ffhmWdg5EglBcld5X7Azd0/Bz4nlMAAQuE84KAqiEskEjNmwNVXH05hIfzxj+E2Ul5e3FGJxKu8PYaepe304MMUxiMSiTVr4NZboUULWLmyBi+/DC+9pKQgAuWvldSvqgMRicrbb4eid4sWwWWXwZlnTuWMM46JOyyRtKFZSZIzfvgBOnWCY4+FdetgzBjo1w922UXzJ0SSKTFITnj99fA8Qu/ecP31YWW1Nm3ijkokPVW4uqpIJvnmG7jhhlDfqEEDeOcdaNky7qhE0pt6DJKV3GHQoFDO4oUX4M47YeZMJQWR8og8MZhZWzP7wMwKzazLFo47y8zczPKjjE8y3/LlYb3lc84JC+hMmwZ33w3bbRd3ZCKZIdLEYGbVgF7AKUBDoIOZNSzluF2A64B3o4xPMps79O8fbhmNGgX33QdTpkDTpnFHJpJZou4xtAAK3X2xu68DBgLtSznub8C9wNoog5PMtXgxnHRSmIbatCnMmQM33wzbahRNpMKi/mezH/BZ0vZS4MjkA8zscKCuu48ws5s3d6HESnIdAfLy8igoKKhwMEVFRZU6L5tkehts2ACvvFKHfv0OYJttnBtuWMzpp3/OsmWwbFn5rpHpbZAKagO1QbK0+n3KzLYBHgQuKetYd+8L9AXIz8/3Vq1aVfjzCgoKqMx52SST22DBgtBDmDIFTj01FL2rW/dg4OAKXSeT2yBV1AZqg2RR30paBtRN2q6T2LfJLkBjoMDMPgGOAoZpAFqSrVsHf/sbHHYYfPghPPtsWEinbt2yzxWRskXdY5gK1DezAwgJ4Tzg/E1vuvt3QO1N22ZWANzk7tMijlPS1LRpoZcwZw6cdx48/DDstVfcUYlkl0h7DIm1GzoBo4GFwEvuPt/M7jazdlHGIpll9Wq45RY48kj4+mt49dXwfIKSgkjqRT7G4O4jgZEl9t21mWNbRRGTpLfx48MCOoWFcMUVYRpqrVpxRyWSvfTks6St77+Hq6+GVq3Ckptjx0LfvkoKIlVNiUHS0ogRoehd377QuXMoenfCCXFHJZIblBgkrXz9NVxwAZx+OtSsGdZffuAB2HHHuCMTyR1KDJIW3GHgwFDO4qWXoGvXsOzmkUeWfa6IpFZaPeAmuWnZMrjmGhg2DI44Iiye06RJ3FGJ5C71GCQ27vDEE6E09pgx4ZbR5MlKCiJxU49BYvHRR2Hq6bhx0Lp1SBC//nXcUYkIqMcgEduwAR58MPQKpk8Ps47GjlVSEEkn6jFIZObNC+Us3nsPzjgDHnsM9tsv7qhEpCT1GKTKrVsH3bvD4YeHdRNeeCGUtFBSEElP6jFIlXrvvdBLmDcP/vQneOghqF27zNNEJEbqMUiVWL0abrwRWraEVatCWexnn1VSEMkE6jFIyo0bF4reLV4MV10F994Lu+4ad1QiUl7qMUjKfPcddOwYahptsw0UFIQBZiUFkcyixCApMXx4eFCtXz+4+WaYPRuOPz7uqESkMpQYZKusWAEdOkC7drDHHvDuu2G9BBW9E8lcSgxSKe7w3HOh6N2QIXD33WHZzXytzi2S8TT4LBX22WdhAZ0RI+Coo8Lto4YN445KRFJFPQYpt40boU+fsIDOuHHhmYS331ZSEMk26jFIuXz4YSh6N348nHhiqHF04IFxRyUiVUE9Btmi9euhRw849FCYNSvcNhozRklBJJupxyCbNXt2KGcxfTq0bw+9e8O++8YdlYhUNfUY5Bd++gnuvDPMMPrss7DU5iuvKCmI5Ar1GKSYyZNDL2HhQrjoorB2wh57xB2ViERJPQYB4Mcf4frr4eijoagIRo6Ep59WUhDJReoxCG++GWYcffIJXHst3HMP7LJL3FGJSFyUGHLYqlVw332H8PrrUL8+TJgAxx4bd1QiEjfdSspRQ4eGB9NGj96bLl3CDCQlBREBJYac8+WXcM458PvfQ14ePPbYdO65B3bYIe7IRCRdRJ4YzKytmX1gZoVm1qWU9zub2QIzm2NmY81s/6hjzEbu8Mwzoejdq6/CP/4Rlt08+OCiuEMTkTQTaWIws2pAL+AUoCHQwcxKVtqZCeS7+6HAYOC+KGPMRkuWwKmnwsUXh8QwezbcdhtUrx53ZCKSjqLuMbQACt19sbuvAwYC7ZMPcPdx7r46sTkFqBNxjFlj40bo1SsUvZs4ER55JPz3N7+JOzIRSWfm7tF9mNnZQFt3vzyxfSFwpLt32szxPYEv3P3vpbzXEegIkJeX13zgwIEVjqeoqIidd965wudlgiVLduD++w9h7txa5Oev5MYbF7H33mt/cVw2t0F5qQ3UBpCbbdC6devp7v6LVVTSdrqqmV0A5AOlLhDp7n2BvgD5+fneqlWrCn9GQUEBlTkvnf38MzzwAHTrFlZRe+opuOii3TE7qtTjs7ENKkptoDYAtUGyqBPDMqBu0nadxL5izKwNcDtwvLv/FFFsGW/mzFDOYuZMOOss6NkT9t477qhEJNNEPcYwFahvZgeYWQ3gPGBY8gFmdhjwONDO3b+KOL6MtHYt3H47HHEEfP45DB4cfpQURKQyIu0xuPt6M+sEjAaqAf3dfb6Z3Q1Mc/dhQA9gZ2CQmQEscfd2UcaZSd55J/QSPvgALrkk3Ebaffe4oxKRTBb5GIO7jwRGlth3V9LrNlHHlImKisKU05494Ve/gtGj4Xe/izsqEckGevI5A73xBjRuHJJCp04wb56SgoikjhJDBlm5Ei69FE4+OZSw2PRsQo7NsBORKqbEkCGGDAlF7wYMCAPNM2eGtRNERFItbZ9jkGD58nC76OWX4bDDYNQoaNYs7qhEJJupx5Cm3MPDaQ0bwogR8K9/haJ3SgoiUtXUY0hDn3wCHTvCmDFwzDHw5JNwyCFxRyUiuUI9hjSycSM8+miYcTR5ciiAN368koKIREs9hjSxcCFcfjlMmgRt28Ljj4fnE0REoqYeQ8x+/jksmtOsGbz/flhMZ+RIJQURiY96DDGaMQMuuywsnHPOOeGZhLy8uKMSkVynHkMM1qyBLl2gRYuwBvMrr8CLLyopiEh6UI8hYhMnhrGERYtC8bsePWC33eKOSkTkf9RjiMgPP8C118Jxx4VxhTFjwjRUJQURSTdKDBF4/fWw7vJjj8H118PcudBGNWRFJE3pVlIV+uYbuOGGUN+oYcMwFfWo0lfYFBFJG+oxVAF3eOklaNAAXngB7rwzzEBSUhCRTKAeQ4p9/nkYSxg6FJo3hzffhEMPjTsqEZHyU48hRdyhX79wy2jUqDDbaMoUJQURyTzqMaTA4sWh6N3YsXD88WG20UEHxR2ViEjlqMewFTZsgIcegiZNQknsPn3grbeUFEQks6nHUEnz54cH1N59F047LSSFOnXijkpEZOupx1BB69bB3/4WVlMrLITnnoPhw5UURCR7qMdQAVOnhl7C3Llw3nmh6N2ee8YdlYhIaqnHUA6rV8Mtt4TnEL75Bl59NTyfoKQgItlIiWFLunWjoACaNg3TTy+/HBYsgHbt4g5MRKTqKDFsxnffwVXd96Z16/CMwltvhVXVataMOzIRkaqlxFCKESPCustPcAU33ghz5kDr1nFHJSISDSUGgG7dAFixAv506BxOPx1qLp3HZFpy/wPGjjvZf48REcl2kScGM2trZh+YWaGZdSnl/e3M7MXE+++aWb0qDahbN7x7dwYODOUsBr1/KN26wYyfGtOCqeE+krsSg4jkjEgTg5lVA3oBpwANgQ5m1rDEYX8GvnX3g4B/A/dWZUzLuj9Be16lQwc48MBQBbVrV6hRoyo/VUQkfUXdY2gBFLr7YndfBwwE2pc4pj3wdOL1YOBEM7OUR9KtGx+1HkBDFvAmbXiAzkx6rxqNB3f73zFdu6b8Y0VE0l3UiWE/4LOk7aWJfaUe4+7rge+APVIaRbdu0L0782lEc6YzlyZ05t9U63pn8VtGun0kIjkoY598NrOOQEeAvLw8CgoKyn9yq1bQqhXHfvsjD/zhRAwoGDcuvFeR62SBoqKiirVdFlIbqA1AbZAs6sSwDKibtF0nsa+0Y5aa2bZATeCbkhdy975AX4D8/Hxv1apVhYMpKCjAALp2pTLnZ4OCgoKc/bNvojZQG4DaIFnUt5KmAvXN7AAzqwGcBwwrccww4OLE67OBt9zdqyyirl11y0hEJEmkPQZ3X29mnYDRQDWgv7vPN7O7gWnuPgzoBwwws0JgJSF5VB0lBRGRYiIfY3D3kcDIEvvuSnq9Fvhj1HGJiEigJ59FRKQYJQYRESlGiUFERIpRYhARkWKsKmeCRsXMVgCfVuLU2sDXKQ4n06gN1AagNoDcbIP93f0Xa1FmRWKoLDOb5u75cccRJ7WB2gDUBqA2SKZbSSIiUowSg4iIFJPriaFv3AGkAbWB2gDUBqA2+K+cHmMQEZFfyvUeg4iIlKDEICIixeRsYjCztmb2gZkVmlmXuOOJmpn1N7OvzGxe3LHEwczqmtk4M1tgZvPN7Lq4Y4qamW1vZu+Z2exEG3SPO6a4mFk1M5tpZq/FHUs6yMnEYGbVgF7AKUBDoIOZNYw3qsg9BbSNO4gYrQdudPeGwFHAtTn4d+An4AR3bwo0A9qa2VHxhhSb64CFcQeRLnIyMQAtgEJ3X+zu64CBQPuYY4qUu08grHeRk9x9ubvPSLz+gfClUHL98azmQVFis3riJ+dmo5hZHeA04Mm4Y0kXuZoY9gM+S9peSo59Kcj/mFk94DDg3ZhDiVziFsos4CtgjLvnXBsADwG3ABtjjiNt5GpiEAHAzHYGhgDXu/v3cccTNXff4O7NCOuvtzCzxjGHFCkzOx34yt2nxx1LOsnVxLAMqJu0XSexT3KImVUnJIXn3P3luOOJk7uvAsaRe+NORwPtzOwTwi3lE8zs2XhDil+uJoapQH0zO8DMahDWlR4Wc0wSITMzwvriC939wbjjiYOZ7WlmtRKvdwBOAt6PNaiIufut7l7H3esRvgfecvcLYg4rdjmZGNx9PdAJGE0YdHzJ3efHG1W0zOwFYDJwiJktNbM/xx1TxI4GLiT8hjgr8XNq3EFFbB9gnJnNIfyyNMbdNV1TVBJDRESKy8keg4iIbJ4Sg4iIFKPEICIixSgxiIhIMUoMIiJpqiLFLs3sKjObm5hh9/am2l9mVs/M1iTNvutT5rU0K0lEJD2Z2XFAEfCMu2/xqXQz23XT0/tm1g64xt3bJkq+vFbW+cnUYxCpBDNrYGZ9zGywmV2dLteS7FJasUsz+7WZjTKz6WY20cx+kzg2uaTLTmxFQUQlBpFKcPeF7n4VcA7hYTkgPEFsZuMTpd0rfS0zq2FmE8xs21THLhmvL/AXd28O3AT03vSGmV1rZh8B9wF/TTrngMR6E+PN7NiyPkCJQaSSEt31EcDIpN2XAS+7+4atuVaiHPxY4NzURCvZIFH08bfAoERV3McJT7AD4O693P3XwP8BdyR2Lwd+5e6HAZ2B581s1y1+jsYYRLbMzM4m/Ga2A/AD8Ht3X5H0/gh3Py3xehJwvrt/ktgeBHxJWAinLvAn4ErgSGCiuxcrRVLiWk2Be9w910p1SJLkMYLEF/oH7r5PGedsA3zr7jVLea8AuMndp23ufHVTRco2zt0HA5hZV+AcM5sP/AHYjsRv+YmCjAduSgoJTYDJ7t7JzG4jFO5rBawAlprZdkDLktdKmAccUYV/Lskw7v69mX1sZn9090GJYpCHuvtsM6vv7h8mDj0N+BBCsURgpbtvMLMDgfrA4i19jhKDSNkuMbNzCV/cewO3uXsBUFDiuNrAqk0bZrY9UIuwEAyEwcB+7r488f4GYN1mrkXiH/I6M9slscqc5JhEsctWQG0zWwp0JfQ6HzOzOwir7g0EZgOdzKwN8DPwLXBx4jLHAXeb2c+ExYiucvctrt6oxCCyBWZ2EWEp2BPcvcjMJgCbq8S7Btg+absRMMPdN60M1hR4LHHdOsDnXva93O2AtZWNXzKbu3fYzFu/WDfD3a/bzDWGENYdKTcNPotsWRNgUiIpnEUY+Jtb2oHu/i1QLdFT2HTu7KRDDgXmJF43TXpdKjPbA/ja3X/eivhFKkyJQWTLngKuMbP3COtCL3b3H7dw/BvAMYnXTYBZ8N/bSjskkgcUTxKb05owU0kkUpqVJJJCZnY4cIO7X5iCa70MdHH3RVsfmUj5qccgkkLuPoOwKlq5H3ArTWKG01AlBYmDegwiIlKMegwiIlKMEoOIiBSjxCAiIsUoMYiISDFKDCIiUowSg4iIFKPEICIixSgxiIhIMf8PruoChCJFzIcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(acube, Tcarre, \"+\", color = \"red\")\n",
    "plt.grid()\n",
    "plt.xlabel('$a^3 (m^3)$')\n",
    "plt.ylabel('$T^2 (s^2)$')\n",
    "print (\"La masse du Soleil est \",\"%0.3e\" % M, \"kg\")\n",
    "print(\"Coeff directeur : \",\"%0.3e\" % droite[0])\n",
    "print(\"Ordonnee à l'origine : \",\"%0.3e\" % droite[1])\n",
    "Tcarremodelise = droite[0]*acube+droite[1]\n",
    "plt.plot(acube, Tcarremodelise, color = \"blue\")\n",
    "plt.title('$T^2 = f(a^3)$')\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}