natronics/coef_drag_exploration.ipynb

## coef_drag_exploration.ipynb
{
 "cells": [
  {
   "cell_type": "raw",
   "metadata": {},
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Explore Drag Coefficient Space\n",
    "\n",
    "We Define a Drag/Velocity Curve as \"True\"."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Imports\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "fineness = 40\n",
    "mach_range = (0,5)\n",
    "\n",
    "# List of Mach numbers from 0 to Max\n",
    "M = [i/fineness for i in range(fineness * mach_range[1] + 1)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# First part of the curve, subsonic incompressable flow\n",
    "def low_speed(m):\n",
    "    if m == 0: return 1.5\n",
    "    return (1 / (m * 160)) + 0.6\n",
    "\n",
    "a = [low_speed(m) for m in M[0:int(fineness * 0.8 + 1)]]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4AAAAEaCAYAAABAVjlbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8XGV97/HPLzeSEJBLOIWdRAFFEK0RRMAjKlZWGy4V\nz6nHAmrFgqEarFarckQxIPZ4OV6qohKBIlpBrJemisKy4uGlFQTkooBgBJQQNEK4h5CE/M4fa20z\nTGbfkp2Z2TOf9+u1XjNrrWfW/PbmEfnu51nPisxEkiRJktT7JnW6AEmSJElSexgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJGqOI+FBE/KLTdXSjiHheRGREPGeYNjvU\nbY5qZ22SJAOgJPWdiNgmIm6LiH0jYnZE/D4idoyIxfV/lA+3Hd+G+qbX33VKROxTvz9mmPZ/V7e5\nKyImNZ2bExHr6/MHbO3aR9Lwsw1uqyPijoi4OCKKLqqr1fbLTtUnSRo/BkBJ6j/vBKZl5s3Ah4Bf\nZub9wP8FdmvYfgJ8penYV5svFpUp41jfnsAG4PvAHsCjwI9H+MwaYBqwoOn4CcDycaxtvLyR6vf5\nLOD1wH3A9yLiQ8N9KCKmbY1iMnMNT/7nfFx96tkNxw7ZGt8tSWovA6Ak9Z/DgDsb3t8BkJmPZObv\nBjdgLfBY47HMfKwecXskIv4iIm6o27241bTIiDisHj3ateHYQRHxnxHxaD36eHFEzG342FHAFzLz\nGuAI4MzMvGuEn+kJ4ItUwWrweyZRBcBzmxtHxEcj4pf1CNxvI+LTETGrqc1BEVFGxMP1dmVE7N/U\n5n/Vo6mP1D/THiPUOeiB+vf5m8y8IjPfRBXM3x0RB9fXHhz9/OuIuCwiVgOnRsTUiDg3Im6PiMci\n4tcRcXpETG2q7V0RsaL+Gb8TEW+orze7VUFN/+zvrw+vbDh+b33dEyPimvp3sjIivhkRT2txyX0i\n4kcRsSYibo2IvxzuF1KPQp8dEb+r+8ZPI+IvRvn7lCSNkgFQkvpMZh6amYfW73fPzOM34zLTgTOA\ntwD7ADeO5kMR8TzgcuAHwP7An1ON3F06GGAy8yOZ+Xf1+7dk5rCjYg2+ABzZEDYXADOAf2/R9hHg\nRGDf+vVw4GMNde4H/BC4BzgUeD7waWBywzWeBhwP/DXwYmAXYMkoa23lM8Bq4FVNxz8CnEc1Gvcv\ndQ13A8dQjSD+I/Dm+nWw/uOAM4EPAvOBbwL/ZwtqazQFeG993SOB7YBvNU+/pfp9fgJ4HvA94OsR\n8YxWF6xHkEvgqcArgedSjT5/OyJeME51S5Ko/iUuSdJYTQZOzsyrBw9ExGg+dwpwcWZ+sOFzrwEe\nAF5OFRQ2S2b+KiJ+BLyBKuwsBM6nGqFsbnt6w+6dEfE+qvB2Un3sPcAvgNdnZtbHbmu6zDTgtfX0\nWSLiY8A5ETE5M5/YjPrXRsSvqKbANvpMZl7UdOy0pvr3oZq2ORjy3gF8MTPPqvd/FdWiLG8da10t\n6vx8w+7tEfFG4HaqMHpTw7lPZebXASLibVSjzW+l+qNBs6OAZwCH1NNRAT4ZES8BFlEFbUnSOHAE\nUJK0OZ4AfrYZn3sBcGw9ZfKRiHgE+D1VoNxrHOpaApwYEXOopo9+oVWjelrljyLinrqG84BZEbFT\n3eT5QNkQ/lr5zWD4q62g+sPqzvV3/Lrh57x2lPUH0PydP21R/5sj4up6CuYjwPupRiQHPQu4sulj\nPxllDcMXGHFwRHw7In4TEQ8DP69PNU8D/eP31b/H/6IaxWzlBVQjifc29Y2jGJ9+IUmqOQIoSdoc\na1qMcm2gCjCNpjbtTwLOoZoa2OzecajrG1RTNb8I/LgeFdynsUE9qvQVqimsb6cafXwpVXgcyyIr\nzSOLg8Ft8I+rL2fj/89uMgrZrF7g5RnAZU2nHm1q9zrg48C7qBbHeQh4LfDuIeoZN/X9g5cC/1F/\n5x+ogttPGdvvrtkkqsV6Xt7i3JoWxyRJm8kAKEkaLyuBXSNiUmZuqI/t39TmGuC5mblsaxRQT6M8\nn+p+uNcM0ezFwPLGaaAR8dqmNtcCRUScOsIo4HC13DnGj5wMzAS+NkK7lwBXZeanBg+0WHzmFuCF\nVCObgw4eYz2tzAe2B96ZmffU333EEG0Pphr1I6r5wS+kuq+ylWuo/pmt34zfmyRpDJwCKkkaLz8A\ndgTeFxFPj4hjqe7Da3QmsH9E/EtEHBARe0bEyyPiM00rgW6JU6kWZGm+b27QrcCciHhd/f1/S7UQ\nTKMPUS1E8sWIeH5EPCMijhnHBUl2iIhdI+KpEfGSiPgs8FHgw5m5yZTPFvXvHxFH1nX9I9VUyUYf\nA/4mqhVbn1H/jIPPUtySkcFfU03/fVv9uzucoReX+fuI+J8RsXddz97APw/R9t+pphT/R/1z7RER\nL4iId9T3iEqSxokBUJI0LjLzRqrVKN9AdV/YcVRhrLHNDVTPk9uF6jl/NwGfp5qR8tA41bE2M+9t\nGIVs9nWqQPLxus5X0jR9MjOvBV4GzAWuAK4D/h5YPx41Ut2beA/VwjIXALOBBZl5yig++2mqUcIv\nU41UPpcqWDfW/xXgfVT3Bt4I/FVDm82eUlmPzp0AvBq4GTgdeNsQzd9JNap3A/CXwKsy81dDXHc9\n1SIx/wl8jirkLqVagfWOza1XkrSp2MyZLZIkaQKJiH+iWtV0TqdrkSR1jvcASpLUYyJiJtVo7KXA\nY2x8BMNHO1mXJKnzHAGUJKnHRMQMqvvq9gNmUU2jPA/4xOY8o1CS1DsMgJIkSZLUJ1wERpIkSZL6\nRE/cA/j2v/3b3PnpT+90GVJL69evZ8qUnvifmnqMfVPdzP6pbmXfVLf67W9/e+/ZZ5+9y0jteqL3\n7vT445z6D/8AM2d2uhRpEytWrGBgYKDTZUibsG+qm9k/1a3sm+pWJ5100m9G0653poA+9linK5Ak\nSZKkrtY7AXD16k5XIEmSJEldrXcCoCOAkiRJkjSs3gmAa9Z0ugJJkiRJ6mq9EwCdAipJkiRJw+qd\nAOgUUEmSJEkalgFQkiRJkvqEAVCSJEmS+kRvBMBJkwyAkiRJkjSC3giAES4CI0mSJEkjaHsAjIjz\nImJlRPxihHYviIj1EfGqES86aZKPgZAkSZKkEXRiBPB8YMFwDSJiMvBh4LJRXdEpoJIkSZI0orYH\nwMy8Alg1QrO3AF8HVo7qmk4BlSRJkqQRTel0Ac0iYg7wP4CXAS8Ypt1CYCHAe17zGh75wx94aMWK\n9hQpjcGqVSP9vUPqDPumupn9U93KvqmJrusCIPBJ4N2ZuSEihmyUmUuAJQAfOOmknDVlCrMGBtpU\nojQ2A/ZNdSn7prqZ/VPdyr6piawbA+ABwEV1+JsNHBER6zPzW0N+wimgkiRJkjSirguAmbnH4PuI\nOB/49rDhD1wERpIkSZJGoe0BMCIuBA4FZkfEcuD9wFSAzPz85lwzIwyAkiRJkjSCtgfAzDx2DG2P\nH1VDRwAlSZIkaUSdeA7g+DMASpIkSdKIeiMADi4Ck9npSiRJkiSpa/VGAJw0CTZsgHXrOl2JJEmS\nJHWtngiAOan+Mdas6WwhkiRJktTFeiIAMvjAeJ8FKEmSJElD6o0AODgC6EIwkiRJkjQkA6AkSZIk\n9QkDoCRJkiT1id4IgIP3ABoAJUmSJGlIPREA0xFASZIkSRpRTwRARwAlSZIkaWS9EwCnTDEASpIk\nSdIweiMAAsyY4XMAJUmSJGkYvRUAHQGUJEmSpCH1VgBcs6bTVUiSJElS1+qtAOgUUEmSJEkaUu8E\nwJkznQIqSZIkScPonQA4fboBUJIkSZKG0TsB0BFASZIkSRpW7wRAVwGVJEmSpGG1NQBGxHkRsTIi\nfjHE+ddExI0R8fOI+K+ImD/qi7sIjCRJkiQNq90jgOcDC4Y5fwfw0sz8U+ADwJJRX9nHQEiSJEnS\nsKa088sy84qI2H2Y8//VsHslMHfUF3cKqCRJkiQNq60BcIxOAL471MmIWAgsBFi0aBH3z57Ntg88\nwL0rVrSrPmlUVq1a1ekSpJbsm+pm9k91K/umJrquDIAR8TKqAHjIUG0ycwn1FNHFixfnjgMDkMnA\nrrvCpN5Z20a9YWBgoNMlSC3ZN9XN7J/qVvZNTWRdl5Qi4rnAOcDRmXnfqD84Y0b16n2AkiRJktRS\nVwXAiHgq8A3gdZl525g+PHNm9ep9gJIkSZLUUlungEbEhcChwOyIWA68H5gKkJmfB04DdgY+GxEA\n6zPzgFFdfHAE0AAoSZIkSS21exXQY0c4fyJw4mZdfPr06tUAKEmSJEktddUU0C3iFFBJkiRJGlbv\nBECngEqSJEnSsHovAK5e3dk6JEmSJKlL9V4AdARQkiRJklrqnQC4zTbVq88BlCRJkqSWeicATppU\nrQTqFFBJkiRJaql3AiBUK4E6BVSSJEmSWuqtADh9ugFQkiRJkobQWwHQEUBJkiRJGlJvBcAZMwyA\nkiRJkjQEA6AkSZIk9QkDoCRJkiT1CQOgJEmSJPWJ3guAPgdQkiRJklrqvQDoCKAkSZIktdR7AXDN\nmk5XIUmSJEldqbcC4MyZTgGVJEmSpCH0VgB0CqgkSZIkDam3AuD06fDEE7BuXacrkSRJkqSu01sB\ncObM6tVRQEmSJEnaRNsDYEScFxErI+IXQ5yPiPhURCyLiBsjYv9RX3zGjOrVAChJkiRJm+jECOD5\nwIJhzh8O7FVvC4HPjfrKgwHQhWAkSZIkaRNtD4CZeQWwapgmRwMXZOVKYIeI2G1UF3cEUJIkSZKG\nNKXTBbQwB7irYX95feyexkYRsZBqhJBFixaxYsUKJj/0EDuuW8eDd93Fuh12aFvB0nBWrRru7x1S\n59g31c3sn+pW9k1NdN0YAEclM5cASwAWL16cAwMDMGUKTJ3KLttuCwMDHa5Q2mjA/qguZd9UN7N/\nqlvZNzWRdeMqoHcD8xr259bHRjZ9evXqFFBJkiRJ2kQ3BsClwN/Uq4EeDDyYmfeM9CHAewAlSZIk\naRhtnwIaERcChwKzI2I58H5gKkBmfh64BDgCWAasBt4w6otPnVpNAzUASpIkSdIm2h4AM/PYEc4n\nsGizv2DGDAOgJEmSJLXQjVNAt4wBUJIkSZJaMgBKkiRJUp8wAEqSJElSn+jNALh6daerkCRJkqSu\n05sB0BFASZIkSdpEbwbANWs6XYUkSZIkdZ3eC4AzZzoFVJIkSZJa6L0A6BRQSZIkSWqp9wLg9OkG\nQEmSJElqofcC4MyZBkBJkiRJaqH3AqBTQCVJkiSppd4MgGvWwIYNna5EkiRJkrpKbwZA8FEQkiRJ\nktTEAChJkiRJfaL3AuDMmdWrzwKUJEmSpCfpvQC47bbV60MPdbYOSZIkSeoyvRcA99wTJk2Cm27q\ndCWSJEmS1FV6LwDOmAH77AM33NDpSiRJkiSpq/ReAASYP98AKEmSJElNejcA3nMPrFzZ6UokSZIk\nqWu0PQBGxIKIuDUilkXEKS3OPzUiLo+I6yLixog4YsxfMn9+9eoooCRJkiT9UVsDYERMBs4CDgf2\nBY6NiH2bmr0XuDgz9wOOAT475i/aZRcYGDAASpIkSVKDdo8AHggsy8zbM3MtcBFwdFObBLav3z8F\nWLFZ3zR/Plx//ebWKUmSJEk9Z0qbv28OcFfD/nLgoKY2i4HLIuItwLbAYa0uFBELgYUAixYtYsWK\nJ+fE6U99KrO+/W3u+/WvyRkzxqd6aTOsWrWq0yVILdk31c3sn+pW9k1NdJsVAIuimFSW5YbxLqZ2\nLHB+Zn4sIl4IfCkinpOZT/q+zFwCLAFYvHhxDgwMPPkqL3sZnH02u913Hxx44FYqVRqdTfqn1CXs\nm+pm9k91K/umJrIxBcCiKPYBzgN2LopiJnAT8E9lWV4xykvcDcxr2J9bH2t0ArAAIDN/EhHTgdnA\n2Jb03HNPmDWrug/QAChJkiRJY74H8BPAyWVZ7l2W5TzgncCZRVH8+Sg/fzWwV0TsERHTqBZ5WdrU\n5rfAywEi4lnAdOAPY6wTJk2C5z7XhWAkSZIkqTbWADi9LMufDe6UZflz4FXAJo9zaCUz1wMnA5cC\nt1Ct9nlTRJwREa+om70DeGNE3ABcCByfmTnGOivz58ONN8KGrTVbVZIkSZImjrHeA7hJ+7IsVxZF\nMerrZOYlwCVNx05reH8z8KIx1tXa/PmwejUsWwbPfOa4XFKSJEmSJqqxBsAVRVFcC6wHbqO6B/Am\noDuX2Xz2s6upoDfcYACUJEmS1PeGnQJaFMXORVHsP7hfluVfl2X5fOAQ4CPAb4CD2dxn9W1tM2bA\nPvv4PEBJkiRJYuQRwPcDU4E3NR4sy3JdURR31/cAdrf58+HyyztdhSRJkiR13EiLwBwJfGCIc6cV\nRXHmONcz/ubPh9/9DlaO7SkSkiRJktRrRgqAs8uyHGp652eBo8e5nvE3f371+sMfdrQMSZIkSeq0\nkQLgyqIo9mx1oizL24A541/SONtlFzjiCPjMZ2BFd96qKEmSJEntMFIA/BrVYi+bKIpiO2DtuFe0\nNbzrXbDddvD+9/tMQEmSJEl9a6QA+E/A04ui+H5RFAc2nTsV+MnWKWuczZoFp58O110HX/5yp6uR\nJEmSpI4YdhXQsiwfKYriUOCTwI+LolgJ/JaNUz//bOuWN44OOACOOw4+9zl44Qthr706XZEkSZIk\ntdVII4CUZflgWZZvAPagGvX7FvBu4NllWS7byvWNr5NPhnnz4H3vg7UTY/aqJEmSJI2XkZ4D+Edl\nWS4Hzt96pbTBtGnwgQ/A618Pp54Kp51W3RsoSZIkSX1gxBHAnrP33nDmmXDVVdWU0Btv7HRFkiRJ\nktQW/RcAAQ47DC68EHbeGU48Eb7wBVcHlSRJktTz+jMAAsyZA+ecA8cfD0uWwGtfC9/5jvcGSpIk\nSepZ/RsAAaZMgTe/uQqAO+5YPSfwqKPg7LPhD3/odHWSJEmSNK5GvQhMT9tvPzjrLLj9drjoIrjg\ngmp08FnPgkMOqbZ99oFJ/Z2XJUmSJE1sBsBGe+4J73lP9biIyy+HH/+4enD8kiXwlKfAc54D++67\ncdt5505XLEmSJEmjZgBsZfvt4eijq23dumql0Kuugptvhq9+FR56aGO73XffuM2ZA7vtBrvuWk0p\njejgDyFJkiRJT2YAHMnUqfD851cbQCbcfTfccgvceWe13XorXHYZrFmz8XPTpsGf/Ek1Sjh7drXt\ntBPssEO17bhjNaq4/fbVNm1aJ346SZIkSX3EADhWETB3brU1yoT774ff/a7a7rkHVq6Ee++ttmXL\n4L774OGHW1932rTqofSzZsG22258nTmz2mbM2Pg6YwZMn77pts021TZt2pPfe++iJEmSJDoQACNi\nAfDPwGTgnMz8UIs2rwYWAwnckJnHtbXIzRFRjfDttFN1f+BQnngCHnwQHnig2h5+uNoeeqh6ffRR\neOSRja/33QerV2/cHnsMHn98bLVNnlwFwalTq9fB94P7U6Zs3B98P2XKxm3q1OoajccmT954rPG1\ncWs8NmnSpu+bjw1ujfuN7SI2bTtpUnV88HWwXeN5SZIkSUCbA2BETAbOAgpgOXB1RCzNzJsb2uwF\n/G/gRZl5f0T8t3bWuNVNnrwxKG6uDRuqEPjYY9W008cf33Rbu3bj+3Xrqv3Bbd26jdvg/vr11bZ2\nbXXdwf3166vzTzzx5GNPPPHkY4373WKHHeD73+90FZIkSVLXaPcI4IHAssy8HSAiLgKOBm5uaPNG\n4KzMvB8gM1e2ucbuN2nSxqmg3Saz2gYDYeO2YUPr/cbX5mOD1xo8N9Q22C5z4/4UZzhLkiRJjdr9\nX8hzgLsa9pcDBzW1eSZARPyYapro4sz8XvOFImIhsBBg0aJFrFixYqsUrK1kcHrm1g5pXdAvVq1a\n1ekSpJbsm+pm9k91K/umJrpuHCKZAuwFHArMBa6IiD/NzAcaG2XmEmAJwOLFi3NgYKDddUqjZv9U\nt7JvqpvZP9Wt7JuayNq9QsbdwLyG/bn1sUbLgaWZuS4z7wBuowqEkiRJkqQt0O4AeDWwV0TsERHT\ngGOApU1tvkU1+kdEzKaaEnp7O4uUJEmSpF7U1gCYmeuBk4FLgVuAizPzpog4IyJeUTe7FLgvIm4G\nLgfemZn3tbNOSZIkSepFbb8HMDMvAS5pOnZaw/sE3l5vkiRJkqRx4lOyJUmSJKlPGAAlSZIkqU8Y\nACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgA\nJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU+0PQBGxIKIuDUilkXEKcO0+6uIyIg4oJ31SZIk\nSVKvamsAjIjJwFnA4cC+wLERsW+LdtsBbwWuamd9kiRJktTL2j0CeCCwLDNvz8y1wEXA0S3afQD4\nMLCmncVJkiRJUi+b0ubvmwPc1bC/HDiosUFE7A/My8zvRMQ7h7pQRCwEFgIsWrSIFStWbIVypS23\natWqTpcgtWTfVDezf6pb2Tc10bU7AA4rIiYBHweOH6ltZi4BlgAsXrw4BwYGtm5x0hawf6pb2TfV\nzeyf6lb2TU1k7Z4Cejcwr2F/bn1s0HbAc4AfRsSdwMHAUheCkSRJkqQt1+4AeDWwV0TsERHTgGOA\npYMnM/PBzJydmbtn5u7AlcArMvOaNtcpSZIkST2nrQEwM9cDJwOXArcAF2fmTRFxRkS8op21SJIk\nSVK/afs9gJl5CXBJ07HThmh7aDtqkiRJkqR+0PYHwUuSJEmSOsMAKEmSJEl9wgAoSZIkSX3CAChJ\nkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmS\nJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIk\nSX3CAChJkiRJfcIAKEmSJEl9ou0BMCIWRMStEbEsIk5pcf7tEXFzRNwYEf8ZEU9rd42SJEmS1Iva\nGgAjYjJwFnA4sC9wbETs29TsOuCAzHwu8G/AR9pZoyRJkiT1qnaPAB4ILMvM2zNzLXARcHRjg8y8\nPDNX17tXAnPbXKMkSZIk9aQpbf6+OcBdDfvLgYOGaX8C8N1WJyJiIbAQYNGiRaxYsWK8apTG1apV\nqzpdgtSSfVPdzP6pbmXf1ETX7gA4ahHxWuAA4KWtzmfmEmAJwOLFi3NgYKCN1UljY/9Ut7JvqpvZ\nP9Wt7JuayNodAO8G5jXsz62PPUlEHAacCrw0Mx9vU22SJEmS1NPafQ/g1cBeEbFHREwDjgGWNjaI\niP2As4FXZObKNtcnSZIkST2rrQEwM9cDJwOXArcAF2fmTRFxRkS8om72UWAW8LWIuD4ilg5xOUmS\nJEnSGLT9HsDMvAS4pOnYaQ3vD2t3TZIkSZLUD9r+IHhJkiRJUmcYACVJkiSpTxgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJ\nkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmS\nJKlPGAAlSZIkqU8YACVJkiSpT7Q9AEbEgoi4NSKWRcQpLc5vExFfrc9fFRG7t7tGSZIkSepFbQ2A\nETEZOAs4HNgXODYi9m1qdgJwf2Y+A/gE8OF21ihJkiRJvardI4AHAssy8/bMXAtcBBzd1OZo4Iv1\n+38DXh4R0cYaJUmSJKknTWnz980B7mrYXw4cNFSbzFwfEQ8COwP3NjaKiIXAQoDjjz/+8ZNOOukX\nW6toaUusXr169syZM+8duaXUXvZNdTP7p7qVfVPdav369XuPpl27A+C4ycwlwBKAiLgmMw/ocElS\nS/ZPdSv7prqZ/VPdyr6pbhUR15x77rkjtmv3FNC7gXkN+3PrYy3bRMQU4CnAfW2pTpIkSZJ6WLsD\n4NXAXhGxR0RMA44Blja1WQq8vn7/KuAHmZltrFGSJEmSelJbp4DW9/SdDFwKTAbOy8ybIuIM4JrM\nXAqcC3wpIpYBq6hC4kiWbLWipS1n/1S3sm+qm9k/1a3sm+pWo+qb4eCaJEmSJPWHtj8IXpIkSZLU\nGQZASZIkSeoTEz4ARsSCiLg1IpZFxCmdrkcaFBHnRcTKiPAZleoqETEvIi6PiJsj4qaIeGuna5IA\nImJ6RPw0Im6o++bpna5JahQRkyPiuoj4dqdrkRpFxJ0R8fOIuD4irhm27US+BzAiJgO3AQXVQ+Wv\nBo7NzJs7WpgERMRLgEeACzLzOZ2uRxoUEbsBu2XmzyJiO+Ba4JX+u1OdFhEBbJuZj0TEVOBHwFsz\n88oOlyYBEBFvBw4Ats/MozpdjzQoIu4EDsjMe0dqO9FHAA8ElmXm7Zm5FrgIOLrDNUkAZOYVVCvZ\nSl0lM+/JzJ/V7x8GbgHmdLYqCbLySL07td4m7l+q1VMiYi5wJHBOp2uRtsRED4BzgLsa9pfjf8RI\n0qhFxO7AfsBVna1EqtRT7K4HVgJlZto31S0+CbwL2NDpQqQWErgsIq6NiIXDNZzoAVCStJkiYhbw\ndeBtmflQp+uRADLzicx8HjAXODAinEKvjouIo4CVmXltp2uRhnBIZu4PHA4sqm9FammiB8C7gXkN\n+3PrY5KkYdT3V30d+NfM/Ean65GaZeYDwOXAgk7XIgEvAl5R32d1EfBnEfHlzpYkbZSZd9evK4Fv\nUt0q19JED4BXA3tFxB4RMQ04Blja4ZokqavVC22cC9ySmR/vdD3SoIjYJSJ2qN/PoFrk7ZedrUqC\nzPzfmTk3M3en+u/NH2TmaztclgRARGxbL+pGRGwL/Dkw5Cr0EzoAZuZ64GTgUqpFDC7OzJs6W5VU\niYgLgZ8+lpvJAAAEbUlEQVQAe0fE8og4odM1SbUXAa+j+gv29fV2RKeLkoDdgMsj4kaqP/KWmely\n+5I0vD8BfhQRNwA/Bb6Tmd8bqvGEfgyEJEmSJGn0JvQIoCRJkiRp9AyAkiRJktQnDICSJEmS1CcM\ngJIkSZLUJwyAkiRJktQnDICSpL5XFMX5RVGc2ek6hlIUxZ1FURzW6TokSRPflE4XIEnSSIqiuBMY\nAAbKsry34fh1wPOAPcqyvLNNtewO3AF8tyzLIxqOfxlYVpbl4nbUIUnS5nAEUJI0UdwBHDu4UxTF\nnwIzO1cOBxVF8d87+P1jVhSFf/iVpD7n/xFIkiaKLwF/A3y63n89cAHwx6mbRVEcWe8/HXgQOLdx\nRK4oikOAjwD7Ag8D7yvL8vz69I5FUXwHeAlwM3BcWZa/HqaejwAfBF7WfKIoiuOBE8uyPKThWAJ7\nlWW5rCiK84HVwB7Ai4EbgL8CTql/rt8Dx5ZleV3DZV9QFMWngN2AbwFvKstyTX3to+qfe/e69r8r\ny/LG+tydwOeA1wB7F0WxbVmW64f5uSRJPcwRQEnSRHElsH1RFM8qimIycAzw5aY2j1KFxB2AI4E3\nFUXxSoCiKJ4GfJcqQO5CNXX0+obPHgOcDuwILKMKd8P5LPDMLbg379XAe4HZwOPAT4Cf1fv/Bny8\nqf1rgL+gCrfPrD9LURT7AecBJwE7A2cDS4ui2Kbhs8dS/T52MPxJUn9zBFCSNJEMjgL+P+AW4O7G\nk2VZ/rBh98aiKC4EXko1YnYc8P2yLC+sz99Xb4O+WZblTwGKovhXNg1gzR6jColnAt/fjJ/lm2VZ\nXlt/3zeBN5dleUG9/1Xg5Kb2nynL8q76/Aepgux7gYXA2WVZXlW3+2JRFO8BDqb6PQF8avCzkqT+\nZgCUJE0kXwKuoJo6eUHzyaIoDgI+BDwHmAZsA3ytPj0PGG5K5+8a3q8GZo2innOAdxZF8ZejaNvs\n9w3vH2ux3/z9jQHuN1SL4gA8DXh9URRvaTg/reF882clSX3MAChJmjDKsvxNURR3AEcAJ7Ro8hXg\nM8DhZVmuKYrik1RTKqEKQQeOcz1ri6I4HfgAcFPDqUdpWKCmKIpdx+Hr5jW8fyqwon5/F/DBsiyH\nm7Ka4/D9kqQeYACUJE00JwA7lmX5aItVLbcDVtXh70CqaZ+X1ef+FXhPURSvBr4BPAWYV5bl9WyZ\nL1Et3rIA+FV97Abg2UVRPA/4JbB4C78DYFFRFN+mGp08FfhqffwLwDeLovg+8FOq4HkocEVZlg+P\nw/dKknqIi8BIkiaUsix/XZblNUOcfjNwRlEUDwOnARc3fO63VCOH7wBWUS0AM38c6nmi/q6dGo7d\nBpxBdW/gr4Afben3UI1uXgbcTjWV9cz6u64B3kg18nk/1QI2x4/D90mSelBkOitEkiRJkvqBI4CS\nJEmS1CcMgJIkSZLUJwyAkiRJktQnDICSJEmS1CcMgJIkSZLUJwyAkiRJktQnDICSJEmS1CcMgJIk\nSZLUJ/4/Rbn2HIrMChsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f5b148c7c88>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Chart\n",
    "fig, ax1 = plt.subplots(figsize=(15,4))\n",
    "plt.title(r\"“True” Mach-Drag Table\")\n",
    "plt.ylabel(r\"$C_D$\")\n",
    "plt.xlabel(r\"Mach Number\")\n",
    "\n",
    "plt.plot(M[:33], a, 'r-', lw=1.4, alpha=0.8, label=\"True Value\")\n",
    "\n",
    "plt.ylim([0, 1.5])\n",
    "plt.xlim([0, 5.0])\n",
    "#plt.legend(loc=1)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": []
  }
 ],
 "metadata": {
  "celltoolbar": "Edit 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.5.3+"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
	{
	"cells": [
	{
	"cell_type": "raw",
	"metadata": {},
	"source": []
	},
	{
	"cell_type": "markdown",
	"metadata": {},
	"source": [
	"# Explore Drag Coefficient Space\n",
	"\n",
	"We Define a Drag/Velocity Curve as \"True\"."
	]
	},
	{
	"cell_type": "code",
	"execution_count": 1,
	"metadata": {
	"collapsed": true
	},
	"outputs": [],
	"source": [
	"# Imports\n",
	"import matplotlib.pyplot as plt\n",
	"%matplotlib inline"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 2,
	"metadata": {},
	"outputs": [],
	"source": [
	"fineness = 40\n",
	"mach_range = (0,5)\n",
	"\n",
	"# List of Mach numbers from 0 to Max\n",
	"M = [i/fineness for i in range(fineness * mach_range[1] + 1)]"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 3,
	"metadata": {},
	"outputs": [],
	"source": [
	"# First part of the curve, subsonic incompressable flow\n",
	"def low_speed(m):\n",
	" if m == 0: return 1.5\n",
	" return (1 / (m * 160)) + 0.6\n",
	"\n",
	"a = [low_speed(m) for m in M[0:int(fineness * 0.8 + 1)]]"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 4,
	"metadata": {},
	"outputs": [
	{
	"data": {
	"image/png": "iVBORw0KGgoAAAANSUhEUgAAA4AAAAEaCAYAAABAVjlbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8XGV97/HPLzeSEJBLOIWdRAFFEK0RRMAjKlZWGy4V\nz6nHAmrFgqEarFarckQxIPZ4OV6qohKBIlpBrJemisKy4uGlFQTkooBgBJQQNEK4h5CE/M4fa20z\nTGbfkp2Z2TOf9+u1XjNrrWfW/PbmEfnu51nPisxEkiRJktT7JnW6AEmSJElSexgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJGqOI+FBE/KLTdXSjiHheRGREPGeYNjvU\nbY5qZ22SJAOgJPWdiNgmIm6LiH0jYnZE/D4idoyIxfV/lA+3Hd+G+qbX33VKROxTvz9mmPZ/V7e5\nKyImNZ2bExHr6/MHbO3aR9Lwsw1uqyPijoi4OCKKLqqr1fbLTtUnSRo/BkBJ6j/vBKZl5s3Ah4Bf\nZub9wP8FdmvYfgJ8penYV5svFpUp41jfnsAG4PvAHsCjwI9H+MwaYBqwoOn4CcDycaxtvLyR6vf5\nLOD1wH3A9yLiQ8N9KCKmbY1iMnMNT/7nfFx96tkNxw7ZGt8tSWovA6Ak9Z/DgDsb3t8BkJmPZObv\nBjdgLfBY47HMfKwecXskIv4iIm6o27241bTIiDisHj3ateHYQRHxnxHxaD36eHFEzG342FHAFzLz\nGuAI4MzMvGuEn+kJ4ItUwWrweyZRBcBzmxtHxEcj4pf1CNxvI+LTETGrqc1BEVFGxMP1dmVE7N/U\n5n/Vo6mP1D/THiPUOeiB+vf5m8y8IjPfRBXM3x0RB9fXHhz9/OuIuCwiVgOnRsTUiDg3Im6PiMci\n4tcRcXpETG2q7V0RsaL+Gb8TEW+orze7VUFN/+zvrw+vbDh+b33dEyPimvp3sjIivhkRT2txyX0i\n4kcRsSYibo2IvxzuF1KPQp8dEb+r+8ZPI+IvRvn7lCSNkgFQkvpMZh6amYfW73fPzOM34zLTgTOA\ntwD7ADeO5kMR8TzgcuAHwP7An1ON3F06GGAy8yOZ+Xf1+7dk5rCjYg2+ABzZEDYXADOAf2/R9hHg\nRGDf+vVw4GMNde4H/BC4BzgUeD7waWBywzWeBhwP/DXwYmAXYMkoa23lM8Bq4FVNxz8CnEc1Gvcv\ndQ13A8dQjSD+I/Dm+nWw/uOAM4EPAvOBbwL/ZwtqazQFeG993SOB7YBvNU+/pfp9fgJ4HvA94OsR\n8YxWF6xHkEvgqcArgedSjT5/OyJeME51S5Ko/iUuSdJYTQZOzsyrBw9ExGg+dwpwcWZ+sOFzrwEe\nAF5OFRQ2S2b+KiJ+BLyBKuwsBM6nGqFsbnt6w+6dEfE+qvB2Un3sPcAvgNdnZtbHbmu6zDTgtfX0\nWSLiY8A5ETE5M5/YjPrXRsSvqKbANvpMZl7UdOy0pvr3oZq2ORjy3gF8MTPPqvd/FdWiLG8da10t\n6vx8w+7tEfFG4HaqMHpTw7lPZebXASLibVSjzW+l+qNBs6OAZwCH1NNRAT4ZES8BFlEFbUnSOHAE\nUJK0OZ4AfrYZn3sBcGw9ZfKRiHgE+D1VoNxrHOpaApwYEXOopo9+oVWjelrljyLinrqG84BZEbFT\n3eT5QNkQ/lr5zWD4q62g+sPqzvV3/Lrh57x2lPUH0PydP21R/5sj4up6CuYjwPupRiQHPQu4sulj\nPxllDcMXGHFwRHw7In4TEQ8DP69PNU8D/eP31b/H/6IaxWzlBVQjifc29Y2jGJ9+IUmqOQIoSdoc\na1qMcm2gCjCNpjbtTwLOoZoa2OzecajrG1RTNb8I/LgeFdynsUE9qvQVqimsb6cafXwpVXgcyyIr\nzSOLg8Ft8I+rL2fj/89uMgrZrF7g5RnAZU2nHm1q9zrg48C7qBbHeQh4LfDuIeoZN/X9g5cC/1F/\n5x+ogttPGdvvrtkkqsV6Xt7i3JoWxyRJm8kAKEkaLyuBXSNiUmZuqI/t39TmGuC5mblsaxRQT6M8\nn+p+uNcM0ezFwPLGaaAR8dqmNtcCRUScOsIo4HC13DnGj5wMzAS+NkK7lwBXZeanBg+0WHzmFuCF\nVCObgw4eYz2tzAe2B96ZmffU333EEG0Pphr1I6r5wS+kuq+ylWuo/pmt34zfmyRpDJwCKkkaLz8A\ndgTeFxFPj4hjqe7Da3QmsH9E/EtEHBARe0bEyyPiM00rgW6JU6kWZGm+b27QrcCciHhd/f1/S7UQ\nTKMPUS1E8sWIeH5EPCMijhnHBUl2iIhdI+KpEfGSiPgs8FHgw5m5yZTPFvXvHxFH1nX9I9VUyUYf\nA/4mqhVbn1H/jIPPUtySkcFfU03/fVv9uzucoReX+fuI+J8RsXddz97APw/R9t+pphT/R/1z7RER\nL4iId9T3iEqSxokBUJI0LjLzRqrVKN9AdV/YcVRhrLHNDVTPk9uF6jl/NwGfp5qR8tA41bE2M+9t\nGIVs9nWqQPLxus5X0jR9MjOvBV4GzAWuAK4D/h5YPx41Ut2beA/VwjIXALOBBZl5yig++2mqUcIv\nU41UPpcqWDfW/xXgfVT3Bt4I/FVDm82eUlmPzp0AvBq4GTgdeNsQzd9JNap3A/CXwKsy81dDXHc9\n1SIx/wl8jirkLqVagfWOza1XkrSp2MyZLZIkaQKJiH+iWtV0TqdrkSR1jvcASpLUYyJiJtVo7KXA\nY2x8BMNHO1mXJKnzHAGUJKnHRMQMqvvq9gNmUU2jPA/4xOY8o1CS1DsMgJIkSZLUJ1wERpIkSZL6\nRE/cA/j2v/3b3PnpT+90GVJL69evZ8qUnvifmnqMfVPdzP6pbmXfVLf67W9/e+/ZZ5+9y0jteqL3\n7vT445z6D/8AM2d2uhRpEytWrGBgYKDTZUibsG+qm9k/1a3sm+pWJ5100m9G0653poA+9linK5Ak\nSZKkrtY7AXD16k5XIEmSJEldrXcCoCOAkiRJkjSs3gmAa9Z0ugJJkiRJ6mq9EwCdAipJkiRJw+qd\nAOgUUEmSJEkalgFQkiRJkvqEAVCSJEmS+kRvBMBJkwyAkiRJkjSC3giAES4CI0mSJEkjaHsAjIjz\nImJlRPxihHYviIj1EfGqES86aZKPgZAkSZKkEXRiBPB8YMFwDSJiMvBh4LJRXdEpoJIkSZI0orYH\nwMy8Alg1QrO3AF8HVo7qmk4BlSRJkqQRTel0Ac0iYg7wP4CXAS8Ypt1CYCHAe17zGh75wx94aMWK\n9hQpjcGqVSP9vUPqDPumupn9U93KvqmJrusCIPBJ4N2ZuSEihmyUmUuAJQAfOOmknDVlCrMGBtpU\nojQ2A/ZNdSn7prqZ/VPdyr6piawbA+ABwEV1+JsNHBER6zPzW0N+wimgkiRJkjSirguAmbnH4PuI\nOB/49rDhD1wERpIkSZJGoe0BMCIuBA4FZkfEcuD9wFSAzPz85lwzIwyAkiRJkjSCtgfAzDx2DG2P\nH1VDRwAlSZIkaUSdeA7g+DMASpIkSdKIeiMADi4Ck9npSiRJkiSpa/VGAJw0CTZsgHXrOl2JJEmS\nJHWtngiAOan+Mdas6WwhkiRJktTFeiIAMvjAeJ8FKEmSJElD6o0AODgC6EIwkiRJkjQkA6AkSZIk\n9QkDoCRJkiT1id4IgIP3ABoAJUmSJGlIPREA0xFASZIkSRpRTwRARwAlSZIkaWS9EwCnTDEASpIk\nSdIweiMAAsyY4XMAJUmSJGkYvRUAHQGUJEmSpCH1VgBcs6bTVUiSJElS1+qtAOgUUEmSJEkaUu8E\nwJkznQIqSZIkScPonQA4fboBUJIkSZKG0TsB0BFASZIkSRpW7wRAVwGVJEmSpGG1NQBGxHkRsTIi\nfjHE+ddExI0R8fOI+K+ImD/qi7sIjCRJkiQNq90jgOcDC4Y5fwfw0sz8U+ADwJJRX9nHQEiSJEnS\nsKa088sy84qI2H2Y8//VsHslMHfUF3cKqCRJkiQNq60BcIxOAL471MmIWAgsBFi0aBH3z57Ntg88\nwL0rVrSrPmlUVq1a1ekSpJbsm+pm9k91K/umJrquDIAR8TKqAHjIUG0ycwn1FNHFixfnjgMDkMnA\nrrvCpN5Z20a9YWBgoNMlSC3ZN9XN7J/qVvZNTWRdl5Qi4rnAOcDRmXnfqD84Y0b16n2AkiRJktRS\nVwXAiHgq8A3gdZl525g+PHNm9ep9gJIkSZLUUlungEbEhcChwOyIWA68H5gKkJmfB04DdgY+GxEA\n6zPzgFFdfHAE0AAoSZIkSS21exXQY0c4fyJw4mZdfPr06tUAKEmSJEktddUU0C3iFFBJkiRJGlbv\nBECngEqSJEnSsHovAK5e3dk6JEmSJKlL9V4AdARQkiRJklrqnQC4zTbVq88BlCRJkqSWeicATppU\nrQTqFFBJkiRJaql3AiBUK4E6BVSSJEmSWuqtADh9ugFQkiRJkobQWwHQEUBJkiRJGlJvBcAZMwyA\nkiRJkjQEA6AkSZIk9QkDoCRJkiT1CQOgJEmSJPWJ3guAPgdQkiRJklrqvQDoCKAkSZIktdR7AXDN\nmk5XIUmSJEldqbcC4MyZTgGVJEmSpCH0VgB0CqgkSZIkDam3AuD06fDEE7BuXacrkSRJkqSu01sB\ncObM6tVRQEmSJEnaRNsDYEScFxErI+IXQ5yPiPhURCyLiBsjYv9RX3zGjOrVAChJkiRJm+jECOD5\nwIJhzh8O7FVvC4HPjfrKgwHQhWAkSZIkaRNtD4CZeQWwapgmRwMXZOVKYIeI2G1UF3cEUJIkSZKG\nNKXTBbQwB7irYX95feyexkYRsZBqhJBFixaxYsUKJj/0EDuuW8eDd93Fuh12aFvB0nBWrRru7x1S\n59g31c3sn+pW9k1NdN0YAEclM5cASwAWL16cAwMDMGUKTJ3KLttuCwMDHa5Q2mjA/qguZd9UN7N/\nqlvZNzWRdeMqoHcD8xr259bHRjZ9evXqFFBJkiRJ2kQ3BsClwN/Uq4EeDDyYmfeM9CHAewAlSZIk\naRhtnwIaERcChwKzI2I58H5gKkBmfh64BDgCWAasBt4w6otPnVpNAzUASpIkSdIm2h4AM/PYEc4n\nsGizv2DGDAOgJEmSJLXQjVNAt4wBUJIkSZJaMgBKkiRJUp8wAEqSJElSn+jNALh6daerkCRJkqSu\n05sB0BFASZIkSdpEbwbANWs6XYUkSZIkdZ3eC4AzZzoFVJIkSZJa6L0A6BRQSZIkSWqp9wLg9OkG\nQEmSJElqofcC4MyZBkBJkiRJaqH3AqBTQCVJkiSppd4MgGvWwIYNna5EkiRJkrpKbwZA8FEQkiRJ\nktTEAChJkiRJfaL3AuDMmdWrzwKUJEmSpCfpvQC47bbV60MPdbYOSZIkSeoyvRcA99wTJk2Cm27q\ndCWSJEmS1FV6LwDOmAH77AM33NDpSiRJkiSpq/ReAASYP98AKEmSJElNejcA3nMPrFzZ6UokSZIk\nqWu0PQBGxIKIuDUilkXEKS3OPzUiLo+I6yLixog4YsxfMn9+9eoooCRJkiT9UVsDYERMBs4CDgf2\nBY6NiH2bmr0XuDgz9wOOAT475i/aZRcYGDAASpIkSVKDdo8AHggsy8zbM3MtcBFwdFObBLav3z8F\nWLFZ3zR/Plx//ebWKUmSJEk9Z0qbv28OcFfD/nLgoKY2i4HLIuItwLbAYa0uFBELgYUAixYtYsWK\nJ+fE6U99KrO+/W3u+/WvyRkzxqd6aTOsWrWq0yVILdk31c3sn+pW9k1NdJsVAIuimFSW5YbxLqZ2\nLHB+Zn4sIl4IfCkinpOZT/q+zFwCLAFYvHhxDgwMPPkqL3sZnH02u913Hxx44FYqVRqdTfqn1CXs\nm+pm9k91K/umJrIxBcCiKPYBzgN2LopiJnAT8E9lWV4xykvcDcxr2J9bH2t0ArAAIDN/EhHTgdnA\n2Jb03HNPmDWrug/QAChJkiRJY74H8BPAyWVZ7l2W5TzgncCZRVH8+Sg/fzWwV0TsERHTqBZ5WdrU\n5rfAywEi4lnAdOAPY6wTJk2C5z7XhWAkSZIkqTbWADi9LMufDe6UZflz4FXAJo9zaCUz1wMnA5cC\nt1Ct9nlTRJwREa+om70DeGNE3ABcCByfmTnGOivz58ONN8KGrTVbVZIkSZImjrHeA7hJ+7IsVxZF\nMerrZOYlwCVNx05reH8z8KIx1tXa/PmwejUsWwbPfOa4XFKSJEmSJqqxBsAVRVFcC6wHbqO6B/Am\noDuX2Xz2s6upoDfcYACUJEmS1PeGnQJaFMXORVHsP7hfluVfl2X5fOAQ4CPAb4CD2dxn9W1tM2bA\nPvv4PEBJkiRJYuQRwPcDU4E3NR4sy3JdURR31/cAdrf58+HyyztdhSRJkiR13EiLwBwJfGCIc6cV\nRXHmONcz/ubPh9/9DlaO7SkSkiRJktRrRgqAs8uyHGp652eBo8e5nvE3f371+sMfdrQMSZIkSeq0\nkQLgyqIo9mx1oizL24A541/SONtlFzjiCPjMZ2BFd96qKEmSJEntMFIA/BrVYi+bKIpiO2DtuFe0\nNbzrXbDddvD+9/tMQEmSJEl9a6QA+E/A04ui+H5RFAc2nTsV+MnWKWuczZoFp58O110HX/5yp6uR\nJEmSpI4YdhXQsiwfKYriUOCTwI+LolgJ/JaNUz//bOuWN44OOACOOw4+9zl44Qthr706XZEkSZIk\ntdVII4CUZflgWZZvAPagGvX7FvBu4NllWS7byvWNr5NPhnnz4H3vg7UTY/aqJEmSJI2XkZ4D+Edl\nWS4Hzt96pbTBtGnwgQ/A618Pp54Kp51W3RsoSZIkSX1gxBHAnrP33nDmmXDVVdWU0Btv7HRFkiRJ\nktQW/RcAAQ47DC68EHbeGU48Eb7wBVcHlSRJktTz+jMAAsyZA+ecA8cfD0uWwGtfC9/5jvcGSpIk\nSepZ/RsAAaZMgTe/uQqAO+5YPSfwqKPg7LPhD3/odHWSJEmSNK5GvQhMT9tvPzjrLLj9drjoIrjg\ngmp08FnPgkMOqbZ99oFJ/Z2XJUmSJE1sBsBGe+4J73lP9biIyy+HH/+4enD8kiXwlKfAc54D++67\ncdt5505XLEmSJEmjZgBsZfvt4eijq23dumql0Kuugptvhq9+FR56aGO73XffuM2ZA7vtBrvuWk0p\njejgDyFJkiRJT2YAHMnUqfD851cbQCbcfTfccgvceWe13XorXHYZrFmz8XPTpsGf/Ek1Sjh7drXt\ntBPssEO17bhjNaq4/fbVNm1aJ346SZIkSX3EADhWETB3brU1yoT774ff/a7a7rkHVq6Ee++ttmXL\n4L774OGHW1932rTqofSzZsG22258nTmz2mbM2Pg6YwZMn77pts021TZt2pPfe++iJEmSJDoQACNi\nAfDPwGTgnMz8UIs2rwYWAwnckJnHtbXIzRFRjfDttFN1f+BQnngCHnwQHnig2h5+uNoeeqh6ffRR\neOSRja/33QerV2/cHnsMHn98bLVNnlwFwalTq9fB94P7U6Zs3B98P2XKxm3q1OoajccmT954rPG1\ncWs8NmnSpu+bjw1ujfuN7SI2bTtpUnV88HWwXeN5SZIkSUCbA2BETAbOAgpgOXB1RCzNzJsb2uwF\n/G/gRZl5f0T8t3bWuNVNnrwxKG6uDRuqEPjYY9W008cf33Rbu3bj+3Xrqv3Bbd26jdvg/vr11bZ2\nbXXdwf3166vzTzzx5GNPPPHkY4373WKHHeD73+90FZIkSVLXaPcI4IHAssy8HSAiLgKOBm5uaPNG\n4KzMvB8gM1e2ucbuN2nSxqmg3Saz2gYDYeO2YUPr/cbX5mOD1xo8N9Q22C5z4/4UZzhLkiRJjdr9\nX8hzgLsa9pcDBzW1eSZARPyYapro4sz8XvOFImIhsBBg0aJFrFixYqsUrK1kcHrm1g5pXdAvVq1a\n1ekSpJbsm+pm9k91K/umJrpuHCKZAuwFHArMBa6IiD/NzAcaG2XmEmAJwOLFi3NgYKDddUqjZv9U\nt7JvqpvZP9Wt7JuayNq9QsbdwLyG/bn1sUbLgaWZuS4z7wBuowqEkiRJkqQt0O4AeDWwV0TsERHT\ngGOApU1tvkU1+kdEzKaaEnp7O4uUJEmSpF7U1gCYmeuBk4FLgVuAizPzpog4IyJeUTe7FLgvIm4G\nLgfemZn3tbNOSZIkSepFbb8HMDMvAS5pOnZaw/sE3l5vkiRJkqRx4lOyJUmSJKlPGAAlSZIkqU8Y\nACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgA\nJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU+0PQBGxIKIuDUilkXEKcO0+6uIyIg4oJ31SZIk\nSVKvamsAjIjJwFnA4cC+wLERsW+LdtsBbwWuamd9kiRJktTL2j0CeCCwLDNvz8y1wEXA0S3afQD4\nMLCmncVJkiRJUi+b0ubvmwPc1bC/HDiosUFE7A/My8zvRMQ7h7pQRCwEFgIsWrSIFStWbIVypS23\natWqTpcgtWTfVDezf6pb2Tc10bU7AA4rIiYBHweOH6ltZi4BlgAsXrw4BwYGtm5x0hawf6pb2TfV\nzeyf6lb2TU1k7Z4Cejcwr2F/bn1s0HbAc4AfRsSdwMHAUheCkSRJkqQt1+4AeDWwV0TsERHTgGOA\npYMnM/PBzJydmbtn5u7AlcArMvOaNtcpSZIkST2nrQEwM9cDJwOXArcAF2fmTRFxRkS8op21SJIk\nSVK/afs9gJl5CXBJ07HThmh7aDtqkiRJkqR+0PYHwUuSJEmSOsMAKEmSJEl9wgAoSZIkSX3CAChJ\nkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmS\nJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIkSX3CAChJkiRJfcIAKEmSJEl9wgAoSZIk\nSX3CAChJkiRJfcIAKEmSJEl9ou0BMCIWRMStEbEsIk5pcf7tEXFzRNwYEf8ZEU9rd42SJEmS1Iva\nGgAjYjJwFnA4sC9wbETs29TsOuCAzHwu8G/AR9pZoyRJkiT1qnaPAB4ILMvM2zNzLXARcHRjg8y8\nPDNX17tXAnPbXKMkSZIk9aQpbf6+OcBdDfvLgYOGaX8C8N1WJyJiIbAQYNGiRaxYsWK8apTG1apV\nqzpdgtSSfVPdzP6pbmXf1ETX7gA4ahHxWuAA4KWtzmfmEmAJwOLFi3NgYKCN1UljY/9Ut7JvqpvZ\nP9Wt7JuayNodAO8G5jXsz62PPUlEHAacCrw0Mx9vU22SJEmS1NPafQ/g1cBeEbFHREwDjgGWNjaI\niP2As4FXZObKNtcnSZIkST2rrQEwM9cDJwOXArcAF2fmTRFxRkS8om72UWAW8LWIuD4ilg5xOUmS\nJEnSGLT9HsDMvAS4pOnYaQ3vD2t3TZIkSZLUD9r+IHhJkiRJUmcYACVJkiSpTxgAJUmSJKlPGAAl\nSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJ\nkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmSJKlPGAAlSZIkqU8YACVJkiSpTxgAJUmS\nJKlPGAAlSZIkqU8YACVJkiSpT7Q9AEbEgoi4NSKWRcQpLc5vExFfrc9fFRG7t7tGSZIkSepFbQ2A\nETEZOAs4HNgXODYi9m1qdgJwf2Y+A/gE8OF21ihJkiRJvardI4AHAssy8/bMXAtcBBzd1OZo4Iv1\n+38DXh4R0cYaJUmSJKknTWnz980B7mrYXw4cNFSbzFwfEQ8COwP3NjaKiIXAQoDjjz/+8ZNOOukX\nW6toaUusXr169syZM+8duaXUXvZNdTP7p7qVfVPdav369XuPpl27A+C4ycwlwBKAiLgmMw/ocElS\nS/ZPdSv7prqZ/VPdyr6pbhUR15x77rkjtmv3FNC7gXkN+3PrYy3bRMQU4CnAfW2pTpIkSZJ6WLsD\n4NXAXhGxR0RMA44Blja1WQq8vn7/KuAHmZltrFGSJEmSelJbp4DW9/SdDFwKTAbOy8ybIuIM4JrM\nXAqcC3wpIpYBq6hC4kiWbLWipS1n/1S3sm+qm9k/1a3sm+pWo+qb4eCaJEmSJPWHtj8IXpIkSZLU\nGQZASZIkSeoTEz4ARsSCiLg1IpZFxCmdrkcaFBHnRcTKiPAZleoqETEvIi6PiJsj4qaIeGuna5IA\nImJ6RPw0Im6o++bpna5JahQRkyPiuoj4dqdrkRpFxJ0R8fOIuD4irhm27US+BzAiJgO3AQXVQ+Wv\nBo7NzJs7WpgERMRLgEeACzLzOZ2uRxoUEbsBu2XmzyJiO+Ba4JX+u1OdFhEBbJuZj0TEVOBHwFsz\n88oOlyYBEBFvBw4Ats/MozpdjzQoIu4EDsjMe0dqO9FHAA8ElmXm7Zm5FrgIOLrDNUkAZOYVVCvZ\nSl0lM+/JzJ/V7x8GbgHmdLYqCbLySL07td4m7l+q1VMiYi5wJHBOp2uRtsRED4BzgLsa9pfjf8RI\n0qhFxO7AfsBVna1EqtRT7K4HVgJlZto31S0+CbwL2NDpQqQWErgsIq6NiIXDNZzoAVCStJkiYhbw\ndeBtmflQp+uRADLzicx8HjAXODAinEKvjouIo4CVmXltp2uRhnBIZu4PHA4sqm9FammiB8C7gXkN\n+3PrY5KkYdT3V30d+NfM/Ean65GaZeYDwOXAgk7XIgEvAl5R32d1EfBnEfHlzpYkbZSZd9evK4Fv\nUt0q19JED4BXA3tFxB4RMQ04Blja4ZokqavVC22cC9ySmR/vdD3SoIjYJSJ2qN/PoFrk7ZedrUqC\nzPzfmTk3M3en+u/NH2TmaztclgRARGxbL+pGRGwL/Dkw5Cr0EzoAZuZ64GTgUqpFDC7OzJs6W5VU\niYgLgZ8+lpvJAAAEbUlEQVQAe0fE8og4odM1SbUXAa+j+gv29fV2RKeLkoDdgMsj4kaqP/KWmely\n+5I0vD8BfhQRNwA/Bb6Tmd8bqvGEfgyEJEmSJGn0JvQIoCRJkiRp9AyAkiRJktQnDICSJEmS1CcM\ngJIkSZLUJwyAkiRJktQnDICSpL5XFMX5RVGc2ek6hlIUxZ1FURzW6TokSRPflE4XIEnSSIqiuBMY\nAAbKsry34fh1wPOAPcqyvLNNtewO3AF8tyzLIxqOfxlYVpbl4nbUIUnS5nAEUJI0UdwBHDu4UxTF\nnwIzO1cOBxVF8d87+P1jVhSFf/iVpD7n/xFIkiaKLwF/A3y63n89cAHwx6mbRVEcWe8/HXgQOLdx\nRK4oikOAjwD7Ag8D7yvL8vz69I5FUXwHeAlwM3BcWZa/HqaejwAfBF7WfKIoiuOBE8uyPKThWAJ7\nlWW5rCiK84HVwB7Ai4EbgL8CTql/rt8Dx5ZleV3DZV9QFMWngN2AbwFvKstyTX3to+qfe/e69r8r\ny/LG+tydwOeA1wB7F0WxbVmW64f5uSRJPcwRQEnSRHElsH1RFM8qimIycAzw5aY2j1KFxB2AI4E3\nFUXxSoCiKJ4GfJcqQO5CNXX0+obPHgOcDuwILKMKd8P5LPDMLbg379XAe4HZwOPAT4Cf1fv/Bny8\nqf1rgL+gCrfPrD9LURT7AecBJwE7A2cDS4ui2Kbhs8dS/T52MPxJUn9zBFCSNJEMjgL+P+AW4O7G\nk2VZ/rBh98aiKC4EXko1YnYc8P2yLC+sz99Xb4O+WZblTwGKovhXNg1gzR6jColnAt/fjJ/lm2VZ\nXlt/3zeBN5dleUG9/1Xg5Kb2nynL8q76/Aepgux7gYXA2WVZXlW3+2JRFO8BDqb6PQF8avCzkqT+\nZgCUJE0kXwKuoJo6eUHzyaIoDgI+BDwHmAZsA3ytPj0PGG5K5+8a3q8GZo2innOAdxZF8ZejaNvs\n9w3vH2ux3/z9jQHuN1SL4gA8DXh9URRvaTg/reF882clSX3MAChJmjDKsvxNURR3AEcAJ7Ro8hXg\nM8DhZVmuKYrik1RTKqEKQQeOcz1ri6I4HfgAcFPDqUdpWKCmKIpdx+Hr5jW8fyqwon5/F/DBsiyH\nm7Ka4/D9kqQeYACUJE00JwA7lmX5aItVLbcDVtXh70CqaZ+X1ef+FXhPURSvBr4BPAWYV5bl9WyZ\nL1Et3rIA+FV97Abg2UVRPA/4JbB4C78DYFFRFN+mGp08FfhqffwLwDeLovg+8FOq4HkocEVZlg+P\nw/dKknqIi8BIkiaUsix/XZblNUOcfjNwRlEUDwOnARc3fO63VCOH7wBWUS0AM38c6nmi/q6dGo7d\nBpxBdW/gr4Afben3UI1uXgbcTjWV9cz6u64B3kg18nk/1QI2x4/D90mSelBkOitEkiRJkvqBI4CS\nJEmS1CcMgJIkSZLUJwyAkiRJktQnDICSJEmS1CcMgJIkSZLUJwyAkiRJktQnDICSJEmS1CcMgJIk\nSZLUJ/4/Rbn2HIrMChsAAAAASUVORK5CYII=\n",
	"text/plain": [
	"<matplotlib.figure.Figure at 0x7f5b148c7c88>"
	]
	},
	"metadata": {},
	"output_type": "display_data"
	}
	],
	"source": [
	"# Chart\n",
	"fig, ax1 = plt.subplots(figsize=(15,4))\n",
	"plt.title(r\"“True” Mach-Drag Table\")\n",
	"plt.ylabel(r\"$C_D$\")\n",
	"plt.xlabel(r\"Mach Number\")\n",
	"\n",
	"plt.plot(M[:33], a, 'r-', lw=1.4, alpha=0.8, label=\"True Value\")\n",
	"\n",
	"plt.ylim([0, 1.5])\n",
	"plt.xlim([0, 5.0])\n",
	"#plt.legend(loc=1)\n",
	"plt.show()"
	]
	},
	{
	"cell_type": "raw",
	"metadata": {},
	"source": []
	}
	],
	"metadata": {
	"celltoolbar": "Edit 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.5.3+"
	}
	},
	"nbformat": 4,
	"nbformat_minor": 2
	}