Skip to content
Permalink
master
Switch branches/tags

Name already in use

A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?

Comp-Mech-Project-2/02_Analyze-data_project.ipynb

Go to file
 
 
Cannot retrieve contributors at this time
executable file 338 lines (338 sloc) 61.8 KB
Raw Blame
Open in GitHub Desktop
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Computational Mechanics Project #02 - Create specifications for a projectile robot\n",
"\n",
"On the first day of class, we threw $2\"\\times~2\"$ dampened paper (spitballs) at a target on the whiteboard. Now, we are going to analyze the accuracy of the class with some cool Python tools and design a robot that has the same accuracy and precision as the class, but we will have the robot move farther away from the target and use a simpler projectile i.e. a tennis ball so we don't need to worry about knuckle-ball physics. \n",
"\n",
"The goal of this project is to determine the precision of necessary components for a robot that can reproduce the class throwing distibution. We have generated pseudo random numbers using `numpy.random`, but the class target practice is an example of truly random distributions. If we repeated the exercise, there is a vanishingly small probability that we would hit the same points on the target, and there are no deterministic models that could take into account all of the factors that affected each hit on the board. \n",
"\n",
"<img src=\"../images/robot_design.png\" style=\"height: 250px;\"/>\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now, we ask ourselves some questions:\n",
"\n",
"1. How do we quantify the class accuracy and precision?\n",
"\n",
"2. If we design a robot, what design components can we control?\n",
"\n",
"3. How can we relate the controlled components to the class accuracy, and specify the component precision?\n",
"\n",
"The first question, we have some experience from our work in [02_Seeing_Stats](../notebooks/02_Seeing_Stats.ipynb). We can define the mean, standard deviation, measure the first, second, and third quartiles, etc. \n",
"\n",
"The second question is a physical question. We cannot control the placement of the robot or the target those are chosen for us. We cannot control temperature, mechanical vibrations, etc. We *can* control the desired initial velocity. The initial velocity will have some speed and direction, and both will be subject to random noise. Once the speed and direction are set, the location on the target is determined by kinematic equations for an object in freefall, as such\n",
"\n",
"$x_{impact} = \\frac{v_x}{v_y}d + x(0)~~~~~~~~~~~~~~~~~~~~(1.a)$\n",
"\n",
"$z_{impact} = d\\left(\\frac{v_z(0)}{v_y}-\\frac{g}{2v_y^2}d\\right)+ z(0)~~~~~(1.b)$.\n",
"\n",
"Where the location of impact is at a $y$-distance of $d$ at a point on the target with coordinates $(x_{impact},~z_{impact})$, and the initial velocity is $\\bar{v}=v_x\\hat{i}+v_y\\hat{j}+v_z(0)\\hat{k}$, the object is released at an initial location $\\bar{r}(0)=x(0)\\hat{i}+0\\hat{j}+z(0)\\hat{k}$, and the only acceleration is due to gravity, $\\bar{a}=-g\\hat{k}$. Equation (1) becomes much easier to evaluate if we assume that $v_x=0$, resulting in an evalution of the accuracy of the height of the impact, $z_{impact}$, as such\n",
"\n",
"$x_{impact} = x(0)~~~~~~~~~~~~~~~~~~~~(2.a)$\n",
"\n",
"$z_{impact} = \\frac{d}{\\cos{\\theta}}\\left(\\sin{\\theta}-\\frac{g}{2v_0^2\\cos{\\theta}}d\\right)+ z(0)~~~~~(2.b)$.\n",
"\n",
"Where $\\theta$ is the angle of the initial velocity and $v_0$ is the initial speed. Equation (2) restricts the analysis to height accuracy. You can incorporate the 2D impact analysis if you finish the 1D analysis. "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The third question, is how we can relate equation (2) to the measured points of impact? For this, we can use Monte Carlo methods *(There are other methods, but Monte Carlo is one of the most straight-forward)*. Our Monte Carlo approach is as such, if we have a desired initial speed, $v_0$, and desired angle, $\\theta$, we can propagate the uncertainty of our actual speeds and angles into the $z_{impact}$ locations. Then, we can choose distributions in speed and angles that match the distributions in $z_{impact}$ locations. Here are the steps:\n",
"\n",
"1. Generate random $\\theta_i$ and $v_{0~i}$ variables\n",
"\n",
"2. Plug into eqn 2 for random $z_{impact~i}$ locations\n",
"\n",
"3. Compare to our measured $z_{impact}$ location statistics\n",
"\n",
"4. Repeat 1-3 until the predicted uncertainty matches the desired uncertainty, we can use a number of comparison metrics:\n",
" \n",
" - standard deviation\n",
" \n",
" - first, second, and third quartiles\n",
" \n",
" - visually, with box plots and histograms"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Project Deliverables\n",
"\n",
"1. Statistical analysis of class accuracy and precision (x- and z-locations) data is in the csv file [../data/target_data.csv](../data/target_data.csv) _Note: if you want to see how I turned the images into data check out the jupyter notebook [process_target_practice](./process_target_practice.ipynb)\n",
"\n",
"2. A Monte Carlo model to generate impact heights based upon uncertainty in $\\theta_0$ and $v_0$. \n",
"\n",
"3. The precision required to recreate the class accuracy and precision with a robot. \n",
"**You must show some validation of your work**\n",
"\n",
"4. [BONUS] Repeat 2-3 taking into account the variation in $x_{impact}$ due to misalignment. \n",
"\n",
"Given constants and constraints:\n",
"\n",
"- $d=$3 m, distance to target\n",
"\n",
"- $g=$9.81 m/s$^2$, acceleration due to gravity\n",
"\n",
"- $z(0)=$0.3 m, the initial height is 0.3 m above the bull's eye\n",
"\n",
"- 4 m/s$<v_0<$12 m/s, the initial velocity is always higher than 9 mph and less than 27 mph"
]
},
{
"cell_type": "code",
"execution_count": 97,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import pandas as pd\n",
"\n",
"#Import rcParams to set font styles\n",
"from matplotlib import rcParams\n",
"\n",
"#Set font style and size \n",
"rcParams['font.family'] = 'sans'\n",
"rcParams['font.size'] = 16\n",
"rcParams['lines.linewidth'] = 3"
]
},
{
"cell_type": "code",
"execution_count": 216,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f6af029feb8>]"
]
},
"execution_count": 216,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAD9CAYAAAB0i+q4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3dfbAddZ3n8feH5+Cum7CGUTKEhAeRMHFAr1vUoC5GMIwiIMwoKsrOKLF0VdCFkgyOssAWYbEAHedBLFdnCCUI8jgwRCQ4a4G4c0NACPIQTICJMgZDeAzhYb77R/eFk07fc/qc231On9OfV9Wtm3T3r7t/3ef29/weWxGBmZlZGbYZ9AmYmdnocFAxM7PSOKiYmVlpHFTMzKw0DipmZlaa7QZ9AoP0ute9LubMmTPo0zAzGyorVqx4PCJm5q1rdFCZM2cO4+Pjgz4NM7OhIunhyda5+svMzErjoGJmZqVxUDEzs9I4qJiZWWkcVMzMrDQOKmZmVhoHFTMzK42DipmZlabRgx/Nrl65jvOW3c+vN25it+nTOHXhvhx94KxBn5bZ0HJQsca6euU6Fl95N5tefBmAdRs3sfjKuwEcWMx65Oova6zzlt3/SkCZsOnFlzlv2f0DOiOz4eegYo31642bulpuZp05qFhj7TZ9WlfLzawzBxVrrFMX7su07bfdYtm07bfl1IX7DuiMzIafG+qtsSYa4937y6w8DirWaEcfOMtBxKxElVZ/Sdpd0hWSnpT0lKQrJc0umDYm+Tkgs902khZLWivpeUl3STq2mhyZmVk7lZVUJO0MLAc2AycAAZwN3CLpzRHxbIHdfA/4VmbZA5n/nwWcApwOrACOAy6XdERE3NB7DszMrFtVVn+dCOwJ7BsRqwEk/QJ4EPgUcH6BfayLiNsnWylpV5KAsiQivpYuvkXS3sASwEHFzKyPqqz+OhK4fSKgAETEGuBW4KiSjrEQ2AFYmlm+FJgvaW5JxzEzswKqDCr7A/fkLF8FzCu4j09L2izpOUnLJb0j5xibgdWZ5avS30WPY2ZmJagyqOwCPJGzfAMwo0D6pcBngEOBRcB/BpZLOiRzjI0RETnHmFi/BUmLJI1LGl+/fn2B0zAzs6Kq7lKcfdgDqFDCiI+1/Penkq4hKfmcDby9ZV9dHSMiLgIuAhgbG8tLa2ZmPaqypPIEOSUFklJKXgmmrYh4GrgeeFvL4g3ADEnZIDKjZb2ZmfVJlUFlFUmbR9Y84N4e95ktmawCdgT2yjkGUziOmZn1oMqgci1wkKQ9JxZImgMcnK7riqTXAu8Dft6y+EbgBeCjmc2PB+5Je5uZmVmfVNmm8m3gs8A1kr5MUsI4C3iUlgGNkvYAHgLOjIgz02WnAPsCtwC/BvYgGY/yeloCSET8VtIFwGJJTwN3AB8CFlBet2UzMyuosqASEc9KWgBcAFxMUnV1M3ByRDzTsqmAbdmy1HQ/8IH05z8BT5GMb/lERPy/zKFOB54BTiIJOvcDH4yI60rPlJmZtaWte+M2x9jYWIyPjw/6NMzMhoqkFRExlrfO71MxM7PSOKiYmVlpHFTMzKw0DipmZlYaBxUzMyuNXydsI+/qlev8HnqzPnFQsZF+6F69ch2Lr7ybTS++DMC6jZtYfOXdACOTR7M6cfVXw008dNdt3ETw6kP36pXrBn1qpThv2f2vBJQJm158mfOW3T+gMzIbbQ4qDTfqD91fb9zU1XIzmxoHlYYb9YfubtOndbXczKbGQaXhhuGhe/XKdRy8ZDlzT7ueg5cs76pq7tSF+zJt+223WDZt+205deG+ZZ+mmeGg0nh1f+hOtc3n6ANncc4x85k1fRoCZk2fxjnHzHcjvVlF3Pur4SYernXt/dWuzafoOR594Ky+5WeUe9KZFeGgYn196HZrmNp83H3ZzNVfVnPD0OYzYdR70pkV4aBitVb3Np9Ww1SqMquKq7+s1ure5tNqt+nTWJcTQOpYqiqD248sj4NKj/wH1T91bvNpderCfbdoU4H6lqqmyu1HNhlXf/Vg1Kc2sd40qfuy249sMi6p9KCMbq42moalVDVVbj+yybik0gP/QVnTDVOvPOsvB5Ue+A/Kmm6YeuVZfzmo9MB/UNZ0TWo/su64TaUHw9TN1awqTWk/su44qPTIf1DDyV3BzarloGKN4bEVZtVzm4o1hsdWmFXPQcUaw13BzarnoGKN4a7gZtWrNKhI2l3SFZKelPSUpCslzS6QbkzSRZLuk/ScpEckXSJpbs62ayVFzs/R1eTKhpW7gptVr7KGekk7A8uBzcAJQABnA7dIenNEPNsm+XHA/sA3gFXALOAvgXFJB0TEo5ntlwFnZJa5oty24K7gZtWrsvfXicCewL4RsRpA0i+AB4FPAee3SXtuRKxvXSDpVmBNut+vZLZ/PCJuL+vEbXS5K7hZtaqs/joSuH0ioABExBrgVuCodgmzASVd9jCwnqTUYmZmNVRlUNkfuCdn+SpgXrc7k7QfsCvwy5zV70/bXjZLut3tKWZmg1FlUNkFeCJn+QZgRjc7krQd8HckJZXvZFZfB3wOWAh8FHgeuErS8ZPsa5GkcUnj69dvVSAyM7MpqHpEfeQsUw/7+SbwR8D7ImKLQBURn9ti59JVwO3AOcDSrU4o4iLgIoCxsbG88zOzHJ7ixoqosqTyBElpJWsG+SWYXJLOARYBfx4RP+q0fUS8DFwO/L6kNxQ9jplNzm87taKqDCqrSNpVsuYB9xbZgaTTgdOAkyLi4i6OPVEacknErASe4saKqjKoXAscJGnPiQWS5gAHp+vakvR5knEtp0fEXxU9aNr+8qfAIxHxWJfnbGY5PMWNFVVlUPk2sBa4RtJRko4ErgEeBb41sZGkPSS9JOkrLcuOAy4EbgSWSzqo5Wdey3YflnSppI9Lelea7hbgrcCXKsybWaN4ihsrqrKgko6YXwA8AFwMXEIyeHFBRDzTsqmAbTPncni6/HDgZ5mfv2nZbg1JN+PzgB+RBKvNwOERcWn5uTJrJk9xY0UpornNDmNjYzE+Pj7o07Ah1MSeUE3Ms+WTtCIixvLW+SVdZl1q6su+PMWNFeGp78265J5QZpNzScWsS+4JVU+unqsHl1TMuuSeUPXjwZn14ZKKWQfZb8DvetNMfrhi3RZVYO4JNVjtqiRbSysuzVTPJRWzNvK+Af9wxTqOfessZk2fhoBZ06dxzjHz/XAaoCJVki7N9IdLKmZtTPYN+Jb71nPraQsGdFaWtdv0aazLCSytVZJFSzM2NS6p2FC4euU6Dl6ynLmnXc/BS5b37dulG+WHQ5HBmb6X/eGgYrU3yGoLN8oPh6MPnMU5x8xvWyXpe9kfrv6y2htktcWpC/fdYqAjuFG+rjoNzvS97A8HFau9QVZbTDyk3GNo+Ple9oeDitVekUbYKnl6ktHhe1k9t6lY7XmGXLPh4ZKK1Z6rLcyGh4OKDQVXWwwXj1xvLgcVMytVt68GcAAaLQ4qZiNm0A/pbrqAN/XdNKPMDfVmI6QO81t10wXc76YZPQ4qZiPigpseqMVDupuR6546ZfQ4qJiNiK/f/GAtHtLddAH31CmdDWreu145qJiNkDo8pIvMwzXBY5Daq0N1ZrfcUG82Quoyv1XRLuAeg9TeME7X76BiNkKG8SHtMUiTq0N1ZrccVMxGjB/So2PQ8971wm0qZmY1NYxtTi6pmJnV1DBWZzqomJnV2LBVZ7r6y8zMSlNpUJG0u6QrJD0p6SlJV0qaXTDtTpLOk/QbSZsk/UzSO3O220bSYklrJT0v6S5Jx5afGzMz66SyoCJpZ2A58CbgBOBjwD7ALZJeU2AX3wFOBL4CHAH8Blgm6YDMdmcBZwDfBP4YuB24XNJ7S8iGmZl1oco2lROBPYF9I2I1gKRfAA8CnwLOnyyhpD8EPgL8eUR8N132z8Aq4EzgyHTZrsApwJKI+Fqa/BZJewNLgBsqyJeZmU2iyuqvI4HbJwIKQESsAW4FjiqQ9kXgspa0LwGXAgsl7ZguXgjsACzNpF8KzJc0d0o5MBsBwzZ3lA23KoPK/sA9OctXAfMKpF0TEc/lpN0B2Ltlu83A6pztKHAcs5E2jHNH2XCrsvprF+CJnOUbgBlTSDuxfuL3xoiIDtu9QtIiYBHA7NmF+gyY9c0FNz3A129+sOf0c067vuM2m158mZMvu5OTL7vzlWUnvXsfvnDYG3s+bt0M+kVlTVb1OJXswx5ABdKpYNqi2716QhEXARcBjI2N5aU1G5gvHPbGnh/uc067nrVL3rfFsrmnXT/pH8iazLZlGuRD3W+THKwqq7+eIKekQFJKySuFtNrQJu3E+onfMyRlg0h2O7Na6Vc7xyCmwh90lVsdXlTWZFUGlVUkbR5Z84B7C6Sdm3ZLzqZ9gVfbUFYBOwJ75WxHgeOY9V0/H7plzh1VNBAO+qE+jDP7jpIqg8q1wEGS9pxYIGkOcHC6rlPa7YE/bUm7HfAh4EcRsTldfCNJkPloJv3xwD1pbzOzWunnQ7ebF2a1000gHPRDvQ4vKmuyKttUvg18FrhG0pdJ2j7OAh4FvjWxkaQ9gIeAMyPiTICIuFPSZcCFkrYH1gCfBubSEkAi4reSLgAWS3oauIMk8Cygc7dls1xltAe020e/H7plzB3VzcuiBj1de11eVNZUlZVUIuJZkof7A8DFwCUkwWFBRDzTsqmAbXPO5c+A7wJnA9cDuwOHR8Qdme1OT7c5CVhGUhL6YERcV2qGrBHKqJrqtI9h/CbdTSAc9HTtZZXOrDeV9v6KiEeAtvNwRcRacnprRcQm4IvpT7v0L5MElbN7PlGzVBmvb+20jzp+k+5UOuum9FGH6dqHbWbfUeKp781alFE11WkfdXjotirSBbfbQOiHenM5qJi1KKM9oMg+6vTQLVI6KzMQemDiaHNQsUbp9EAro2qqjtVb7RQtnWUD4UQX426CQ9kDEx2g6scv6bLGKNIIX0Yj77A1FPfScaDXDg1ldqce9CBLy+eSijVG0Ub4Mqqm6lS91UkvJateOzSU2Z26jE4VVj6XVKwxBj0or656KVn1ei3L7E7t+1lPLqlYYwx6UF6ddVuy6vValtne1O4c3NYyOC6pWGMMelBenjInljzp3fuUeGbt9Xoty2xvmuwc3vWmmW5rGSBt/SqS5hgbG4vx8fFBn4b1UZ2+wWZ7QkHyUKxzo36rOlzLvHM4b9n9uSWYWdOncetpC/p6fqNK0oqIGMtd56DioGKDcfCS5X74VWBQ75BpknZBxW0qNpTq8C15qtzQXA23nQ2W21Rs6IzK+IQ6TyzZr5eIVaGObWdN4qBiQ2fQL4EqS10ffsMetIdt8OmocfWXDZ1RqTaq28SSE0ZhUOEwDT4dNQ4qNnRGqc68zIdfWe1MoxK0bTBc/WVDp67VRoNUZpVVndt6rP4cVGzoVFlnPqwN1GW2MxUN2sN6raxarv6yoVRFnXnZ07L3U5lVVkXaeob5Wlm1HFTMUsPcQF12O1OnoF3kWo3CWCLrnqu/zFJVN1BXWV3U73amTtdq2LslW+8cVMxSVTZQV/2Q7ffYjE7XalTGEln3XP1llqryNcC9VK11W33Uz7EZna6VuyU3l4OKWarKwYjdPmTr3hDe6VqN0lgi646DilmLqr7td/uQHYZOA+2uVZWlPqs3t6mY9UG3DenDXn3k+beayyUVsz7otmptFKqPPP9WMzmomPVJNw9ZVx/ZsHJQMStZGYP+6jqDsVknDipmJSqz15arj2wYVdZQL2kbSYslrZX0vKS7JB1bIN1rJX1F0m2SfidpY/rvo3O2PUNS5PxcXU2uzNrzoD9ruipLKmcBpwCnAyuA44DLJR0RETe0STcb+Azw3XQf/w58GLhK0mcj4q9z0rwdaP1L3lDC+Zt1bdh7bdWB5wwbbpUEFUm7kgSUJRHxtXTxLZL2BpYA7YLKGmDPiHiuZdkySbsDXwLygsrPI+KlEk7dbEpGodfWINV90Kd1VlX110JgB2BpZvlSYL6kuZMljIhnMwFlwjiwW3mnaFa+TuNR/A6S9lx9OPyqCir7A5uB1Znlq9Lf83rY5zuB+yZZ96iklyU9LOlcSf5aaAPRbtCfZ+7tzNWHw6+qNpVdgI0REZnlG1rWFyZpEXAQcHxm1WrgNGAlEMB7gC8AbwEOa7OvRQCzZ8/u5jTMCpms19YwTL0yaK4+HH6FSiqSDp2kl1X25ycTSUge8lvtqtsTlHQI8A3g4oi4pHVdRCyNiHMj4kcRcVNEnAqcChwq6dC8/UXERRExFhFjM2fO7PZ0zHrmb+Gd9fu9MFa+oiWV24D9Cmw30RayAZghSZnSyoyW9R1JehtwLbAc+ETBc/0+cCHwNuDHBdOYVc7fwjvzoM/hVyiopA3nk7Vn5FkF7AjsxZbtKhNtKfd22oGk+cAy4E7g2Ih4sYvjQ35JyWxgPPVKMR70Odyqaqi/EXgB+Ghm+fHAPRGxpl1iSfsANwG/Ao6IiG7qByaO+fMu0phVzjP3WhNU0lAfEb+VdAGwWNLTwB3Ah4AFwFGt20q6GdgjIvZO/78rSUDZAfgqME/aoilmZURsTrddCfwDcD9JyeQw4HPAjRFxSxV5M5sKfwu3UVfliPrTgWeAk4DXkzz4PxgR12W22zZzHvOAPdJ//2POfucCa9N/3w98FnhDup+HgDOB/z310282j2ruja+bNZ227vXbHGNjYzE+Pj7o06id7KhmSOr+XVXTnq/bqxxcR5ukFRExlrfOsxTbVjyeojd1vm79fMh7qpV66tdnwK8Ttq14PEVv6nrd+j2S31Ot1E8/PwMOKraVycZNeDxFe3W9bv1+yNc1uDZZPz8DDiq2FY9q7k1dr1u/H/JVBFdPxDk1/fwMOKjYVjyeojd1vW79LkGVHVw9EefU9fMz4N5f7v1lI24QvdLKbBQ+eMny3OltZk2fxq2nLZjqqTZC2Z8B9/4ya7BBzKdV5iBPt9FMXT8/Aw4qZg0wzCP5PRFnOfr1GXCbipnVWl07QFg+l1TMrNaGcTr8Js8o4KBiZrU3TNV3TZ9RwNVfZmYlavqMAg4qZmYlanpvNQcVM7MS1XW6nn5xUDEzK1HTe6u5od7MrETD2FutTA4qZmYlG6beamVz9ZeZmZXGQcXMzErjoGJmZqVxUDEzs9I4qJiZWWnc+8sar8mT/5mVzUHFGq3fk/85gNmoc/WXNVo/J//zu9atCRxUrNH6Oflf02evtWZwULFG6+fkf02fvdaawUHFGq2fk/81ffZaa4bKgoqkbSQtlrRW0vOS7pJ0bMG035MUOT8X5mz7dkm3Sdok6TFJ50vyX6kVcvSBszjnmPnMmj4NAbOmT+OcY+ZX0nje9NlrrRmq7P11FnAKcDqwAjgOuFzSERFxQ4H064EjM8t+0/ofSW8GbgKWAUcAc4HzgFnAh6Z09tYY/Zr8r+mz11ozKCLK36m0K/AosCQivtqy/GZgZkS8uUP67wGHRsTvd9juKuAPgHkR8WK67OPA3wNvjYg72qUfGxuL8fHxAjkyM7MJklZExFjeuqqqvxYCOwBLM8uXAvMlzZ3qASRtDxwO/GAioKR+ALwAHDXVY5iZWXeqqv7aH9gMrM4sX5X+nges6bCPXSU9DkwHfgV8B/haREz0ydwL2Am4pzVRRDwv6aH0GGbWRx7caVUFlV2AjbF13dqGlvXt3EnSDrOKJHB8ADgH2Af4ZGYfT+Sk3zDZMSQtAhYBzJ49u8NpmFlR/Z6dwOqpUPWXpEMn6Y2V/fnJRBIgr7FGRY4XERdGxF9FxPKIuCEiTgS+DnxC0j6ZfXV1nIi4KCLGImJs5syZRU7HzArw4E6D4iWV24D9Cmz3XPp7AzBDkjKllRkt67v1feBkYAx4kPalnhm8WtVmZn3gwZ0GBYNKRDwH3NfFflcBO5K0e7S2q0y0c9zbxb4mZEsmD5G02+y/xUbSTsCewOU9HMPMerTb9GmsywkgHtzZLFX1/rqRpAfWRzPLjwfuiYhOjfR5PkISUP4FICJeSI/zQUmtwfFPSALatT0cw8x65MGdBhU11EfEbyVdACyW9DRwB8lgxAVkuvqmY1f2iIi90//vAVwMXEpSytmRpKH+vwHfioiHWpKfAfwM+IGkvwbmkAx+vCIiVlSRNzPL58GdBtWOqD8deAY4CXg9cD/wwYi4LrPdtpnzeJqkveRLwO+RlE5+CXwe+JvWhBFxp6SFwLnA9cCTwD8Af1F2Zsyss37NTmD1VcmI+mHhEfVmZt0bxIh6MzNrIAcVMzMrjYOKmZmVxkHFzMxK46BiZmalcVAxM7PSVDlOxcxsqHjq/qlzUDEzw1P3l8XVX2ZmeOr+sjiomJnhqfvL4qBiZsbkU/R76v7uOKiYmeGp+8vihnozMzx1f1kcVMzMUp66f+pc/WVmZqVxUDEzs9I4qJiZWWkcVMzMrDQOKmZmVhoHFTMzK42DipmZlUYRMehzGBhJ64GHO2z2OuDxPpzOIDmPo8F5HA3DkMc9ImJm3opGB5UiJI1HxNigz6NKzuNocB5Hw7Dn0dVfZmZWGgcVMzMrjYNKZxcN+gT6wHkcDc7jaBjqPLpNxczMSuOSipmZlcZBxczMStPYoCJpG0mLJa2V9LykuyQdWyDdHEnR5ue4lm3PmGSbq6vN3SvH7ymPadrvTXLuF+Zs+3ZJt0naJOkxSedL6ss7WKdwH18r6Svpef9O0sb030fnbNuX+yhpd0lXSHpS0lOSrpQ0u2DanSSdJ+k36X34maR35mzX82eiDL3mUdKYpIsk3SfpOUmPSLpE0tycbddOcr+2urdVmOJ9nOy5ckBmu4Hex7YiopE/wP8CNgOnAO8CvgX8O/DeDul2BA7K+fkx8DywS8u2ZwABHJzZ9o11zmOa9nvAb3PyuUdmuzcDm4CrgXcDnwSeAC6rcx6BPwAeA84B/hhYmOY5gP+e2bby+wjsDDwI3AMcDRwF3A08BLymQPpLgI3Aiel9uDK9LweU9ZkYZB6BrwG3Ap8B/ivwEeCXwO+A3TPbrgVuzPnszqhzHtP0AXw359x3rst97JiHQZ/AQDINu6Y35H9mlt8M/KLHD9JTwOWZ5RMPo+2GLY/pA/ZfC2x3VfpHtH3Lso+n+X5LXfMIvCb7h9qS9pF+30fgJOBlYO+WZXOBl4Avdkj7h+n5/VnLsu2A+4Fry/pMDDiPM3OW7ZE+SM/MLF8LLK06P2XnMd02gLM7bDPQ+9jpp6nVXwuBHYClmeVLgfl5ReoOjgH+I/D3JZxbWcrO41YkbQ8cDvwgIl5sWfUD4AWSb2lV6jmPEfFsRDyXs2oc2K28UyzsSOD2iFg9sSAi1pB8O+90HY8EXgQua0n7EnApsFDSjuniyj8TBc6zpzxGxPqcZQ8D64E6vf93KvexqEHfx7aaGlT2J4n0qzPLV6W/53W5vxNIqopunGT9o5JelvSwpHP71N5QRh53lfS4pJckPSDpS5K2bVm/F7ATSVH/FRHxPElxv9vr2K2y7yPAO4H7JllX5X3cn8x1TK2icz72B9bkBMlVJA+fvVu2K/t6dWMqedyKpP1IvrX/Mmf1+9O2l82Sbu9Xewrl5PHT6Xk/J2m5pHfkHGOQ97Gt7QZ58AHaBdgYaZmxxYaW9YVImgUsAL6efjtstRo4DVhJUqx9D/AF4C3AYT2cdzemmsc7gRUkH9SdgA+QtD/sQ9Ju0rqPJ3LSbyhwjKkq7T4CSFpEUn99fGZVP+7jLkx+HWdMIe3E+onfpV2vHkwlj1uQtB3wdyQlle9kVl8H/AuwBvg94LPAVZI+FhHZb/dlm2oelwL/CPyapHrvVGC5pMMi4ictxxjkfWxrJIKKpEOBmwps+s8RcQggkofDVrvq4fAfIynxbVX1lfMBvknSvwIXSjo0In5c9CD9zmNEZHt53SDpGeBkSedGxIMt+yrlWg7yPko6BPgGcHFEXNK6rsz72EGveSl6Hcr83PeqrON/E/gj4H0RscVDPCI+t8XOpauA20m+FFUdVGBqf3cfa/nvTyVdQ1LyORt4e8u+Bn0fJzUSQQW4DdivwHYT1QMbgBmSlIn2M1rWF/Vx4M6IuKvg9t8HLgTeRtJjrKhB5nHC94GTgTGSxvl234xm8GpxvKiB5FHS24BrgeXAJwqea6/3cTJPMPl1zPvm22oDkNdlNXsdqvhMdGMqeXyFpHOARcAJEfGjTttHxMuSLgfOlfSGiPhN0WP1oJQ8ToiIpyVdz5afy0Hfx7ZGIqikdcmT1YPnWUXSNXgvtqyXnKiLvLfITtKH0X4kVSHd6mp+nEHlMSNbMnmIpG53/y02knYC9gQu72bng8ijpPnAMpLqvmMzHQ6KKGueo1VkrmNqHp3zsQr4gKSdM+0q80g6TKxu2a7sz0Q3ppJHACSdTlIV+fmIuLiLY7crVZdpynnMkS2ZDPo+tjfo7meD+OHVLnlfzSz/MXB3F/v5Jkmvm127SPMFkg/Iu4Yhj5m0XyfpwrlXy7KrgQdo6W5L0iYRwFvrnEeS9qHHSHp8vbbLY5d6H0lKgC8Be7Ysm5N+vv5Hh7QHpOdyQsuy7UgasK+r8jPRrzym234+zedfdHnc7dJ7/HDd85izv9cCj5BU+dbiPnY850GfwMAyDktIBit+ETgE+Nv0gfn+zHY3A6tz0m9P8na2a9scY2X68HkvyQC789MP1z/VOY8kDYT/l2Sg2XuA9wP/J037t5m0B5AMsruSZNDdJ0iK35dXnb8p5nFXkvEMG4D3sfVgsx37eR9Jxs2sJhkodxRJ19S7gF8B/yFzb14CvpJJfylJ9con0/twRXpd3pLZrtD1quhe9ZxH4Lj0PP8p517Na9nuw+m1+DjJoMDjgJ+SBKPjap7HU4BvkwzsPISkV+ndJKXNd9TlPna8BoM+gYFlHLYFvkzyOuHNwC+AP8nZ7ifA2pzlH0g/qMe2OcalJFVEz6XHuBf4y9YHVh3zSFInfHWa7nmSoHEHSS+abXLSvxP4Wbrtv5G0NWw1sLBmeTwkvX+T/czp930kaRf5IclA2qfTezAns82c9PzOyCyfRhLsHkvvw8+BQ3q9XhXer57yyKuzHeT9/KRlu4NI2sb+jSTwP0nyDX7hEOTx/STjWSi86OIAAABMSURBVB5Pz/13JG19/6Vu97Hdj6e+NzOz0jR18KOZmVXAQcXMzErjoGJmZqVxUDEzs9I4qJiZWWkcVMzMrDQOKmZmVhoHFTMzK83/B/TMTg0zCudGAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"student = pd.read_csv('../data/target_data.csv')\n",
"student.columns\n",
"s_xpos = student[' x position (m)']\n",
"s_ypos = student[' y position (m)']\n",
"s_xmean = np.mean(s_xpos)\n",
"s_ymean = np.mean(s_ypos)\n",
"s_xstd = np.std(s_xpos)\n",
"s_ystd = np.std(s_ypos)\n",
"\n",
"Q1_x = np.percentile(s_xpos, q=25)\n",
"Q2_x = np.percentile(s_xpos, q=50)\n",
"Q3_x = np.percentile(s_xpos, q=75)\n",
"\n",
"Q1_y = np.percentile(s_ypos, q=25)\n",
"Q2_y = np.percentile(s_ypos, q=50)\n",
"Q3_y = np.percentile(s_ypos, q=75)\n",
"\n",
"plt.scatter(s_xpos, s_ypos)\n",
"plt.plot(0,0,'+',markersize=30)"
]
},
{
"cell_type": "code",
"execution_count": 261,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Robot mean is -0.04556045690004387\n",
"Student mean is -0.047397370492807414\n",
"Robot standard deviation is 0.28715699506722125\n",
"Student standard deviation is 0.2548611138551949\n"
]
},
{
"data": {
"text/plain": [
"{'whiskers': [<matplotlib.lines.Line2D at 0x7f6aef6f7cc0>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef6f7da0>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef70e438>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef70e780>],\n",
" 'caps': [<matplotlib.lines.Line2D at 0x7f6aef704390>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef7046d8>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef70eac8>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef70ee10>],\n",
" 'boxes': [<matplotlib.lines.Line2D at 0x7f6aef6f78d0>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef704e10>],\n",
" 'medians': [<matplotlib.lines.Line2D at 0x7f6aef704a20>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef70eef0>],\n",
" 'fliers': [<matplotlib.lines.Line2D at 0x7f6aef704d68>,\n",
" <matplotlib.lines.Line2D at 0x7f6aef69a4e0>],\n",
" 'means': []}"
]
},
"execution_count": 261,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAD9CAYAAAB0i+q4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAASv0lEQVR4nO3df7DldV3H8ecLEMzK2k1JM/eXOOYymFPHhkkrIgxMghoUSSjTch0bf9AMM0ZbhmAjpiOMZuU6jhbr+APTgGRAgx3GUde8i2IsSkC7hKXj2i7ySxaxd3+cc/Nw9ty75977+d673Pt8zJw5536+n8/38/nOnLmv+Xw/3+/3pKqQJKmFw5Z6AJKk5cNQkSQ1Y6hIkpoxVCRJzRgqkqRmjljqASylJzzhCbVu3bqlHoYkPars2LHj21X1xHHbVnSorFu3jqmpqaUehiQ9qiS5c6Ztnv6SJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1IyhIklqZkXf/ChpeUky5zb+plRbhoqkZWOmgEhieCwST39JkpoxVCRJzRgqkqRmDBVJUjOdhkqSpyb5WJLvJLknyceTrJmwbc3wevZIvcOSnJ9kd5IHk9yU5IxujkiSNJvOrv5K8jjgemA/8DKggDcD25I8q6run2A3HwDeM1L27yN/XwScB2wGdgBnAZcnObWqrp7/EUiS5qrLS4pfCWwAnlFVtwMk+QpwG/Aq4B0T7OO/qmr7TBuTHE0/UC6uqrcPirclOQa4GDBUJGkRdXn66zRg+3SgAFTVLuCzwOmN+jgZOBLYOlK+FTguyfpG/UiSJtBlqBwL3DymfCewccJ9vDrJ/iQPJLk+yS+N6WM/cPtI+c7B+6T9SJIa6DJUVgP7xpTvBVZN0H4r8EfAScAm4CeA65OcMNLH3XXgrbJ7h7Y/QpJNSaaSTO3Zs2eCYUiSJtX1Y1rGPRdhoofzVNXvDv35mSRX0J/5vBl43tC+5tRHVW0BtgD0ej2f2yBJDXU5U9nHmJkC/VnKuBnMrKrqXuCTwHOGivcCq3LgU+RWDW2XJC2SLkNlJ/01j1EbgVvmuc/RmclO4CjgaWP6YAH9SJLmoctQuRI4PsmG6YIk64DnDrbNSZLHAy8EvjBUfA3wEHD2SPVzgJsHV5tJkhZJl2sq7wVeA1yR5M/ozzAuAu5i6IbGJGuBO4ALq+rCQdl5wDOAbcB/A2vp34/yJIYCpKq+leQS4Pwk9wI3Ai8BTqTdZcuSpAl1FipVdX+SE4FLgMvon7q6Dji3qu4bqhrgcB45a7oV+O3B68eAe+jf3/IHVfWvI11tBu4DXk8/dG4Fzqyqq5oflCRpVlnJP1zT6/VqampqqYchqWP+SFdbSXZUVW/cNp9SLElqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1IyhIklqxlCRJDVjqEiSmun6lx+1zBz4e2iT8blL0spgqGhOZgsHH9onydNfkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkprpNFSSPDXJx5J8J8k9ST6eZM0E7XpJtiT5WpIHkvxnkg8mWT+m7u4kNeb1W90clSRpJp39SFeSxwHXA/uBlwEFvBnYluRZVXX/LM3PAo4F3gnsBJ4C/DkwleTZVXXXSP1rgQtGym5d8EFIkuaky19+fCWwAXhGVd0OkOQrwG3Aq4B3zNL2rVW1Z7ggyWeBXYP9vnGk/reranurgUuS5qfL01+nAdunAwWgqnYBnwVOn63haKAMyu4E9tCftUiSDkFdhsqxwM1jyncCG+e6syTPBI4Gvjpm828O1l72J9nueookLY0uQ2U1sG9M+V5g1Vx2lOQI4O/oz1TeN7L5KuC1wMnA2cCDwCeSnDPDvjYlmUoytWfPARMiSdICdLmmAv3F+VGZx37+GvhF4IVV9YigqqrXPmLnySeA7cBbgK0HDKhqC7AFoNfrjRufJGmeupyp7KM/Wxm1ivEzmLGSvAXYBLyiqj51sPpV9X3gcuCnkzx50n4kSQvX5UxlJ/11lVEbgVsm2UGSzcCfAK+rqsvm0Pf0bMiZiCQtoi5nKlcCxyfZMF2QZB3w3MG2WSV5Hf37WjZX1bsm7XSw/vJi4D+r6ptzHLMkaQG6DJX3AruBK5KcnuQ04ArgLuA905WSrE3ycJI3DpWdBVwKXANcn+T4odfGoXq/k+TDSX4vya8O2m0Dfh54Q4fHtuytXr2aJHN6AXOqv3r1uLOjkh7NOjv9VVX3JzkRuAS4jP4pqeuAc6vqvqGqAQ7nkQF3yqD8lMFr2A3ACYPPu+hfZvw2+us3DwBfBE6pqmtbHs9Ks2/fPqq6PXs4HUSSlo90/Y/jUNbr9Wpqamqph3FISrIoobKSv39aPH7X2kqyo6p647b5lGJJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSLpUWeuT9GGuT1B26doz1/XPycsSc35FO1DlzMVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc10GipJnprkY0m+k+SeJB9PsmbCto9N8rYk30jy3SSfT/LLY+odluT8JLuTPJjkpiRntD8aSdLBdBYqSR4HXA/8DPAy4HeBpwPbkvzwBLt4H/BK4I3AqcA3gGuTPHuk3kXABcBfAy8AtgOXJ/mNBochSZqDLn/58ZXABuAZVXU7QJKvALcBrwLeMVPDJD8LvBR4RVW9f1B2A7ATuBA4bVB2NHAecHFVvX3QfFuSY4CLgas7OC5J0gy6PP11GrB9OlAAqmoX8Fng9Anafg/4yFDbh4EPAycnOWpQfDJwJLB1pP1W4Lgk6xd0BJKkOekyVI4Fbh5TvhPYOEHbXVX1wJi2RwLHDNXbD9w+ph4T9CNJaqjL01+rgX1jyvcCqxbQdnr79PvdVVUHqff/kmwCNgGsWTPRNQMrUv3F4+GCH+u+D0nLSpehAjD6zx4gE7TLhG0nrfeDAVVtAbYA9Hq9cW0F5E33cGBWN+4joS7otAtJi6zL01/7GDNToD9LGTcLGbZ3lrbT26ffVyUZDZHRepKkRdBlqOykv+YxaiNwywRt1w8uSx5t+xA/WEPZCRwFPG1MPSboR5LUUJehciVwfJIN0wVJ1gHPHWw7WNvHAC8eansE8BLgU1W1f1B8Df2QOXuk/TnAzYOrzSRJi6TLNZX3Aq8BrkjyZ/TXPi4C7gLeM10pyVrgDuDCqroQoKq+nOQjwKVJHgPsAl4NrGcoQKrqW0kuAc5Pci9wI/3gOZGDX7YsSWqss1CpqvuTnAhcAlxGf/H8OuDcqrpvqGqAwzlw1vRy4C+BNwM/DtwEnFJVN47U2wzcB7weeBJwK3BmVV3V9ogkSQeTrq/wOZT1er2amppa6mEckpIsztVfK/j7p/nz+7m0kuyoqt64bT6lWJLUjKEiSWrGUJEkNWOoSJKa6foxLZLUnM+mO3QZKpIedXw23aHL01+SpGYMFUlSM4aKJKkZQ0WS1IyhIklqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1IyhIklqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1IyhIklqxlCRJDXTWagkOSzJ+Ul2J3kwyU1Jzpig3eOTvDHJ55L8T5K7B59/a0zdC5LUmNc/dXNUkg4VSTp9rVq1aqkP8VHpiA73fRFwHrAZ2AGcBVye5NSqunqWdmuAPwLeP9jH/wK/A3wiyWuq6t1j2jwP+P7Q33sbjF/SIaqq5lQ/yZzbaH46CZUkR9MPlIur6u2D4m1JjgEuBmYLlV3Ahqp6YKjs2iRPBd4AjAuVL1TVww2GLklagK5Of50MHAlsHSnfChyXZP1MDavq/pFAmTYF/FS7IUqSWusqVI4F9gO3j5TvHLxvnMc+fxn42gzb7kry/SR3Jnlrkh+ax/4lSQvU1ZrKauDuOvAk5t6h7RNLsgk4HjhnZNPtwJ8AXwIK+HXgj4GfA54/y742AaxZs2Yuw5AkHcREoZLkJODTE1S9oapOAEL/n/wBu5p8aP/f9wnAO4HLquqDw9uqavT02qeTfB24NMlJVfUvo/urqi3AFoBer+fKnSQ1NOlM5XPAMyeoN70WshdYlSQjs5VVQ9sPKslzgCuB64E/mHCsHwIuBZ4DHBAqkqTuTBQqg4XzmdYzxtkJHAU8jUeuq0yvpdxysB0kOQ64FvgycEZVfW8O/cP4mZIkqUNdLdRfAzwEnD1Sfg5wc1Xtmq1xkqfTP932H8CpVfXdOfQ93ecX5tBGktRAJwv1VfWtJJcA5ye5F7gReAlwInD6cN0k1wFrq+qYwd9H0w+UI4G/ADYmj1iK+VJV7R/U/RLwD8Ct9GcmzwdeC1xTVdu6ODZJ0sy6vKN+M3Af8HrgSfT/8Z9ZVVeN1Dt8ZBwbgbWDz/88Zr/rgd2Dz7cCrwGePNjPHcCFwF8tfPgaCfPmfAyGtPxkJT+6oNfr1dTU1FIPY9nwURg6VPndbCvJjqrqjdvmU4olSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqZnOQiXJYUnOT7I7yYNJbkpyxoRtP5CkxrwuHVP3eUk+l+S7Sb6Z5B1Jfqj9EUmSDuaIDvd9EXAesBnYAZwFXJ7k1Kq6eoL2e4DTRsq+MfxHkmcBnwauBU4F1gNvA54CvGRBo5ckzVknoZLkaPqBcnFVvX1QvC3JMcDFwCSh8lBVbT9InTcBXwdeXFXfG/T9EPD3Sd5aVTfO7wgkSfPR1emvk4Ejga0j5VuB45KsX2gHSR4DnAJ8dDpQBj4KPAScvtA+JElz01WoHAvsB24fKd85eN84wT6OTvLtJA8n+fckb0hy+ND2pwGPBW4eblRVDwJ3TNiHJKmhrtZUVgN3V1WNlO8d2j6bL9Nfh9lJPzh+G3gL8HTgD0f2sW9M+70z9ZFkE7AJYM2aNQcZhiRpLiYKlSQn0V8QP5gbquoEIMBooDAoP6iqGr3K6+ok9wHnDtZKbhva15z6qaotwBaAXq83rq0kaZ4mnal8DnjmBPUeGLzvBVYlychsZdXQ9rn6EHAu0ANuY/ZZzyp+cKpNkrRIJgqVqnoA+Noc9rsTOIr+usfwusr0Osctc9jXtNGZyR30122OfUSl5LHABuDyefQhSVqArhbqr6F/BdbZI+XnADdX1a557POl9APliwBV9dCgnzOTDIfji+gH2pXz6EOStACdLNRX1beSXAKcn+Re4Eb6NyOeyMilvkmuA9ZW1TGDv9cClwEfpj/LOYr+Qv3vA++pqjuGml8AfB74aJJ3A+vo3/z4sara0cWxSZJm1uUd9ZuB+4DXA08CbgXOrKqrRuodPjKOe+mvl7wB+En6s5OvAq8D/ma4YVV9OcnJwFuBTwLfAf4B+NPWByNJOrgceNXvytHr9Wpqamqph7FsJGElf5906PK72VaSHVXVG7fNpxRLkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ10+Ud9ZK0qJKZf11jpm3eFNmWoSJp2TAglp6nvyRJzRgqkqRmDBVJUjOGiiSpGUNFktSMV39pTma7ZHO27V6VI60MhormxHCQNBtPf0mSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1IyhIklqxlCRJDWTlXwzW5I9wJ1LPY5l5AnAt5d6ENIYfjfbWltVTxy3YUWHitpKMlVVvaUehzTK7+bi8fSXJKkZQ0WS1Iyhopa2LPUApBn43VwkrqlIkppxpiJJasZQkSQ1Y6ho3pL8dJJ3Jfl8kgeSVJJ1Sz0uKcmLkvxjkjuTfDfJrUnekuRHl3psy52hooU4BjgT2Ad8ZonHIg07D/g+8KfAKcDfAq8GPp3E/3sdcqFe85bksKr638HnPwTeC6yvqt1LOjCteEmeWFV7Rsp+D/h74Neq6vqlGdnyZ2Jr3qYDRTrUjAbKwBcH709ZzLGsNIaKpJXiVwbvX13SUSxzhoqkZS/JU4ALgX+pqqmlHs9yZqhIWtaS/AhwBfAw8PIlHs6yd8RSD0CSupLkscCVwAbgV6rq60s8pGXPUJG0LCV5DPCPwC8AJ1XVvy3xkFYEQ0XSsjO4F+WDwK8BL6yq7Us8pBXDUNGCJHnR4OPPD95fMPhFzT1VdcMSDUt6N/Bi4C+B+5McP7Tt654G6443P2pBksz0Bbqhqk5YzLFI05LsBtbOsPlNVXXB4o1mZTFUJEnNeEmxJKkZQ0WS1IyhIklqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM/8HYTSh5/s5BdgAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"ang_l = np.pi/15\n",
"ang_h = np.pi/4\n",
"vel_h = 5.5\n",
"vel_l = 5.4\n",
"theta = ang_l + (ang_h - ang_l)*np.random.rand(10000)\n",
"v_0 = vel_l + (vel_h - vel_l)*np.random.rand(10000)\n",
"d=3\n",
"g=9.81\n",
"z_0 = 0.3\n",
"\n",
"z_r = d/np.cos(theta) * (np.sin(theta)-(g*d/(2*v_0**2*np.cos(theta)))) + z_0\n",
"z_rmean = np.mean(z_r)\n",
"z_rstd = np.std(z_r)\n",
"\n",
"print('Robot mean is',z_rmean)\n",
"print('Student mean is', s_ymean)\n",
"print('Robot standard deviation is',z_rstd)\n",
"print('Student standard deviation is', s_ystd)\n",
"\n",
"plt.boxplot([z_r, s_ypos])"
]
},
{
"cell_type": "code",
"execution_count": 262,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEoCAYAAACgg32xAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3de7xVdZ3/8ddbQRAvCUn1y0tgdhEinDrT9CsTtWagQrQxMy+VWlrWL7uM1qj91EzHUhvtopmXn0xhSWoZTiXewKaipoNpQmqReCmxOQpeEESNz++P9d253Kxz2N+z1zn7AO/n47Efm/1dn7X257s37M/6ru9aC0UEZmZmOTbrdAJmZrbhcfEwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYUOWpMMlhaTDO53LxkrSfElZ5+un72T+AKVkGwgXD+sXSePSj0j58bSk+yV9W9Junc5xfSTNTHmP63QuQ4mkU9Pnslenc7Gha1inE7AN3t3AFenP2wJvAg4D9pf0DxHxu45lZq14PzAqc53dgFUDkIttQFw8rF13RcSp5QZJ5wMfBf6V4sfJhqiIuL8f69w1ELnYhsWHrWwgzEzPr69aKGl/Sf8l6QlJT0r6taQP9rVBSQdKulXSakkPSvqypK0q4oZLOk7SHSl2haTrJL2lKe5e4APp5dLSobeZzdssrbOVpMcl3dbL8m0krZT02776kmL/dmhI0jGS7pT0lKSlkj4naZ0dO0lbSzpd0u8lrZHUI+lqSZMqYneQdL6kJelzeDh9fmc3xT1vziPNZZySXs4rfS7zSzGVcx6SXpvy6Un5/T7lu1VTXOOQ50xJu0q6RtJj6e/DtZJevr7PzzrPIw8bSM80N0j6F+AcoAf4D+Bp4F3AJZImR8SxFds5EHgrMBu4HvhH4NPA30l6W0SsTdsWcBUwA/gd8DVgNHAQxQ/hwRFxZdrmecDhwGTgK8Cjqb2yMABExJOSvgN8WNLrI2JhU8hBwFbAJb1+Ius6DngLxaG//0y5fwF4Bc8VNySNBG4G/h74FXA1sCPwHmCapKkR8bMUuxXwc+ClaZtXAVsDrwL+D3B8H/nMTM9TKL6fe9Preyti/yYV57nA5sD3gD8BewEnAf8oaUpEPNW02jjgl8DtwMXAJGA6MFHSxIhY3dd7WodFhB9+ZD8o/uEHcE3FsgvTsgua2l9OUVD+DLyk1L41xQ9IAHuW2g9PbQFMKbVvBlyT2o8stX8gtc0FhpXadwOepCgQ25TaZ6b4cRn9fn1V39KyBcBTwJgWtnNq2s4q4NWl9pHAf6dl+5TaT0ltlwIqtU8B1gJ/ADZLbTNS7LEV7/vCptfzi5+Bytz26iX3AOY3fR9LUh5TmmIvTfEnV/zdCeBfmuIvS+0Hd/rvuB99P3zYytr16nQI5tR0KOlXwIeBe4DTm2IPpRjtfikiHmo0RsRKnjtU8gHWNTcibinFrwVOTC8PK8U11j0+Ip4txd9JsWf7AmD/3A6WRTHa+A1wsKQtG+3p7LI3UhTT5Rmb/FaU5hCi2Ds/Nb1s7ttTwImRfmVT/C3AHGBX4M1N215nzz0iHsnIrVV7UOwYXFP+npITgTVUf6/3AOc2tc1Mz111Jmj1c/Gwdr2K4of/FIpDSW8AlgJvjogHm2J3T8/zK7Yzrymm7GfNDVGcxbWC4rBTefsrIqJqzmF+KaZdFwPbAf9camvM2Vyaua11+lZqmwwgaVtgPHBnRPylIn5+em707RbgIeB8Sd+TdKSkV2bmlaPX7zXlexewi6Rtmhb/Nu0IlP05PW9Xa4ZWOxcPa9cPI0IRIeAlFMfrxwPfq5j03TY9r/MDGBGPUeyhbtu8jGJ+pMpfmuK3rdp28lAppl2XUxxu+iAUk/TA+4D7gBszt7VO3yLicYpRRiPXXj+35Hl9S5/lmyjmHv6JoqDdnSbl35WZXyuy8it5rCK2MWLcvN2kbGC5eFhtIuIvEXEyxZ75W4BPNIU8np5f3LyupBcAI0oxZWN7ecsXN8U/XrXtpves2n6W9OP+PWAvSbtQTPK+CPh/5UNKLVqnb2kPfWQp114/t6b2v/UtIpZGxGHA9hSjwVNSjldKqvuQUHZ+tuFz8bCBcCKwEjgpHXJpaJzJtGfFOlOaYsr2aG6QNIHiTKrbm7Y/WtJrWtz+X9Nzf/ZyLwYEHEExAlnLc8frc6zTN4rCC6lvqVgtBXaTVFVIe/3sIuLZiPh1RJxGUcw3B965npxyP5dev1dJLwJeDdwTEU+0uD3bALh4WO0i4mHgfIof9/Lo4zsUP0zHl38E06mln08vv1WxyamSppTiNwPOSC9nleIa654pafNS/CuBoykOk/ywFN+Y2N6htZ49JyJ+ASxO250G3BD9uOAOeL+kV5dyHclzE+bNfRsJnFZeWdIeFCcB/JHi9FwkvUbSThXv1RgBrO8U2NzP5WcUk9/vkvSmpmWnp7yrvlfbgPk6DxsoX6a4puBTkr4SEY9HxBJJJwJfAu6QdCXPXecxnuL01+azdQB+DMyVNBt4kOI4/usoJmhnluK+Bbyb4jDSbyT9hOeu89gSODTtxTfMo7jO4kJJV1HMY9wRET9qsY8XU1wvAnnXdpTdDPxK0hUUh3VmAK+kOAvr5lLclyj69RFJrwV+SvHjfhDF/MiRpcnntwFflvRTitvHrKA4XXk6xRzLt9eT0y0Up8ueIelVwBPA/RHxnargiFgr6UjgJ8DNkr5HMfE9BfjfwELgrBY/D9tQdPpcYT82zAd9XOdRijk7xXyuqf2fKfZWV1L8YC8EjqpY//C0/uEUFwr+hmKveRnw78BWFesMBz5LMSp4iuLajrk0XX9Qij+BYq/9mfReMzM+gxdSHK7qAbbI/PxOTe+3F8WtXO6iOGHgXuD/UrpOpbTONsC/UVxT8TTwCPB9YHJT3G4UFz7+hmIUsYriOpCvATs1xc6n6TqP1P7B9BmuYd3rOp73utS+e8rnkZTfkpTv1r383Vnns+5rmR9D66H0hZlZJkn7ADcB50bEpzPXPZViEnvviJhff3ZmA8tzHmb990mKveSLOp2I2WDznIdZBkk7A4dQXMC3LzA7fJdZ2wS5eJjl2QU4k2IS+WrgY51Nx6wzPOdhZmbZPOdhZmbZXDzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZXDzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZXDzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsg3rdAIDafvtt49x48Z1Og0zsw3KwoULH46IsX3FbNTFY9y4cXR3d3c6DTOzDYqk+9YX48NWZmaWzcXDzMyyuXiYmVk2Fw8zM8vWUvGQtKOkr0laIGmVpJA0riJutKRLJD0s6UlJN0qaVBE3UtLZkpZJWp22u2dF3GaSTpB0r6SnJN0u6YD+dNTMzOrT6shjV+A9wArgv6oCJAmYA0wDPg4cAAwH5knasSn8UuAo4GRgOrAMmCtp96a4LwCnAl8H3g78ErhS0jtazNvMzAaAImL9QdJmEbE2/flDwMXA+Ii4txSzH3ANsE9EzEttLwCWArMi4tjUNhm4DTgyIi5LbcOAxcDdETEjtb0IeAD4YkScUnqfm4CxEfHa9eXd1dUVPlXXzCyPpIUR0dVXTEsjj0bhWI8ZwIONwpHWewy4FtivKe4ZYHYp7lngCmCqpBGpeSqwBTCr6X1mAZMkjW8ldzMzq1+dE+YTgUUV7YuBnSVtXYpbGhGrKuK2oDhE1ohbAyypiAOY0HbGZmbWL3UWjzEUcyLNlqfn0S3GjSk9PxrrHldrjnseSUdL6pbU3dPT01LiZmYDQer/Y6irs3gIqJpAaf4Y6o57noi4KCK6IqJr7Ng+b81iZmb9VGfxWE71aKAx4ljRYtzy0vPodBZXX3FmZjbI6iweiynmKZpNAO6PiJWluPGSRlXEPc1zcxyLgRHAyyviAH7XdsZmZtYvdRaPOcAOkqY0GiRtC+yblpXjhgMHluKGAQcB10fEmtR8HUUxObTpfQ4DFkXE0hpzNzOzDC3fkl3Su9MfX5+e3y6pB+iJiFsoisICYJak4ykOU51AMUdxVmM7EXGbpNnAeZKGU1wHcgwwnlKhiIj/kXQucIKkJ4BbKQrMPjz/1F8zMxtkOf+fx5VNry9Iz7cAe0XEWknTgXPSspEUxWTviHigad0jgDOA04HtgNuBaRFxa1PcScBK4BPAS4C7gfdExLUZeZuZWc1ausJ8Q+UrzM2sk9o55baTP821XWFuZmZW5uJhZmbZXDzMzCybi4eZmWXLOdvKzGyTsyHcZ6oTPPIwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsW63FQ9KbJV0v6X8kPS7pVklHNsWMlnSJpIclPSnpRkmTKrY1UtLZkpZJWi1pgaQ968zXzMz6p7biIem1wI3AcOAo4ADg18Clko5JMQLmANOAj6eY4cA8STs2bfLStJ2TgenAMmCupN3rytnMzPpnWI3bei+wObBvRKxMbTdImgy8H/gGMAPYA9gnIuYBSFoALAU+Axyb2iYDhwBHRsRlqe0WYDFwWtqOmZl1SJ2HrbYAngFWN7U/WnqfGcCDjcIBEBGPAdcC+5XWmZG2NbsU9yxwBTBV0oga8zYzs0x1Fo+Z6fmrkl4qaTtJRwFvBc5NyyYCiyrWXQzsLGnrUtzSiFhVEbcFsGuNeZuZWabaikdELAL2ohhB/BlYAZwPfCQirkhhY1J7s+XpeXSLcWN6y0PS0ZK6JXX39PRk9cHMzFpT54T5K4CrKUYH+wJvAy4ELpR0aCMMiKrVK163EreOiLgoIroiomvs2LGtpm9mZhnqnDD/N4p5iukR8Uxqu0nSC4GvSPouxcihatTQGHE0RhvLgZ37iFtesczMzAZJnXMek4DbS4Wj4b+BFwIvohiVTKxYdwJwf+ksrcXAeEmjKuKeBpbUlrWZmWWrs3g8BOwuaYum9n8AnqIYLcwBdpA0pbFQ0rYUh7nmlNaZQ3H9x4GluGHAQcD1EbGmxrzNzCxTnYetvg5cCVwr6QKKU3ZnAAcD50bE05LmAAuAWZKOpzhMdQLFXMZZjQ1FxG2SZgPnSRpOcR3IMcB44FDMzKyj6jzb6irgHcAI4BKKyfM9gI8Bx6eYtRRXi98AXAD8APgrsHdEPNC0ySOAy4DTgR8BOwHTIuLWunI2M7P+UUTVSU0bh66uruju7u50Gma2AdN6z/EcGJ38aZa0MCK6+orxXXXNzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZXDzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZXDzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZhnU6ATOzgSR1OoONk0ceZmaWzcXDzMyyuXiYmVk2Fw8zM8tWe/GQ9A5JP5W0UtLjkrol7VNaPlrSJZIelvSkpBslTarYzkhJZ0taJmm1pAWS9qw7XzMzy1dr8ZD0YeCHwELgXcCBwJXAqLRcwBxgGvBx4ABgODBP0o5Nm7sUOAo4GZgOLAPmStq9zpzNbOiT+v+wgVHbqbqSxgHnAcdHxHmlRXNLf54B7AHsExHz0noLgKXAZ4BjU9tk4BDgyIi4LLXdAiwGTkvbMTOzDqlz5HEksBa4sI+YGcCDjcIBEBGPAdcC+zXFPQPMLsU9C1wBTJU0osa8zcwsU53FYw/gLuC9kv4o6VlJSyR9rBQzEVhUse5iYGdJW5filkbEqoq4LYBda8zbzMwy1Vk8Xgq8Ajgb+CLwT8ANwNclfSLFjAFWVKy7PD2PbjFuTG9JSDo6TdJ39/T05PXAzMxaUmfx2AzYBvhwRFwcETdHxDHAdcAJabJcQFSs2zyt1WrcOiLioojoioiusWPH5vXAzMxaUmfxeCQ939DUfj3wYuB/UYwcqkYNjRFHY7SxvrjlFcvMzGyQ1Fk8FvfS3hgtrE0xEytiJgD3R8TK0rbGSxpVEfc0sKTNXM3MrA11Fo8fpOepTe1TgT9FxEMU13jsIGlKY6GkbYF907KGORTXfxxYihsGHARcHxFraszbzMwy1XlL9h8D84BvStoeuAd4N8XE+REpZg6wAJgl6XiKw1QnUIxOzmpsKCJukzQbOE/ScIrrQI4BxgOH1pizmZn1Q23FIyJC0v7AmcDnKeYn7gIOjYjvpJi1kqYD5wAXACMpisneEfFA0yaPAM4ATge2A24HpkXErXXlbGaDx1d7b1wUUXVS08ahq6sruru7O52GmeHikauTP82SFkZEV18xvquumZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy1bnRYJmtpHz6bbW4JGHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTafbWW2ifEZU1YHjzzMzCybi4eZmWVz8TAzs2wuHmZmls3Fw8zMsrl4mJlZNhcPMzPL5uJhZmbZfJGg2QbGF/nZUODiYWY2BLWzkxBRXx698WErMzPL5pGHWT/58JFtyjzyMDOzbB552CbNowez/nHxsI4b6hODZrauAT1sJek6SSHp9Kb20ZIukfSwpCcl3ShpUsX6IyWdLWmZpNWSFkjacyBztg2L1N7DzPpnwIqHpIOByRXtAuYA04CPAwcAw4F5knZsCr8UOAo4GZgOLAPmStp9oPI2M7P1G5DiIWk74Fzg0xWLZwB7AO+LiO9GxHWpbTPgM6VtTAYOAT4VERdHxE3Ae4D7gdMGIm8zM2vNQI08zgIWR8R3K5bNAB6MiHmNhoh4DLgW2K8p7hlgdinuWeAKYKqkEQORuJmZrV/txUPSHsD7gY/2EjIRWFTRvhjYWdLWpbilEbGqIm4LYNca0t3o+Pi/mQ2GWouHpOHAN4FzIuLuXsLGACsq2pen59Etxo3pJYejJXVL6u7p6WktcTMzy1L3yOOzwJbAGX3ECKg6wbJ537fVuOeJiIsioisiusaOHdtXqJmZ9VNt13lI2hk4CfgQMKJpTmJEmkR/gmLkUDVqaIw4GqON5cDOfcQtr1hmZmaDoM6Rxy7ASGAWRQFoPACOS3+eRDFnMbFi/QnA/RGxMr1eDIyXNKoi7mlgSY25m5lZhjqLx23A3hUPKArK3hQ/+HOAHSRNaawoaVtg37SsYQ7F9R8HluKGAQcB10fEmhpztzZ5ot5s01LbYauIeBSY39xeXBPIfRExP72eAywAZkk6nmJEcgLFXMZZpe3dJmk2cF6aiF8KHAOMBw6tK297jn/IzaxVg35vq4hYK2k6cA5wAcWhrgXA3hHxQFP4ERST76cD2wG3A9Mi4tZBTDmb79VkZhs7xUb8a9XV1RXd3d2D/r7egzezTmr3Z13Swojo6ivG/5+HmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2Vw8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLLVVjwkvVvS1ZLuk7Ra0t2SzpS0TVPcaEmXSHpY0pOSbpQ0qWJ7IyWdLWlZ2t4CSXvWla+ZmfVfnSOP44C/AicC04BvAMcAN0jaDECSgDlp+ceBA4DhwDxJOzZt71LgKOBkYDqwDJgrafcaczYzs34YVuO29o2IntLrWyQtB/4D2Au4GZgB7AHsExHzACQtAJYCnwGOTW2TgUOAIyPistR2C7AYOC1tx8zMOqS2kUdT4Wj4dXreIT3PAB5sFI603mPAtcB+pfVmAM8As0txzwJXAFMljagrbzMzyzfQE+ZT0vOd6XkisKgibjGws6StS3FLI2JVRdwWwK51J2pmZq0bsOIhaQeKQ0w3RkR3ah4DrKgIX56eR7cYN6aP9z1aUrek7p6eqsGQmZm1a0CKRxpB/BB4FjiivAiIqlUqXrcSt46IuCgiuiKia+zYsS1mXPFG6v/DzGxjV3vxkDSS4oyqXYCpEfGn0uLlVI8aGiOOFS3GLa9YZmZmg6TW4iFpOHA18AbgHRFxR1PIYor5jGYTgPsjYmUpbrykURVxTwNL6svazMxy1XmR4GbA5cBbgf0i4pcVYXOAHSRNKa23LbBvWlaOGw4cWIobBhwEXB8Ra+rK28zM8tV5ncf5FD/2ZwBPSnpjadmf0uGrOcACYJak4ykOU51AMZdxViM4Im6TNBs4L41mllJccDgeOLTGnM3MrB/qPGz19vR8EkWBKD8+BBARaymuFr8BuAD4AcVV6XtHxANN2zsCuAw4HfgRsBMwLSJurTFnMzPrB0VUndS0cejq6oru7u71B1bwWVNmtqFq92dd0sKI6OorxnfVNTOzbC4eZmaWzcXDzMyyuXiYmVk2Fw8zM8vm4mFmZtlcPMzMLJuLh5mZZXPxMDOzbC4eZmaWzcXDzMyyuXiYmVk2Fw8zM8vm4mFmZtlcPMzMLJuLh5mZZXPxMDOzbC4eZmaWzcXDzMyyuXiYmVk2Fw8zM8vm4mFmZtlcPMzMLJuLh5mZZXPxMDOzbC4eZmaWzcXDzMyyuXiYmVk2Fw8zM8s25IuHpJ0kXSXpMUmPS/q+pJ07nZeZ2aZsSBcPSaOAm4FXAx8A3ge8ApgnaatO5mZmtikb1ukE1uMoYBfgVRGxBEDSb4E/AB8G/r2DuZmZbbKG9MgDmAH8slE4ACJiKfBzYL+OZWVmtokb6sVjIrCoon0xMGGQczEzs2SoH7YaA6yoaF8OjK5aQdLRwNHp5UpJdw9QbmXbAw8PwvsMFe7vxm1T6u9G2Vep10Wt9vdl6wsY6sUDICraev1oIuIi4KKBS2ddkrojomsw37OT3N+N26bU302pr1Bvf4f6YasVFKOPZqOpHpGYmdkgGOrFYzHFvEezCcDvBjkXMzNLhnrxmAO8UdIujQZJ44A3p2VDxaAeJhsC3N+N26bU302pr1BjfxVRNaUwNKQLAW8HVgOfo5j/+AKwDfDaiFjZwfTMzDZZQ3rkERFPAvsAvwe+DVwOLAX2ceEwM+ucIT3yMDOzoWlIjzyGKkmbSTpB0r2SnpJ0u6QDMtbfUtKpkv4gaY2kv0j6T0lbDGTe/dVuf0vb2UXSKkkhadeByLVd/e2rpG0lnSzpF5IekfRo+vP+g5H3+rRzg1FJIyWdLWmZpNWSFkjac6Bzbkd/+yupS9JFku5Kf1fvl3S5pPGDkXd/1HXz2PT3PiT9rKUVIsKPzAdwBrAGOA7YG/gmsBZ4RwvrDgfmA/cBHwH2BA4ALgS27HTf6u5v03auA5ZRzF3t2ul+1dlX4DXAQ8CZwNuBqcDM1NePdbhPoyjuB7cI2J/i1j53AH8Etmph/cuBRynuNfdW4PsU85C7d/r7qru/wDkUtz/6KDAFOAS4E3gE2KnTfav7uy1tZxdgJfAX4GctrdPpzm9oD+BF6cfl803tNwG/bWH9fwUeH4p/EQeiv6X4Q9JfzE8O1eLRTl+BrYBRFe03Afd3uF+fAP5a/syB8cCzwKfXs+7k9H0dUWobBtwNzOn0dzYA/R1b0faytANxWqf7Vmdfm7Yzl2JHaX6rxcOHrfJNBbYAZjW1zwImtTC8/ShwZUQ8MBDJDYB2+4uk0RR3QD6OYg92qOp3XyPiyYhYVbGoG3hpfSn2Szs3GJ0BPAPMLq37LHAFMFXSiPrTbVu/+xsRPRVt9wE9wA4151mHtm8eK+kQ4HXACTlv7OKRbyLF3umSpvbF6bnXGzam45A7AfdIujgdn3xK0k2Sdh+YdNvW7/6WnAXcFRHfrjOxAVBHX5vtCdzVTlI1aOcGoxOBpRWFcTFFoR2Kc1e13lBV0m4Uo9I728xrILTV17Rjdy7wmYhYnvPGLh75xgCPRhrrlSwvLe9NYw/0sxTHGN8LHAyMBeYP0f8hsZ3+ImkP4P0UI66hrq2+Nks36XwjxTxIJ2XfYLTFdRvLh5p2+vs8koZRzEf2AJe2n1rt2u3r2RSXQszMfeNNvnhIels6w2B9j/mNVci8WWNJ4/NeBewbET+OiB8A7wS2BD7Wbn/WZzD7m84e+yZwbkQM+u1kBvm7bX7vvYCvAt+OiMv73Yn69LdftX0mg6yunL8OvAk4LCKG6v30+vvv8y0UO3bHVOwwrdeGcFfdgfYLYLcW4hrD9uXAaElq+sBHl5b35pH0/PPyYYCIeEDSXcDftZhzOwazv5+k2DP6qqTtUtuo9LyNpG0i4okW8+6Pwezr30j6e4rb59wMfLDFXAdSOzcYXQ5UjYizPpNBVssNVSWdSfHfO3wgIq6vKbe6tdPXb1KMpv5U+vc5DNg8vV4dEWt6W3mTLx7pRzznmPRiYATwcp5/bLxxfLGvPex7KE5x7G1PYW1GHv0yyP2dALwE+HPFslspbj0zYHM9g9xXACRNojhz5TbggIh4JuP9B0o7NxhdDLxL0qimeY8JwNOsOz80FLR9Q1VJJ1GcGXnsEJ+ra6evu6XHRyqWrQA+BZzX28qb/GGrfriO4h/NoU3thwGL0pkOldIPyY+At6i4bxfwt4n0VwG/rj/dtvW7v8AXKa6VKD++VFr/Q/Wm2rZ2+oqkVwA3UOwkTI+I1QOSZb52bjA6h+LapANL6w4DDgKu72vPtIPauqGqpGOB04GTIuJrA5RjXdrpa/O/zb0pdugWpT9f1efanT5PeUN8UPwoPgV8GtgL+AbFqGHfpribgCVNbRMoLsaZD+xL8Y9yEcU1EC/udN/q7m/Ftg5niF7n0U5fKc7GuZfiMM47KSbKy48RHezTVhQjhDsoTt+ckX4k7gG2LsW9jOL6gJOb1r+CYk/0QxQXCV6VPqPXdfr7qm8BvLAAAADASURBVLu/FCexrAV+UvEdTuh03+r+biu2Nx9fJDigX9jmFHf5vY/i1M7fAu/u5Yu4t6L9DcA8imPtjwHXDNUf0zr62xQz1ItHv/qaCk308RjX4X7tDFxNcYHqE+nv3LimmHEp11Ob2rekuE7noVQ0fgXs1envaiD6y3N3Bah6zO90v+r+biu21XLx8I0Rzcwsm+c8zMwsm4uHmZllc/EwM7NsLh5mZpbNxcPMzLK5eJiZWTYXDzMzy+biYWZm2f4/omyLb81wVm8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.hist(z_r, bins = 20, color = 'b')\n",
"plt.title('Robot y position \\n');"
]
},
{
"cell_type": "code",
"execution_count": 232,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEoCAYAAAB4oxv+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAZ+UlEQVR4nO3de7hk54Dv8e+PjkRECOkYzHRayLiNCNPIGSGNkBGTkBMOBkkI7bhOEMYwQzghx+VEGCPECTlHQlwiF0MQcnE/I3GLkERCxwSNpOV+kct7/njXlurq2r1r7121a7/d38/z1FO93/WuWu9616pfrXrXqtUppSBJasftJt0ASdLsGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuDUnSQ5JUpKsnHRbNhVJlnd9fsws5jmgm+eA8bVMC83gXuSSbJXkTUl+kOTaJNcn+WWS05O8Ock9+uqXJGdOqLljNckQ6gnN3seNSVYn+UiS+y50m3ratjrJ6kktXwtvyaQboOkl2Rr4JvBXwIXAx4ArgO2BnYFDuum/nVATN0UXAMd3/94aWAk8H9gnySNLKT8b47J/BTwQuHIW85wIfAf4zVhapIkwuBe3g6ih/WHgxaXvZ65JHgxcNYmGbcLOL6UcMvVHkgAfBfYH3ggcMK4Fl1JuAs6f5TxXMrugVwMcKlncdumeP9Af2gCllPNKKf8JkGRlkqk6u/V9pV/Z1Tmm+3t5/2tNNy3J0iRHJ/l9N1TzjSS7bajRSfZNcmaSK7uhne8nWTWg3p/GyZM8J8kPk9yQ5NIkhya5fW/7qAEJ8NGedVs9Q1tOT/LHJNtNM/1rSW7qH3IaVrddPtD9uaLvtXdKckLXdzcmubBbrzsNaMfuSU5LsqanD76UZO+eOuuMcU/9Tf0Gtn3fNj+gqzPt8FKSA5N8t9uuV3V98dQB9YbeTloYHnEvbmu75/sBP5ih7mrgLcCbgUuAY/qmzVqSrYCzqF/PzwK+Bfwl8GXgzGnmeRdwcNeGTwLXAU8EPpTkgaWUVw2Y7RVdnZOBM4C9qUevS4DXd3VOAu4KPLWrN9UfV8ywGkcBjwOeB/yvvrbuCDwGOLGUMp/hpqxXkDwG+BJwe+BTwKXUYZU3Ak9Mslsp5Yau7t8BpwBrqOu2Fvgz4FHUvjhlmuVeQd3mB3V/H9EzbYP7S5J/BV5O3U4fAu4APAM4KclrSimHD5htmO2khVBK8bFIH8A+QKEOh7wLeAJwlxnmKcCZ00w7ppu+fJhpwFu7svf11X1BV16AlT3le3RlJwFb9JRv1pUV4BE95Yd0ZWuB+/WU3w24HLgauENP+QFd/QNm0YebA5cBPxkw7bDu9Z4yxOssn1q3vvL09N1Hu7LbARcBtwK79dU/uqv7pp6yzwI3AksHLPfuA9pwTF+d1cDqadq9Xp8Bu3VlPwS26in/M+o4+k3ADnPdTj7G/3CoZBErpZwIvIF61HYw8BXgiiTnJ3l3kj8fcxOeC1xPParr9VHqSbp+L6O+wVeV7mgS/jQ2+8/dn88cMN/7SikX9dRfSz3K3Aq4/5xbX1/rRuD/Ag9MMjX0RPf1fj/g18AXZ/GSD+iGDg5JcjhwNnV8+w/A27s6uwL3pYb8WX3zv4Ea0vv3ld/UPfrbf/ks2jasqWX/Synlmp5lraEeICwBnjNgvrFtJ82OQyWLXCnlsCQfAPYE/gZ4BHUs9TXAC5M8uZTy7VEvN8mdgfsA3+sPj1JKSfIt1n+zPor67eCl9ZzdOjbrnh8wYHHfH1D2q+75rrNp9zQ+DLwKOJB6hQXAk4F7AW8vpdwyi9e6P3U4CmrQ/hr4CHBoKeUXXfnO3fOZ/TOXUn6b5HzgoUnuXEq5mjqktA/w4yTHd/N9o5Qy0zDQXE3bvp6ynQdMG/d20pAM7gaUemXAJ7oH3Ym2DwD7UscndxrDYu/SPf9umumDxoTvRt2n3jxg2pT1Tswx+KqHm7vneZ/4KqX8NMk3gWcmOaiUci01xAs1dGfj5FLK02aos3X3PN24+RrgoV29q0spn0xyM/Bq6gfMa4Cbk3wOOKiU8stZtnEmWwM3lFIGXZG0pqdOv7FuJw3PoZIGlVJ+Rz3ZdiPwkCR3H3LWW7vnQW+y/jfq1Jt64NUYwKCrMK4CflVKyQYejxuyraP2YeDOwDO6D76nUM8FXDyGZU313XRXqtyjrx6llBNKKY8GtqWegP0M9Sj8lAz4+jKC9m3R/U5gxrZp8TG423Ujt42J9r6xb2X6o5+pr9737i1McjvqEeCfdEdjv6CODd+9r36owzb9/gO4d5K/GGYF5mBqSGMuR3efoq7/gdSx7c2oJwrHYeqKjsf2T+g+NB4A/LwbJllHKeUPpZRTSinPBr5K3S7LZljeLcyuT6ZtH/XEZW8dLUIG9yKWZFWSh08z+eXUk0IXlFIu6ylfS18w9zine96vr/yVwA4D6h8L3JH1hz6ez+CTUf/aPR+d5C79E5PcZ9A15LMwdXnkdOs3rVLK9cBx1BOHr6aG+AnzaMuGfAP4OfXXlP0fcIcCW1BPmAKQ5PFJNu+tlGQJdegJ4AY2bC2wbf9rbMDUst+SZMueZW4HvI46/PHxIV9LE+AY9+K2J/X65wuo11D/hjr2vAvw19Q39Ev75jmDOhzwSeBH1KOxj3fjpCdSr9s9sLsi5VzgYdQTUWdx29HWlHcCTwdekWQnbruOey/qtdxP6q1cSvlCksOAfwIuSvIl6vXLS6nXgu8C/D1zvK6cemLxBuAfuq/5lwFXllKOHHL+D1OvfLkn8G+9V76MUinl1iQvAE4FTk/yKepJvN2A/0L9AH1nzyyHA3+eeo+Z1dSj592pv5o9tsx8jfkZ1BPWJ3dj+TcBny+lnDtN+85MciTwEuoJ0ROp13H/N+rQ2OvGNISkUZn09Yg+pn9Qj2pfB5xGHba4gfqDlguoJyXvP2Cee1HHRy+nDpv0X2t9P+DfgWuo45inADsyzTXe1ND9CDUkr6XeG2Ult13bu3JAG54MfL6b54/U0DqLetJt2556G3qNgdOoP/o4h3qZYmGa65c30Kc/6uZ72CznW86A67hnmGdn6jXal3f9cBH1ksGt+uo9kzqUc3G3fS+nDju9GFgyoA3H9M2/NXXYZw31g/pP120zzbXv1OG1F3V9eV23P3wd+K/DbouZpvkY3yNd50sbve4Sx98AF5ZSphuCkhY9x7i1KVlFvRzxg5NuiDQfHnFro5fk9dTL3F5MvS79AWVM49vSQjC4tdHr7qD3R+B7wEtLKYN+ASg1w+CWpMY4xi1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMQa3JDXG4JakxhjcktQYg1uSGmNwS1JjDG5JaozBLUmNMbglqTEGtyQ1xuCWpMYY3JLUGINbkhqzZNwL2Hbbbcvy5cvHvRhJ2qicc845l5VSlg6aNvbgXr58OWefffa4FyNJG5Ukl0w3zaESSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqzNh/OSmNXTL3eUsZXTukBeIRtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMQa3JDVmqOBO8ugkX07yuyRXJflekheMu3GSpPXNGNxJdgK+AmwGvAjYF/gucHSSl4y3eZKkfsP8n5PPAm4P7FVKuaYrOy3JQ4H9gCPH1ThJ0vqGGSq5A3ATcH1f+RVDzi9JGqFhgveY7vl9Se6V5K5JXgQ8AXjP2FomSRpoxqGSUsqPk6wETgRe2hXfBPz3Usrxg+ZJsgpYBbBs2bLRtFSSBAx3cnJH4ATgPGAvYHfgg8AHkzxn0DyllKNKKStKKSuWLl06yvZK0iZvmJOTb6ceYf9dKeWmruyrSe4OvDfJJ0opt46thZKkdQwzxv0Q4Ic9oT3lP4C7A9uNvFWSpGkNE9xrgJ2T3KGv/FHADcDakbdKkjStYYZK3g98Gvhckg9QLwvcG3g28J5Syh/H2D5JUp8Zj7hLKZ8B9gQ2B/439UTlrsDLgNeOtXWSpPUMc8RNKeVU4NQxt0WSNAR/+ShJjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMQa3JDXG4JakxhjcktQYg1uSGmNwS1JjDG5JaozBLUmNMbglqTEGtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0ZOriT7Jnka0muSXJVkrOTPH6cjZMkrW+o4E7yYuBk4BxgH+AZwKeBLcfXNEnSIEtmqpBkOXAE8NpSyhE9k740pjZJkjZgmCPuFwC3Ah8cc1skSUMYJrh3Bc4HnpXk4iQ3J7koycvG3DZJ0gAzDpUA9+oe7wLeAFxMHeN+f5IlpZT39s+QZBWwCmDZsmWja60Wr2Tu85YyunZIm4CUGd40SS4EdgT2LaV8tqf8VOBhwD3LBl5kxYoV5eyzzx5Rc7VoTTK4/dDQRijJOaWUFYOmDTNUcnn3fFpf+ZeBewD3nEfbJEmzNExwnzdN+dRhzq0jaoskaQjDBPeJ3fMefeV7AJeWUtaMtkmSpA0Z5uTkF4AzgA8l2Rb4OfB04EnA88fYNknSADMGdymlJHkacBjwFmAb6uWBzymlfHzM7ZMk9RnmiJtSylXAy7qHJGmCvDugJDXG4JakxhjcktQYg1uSGmNwS1JjDG5JaozBLUmNMbglqTEGtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjRnqftzSWM3nf2mXNkEecUtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMQa3JDXG4JakxhjcktQYg1uSGmNwS1Jj5hTcSb6YpCQ5dNQNkiRt2KyDO8mzgYeOoS2SpCHMKriT3BV4D/Dq8TRHkjST2R5xvxM4r5TyiXE0RpI0s6H/l/ckuwL74TCJJE3UUMGdZDPgQ8C7SykXDFF/FbAKYNmyZfNq4CYnmd/8pYymHRqv+WznVrfxprjOYzLsUMk/AncE3jZM5VLKUaWUFaWUFUuXLp1z4yRJ65vxiDvJMuCNwAuBzZNs3jN58+6E5dWllFvG1EZJUo9hjrh3ALYAjgX+0PMAOLj790PG0jpJ0nqGGeP+AfC4AeVnUMP8aOCiUTZKkjS9GYO7lHIFcGZ/eeqJhktKKetNkySNj/cqkaTGDH0dd79SyjyvW5MkzYVH3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMQa3JDVmzvfj1iKVedwmvZTRtWNTMJ++lubBI25JaozBLUmNMbglqTEGtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqzIzBneTpSU5IckmS65NckOSwJHdeiAZKktY1zBH3wcAtwBuAvwWOBF4CnJbEI3ZJWmDD/J+Te5VSft/z91lJ1gL/B1gJnD6OhkmSBpvxiLkvtKd8t3u+92ibI0mayVyHOnbrnn86qoZIkoYzzFDJOpLcG3gr8JVSytnT1FkFrAJYtmzZvBo4Mcnc5y1ldO1YSPNZ51a1uM7zbXOL++em+H7cgFkdcSfZCjgZuBl4/nT1SilHlVJWlFJWLF26dJ5NlCT1GvqIO8kWwCnADsBupZRLx9YqSdK0hgruJJsBJwCPBHYvpZw71lZJkqY1Y3B312ofBzwBeEop5Ttjb5UkaVrDHHH/G/AM4G3AtUl26Zl2qUMmkrSwhjk5+eTu+Y3At/seLxxTuyRJ05jxiLuUsnwB2iFJGpL3GpGkxhjcktQYg1uSGmNwS1JjDG5JaozBLUmNMbglqTEGtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxQ/9nwROTTLoFs9dim7XpcP+cnfn0Vymja0cPj7glqTEGtyQ1xuCWpMYY3JLUGINbkhpjcEtSYwxuSWqMwS1JjTG4JakxBrckNcbglqTGGNyS1BiDW5IaY3BLUmMMbklqjMEtSY0xuCWpMUMFd5K/SPKZJFcmuSrJZ5MsG3fjJEnrmzG4k2wJnA48ANgfeB6wI3BGkjuNt3mSpH7D/J+TLwJ2AO5fSrkIIMmPgJ8BLwYOH1/zJEn9hhkq2Rv4zlRoA5RSfgF8E3jquBomSRpsmOB+MPDjAeXnAQ8abXMkSTMZZqjkbsAfBpSvBbYZNEOSVcCq7s9rklwwt+YtmG2ByybdiEXAfrAPpmw8/ZDMdc7598Hclw2w/XQThglugDKgbNoWlVKOAo4a8rUnLsnZpZQVk27HpNkP9sEU+2Fx98EwQyV/oB5199uGwUfikqQxGia4z6OOc/d7EPCT0TZHkjSTYYL7FGCXJDtMFSRZDjy6m7YxaGZYZ8zsB/tgiv2wiPsgpQwavu6pUH9k80PgeuCfqePd/wO4M7BTKeWacTdSknSbGY+4SynXAo8HLgQ+BhwH/AJ4vKEtSQtvxiNuSdLislHfHTDJ7ZL8U5LVSW5I8sMk+w457zFJyoDHEQPq7prkW0muT7ImyeFJ7jj6NZqbufZDkq2TvKlbt8uTXNH9+2kD6h4yTX+dNJ61mrbNc74hWpItkrwryW+6bfntJI8dUG/O+9VCmWs/JFmR5Kgk5ye5LskvkxyX5D4D6q6eZpuvt39Mwjz3hUHrVZLs3FdvMvtCKWWjfQBvA24EDgYeB3wIuBXYc4h5jwF+B+zS99i+r95O1PH/k4AnAC+kXib5yUmv/3z7AfgrYA1wGPBkYI+uXwrwsr66h3Tlj+7rr79cwPXcknoPnR8DT6PekuFc4GLgTkPMfxxwBfX+PE8APttt251HtV8t9n4A3k29ncVLgd2Avwd+ClwO/EVf3dXAFwe8R7ZpuQ+6+Qvw0QHrtuVi2BcmvpONccNt13XoW/rKvwr8aIj5jwEuHaLeid0OsllP2X7dhn94y/0A3Kl/R+2Z95d9ZVPBvWSC6/oPwC3A/XrK7gPcDLx6hnkf2rX/+T1lS4ALgFNGtV810A9LB5Rt34XRW/vKVwPHTnp9R90HXd0CHDpDnYntCxvzUMkewB2AY/vKjwUeMuir32wl2Qz4W+BTpZSbeiZ9Cvgji+MmXHPuh1LKtaWU6wZMOhu41+iaODLzuSHa3sBNwCd75r0ZOB7YI8nmXfHY96sRmHM/lFJ+P6DsEuD3wL1H3M5xWoib401sX9iYg/vB1E/Di/rKz+ueh7lB1nZJLktyc5ILk/xjktv3TL8vsAV9N+EqpdxA/Uq2GG7CNYp+6PdY4Pxppv1nkluSXJLkHQs81j+fG6I9GPjFgA+q86hvzvv11Bt1f47aSG8Ml+SB1KPLnw6YvFc3Fn5jku8slvFtRtMHL+nW67okpyd5zIBlTGRfGPZeJS26G3BF6b679FjbM31DfgCcQ90IWwD7UMd6d6SOY/e+xnQ34ZppGQthvv2wjtQbiO0CPLdv0kXA64HvU79mPgl4FfBw4ImzbPNczfqGaEPOOzV96nlk/Tkm8+mHdSRZAnyQesR9dN/kzwHfpV4efA/g5cCJSZ5XSuk/Cl1o8+2DY4F/B35NHSp6LXB6kieWUs7sWcZE9oVmgjvJ7sBpQ1Q9q5SyknoTrFndHKtXKaX/6pEvJLkGOCjJO0opP+t5rTkvZ7YWuh/6lr0SeB/wsVLKcb3TBrxRT0tyKXBEkt1LKV+Z7fLmaK7rOmw/jaw/x2xUbXw/8DfAU0op6wRhKeUV67x4ciLwHeoBzqSDG+b3/n9ez59fT3Iy9Qj+UGDXnteayL7QTHAD3wIeOES9qa+6a4FtkqTvE3Gbnumz9QngIGAF9YTkhj5Zt+G2r0yjNJF+SPII6i0OTgcOHLKtnwCOAB4BLERwz+eGaGuBQZeK9ffTOParURvJjeGSHEa9PfP+pZQvz1S/lHJLkk8D70hyz1LKb4Zd1hiM9OZ4pZSrk3yedff9ie0LzQR3N/Y43bjqIOcBm1PHoXvHoKbGneZyg6z+I+yLqWNc69yEK8kW1P/u7dNzWMYGTaIfkjwE+BJ1+GjfvhOxw1ioX3nN54Zo5wH7JNmyb5z7QdQTzRf11Bv1fjVq874xXJI3Uoe+XllK+dgslr2hb6ELaRw3x+s/wp7cvjDpy3bG9eC2S3Xe3Ff+FeDcOb7me6mXRd23p+wk6u0AlvSUPZe6gf+69X6gjumvoV5JsvUsl/2qrh8et0DrehD1cq8desqWU68Wec0M8+7ctXX/nrIl1BNynxvnfrWY+qGr+8quL94wy+Uu6faTS1rvgwGvtzXwS+oQ5MT3hYnvZGPeeP8TuAF4NbASOLIL3r366n0VuKjn7+2Br1F/hPAkYC/gI928R/bNuzP1Rxqfpf5o40DqV6RPT3r9R9AP21Gv1V0LPIX1f4yweU/d73dBvSf1xzqHd2+SUxdwPe9EPfI5l3rJ197UG6T9HNiqb/veDLypb/7jqV+jX9hty890/fbwvnpD9ecEt/ec+wF4Vrcupw7Y3g/qqffsrr/2o/7w5FnA16mB/6zG++Bg4MPUHx+tBPbvXuePwGMWw74w8Z1szBvv9tQ7Gl5C/WT8EfD0AfXOBFb3/H036pH0Jd1GuR74HvWs+e0GzP9Y4Ntd3d9Sx3XX++FKg/2wsnsjTvdY3lP3eOrQ0XXdMn4C/As94b5A67oMOAG4Cri6247L++os79p/SF/5HakfOGu6bfn/gJVz7c8Jb/M59QO3/TJ20OPMnnq7UM93/Jb6AX0l9Uhzj0mv+wj6YC/q9d6Xdet2OfX8ziMXy77gTaYkqTEb8w9wJGmjZHBLUmMMbklqjMEtSY0xuCWpMQa3JDXG4JakxhjcktSY/w+wjTqqYcdT3AAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.hist(s_ypos, bins = 20, color = 'r')\n",
"plt.title('Student y Position \\n');"
]
},
{
"cell_type": "code",
"execution_count": 96,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEDCAYAAAAyZm/jAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3dfdhcVXnv8e+PkBeBRJKC9ggkEKNAOCjV1HIOKAGtoELQI76eqkgVUasWBZWi1CKKvKlVC5JaKwcoUkRKEEUChNijCZpSRQKIUELIAZSQkACBhMB9/th7zDCZZ2bP8+yZ2S+/z3XNtZ/Zs9fslbngXmvfe621FRGYmVl9bDPsCpiZ2WA58JuZ1YwDv5lZzTjwm5nVjAO/mVnNOPCbmdWMA7+Z2RhJOlTSDZIelLRR0ipJ/yppdpdycyTNl3SHpA2SVkq6WNIefa2vx/GbmY2NpHcALwNuAh4CpgOfBnYD9o2Ie0codzbwP4CLgeXALsBngecB+0XEfX2prwO/mVn+JO0J3AGcEBHnjHDMzhHxUMu+GcA9wGkRcUo/6uZUj5lZfzycbp8a6YDWoJ/uu5fkqmGXPtXLgd/MLC+SxkmaIOlFwPnAg8B3e/yOvUlSPbf3oYrJOYqS6pkwfvuYNGlq1+OenqSev/vpCaOp0dZiQv9/qwkTNvf9HHmZsu2Tw65CLUwdt2HYVSic5b9+anVE7DyW73jl3Emxds0zWc+3HGj+D35+RMxvPU7SMuDl6du7gHkRkTmAS9oWuB7YG9gzItZmLduLbfvxpaMxadJU/nTOh7set+6FE3v+7kdn9N5YtLNx+qZcvqeT3Xfd6sqvsP78j+8YdhVq46gpNw+7CoWy1/QH2t4s7cXaNc9w+dU7ZT3fkxExJ8Oh7wKmADOBE4CFkg6MiBUZq/UN4H8Cb+hX0AenemwMFj6417CrUBvfW/+yP7ysuCLi9oi4KSIuAV4N7EAyuqcrSacDxwLHRMS1faxmuQL/aHr7ZbNi1ZiuXq0G3ACUQ0Q8QpLumdXtWEknkzQQH4uIC/tdt1IFfise9/qHx1cBxSbp+cBewN1djvsocBpwckR8fRB1K0yO38pr4YN7Od8/ZK3B3/cEBkvSFcDNwC3AeuDFwPHAZuCc9JgZJI3AqRFxarrv7cBXgWuAGyTt3/S16yPitn7Ut/KBP68bu9aZg3+xNDcEbgQGYinwVuATwATgPuBG4PSmG7sCxvHsTMth6f7D0lezxcDcflS28oG/jFas2rlUo3us2Hw10H8RcQZwRpdjVpAE+eZ9RwNH96teIylN4K/Djd2yc6+/HNwQWGkCv5WDg3/5OC1UPw78ZvYHvhqoBwf+gipznt+9/upwQ1BNDvzWFw7+1eS0UDWUIvD7xm45OfhXW7uJY24MyqEUgd/MysGNQTk48BdYmfP8De71mxuD4ql04Pes3WJw8LdWbgyGq9KB34rDwd+6cWMwOIUP/L6xWx2NlTzdAFhWI688evVA61E1o1qWWdI1kkLSaS37p0r6lqTVkh6XdJ2kffOpaj1VcX1+L+VsNlw9B35J7wBe2ma/gAUkK8x9BHgzMB5YJGnXMdbTKsbB32x4egr8knYEvgJ8vM3H84ADgXdFxCURcU26bxvgk2OtqFWPg7/ZcPTa4z8TWJ4+T7LVPOD+iFjU2BER64CrgCNHX0WrMgd/qwJJR0m6XNK9kp6Q9BtJp0ua3OP3nJSm0f9vv+oKPQR+SQcC7wY+NMIh+wC3ttm/HJguaYdeK+cbu4kq5vmbLXxwLzcAVnYnAE8Df0OS7j4P+CCwUFKmOCtpJnAy8Pt+VbIh06geSeOB84GzI+I3Ixw2DVjRZv+adDsVeKzXClp9eMinldgREdE823KxpDXABSRP0bohw3ecB1wM7EmfR1xm7fF/CngO8IUOxwiIEfa3LyAdK2mZpGWbnno8Y1WsytzztzJqCfoNv0i3u3QrL+mdwMuAk/Ks10i6Bn5J00kuPz4LTJS0Y3qTl6b340h69tPafMXUdLu29YOImB8RcyJizoTx24/uX1ATVU/3NHPwt4o4KN3e3ukgSVNJBs18MiLWdDo2L1l6/DOBScBFJMG78YIkr7UW2Jckl79Pm/KzgZUR4TSPZea8v5WZpF2AU4HrImJZl8PPAu4EvtPvejVkySP9Eji4zf5FJI3BPwF3kYzhf6+kgyJiMYCkKcARwL/kU12rG+f9rV/WPr1dh5nBra7eSVJzAJ8fEfPbHZkOZLkS2Ay8t9O3SnolyaCZl0VEu1R5X3QN/BHxCHBj6/5kvhb3RsSN6fsFwBLgIkknklwJnESS4z8ztxrXWBVW6xwNB38rgNURMafbQZImkXSCZwIHRcSqLkXOJ+k8r2pKoW8LjEvfPxERG8dQ77ZGtWRDOxHxDHA4sBA4F7iCZHjTwRFxX17nycorc1aLUz9WdOnox8uBVwCvj4hfZyi2N3Acz06jHwDsn/79wX7UddRDhiJiq8ia3pg4Jn2Nicfwt1fXXn+De/9WROlY/YuBVwNviIilGYu2S6N/FRhHsvTNXfnU8NkKvzqnWSsHfyugfwDeQjLk/XFJ+zd9tioiVkmaAdwNnBoRpwI0UuXNJD0CbNvus7zkluqpg4krJwy7CpZy6scK5nXp9mSSe53Nr/eln4mkJz/0uOsefwnVPd3TzL1/K4KI2D3DMSvoMKG16bi5Y69RZ0NveczGyj1/s9448JdUnWbyZuHUj1l2DvxWKQ7+Zt058JeYe/3tOfibdVbIwO8x/DZWTv2YjayQgd8sL24AzLbmwF9yTvdk4wbAbAsHfqsVNwBmnsDVs4krJ7Bx+qZhV+NZPKGrd83B3xPArG7c47fa81WA1Y0Dv1nKDYDVhQN/Rfgmb37cAFjVFS7wewy/FYUbAKuqwgV+Gz33+vvDDYBVjQO/WUZuAKwqHPgrxr3+/nMDYO1I2lXS1yUtkbRBUkjavYfye0u6TNJqSU9I+o2kj/Wjrg78o+AncRm4AbCtzALeSvKQ9H/vpaCkOcBNwESSJ3a9HjiH5IldufMErgryhK7B8mQwS/0kIp4PIOl9wGuzFEof1H4BcH1EvKnpo0X5VzHhwG+Wo0Yj4AagfiLimVEWnQvMBo7LrzadOdVTUc71D1cjDeRUkGVwYLqdJGmppKck/V7S1yQ9px8ndI/frM+cCiqm9Zsn9dAwX72TpGVNO+ZHxPycqvKCdHsp8A3g08Ac4FRgN+BNI5QbNQf+USriYm2tnOsvHjcCpbU6Iub06bsbmZeLIuKU9O8bJY0DviRpdkTc1o8TmtmAOR1kqYfT7cKW/dem2/3yPqF7/BXnXn85+Eqg1pan22jZr3Q72pvGI3LgNysYNwK18yNgI3AY8IOm/Yem22VblRgjB/4acK+/vNwIlIuko9I/X55uXyfpIeChiFgsaQZwN3BqRJwKEBEPSzod+Kyk9cANJDd3TwEuiIi78q6nA/8YlOEGb4ODf/m5ESiFy1ren5tuF5OM1xfJbNzW+6unAo8CHwJOAB4AzgI+349KOvCblZAbgWKKCHX5fAVbcvfN+wP4cvrqOwf+GnGvv5paRwW5IbBuHPhrxsG/+twQWDeFCvx++tZgOPjXixsCa1WowF9GZbrBawZbNwTgxqBuHPhryr1+a+argnrpumSDpEMl3SDpQUkbJa2S9K+SZrccN1XSt9Knxzwu6TpJ+/av6p1Nvrd1Epy18gqeNpLm5SS8pET1ZOnxTwP+g2Q86kPAdJLV45ZK2jci7pUkYAGwB/ARkifQnAQskrRfRKzqS+1tzNzztyx8RVAtXQN/RFwCXNK8T9LPgTuAo0geDzaPZE3pQyJiUXrMEuAe4JPAR/OttuXJwd96NdJVgBuEchhtjr+xmtxT6XYecH8j6ANExDpJVwFHUvHAX4UbvA7+lgffOC6HzIE/XRt6HDAD+BLwIPDd9ON9gFvbFFsOvFvSDhHx2Bjran3m4G/94KuD4umlx38TWxYeuoskrfP79P00YEWbMmvS7VRgq8Av6VjgWICJE3fsoSrWLw7+Nii+OhieXgL/u4ApwEySRYQWSjqwae2JdsNouq1bMR+YDzBl8q4ehlMQDv42LL46GIzMgT8ibk//vEnSj0h6+J8meTL8GpJef6up6XbtGOpYClXI8zdz8Lci2bpBuHoo9aiKUT16MSIeIUn3zEp3LSfJ87eaDax0fr+cPM7frJpGFfglPR/Yi+SBApCM4d9F0kFNx0wBjkg/s5Jy8DfLRtJukr4naZ2k9ZK+L2l6xrLTJV0gaaWkDZLulHSapO37UdeuqR5JVwA3A7cA64EXA8cDm0nG8EMS3JcAF0k6kS0TuAScmX+1i6lq6Z4Gp33MOpO0HcmTszYC7yG553kaySTWl0TE4x3Kbg9cB4wHPgusBP4U+DvgRcDb8q5vlhz/UuCtwCeACcB9wI3A6emNXSLiGUmHA2eTzPCdRNIQHBwR9+VdaRs8B3+zjt5PMvBlz8ajEiXdAvwW+ACdH7ByAEmAPzQirk33LZI0DThB0nYRsSHPymaZuXsGcEaG49YAx6QvqyAHf7MRzQOWNj8fNyLukfRTkkmsnQL/hHS7vmX/IyTp+I6jI0djVDl+qy/n/M3a6jSJdXab/c2uI7kyOEPSbEk7SDoE+BjwzU5potGq9LLMk+8NHp2Re2PZUVXz/M3c87cq2LRp2146MjtJWtb0fn46D6lhGu2Hra9hy7D2tiLiSUkHApeTNBQN3wL+KmsFe1HpwG/94+BvNbM6IuZ0OabnSawAkiYBlwLPI5kouxJ4BXAKySCaD/ZW1e4c+G3UHPzN/mAtI09i7TaB9S+BucCsiGgMkf+JpHXAfEnfjIhf5VZTnOPvi4krJ3Q/qCKc8zcDOk9iva1L2X2BtU1Bv+Hn6XbvMdZtKw78NmYO/mYsAPaXNLOxQ9LuJEM1u01ifRCYKmlWy/4/S7f/L6c6/oEDv+XCwd9q7h9J1i+7UtKRkuYBV5LMezq/cZCkGZI2Szqlqex3gEeBH0p6j6SD04mwZ5M8/fCneVfWgb9P6pTuaXDwt7pKh1weAtwJXAhcTPIEwkNa1ioTyXNNtmkquwLYH/glyWzfH5JMCJsP/HlEPJN3fX1z13LlG75WVxGxEnhzl2NW0GakT0TcRrJCwkC4x2+5c8/frNgc+PuojumeBgd/s+KqfOCffK8f7DUsDv5mxVT5wG/D5eBvVjyFCvzPvXvjsKtgfbBi1c5uAMwKpFCBv4rqnOdv5QbArBgc+G3gHPzNhsuB34bCvX+z4XHgHwCne0bm4G82eA78NnTu/ZsNlgO/FYYbALPBqEXgL8IkLqd7snPwN+uvWgR+Kx/3/s36x4F/gNzr750bALP8OfBbKTj4W9FJ2k3S9yStk7Re0vclTc9YdpKksyQ9IOkJSUskvapfdXXgt9Jw79+KStJ2wA3AXsB7gHcBLwIWSdo+w1f8E8nDV04BDgceAH4sab9+1NcPYhmwiSsnsHH6pmFXo9Qawd8PfLECeT8wE9gzIu4CkHQL8FvgA8CXRyoo6aXAO4FjIuKf032LSR7gfiowL+/K1qbHX4SRPZYv9/6tQOYBSxtBHyAi7iF5Xu6RGco+BVzaVHYz8F3gUEkT865s4QK/V+i0Xjj9YwWxD3Brm/3LgdkZyt4TERvalJ0AzBp79Z6tcIG/Djy6J39uAKzPdpK0rOl1bMvn04C1bcqtAaZ2+e5OZRuf58o5/iFxrr8//LB3y0qb1EsnbHVEzOlyTLt88lYPVh/hmNGWHRX3+K1y3Pu3IVhL+575VNr35put6VC28XmuHPiHyCmf/mo0AG4EbACWk+TqW80GbstQdo90SGhr2U3AXVsXGZtaBf4ijuxx8B8MNwDWZwuA/SXNbOyQtDtwQPpZt7Ljgbc0ld0WeBtwbUTkPuKlVoHfzFcB1if/CKwArpR0pKR5wJXAfcD5jYMkzZC0WdIpjX0R8UuSoZxflfQ+Sa8mGcq5B/C3/aisA38BuNc/HG4ALC8R8ThwCHAncCFwMXAPcEhEPNZ0qIBxbB173wv8M3AacDWwG3BYRNzcj/p2HdUj6SjgHcAc4HnASuD7wBcj4tGm46YCZwFvBJ4DLAGOj4hf96HeleNRPsPTHPw9IshGKyJWAm/ucswK2ozWiYgngI+nr77L0uM/AXga+BvgMOA84IPAQknbAEgSSZ7qMOAjJP/48STrVOzah3qb9YWvAqwOsozjPyIimrtBiyWtAS4A5pIsTDQPOJDksmYRgKQlJJc6nwQ+mmelq8q9/uLwVYBVWdcef0vQb/hFut0l3c4D7m8E/bTcOuAquq9TsZU6L9vgfH/x+CrAqma0N3cPSre3p9tO61RMl7TDKM+TuyIO6bRycANgVdHzkg2SdiFZKvS6iFiW7p5GMpSpVWPG2VTgsdYP0/UujgWYOHHHXqtSWU75FJvTQFZ2PfX40577lcBmkuFHf/iIUaw1ERHzI2JORMyZMD7LswrqwymfcvBVgJVR5h6/pEkkI3dmAgdFxKqmj7utNdFtrQqzUvNVgJVJpsAvaTxwOfAK4DVtxuYvB17bpuhsYGXLBAbLyCmfcnIjYEXXNdWTjtW/GHg1cGRELG1z2AJgF0kHNZWbAhxB93UqBq5MN3id8im35iUinBKyosjS4/8HksWDvgA8Lmn/ps9WpSmfBSQzdS+SdCJJauckkhz/maOp2HPv3si6F+b+xLFScs+/Onw1YEWQJfC/Lt2enL6a/R3wuYh4RtLhwNnAucAkkobg4Ii4L6/KmlVJ6xWAGwIblK6BPyJ2z/JFEbEGOCZ9Wc7c668+Xw3YoPjRiyXi4F8fbgSsn2ob+CffGzw6o2+PtDTLjVNCljevx18yHuVjHiVUHZK2kXSSpBWSnpT0K0kdl3ZuKjtO0vGSbpX0uKQHJF0h6SXdyhY68Nd5sbZOHPytwcNFS+/zwOeAb5AMpFkKXCbp9RnLng38G8nQ+Y8BLyTDcvi1TfWUnfP91k674O/UUDFJeh7J806+FBFnp7sXSZoFfAn4YZevOBq4NCI+0/Sdt5AsnvkGmh752KrWgd95fqsD3yMorEOBCcBFLfsvAr4taY+IuKdD+QnA+pZ9j6TbjtmcWgf+snOv30bDVwWFsQ+wEbirZf/ydDub5GFWIzkX+LikHwGLgJ1IUj+rSB7ePiIH/pJz8Lc8+Kqgq50kLWt6Pz8i5o/xO6cBj0RE6xoya5o+H1FEnCJpI8kz0Bs9/DuBuem8qhHVPvBXId3j4G95q8NVwbhNPa3btToi5nQ6QNJrgIUZvmtxRMxllMvZN53vg8BngNPY0uP/NHCtpFdGxP0jla194DezbOrQGIzRz4C9Mxy3Id2uAaZKUkuvf2rT521JmgZ8BTgrIv62af8NJA/FOhE4fqTyhQ/8XqwtG/f6bRjcGGwRERuAO3ooshyYSDIEsznPPzvd3tah7IvTsr9o3hkRayTdTZcGqPCBfxCqkO4BB38rhpHmE9S1QejgGmAT8L9JFrxs+Avg1i4jeh5Mt6+gaen79EpgFnBzpxM78FeMg78VlRuEZ4uI30v6CnCSpEdJgvXbgEOAI5uPlXQ9MCMiZqVlV0j6AXCipGeAxcAfAZ8kuRI4r9O5HfjNbKhq3iCcDDxGMuv2j4HfAG+NiKtajhvH1vH6bcAngHek2/UkjceBEbGMDhz4U1VJ94B7/VYNdWgQIuJpklE5p3U5bm6bfRtIlm34fK/nLfRaPTZ6Xs/HqsprEo1dKQL/oBZrK9OzeLNw8DezdkoR+M3MLD8O/BXnXr+ZtXLgb1G1dA84+JvZsznw14SDv5k1lCbwD/JpXFXs9YODv5klShP4LR8O/mbmwD+Cqvb6IQn+bgDM6suBv8Yc/M3qqVSBf5B5fqh2r7/Bwd+sfkoV+K0/HPzN6sWBv4s69PrBeX+zOnHgt2dx8DervtIF/kHn+aE+vf4GB3+zwZD0cUlXSXpAUkj6XI/lnyPpc5J+K2mjpN9J+oGkjv8Tly7w22A4+JsNxPuB5wH/1mtBSeOBHwHvBc4B/hz4ELCK5MEtI/KDWGxEjeDvh7qY9c0+EfGMpG2B43os+wngZel33Ne0//JuBd3jz6hu6Z5m7v2b9UdEPDOG4h8CLmsJ+pk48FsmDv5mxSFpOrAb8F+S/lHSeklPSrpe0n7dypcy8A/jBi/Uu9cPHvJptbaTpGVNr2OHXJ8XpNtPATOBt5M8dH1n4Ma0YRhRphy/pF3TE8wBXgo8B9gjIla0HDcVOAt4Y3rMEuD4iPh1xn+MlYAf5m5VMO7J6KUTuToi5nQ6QNJrgIUZvmtxu4en96jRad8AHJE+eB1Jy4C7gA+TxOy2st7cnQW8FfgP4N+B17YeIEnAAmAP4CPAWuAkYJGk/SJiVcZzFdrke4NHZ2jY1Rg6B3+zrfwM2DvDcRtyONfD6fanjaAPEBH3SboD+JNOhbMG/p9ExPMBJL2PNoEfmAccCBwSEYvSY5cA9wCfBD6a8VxWEg7+ZlukAfiOAZ3uv4AngHb5ZwEdbxpnyvFnvPM8D7i/EfTTcuuAq4Ajs5ynF8PK84Nz/c2c9zcbvIh4CrgaeKWk7Rv709z+nsAvOpXP8+buPsCtbfYvB6ZL2iHHc1nBOPib9U7SHElHAf8r3TVb0lHpa7um466XdFdL8b8FtgeulnSEpLcAPwQeAb7R6bx5Bv5pJHn9VmvS7dTWDyQd27hLvumpx3OsSv+51781B3+znv0VcBlwafr+Len7y0hm9DaMoyU1HxG3AYeQpHsuBb5FcmP3gIj4XaeT5jlzV4ycb2orIuYD8wGmTN7VkbQCPNvXLLuIOBo4OsNxc0fY/3Pg4F7Pm2ePfw1Jr79Vo6ff7mpgTIaZ5wf3+jtx79+suPIM/MtJ8vytZgMrI+KxHM9lJeDgb1ZMeQb+BcAukg5q7JA0BTgi/ayS3OvvzKN+zIonc44/vfMM8PJ0+zpJDwEPRcRikuC+BLhI0olsmcAl4Mz8qmxl5Ny/WXH0cnP3spb356bbxcDcdGnRw4Gz088mkTQEB49m9bisnnv3Rta9cGK/vj4Tz+bNzg2A2fBlDvwR0TWyRcQa4Jj0ZTYiNwBmw1PK1TmLyLn+0fE9ALPBq0TgH/awThs7NwBmg1OJwF8U7vWPnRsAs/5z4LdCcgNg1j8O/Dlzrz9fbgDM8leZwF+kPL+Df/4c/M3yU5nAb9Xn3r9ZPhz4+8S9/v5xA2A2NpUK/EVK91j/uQGwMpP0Ykl/L+kWSY9JekDSAkkvHcV3zZS0QVJImtXt+EoF/qJxr38w3ABYSb2WZC39C0gWs/wQsDNwk6SXdyrYxrnAuqwHO/BbZbgBsJL5LvDSiDgnIhZFxBXAYSQPUf9Y1i+R9E7gT4AzspZx4O8z9/oHzw2AlUFErI6IaNm3DrgT2CXLd0iaCnwZOIHkWbuZVC7wO89vDY0GwI2AlYWkacB/B27PWORM4I6IuLCX8+T5zF0bgZdtHr7m4O8VQa3Avk7yDJOvdjtQ0oHAu0nSPD1x4LfacSNgAHpyExNuz/yokJ0kLWt6Pz8i5j/r+6TXAAszfNfidg9Pl3QS8E7gLyPirk5fIGkCcD7wlYi4LcM5n8WBf0Dc6y8mPxfAMlodEXO6HPMzYO8M37WhdYek44AvAp+JiG9n+I6/BqYBX5O0Y7pvu3Q7WdLkiHh0pMKVDPxFeCqXlYuvAmysImIDcEev5SS9i2Q45jkR8YWMxWYDfwz8vzaf3Qz8CthvpMKVDPxF5V5/ObgRsEGR9Cbgn4FvRcQJPRT9EvCdln2HAZ8C/gL4TafCDvwD5uBfLm4ErF8kvQq4BLgF+I6k/Zs+3hgR/9l07PXAjIiYBRARd9BydSFp9/TPm7rdI6hs4He6x/LmRsBydggwkWRUzk9bPrsX2L3p/ThyjNeVG8dfBp7UVX6eH2BjFRGfiwiN8Nq95di5rfvafN930rIde/vgwD80Dv7V4EliVkaVDvxFn8Xr4F8tbgSsLCqb4y8L3+ytptbg73sCViQO/AXg4F99bgisSCqd6oHip3sanPapl+a0kFNDNmju8ReIe/715SsCG6TK9/jLxj1/A98otv6qRY/fk7mszHw1YHmrReAvG6d8rBM3BDZWDvwF5eBvWbkhsF7VJvCXMd3j4G+j0e6+gBsDa1abwF9WDv6WBzcG1syBvwQc/K0f3BjUV60CfxnTPQ0O/jYIbgzqoVaBv+wc/G0Y3BhUT64TuCTtJul7ktZJWi/p+5Km53mOuvMELyuC1iUnPNmsd5ImS/pXSXdJelzSI5JukvQXGcpOkXSKpJ9Jejgt+zNJb8xy7twCv6TtgBuAvYD3AO8CXgQskrR91/JPDqYHUZa1ezpx8LeicoPQkwnAZuB0YB7wTpLHKV4o6fguZacDHwIWkzxj923AncAVkj7c7cR5pnreD8wE9mw8AUbSLcBvgQ8AX87xXLXntI+VyUjBv84po4h4mCTYN/uhpBcDxwBf6VD8HmBmRGxo2vdjSbuRPHD9HzqdO89UzzxgafNjvyLiHpJnSR6Z43ks5Z6/lZ2vENp6GHiq0wER8XhL0G9YBryg2wny7PHvA1zZZv9y4C1ZvmDC7fexae/dcqxSe2Ue3dPKPX+rok7Bv2pXCZJE8jD15wJvBg4F/nKUX/cqknRRR3kG/mnA2jb71wBT2xWQdCxwbPp24zW/O+9WfpdjjTq5cUDnGZ2dgNXDrkRB+LfYwr/FFnuO9QvWb37ox9f87rydMh4+SdKypvfzI2L+WOuQ+jDw9fTvp4CPRcT/6fVL0ni6P0nOv6O8h3O2yz2M2B1Nf7j5AJKWRcScnOtTSv4ttvBvsYV/iy1agvCoRMRhedSlQdJrgIUZDl0cEXOb3l8KLCVp2OcBX5f0dESc38O55wJfAy6MiIu7HZ9n4F9L0utvNZX2VwJmZlXyM2DvDMc9KzcfEQ8BD6Vvr0lHSJ4t6dsR0THXDyDpT4EFJKMqM6WI8gz8y0ny/K1mA7fleB4zs8JJb7Z2za9nsIxkSPzzgVWdDpS0L/Bj4JfAm7M0FJDvqJ4FwP6SZjZVanfggPSzbvLKl1WBf4st/Fts4d9iiyr/FoQyryoAAAUVSURBVAcBjwG/73SQpBeRpJb+Czg8Ip7IegJF5DMkMJ2k9SvgCeAzJPn+zwOTgZdExGO5nMjMrAIkfYDkZux1JD37PwLeSjIZ69MRcUbTsdcDMyJiVvr+ecDPgSkkk2Ufbvn6/4yIEWer5pbqiYjHJR1CMungQpKbutcDf+2gb2a2lV+TzHE6m+T+6GrgdpLe+9Utx47j2fF6NjAj/fsHbb57D2DFSCfOrcdvZmblkOsibWMl6RpJIem0YddlGCS9XtJPJD2WLnK3LL2KqhVJB0i6VtLv09/hZknHDLte/SZpV0lfl7RE0ob0/4Xd2xw3VdK3JK1OF/e6Lr3JVxlZfgtJr5Z0kaS7JT2Rbs9L0yDWQWECv6R3AC8ddj2GJc33XQn8B/AmktnOlwHbDbNegybpJSQ5z/Ek6z+9GfgF8E+SPjjMug3ALJIc71rg39sdkM7yXAAcBnyE5PcZT7IY4q4DqucgdP0tgONI8uKnkfwejcXOlkraYRCVLK2IGPoL2BF4EHgHyU3h04ZdpwH/+3cnuSn+18Ouy7BfwBeBTcAOLfuXAkuGXb8+/9u3afr7fen/C7u3HHNkuv/gpn3PJZkh/7Vh/xsG/Fvs3Kbcq9Jjjxn2v6HIr6L0+M8ElkfEJcOuyJAcAzwDfHPYFSmACSTT1luHpj1Cga5Q+yEinslw2Dzg/ohY1FRuHXAVFVoMMctvEcnEp1a/SLe75Fujahn6/0iSDgTeTbK2dF0dSDLx4+1pnnJz+nCGrutqV9B30u3XJL1A0o6S3g+8ms7L1NbFPsCtbfYvB6Y7xcFB6fb2odai4Ib66EVJ44HzgbMj4jfDrMuQvSB9nQX8DXA3SY7/G5K2jYi/H2blBikibk3XHbmCLZ2Bp4DjIuK7Q6tYcUyj/TC9Nel2Ksnkn9qRNBn4KknQ/7chV6fQhv3M3U8BzwG+MOR6DNs2JBPdjo6I76f7bkhHMZwk6WuRJjCrLp2NeDlJD/Y4kpTPkcA3JT0ZGRagqjjR42KIdSBpW+ASkhTPARGxechVKrShBf70Wbwnk9y4mSipeYH8iZJ2BB6NiKeHUsHBepjkMZWtK/tdSzJa4b8B9w+6UkPyRZIe/uGxZd2R6yX9EfD3ki7JmAuvqjWMvBgi1HBBREnbABcArwHeEBG3DLlKhTfMHP9MYBJwEcl/rI0XwAnp35Uam9zB8hH2N3pxdQp0+wK/iq0Xm/o5ydC9uo/R7rQY4sqo5yz5b5Isc/D2iLh+2JUpg2EG/l8CB7d5QdIYHAzc1b5o5VyRbg9t2X8osCoiHhxwfYbpQWA/Sa2PYPoz4Em25LLragGwi6TGTUwkTQGOINtiiJUi6RySrMF7I8J5/YyGluqJiEdo8xysZH4K90bEVp9V2A+BRcD5knYiWW3vKOC1wHuHWbEh+AbJxLWrJJ1LkuOfRzLH4ysRUa3n7rWQdFT658vT7eskPQQ8FBGLSYL7EuAiSSeSXBmfRHJ1eOag69tP3X4LSZ8CPg58G/itpP2bij8UEXcPsLqlUri1eiQF8IWI+Myw6zJIaa/tdJKAP5VkeOeXIuJfhlqxIZD0OpIb//uQpAPvJlmG9/yq3/NJ//tv5w9PbZI0jWRhrzeS/D5LgI9HxK8GUskB6fZbSLqRLcM3W10QEUf3pWIVULjAb2Zm/TX0CVxmZjZYDvxmZjXjwG9mVjMO/GZmNePAb2ZWMw78ZmY148BvZlYzDvxmZjXz/wHGPP+Nux7VUQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"A,V = np.meshgrid(np.linspace(0,np.pi/4,10),np.linspace(4,12,10))\n",
"Z = d/np.cos(A)*(np.sin(A)-\\\n",
" (g*d/(2*V**2*np.cos(A)))) +\\\n",
" z_0\n",
"plt.contourf(V,A*180/np.pi,Z)\n",
"plt.colorbar();"
]
},
{
"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.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}