{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Chapter 1\n",
    "\n",
    "An introduction to the Jupyter Notebook and some practice with probability ideas from Chapter 1.\n",
    "## 1.1 Probability\n",
    "### 1.1.1 Moments of Measured Data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The Jupyter Notebook has two primary types of cells \"Markdown\" cells for text (like this one) and \"Code\" cells for running python code. The cell below this one is a code cell that loads the plotting functions into the `plt` namespace and loads several functions from the `numpy` library. The last line requests that all plots show up inline in the notebook (instead of in other windows or as files on your computer)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "from numpy import array, sin, sqrt, dot, outer\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Arrays of data, and their average and sum in Python. We use some definitions from `numpy`. Notice the way these operators can be applied to arrays (with the \".\" operator)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "x = array([1,2,3])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "6"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2.0"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The formal definition (and to make sure we match with the book) is to take the sum and divide by the number of items in the sample:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2.0"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.sum()/len(x)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Higher order moments\n",
    "Operate on the sample, python does this element-by-element, then do the same thing as above. You may be surprised that to raise a power is \"**\" instead of \"^\". This is a difference between python and other languages, just something to keep in mind!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1, 4, 9])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "14"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(x**2).sum() # Note the parenthesis!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4.666666666666667"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(x**2).sum()/len(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "inputHidden": false,
    "outputHidden": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4.666666666666667"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(x**2).mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### It works for functions too\n",
    "`sin(x)` will find the `sin` of each element in the array `x`:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0.84147098, 0.90929743, 0.14112001])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sin(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6306294732311485"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sin(x).sum()/len(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6306294732311485"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sin(x).mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Variance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6666666666666666"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.var()  # Variance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.816496580927726"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.std()  # Standard Deviation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can test if two quanties are equal with the `==` operator. Not the same as `=` since the `=` is an assignment operator. This will trip you up if you are new to programming, but you'll get over it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x.std()**2 == x.var()  # Related by a square root"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example 1.1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "x_m = array([9,5,25,23,10,22,8,8,21,20])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "15.1"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x_m.mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "281.3"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(x_m**2).sum()/len(x_m)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7.300000000000002"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sqrt(281.3 - (15.1)**2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7.3"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x_m.std()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "False"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sqrt(281.3 - (15.1)**2) == x_m.std()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Close enough!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1.1.2 Probability"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example 1.2\n",
    "This is an illustration of how to implement the histogram from Example 1.2 in the text. Note the use of setting the number of bins. The `hist` command will pick for you, and you should try other values to see the impact. There is no one correct value, but the too many bins doesn't illustrate clusters of data, and too-few bins tends to oversimplify the data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAARu0lEQVR4nO3da4xcd33G8e+DYy4tEQG8LW5is7SgSoC4hFUIpUURUBoCStoSKqMWwk0WiAiQqNoEpAB5RVoVKi4iMiQiUASh4VIDiSAUEPAiho3rmARDMSgoblJiEnCIuNX01xdzUg2TmZ2z9szu8uf7kUZ75pz/zDw6c/z4zJnLSVUhSfr1d5/1DiBJmg0LXZIaYaFLUiMsdElqhIUuSY04Yb0eeMuWLbW4uLheDy9Jv5auv/76H1TVwrhl61boi4uLLC8vr9fDS9KvpSTfm7TMQy6S1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpEb0LPcmmJP+R5FNjlt0vyZVJDibZk2RxliElSdOtZg/9NcCBCcteBvywqh4JvA245HiDSZJWp1ehJzkFeA7w3glDzgGu6KavAp6RJMcfT5LUV99viv4z8HfAiROWnwzcAlBVR5McAR4K/GB4UJKdwE6A7du3H0veX0uLF3x6vSP0dvNbnrPeESQdo6l76EmeC9xeVdevNGzMvHudCqmqdlXVUlUtLSyM/SkCSdIx6nPI5anA2UluBj4MPD3Jv4yMOQRsA0hyAvAg4M4Z5pQkTTG10Kvqwqo6paoWgR3A56vqb0aG7QbO66bP7cZ4slJJWkPH/GuLSS4GlqtqN3AZ8IEkBxnsme+YUT5JUk+rKvSq+iLwxW76oqH5PwOeP8tgkqTV8ZuiktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1wkKXpEZY6JLUCAtdkhphoUtSIyx0SWqEhS5JjbDQJakRFrokNcJCl6RG9DlJ9P2TfDXJDUluSvLmMWNenORwkn3d5eXziStJmqTPGYt+Djy9qu5Oshn4SpJrquq6kXFXVtX5s48oSepjaqF3J3u+u7u6ubt4AmhJ2mB6HUNPsinJPuB24Nqq2jNm2POS7E9yVZJtM00pSZqqV6FX1S+r6gnAKcBpSR47MuSTwGJVPQ74HHDFuPtJsjPJcpLlw4cPH09uSdKIVX3Kpap+BHwROHNk/h1V9fPu6nuAJ024/a6qWqqqpYWFhWOIK0mapM+nXBaSnNRNPwB4JvDNkTFbh66eDRyYZUhJ0nR9PuWyFbgiySYG/wF8pKo+leRiYLmqdgOvTnI2cBS4E3jxvAJLksbr8ymX/cATx8y/aGj6QuDC2UaTJK2G3xSVpEZY6JLUCAtdkhphoUtSIyx0SWqEhS5JjbDQJakRFrokNcJCl6RGWOiS1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpERa6JDXCQpekRvQ5p+j9k3w1yQ1Jbkry5jFj7pfkyiQHk+xJsjiPsJKkyfrsof8ceHpVPR54AnBmktNHxrwM+GFVPRJ4G3DJbGNKkqaZWug1cHd3dXN3qZFh5wBXdNNXAc9IkpmllCRNNfUk0QBJNgHXA48E3lVVe0aGnAzcAlBVR5McAR4K/GDkfnYCOwG2b99+fMklzdXiBZ9e7wi93fyW56x3hA2h15uiVfXLqnoCcApwWpLHjgwZtzc+uhdPVe2qqqWqWlpYWFh9WknSRKv6lEtV/Qj4InDmyKJDwDaAJCcADwLunEE+SVJPfT7lspDkpG76AcAzgW+ODNsNnNdNnwt8vqrutYcuSZqfPsfQtwJXdMfR7wN8pKo+leRiYLmqdgOXAR9IcpDBnvmOuSWWJI01tdCraj/wxDHzLxqa/hnw/NlGkyStht8UlaRGWOiS1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1wkKXpEb0OafotiRfSHIgyU1JXjNmzBlJjiTZ110uGndfkqT56XNO0aPA66pqb5ITgeuTXFtV3xgZ9+Wqeu7sI0qS+pi6h15Vt1XV3m76x8AB4OR5B5Mkrc6qjqEnWWRwwug9YxY/JckNSa5J8pgJt9+ZZDnJ8uHDh1cdVpI0We9CT/JA4KPAa6vqrpHFe4GHV9XjgXcAnxh3H1W1q6qWqmppYWHhWDNLksboVehJNjMo8w9W1cdGl1fVXVV1dzd9NbA5yZaZJpUkrajPp1wCXAYcqKq3ThjzsG4cSU7r7veOWQaVJK2sz6dcngq8EPh6kn3dvNcD2wGq6lLgXOCVSY4CPwV2VFXNIa8kaYKphV5VXwEyZcw7gXfOKpQkafX8pqgkNcJCl6RGWOiS1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1os85Rbcl+UKSA0luSvKaMWOS5O1JDibZn+TU+cSVJE3S55yiR4HXVdXeJCcC1ye5tqq+MTTm2cCjusuTgXd3fyVJa2TqHnpV3VZVe7vpHwMHgJNHhp0DvL8GrgNOSrJ15mklSRP12UP/f0kWgScCe0YWnQzcMnT9UDfvtpHb7wR2Amzfvn11SSVpgsULPr3eEVbl5rc8Zy732/tN0SQPBD4KvLaq7hpdPOYmda8ZVbuqaqmqlhYWFlaXVJK0ol6FnmQzgzL/YFV9bMyQQ8C2oeunALcefzxJUl99PuUS4DLgQFW9dcKw3cCLuk+7nA4cqarbJoyVJM1Bn2PoTwVeCHw9yb5u3uuB7QBVdSlwNXAWcBD4CfCS2UeVJK1kaqFX1VcYf4x8eEwBr5pVKEnS6vlNUUlqhIUuSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1wkKXpEZY6JLUCAtdkhphoUtSIyx0SWpEn3OKXp7k9iQ3Tlh+RpIjSfZ1l4tmH1OSNE2fc4q+D3gn8P4Vxny5qp47k0SSpGMydQ+9qr4E3LkGWSRJx2FWx9CfkuSGJNckecykQUl2JllOsnz48OEZPbQkCWZT6HuBh1fV44F3AJ+YNLCqdlXVUlUtLSwszOChJUn3OO5Cr6q7qurubvpqYHOSLcedTJK0Ksdd6EkeliTd9Gndfd5xvPcrSVqdqZ9ySfIh4AxgS5JDwBuBzQBVdSlwLvDKJEeBnwI7qqrmlliSNNbUQq+qF0xZ/k4GH2uUJK0jvykqSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1wkKXpEZY6JLUCAtdkhphoUtSIyx0SWqEhS5JjZha6EkuT3J7khsnLE+Styc5mGR/klNnH1OSNE2fPfT3AWeusPzZwKO6y07g3ccfS5K0WlMLvaq+BNy5wpBzgPfXwHXASUm2ziqgJKmfqSeJ7uFk4Jah64e6ebeNDkyyk8FePNu3bz/mB1y84NPHfFtJatUs3hTNmHk1bmBV7aqqpapaWlhYmMFDS5LuMYtCPwRsG7p+CnDrDO5XkrQKsyj03cCLuk+7nA4cqap7HW6RJM3X1GPoST4EnAFsSXIIeCOwGaCqLgWuBs4CDgI/AV4yr7CSpMmmFnpVvWDK8gJeNbNEkqRj4jdFJakRFrokNcJCl6RGWOiS1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqRG9Cj3JmUm+leRgkgvGLH9xksNJ9nWXl88+qiRpJX3OKboJeBfwp8Ah4GtJdlfVN0aGXllV588hoySphz576KcBB6vqu1X1C+DDwDnzjSVJWq0+hX4ycMvQ9UPdvFHPS7I/yVVJto27oyQ7kywnWT58+PAxxJUkTdKn0DNmXo1c/ySwWFWPAz4HXDHujqpqV1UtVdXSwsLC6pJKklbUp9APAcN73KcAtw4PqKo7qurn3dX3AE+aTTxJUl99Cv1rwKOSPCLJfYEdwO7hAUm2Dl09Gzgwu4iSpD6mfsqlqo4mOR/4DLAJuLyqbkpyMbBcVbuBVyc5GzgK3Am8eI6ZJUljTC10gKq6Grh6ZN5FQ9MXAhfONpokaTX8pqgkNcJCl6RGWOiS1AgLXZIaYaFLUiMsdElqhIUuSY2w0CWpERa6JDXCQpekRljoktQIC12SGmGhS1IjLHRJaoSFLkmNsNAlqREWuiQ1olehJzkzybeSHExywZjl90tyZbd8T5LFWQeVJK1saqEn2QS8C3g28GjgBUkePTLsZcAPq+qRwNuAS2YdVJK0sj576KcBB6vqu1X1C+DDwDkjY84BruimrwKekSSziylJmqbPSaJPBm4Zun4IePKkMVV1NMkR4KHAD4YHJdkJ7Oyu3p3kW8cSGtgyet8bxEbNBT2zZe1fW23UdWau1THXKuSS48r18EkL+hT6uD3tOoYxVNUuYFePx1w5ULJcVUvHez+ztlFzwcbNZq7VMdfq/Kbl6nPI5RCwbej6KcCtk8YkOQF4EHDnLAJKkvrpU+hfAx6V5BFJ7gvsAHaPjNkNnNdNnwt8vqrutYcuSZqfqYdcumPi5wOfATYBl1fVTUkuBparajdwGfCBJAcZ7JnvmGdoZnDYZk42ai7YuNnMtTrmWp3fqFxxR1qS2uA3RSWpERa6JDViQxd6kpuTfD3JviTLY5Ynydu7nxzYn+TUNcj0h12eey53JXntyJgzkhwZGnPRHPNcnuT2JDcOzXtIkmuTfLv7++AJtz2vG/PtJOeNGzPDTP+Y5Jvd8/TxJCdNuO2Kz/mcsr0pyX8NPV9nTbjtij+BMYdcVw5lujnJvgm3ncs6S7ItyReSHEhyU5LXdPPXe/ualGtdt7EVcq3d9lVVG/YC3AxsWWH5WcA1DD4HfzqwZ43zbQL+G3j4yPwzgE+tUYanAacCNw7N+wfggm76AuCSMbd7CPDd7u+Du+kHzzHTs4ATuulLxmXq85zPKdubgL/t8Vx/B/h94L7ADcCj55lrZPk/ARet5ToDtgKndtMnAv/J4Oc/1nv7mpRrXbexFXKt2fa1offQezgHeH8NXAeclGTrGj7+M4DvVNX31vAxf0VVfYl7f+Z/+KcYrgD+fMxN/wy4tqrurKofAtcCZ84rU1V9tqqOdlevY/B9hjU3YX310ecnMOaSK0mAvwI+NKvH65nptqra203/GDjA4Fvh6719jc213tvYCuurj5lsXxu90Av4bJLrM/jZgFHjfpag7wqchR1M/kf2lCQ3JLkmyWPWMBPA71bVbTDYyIDfGTNmPdfdSxm8shpn2nM+L+d3L9Uvn3AIYT3X158A36+qb09YPvd1lsEvqD4R2MMG2r5Gcg1b121sTK412b42eqE/tapOZfBLj69K8rSR5b1+cmAeMviS1dnAv45ZvJfBYZjHA+8APrEWmVZpXdZdkjcAR4EPThgy7Tmfh3cDfwA8AbiNweGNUeu2rQEvYOW987musyQPBD4KvLaq7up7szHzZrq+JuVa721sTK412742dKFX1a3d39uBjzN4WTKsz88SzMuzgb1V9f3RBVV1V1Xd3U1fDWxOsmWNcgF8/55DT93f28eMWfN1170x9lzgr6s7cDiqx3M+c1X1/ar6ZVX9L/CeCY+5LttaBj+l8ZfAlZPGzHOdJdnMoJw+WFUf62av+/Y1Ide6b2Pjcq3l9rVhCz3Jbyc58Z5pBm943DgybDfwogycDhy556XgGpi415TkYd1xT5KcxmA937FGueBXf4rhPODfxoz5DPCsJA/uXgI+q5s3F0nOBP4eOLuqfjJhTJ/nfB7Zht93+YsJj9nnJzDm4ZnAN6vq0LiF81xn3TZ8GXCgqt46tGhdt69JudZ7G1sh19ptX7N+p3dWFwbv9t7QXW4C3tDNfwXwim46DE6+8R3g68DSGmX7LQYF/aChecO5zu8y38DgzZk/mmOWDzF4Gfc/DP6XfxmDny7+d+Db3d+HdGOXgPcO3falwMHu8pI5ZzrI4Bjhvu5yaTf294CrV3rO12B9faDbfvZ3/4i2jmbrrp/F4JML35l1tnG5uvnvu2e7Ghq7JusM+GMGL/v3Dz1vZ22A7WtSrnXdxlbItWbbl1/9l6RGbNhDLpKk1bHQJakRFrokNcJCl6RGWOiS1AgLXZIaYaFLUiP+D8O11B/Db53pAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "n, bins, patches = plt.hist(x_m,bins=7)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The hist function has several possible arguments, we use bins=7 to match the example."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1., 4., 0., 0., 0., 3., 2.])"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# an array of the counts in each bin:\n",
    "n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0.6, 3.6, 0. , 0. , 0. , 6.3, 4.8])"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "n/10.0*array([6,9,12,15,18,21,24])  # counts times each bin-center value"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "15.3"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# sum of the last cell should be the mean:\n",
    "sum(_)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "n/10.0*array([6,9,12,15,18,21,24])**2  # counts times each bin-center value"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# sum of the last cell should be the second moment:\n",
    "sum(_)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Both of these results are close to the previous value, but not exact. Remember, the historgram is a representation of the data and the agreement will improve for larger data sets."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1.2 Linear Algebra\n",
    "### 1.2.1 Vectors and Basis sets\n",
    "We'll use the qutip library later, even though this is all just linear algebra. For now, try using standard-python for vector math:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "rvec = array([1,2])  # A row vector"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "rvec"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cvec = array([[1],[2]])  # A column vector"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cvec"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cvec*rvec # Actually the outer product:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "inputHidden": false,
    "outputHidden": false
   },
   "outputs": [],
   "source": [
    "rvec*cvec # still the outer product... so this simple `*` doesn't respect the rules of linear algebra!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `dot` function properly computes the dot product that we know and love from workshop physics:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "dot(rvec,cvec)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "outer(cvec,rvec)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "dot(cvec,rvec)  # This doesn't work, because `dot` knows what shape the vectors should be"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This probably isn't your first error message, but it's important to look at what Python is telling you. First, it lists the type of error (ValueError). Then it shows where the error occured (in the dot function). This is helpful when you have larger cells with more lines in them. Obviously we already know what line caused this error since there is only one. Finally, the error is explained as follows: the shapes are not aligned, meaning the vectors don't have the right dimensions for a dot product. Some error messages are more helpful than others, but they all look like this. An error and a traceback."
   ]
  }
 ],
 "metadata": {
  "kernel_info": {
   "name": "python3"
  },
  "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.4"
  },
  "nteract": {
   "version": "0.15.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}