{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Statistical Graphics: Histograms and Box Plots\n", "\n", "- [Download the lecture notes](https://philchodrow.github.io/PIC16A/content/np_plt/plt_3.ipynb). \n", "\n", "Quite often, it's most insightful to visualize not \"raw\" data points, but rather some distributional information about the points. Histograms and box plots are three of the most common statistical graphics for visualizing distributional properties of data. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from matplotlib import pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1d Histograms" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "x = np.random.randn(1000) # 1000 random normal variables" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAM4klEQVR4nO3df6idhX3H8fen6tywHVO8SqZhV0YYjWOzENygYzjsalZHY/9wRMYITMgKyixs0NjC3A8CKWPdYMyxDKUZWF3ASoV0nc51uP5h9SquNUbXUFNNFXNbV6oMHNHv/riPeBvv9d57zj2e3G/eL7jcc57znPN8HxLfPnnuOc9NVSFJ6uV90x5AkrT+jLskNWTcJakh4y5JDRl3SWro7GkPAHDhhRfW7OzstMeQpA3l8ccf/35VzSz12GkR99nZWebm5qY9hiRtKEm+u9xjnpaRpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhk6LT6hKp7PZPYemst1j+66dynbVg0fuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSF/QbY2hGn9kmppo/LIXZIaMu6S1JBxl6SGjLskNWTcJamhFeOeZHOSryU5kuRwkluG5RckeTDJt4fv5y96zq1JjiZ5Nsk1k9wBSdI7rebI/STwR1X1QeBXgZuSbAX2AA9V1RbgoeE+w2M7gcuB7cDtSc6axPCSpKWtGPeqeqmqnhhuvwocAS4BdgAHhtUOANcNt3cA91TV61X1HHAUuHK9B5ckLW9N59yTzAIfAr4BXFxVL8HC/wCAi4bVLgFeWPS048MySdJ7ZNVxT/J+4F7gU1X1o3dbdYlltcTr7U4yl2Rufn5+tWNIklZhVXFPcg4LYb+rqr40LH45yabh8U3AiWH5cWDzoqdfCrx46mtW1f6q2lZV22ZmZkadX5K0hNW8WybAHcCRqvr8oofuB3YNt3cBX160fGeSc5NcBmwBHl2/kSVJK1nNhcM+DPwe8K0kTw7LPgPsAw4muRF4HrgeoKoOJzkIPM3CO21uqqo31n1ySdKyVox7VX2dpc+jA1y9zHP2AnvHmEuSNAY/oSpJDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGlrNVSElTcHsnkNT2/axfddObdtaHx65S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqaEV457kziQnkjy1aNmfJvlekieHr48teuzWJEeTPJvkmkkNLkla3mqO3L8AbF9i+V9X1RXD11cAkmwFdgKXD8+5PclZ6zWsJGl1Vox7VT0MvLLK19sB3FNVr1fVc8BR4Mox5pMkjWCcc+43J/nmcNrm/GHZJcALi9Y5PiyTJL2HRo373wM/D1wBvAT81bA8S6xbS71Akt1J5pLMzc/PjziGJGkpI8W9ql6uqjeq6k3gH3n71MtxYPOiVS8FXlzmNfZX1baq2jYzMzPKGJKkZYwU9ySbFt39BPDWO2nuB3YmOTfJZcAW4NHxRpQkrdXZK62Q5G7gKuDCJMeB24CrklzBwimXY8AfAFTV4SQHgaeBk8BNVfXGZEaXJC1nxbhX1Q1LLL7jXdbfC+wdZyhJ0nhWjLu02OyeQ9MeQdIqePkBSWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGlox7knuTHIiyVOLll2Q5MEk3x6+n7/osVuTHE3ybJJrJjW4JGl5qzly/wKw/ZRle4CHqmoL8NBwnyRbgZ3A5cNzbk9y1rpNK0lalRXjXlUPA6+csngHcGC4fQC4btHye6rq9ap6DjgKXLlOs0qSVmnUc+4XV9VLAMP3i4bllwAvLFrv+LDsHZLsTjKXZG5+fn7EMSRJS1nvH6hmiWW11IpVtb+qtlXVtpmZmXUeQ5LObKPG/eUkmwCG7yeG5ceBzYvWuxR4cfTxJEmjGDXu9wO7htu7gC8vWr4zyblJLgO2AI+ON6Ikaa3OXmmFJHcDVwEXJjkO3AbsAw4muRF4HrgeoKoOJzkIPA2cBG6qqjcmNLskaRkrxr2qbljmoauXWX8vsHecoSRJ4/ETqpLUkHGXpIZWPC0j6cwzu+fQVLZ7bN+1U9luRx65S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIb8Bdkb0LR+ebGkjcMjd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoa69oySY4BrwJvACeraluSC4B/BmaBY8DvVNX/jDemJGkt1uPI/Teq6oqq2jbc3wM8VFVbgIeG+5Kk99AkTsvsAA4Mtw8A101gG5KkdzFu3At4IMnjSXYPyy6uqpcAhu8XLfXEJLuTzCWZm5+fH3MMSdJi417P/cNV9WKSi4AHkzyz2idW1X5gP8C2bdtqzDkkSYuMdeReVS8O308A9wFXAi8n2QQwfD8x7pCSpLUZOe5JzkvygbduAx8FngLuB3YNq+0CvjzukJKktRnntMzFwH1J3nqdL1bVV5M8BhxMciPwPHD9+GNKktZi5LhX1XeAX15i+Q+Aq8cZSpI0Hj+hKkkNGXdJasi4S1JDxl2SGjLuktTQuJ9QlaR1M7vn0FS2e2zftVPZ7iR55C5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDXk9dzHMK1rT0vSSjxyl6SGjLskNWTcJakh4y5JDRl3SWrId8tIOuNN851vx/ZdO5HX9chdkhoy7pLUkHGXpIaMuyQ1ZNwlqaEW75bxGi+S9OM8cpekhiYW9yTbkzyb5GiSPZPajiTpnSYS9yRnAX8H/BawFbghydZJbEuS9E6TOnK/EjhaVd+pqv8D7gF2TGhbkqRTTOoHqpcALyy6fxz4lcUrJNkN7B7uvpbk2QnNciHw/Qm99nupy35An33psh/QZ1823H7kc8s+tJp9+bnlHphU3LPEsvqxO1X7gf0T2v7bgyRzVbVt0tuZtC77AX32pct+QJ996bIfMP6+TOq0zHFg86L7lwIvTmhbkqRTTCrujwFbklyW5CeAncD9E9qWJOkUEzktU1Unk9wM/CtwFnBnVR2exLZWYeKnft4jXfYD+uxLl/2APvvSZT9gzH1JVa28liRpQ/ETqpLUkHGXpIbaxz3JXyT5ZpInkzyQ5GenPdOokvxlkmeG/bkvyc9Me6ZRJLk+yeEkbybZkG9b63J5jSR3JjmR5KlpzzKOJJuTfC3JkeHv1i3TnmlUSX4yyaNJ/mvYlz8b6XW6n3NP8tNV9aPh9h8CW6vqk1MeayRJPgr8+/AD688BVNWnpzzWmiX5IPAm8A/AH1fV3JRHWpPh8hr/DfwmC2/7fQy4oaqenupgI0jy68BrwD9V1S9Oe55RJdkEbKqqJ5J8AHgcuG6D/pkEOK+qXktyDvB14JaqemQtr9P+yP2tsA/O45QPU20kVfVAVZ0c7j7CwucHNpyqOlJVk/pE8nuhzeU1quph4JVpzzGuqnqpqp4Ybr8KHGHhk/IbTi14bbh7zvC15m61jztAkr1JXgB+F/iTac+zTn4f+JdpD3GGWuryGhsyJB0lmQU+BHxjupOMLslZSZ4ETgAPVtWa96VF3JP8W5KnlvjaAVBVn62qzcBdwM3TnfbdrbQvwzqfBU6ysD+npdXsxwa24uU1NB1J3g/cC3zqlH+1byhV9UZVXcHCv86vTLLmU2YtfhNTVX1klat+ETgE3DbBccay0r4k2QX8NnB1ncY/MFnDn8lG5OU1TkPD+el7gbuq6kvTnmc9VNUPk/wHsB1Y0w+9Wxy5v5skWxbd/TjwzLRmGVeS7cCngY9X1f9Oe54zmJfXOM0MP4S8AzhSVZ+f9jzjSDLz1jvhkvwU8BFG6NaZ8G6Ze4FfYOHdGd8FPllV35vuVKNJchQ4F/jBsOiRjfjOnySfAP4WmAF+CDxZVddMd6q1SfIx4G94+/Iae6c80kiS3A1cxcLlZV8GbquqO6Y61AiS/Brwn8C3WPhvHeAzVfWV6U01miS/BBxg4e/W+4CDVfXna36d7nGXpDNR+9MyknQmMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWro/wFcTBYf9QMRiQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1)\n", "h = ax.hist(x) # automatically chooses a binwidth, usually pretty smart" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAANlElEQVR4nO3dcaid913H8fdn3exkm9iRtMQ2eIsEWSfaQajCRCqda1ll6f6oZIoELMRBix0oLN3AOiUQEacgToy0LMLWGujKCp1zNU7q/ti621K3pmldWLM2S2iy1bkVcZLu6x/3ibtN782995x7cu755v2CyznP73nOeb4/knzuL7/nOb+TqkKS1Mvrpl2AJGn9Ge6S1JDhLkkNGe6S1JDhLkkNvX7aBQBs2rSp5ubmpl2GJM2Uxx9//NtVtXmpfRsi3Ofm5pifn592GZI0U5J8c7l9TstIUkOGuyQ1ZLhLUkOGuyQ1ZLhLUkOGuyQ1ZLhLUkOGuyQ1ZLhLUkMb4hOq0kY2t+fhqZz32L6bp3Je9eDIXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSHDXZIaMtwlqSG/IFszYVpfUi3NKkfuktSQ4S5JDRnuktSQ4S5JDRnuktTQiuGeZGuSLyQ5kuRwkjuH9rcmeSTJ14fHyxa95q4kR5M8m+TGSXZAkvRaqxm5nwF+v6reBvwScHuSa4A9wKGq2gYcGrYZ9u0E3g7cBHw8ySWTKF6StLQVw72qTlbVE8Pz7wNHgCuBHcCB4bADwC3D8x3A/VX1g6p6DjgKXLfehUuSlremOfckc8A7gC8DV1TVSVj4BQBcPhx2JfDCopcdH9okSRfIqsM9yZuBB4APVtX3znfoEm21xPvtTjKfZP706dOrLUOStAqrCvckb2Ah2D9ZVZ8eml9MsmXYvwU4NbQfB7YuevlVwIlz37Oq9lfV9qravnnz5lHrlyQtYTV3ywS4BzhSVR9btOshYNfwfBfwmUXtO5NcmuRqYBvw2PqVLElayWoWDnsn8NvA15I8ObR9GNgHHExyG/A8cCtAVR1OchB4moU7bW6vqlfWvXJJ0rJWDPeq+iJLz6MD3LDMa/YCe8eoS5I0Bj+hKkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1NBqVoWUNAVzex6e2rmP7bt5aufW+nDkLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1JDhLkkNGe6S1NCK4Z7k3iSnkjy1qO2PknwryZPDz3sW7bsrydEkzya5cVKFS5KWt5qR+yeAm5Zo/4uqunb4+SxAkmuAncDbh9d8PMkl61WsJGl1Vgz3qnoUeGmV77cDuL+qflBVzwFHgevGqE+SNIJx5tzvSPLVYdrmsqHtSuCFRcccH9okSRfQqOH+N8DPANcCJ4E/H9qzxLG11Bsk2Z1kPsn86dOnRyxDkrSUkcK9ql6sqleq6ofA3/GjqZfjwNZFh14FnFjmPfZX1faq2r558+ZRypAkLWOkcE+yZdHm+4Czd9I8BOxMcmmSq4FtwGPjlShJWqvXr3RAkvuA64FNSY4DdwPXJ7mWhSmXY8DvAlTV4SQHgaeBM8DtVfXKZEqXJC1nxXCvqvcv0XzPeY7fC+wdpyhJ0nhWDHdpsbk9D0+7BEmr4PIDktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktSQ4S5JDRnuktTQiuGe5N4kp5I8tajtrUkeSfL14fGyRfvuSnI0ybNJbpxU4ZKk5a1m5P4J4KZz2vYAh6pqG3Bo2CbJNcBO4O3Daz6e5JJ1q1aStCorhntVPQq8dE7zDuDA8PwAcMui9vur6gdV9RxwFLhunWqVJK3SqHPuV1TVSYDh8fKh/UrghUXHHR/aXiPJ7iTzSeZPnz49YhmSpKWs9wXVLNFWSx1YVfurantVbd+8efM6lyFJF7dRw/3FJFsAhsdTQ/txYOui464CToxeniRpFKOG+0PAruH5LuAzi9p3Jrk0ydXANuCx8UqUJK3V61c6IMl9wPXApiTHgbuBfcDBJLcBzwO3AlTV4SQHgaeBM8DtVfXKhGqXJC1jxXCvqvcvs+uGZY7fC+wdpyhJ0nj8hKokNWS4S1JDK07LSLr4zO15eCrnPbbv5qmctyNH7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ35BdkzaFpfXixpdjhyl6SGHLlLYzj2xt887/65//nUBapEejXDXRe9aQb0+c7tLwaNw2kZSWrIkbsuCiuNzqVuHLlLUkOGuyQ15LSMtEF5J47G4chdkhoy3CWpIadlpAma5F063iOv8xkr3JMcA74PvAKcqartSd4K/AMwBxwDfqOq/nO8MqXp8TZKzaL1mJb51aq6tqq2D9t7gENVtQ04NGxLki6gScy57wAODM8PALdM4BySpPMYN9wL+HySx5PsHtquqKqTAMPj5Uu9MMnuJPNJ5k+fPj1mGZKkxca9oPrOqjqR5HLgkSTPrPaFVbUf2A+wffv2GrMOSdIiY4V7VZ0YHk8leRC4DngxyZaqOplkC3BqHeqUzsuLntKrjTwtk+RNSd5y9jnwbuAp4CFg13DYLuAz4xYpSVqbcUbuVwAPJjn7Pp+qqs8l+QpwMMltwPPAreOXKUlai5HDvaq+AfzCEu3fAW4YpyhJ0nhcfkCSGjLcJakh15bRhuFaKdL6Mdx1wYxzu6K3Okpr47SMJDXkyF3ShjG35+GpnPfYvpunct5JMtylhvyKPhnuWjfOi0sbh3PuktSQ4S5JDTktI12EnJPvz5G7JDXkyF1r4kXTi4OfFp59jtwlqSHDXZIaclpG0pp4MXY2OHKXpIYcuUtaV16M3RgcuUtSQ47cG3JOVJIjd0lqyJH7GKa19vS4nBOV+nPkLkkNGe6S1JDhLkkNOeeuV3FhMKkHw30GGcCSVmK4S2phnEHP3J7p3SV2bN/NE3lfw13SBeMH7C4cL6hKUkOGuyQ15LTMBuVFU12MnLZZP47cJamhFiP3jbrGi6NvaX35b2r1WoS7JI1jktNBK/9C+q+R3/t8JjYtk+SmJM8mOZpkz6TOI0l6rYmM3JNcAvw18GvAceArSR6qqqcncb5JLWHrxR1JMJvTQZOalrkOOFpV3wBIcj+wA5hIuJ/PdP+7JUnTMalwvxJ4YdH2ceAXFx+QZDewe9h8OcmzE6plE/Dt5Xf/+rJ7sv61jGOFfsyULn3p0g/o05fZ68dHl02a1fTlp5fbMalwX6raetVG1X5g/4TO/6NCkvmq2j7p80xal35An7506Qf06UuXfsD4fZnUBdXjwNZF21cBJyZ0LknSOSYV7l8BtiW5OsmPATuBhyZ0LknSOSYyLVNVZ5LcAfwTcAlwb1UdnsS5VmHiUz8XSJd+QJ++dOkH9OlLl37AmH1JVa18lCRppri2jCQ1ZLhLUkPtwz3JnyT5apInk3w+yU9Nu6ZRJfmzJM8M/XkwyU9Ou6ZRJLk1yeEkP0wyk7etdVleI8m9SU4leWratYwjydYkX0hyZPi7dee0axpVkjcmeSzJvw99+ehI79N9zj3JT1TV94bnvwdcU1UfmHJZI0nybuBfhgvWfwpQVR+acllrluRtwA+BvwX+oKrmp1zSmgzLa/wHi5bXAN4/qeU1JinJrwAvA39fVT837XpGlWQLsKWqnkjyFuBx4JYZ/TMJ8KaqejnJG4AvAndW1ZfW8j7tR+5ng33wJs75MNUsqarPV9WZYfNLLHx+YOZU1ZGqmtQnki+E/19eo6r+Fzi7vMbMqapHgZemXce4qupkVT0xPP8+cISFT8rPnFrw8rD5huFnzbnVPtwBkuxN8gLwW8AfTruedfI7wD9Ou4iL1FLLa8xkkHSUZA54B/Dl6VYyuiSXJHkSOAU8UlVr7kuLcE/yz0meWuJnB0BVfaSqtgKfBO6YbrXnt1JfhmM+ApxhoT8b0mr6McNWXF5D05HkzcADwAfP+V/7TKmqV6rqWhb+d35dkjVPmbX4so6qetcqD/0U8DBw9wTLGctKfUmyi4XVzm6oDXzBZA1/JrPI5TU2oGF++gHgk1X16WnXsx6q6rtJ/hW4CVjTRe8WI/fzSbJt0eZ7gWemVcu4ktwEfAh4b1X997TruYi5vMYGM1yEvAc4UlUfm3Y940iy+eydcEl+HHgXI+TWxXC3zAPAz7Jwd8Y3gQ9U1bemW9VokhwFLgW+MzR9aRbv/EnyPuCvgM3Ad4Enq+rG6Va1NkneA/wlP1peY++USxpJkvuA61lYXvZF4O6qumeqRY0gyS8D/wZ8jYV/6wAfrqrPTq+q0ST5eeAAC3+3XgccrKo/XvP7dA93SboYtZ+WkaSLkeEuSQ0Z7pLUkOEuSQ0Z7pLUkOEuSQ0Z7pLU0P8BacNsICh0qSYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ax.hist(x, bins = 40) # same data, in finer bins\n", "fig" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "When comparing distributions of several data sets, it's often useful to use *transparency* via the keyword `alpha`. `alpha = 1` generates fully opaque data markings, while `alpha = 0` is fully transparent (and therefore invisible). " ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAARYklEQVR4nO3dfYylZXnH8e9P0IpYAsjsdpXdrqabVWIV6MRqSczWhQbUsPyDwUSzsSSbJtZqY6OrJm363zZtrCYa0w2+jJFaKWLYmGLFpcQ0UeqCqOhK11IFZNxd0S2+pRa9+sc8C4fZmZ0zM+ftnvl+ksnzcs6ZcyWz+5t7rnPfz5OqQpLUnqeNuwBJ0soY4JLUKANckhplgEtSowxwSWrUmaN8swsuuKC2bt06yreUpObdfffdP6yqqfnnRxrgW7du5dChQ6N8S0lqXpLvLXTeFookNcoAl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY0ywCWpUQa4JDVqpCsxtTZ84VtHn9i//KKNY6xEWt8cgUtSowxwSWqUAS5JjTLAJalRSwZ4ku1J7u35eizJ25Kcn+T2JEe67XmjKFiSNGfJAK+q+6vq4qq6GPg94OfAZ4C9wMGq2gYc7I4lSSOy3BbKTuC/qup7wC5gpjs/A1wzyMIkSae33Hng1wGf7PY3VtUsQFXNJtmw0AuS7AH2AGzZsmWlda4v99/25P72q8ZXh6SJ1vcIPMkzgKuBf17OG1TV/qqarqrpqalTbukmSVqh5bRQrgLuqaqTy/COJtkE0G2PDbo4SdLiltNCeT1Ptk8ADgC7gX3d9tYB1qUJ4JJ5abL1NQJP8izgCuCWntP7gCuSHOke2zf48iRJi+lrBF5VPweeM+/co8zNSpEkjYFXI5x0vTNSYOJmpdhmkcbHpfSS1CgDXJIaZQtlnNbYgp311E6586E7F31sx+YdI6tD65sjcElqlAEuSY2yhTIp5s82UbNsr2hUHIFLUqMcgWtdcFSstcgRuCQ1ygCXpEYZ4JLUKANckhplgEtSo5yFoqFYT8vqpXFxBC5JjTLAJalRtlA0dJPeTjndIh9pkjkCl6RG9XtT43OT3Jzk20kOJ3lFkvOT3J7kSLc9b9jFSpKe1O8I/P3A56rqhcBLgcPAXuBgVW0DDnbHkqQRWTLAk5wDvBL4MEBV/bKqTgC7gJnuaTPANcMqUpJ0qn5G4C8AjgMfTfLVJDckORvYWFWzAN12wxDrlCTN00+AnwlcCnyoqi4BfsYy2iVJ9iQ5lOTQ8ePHV1imJGm+fgL8YeDhqrqrO76ZuUA/mmQTQLc9ttCLq2p/VU1X1fTU1NQgapYk0UeAV9UPgIeSbO9O7QS+BRwAdnfndgO3DqVCSdKC+l3I8xbgxiTPAB4A3sRc+N+U5HrgQeDa4ZQoSVpIXwFeVfcC0ws8tHOw5UiS+uVSej1F77J3SZPNpfSS1CgDXJIaZQtFfbG1Ik0eR+CS1CgDXJIaZQulZfff9uT+9qvGV8cE8eYMWk8cgUtSowxwSWqULRQ1xzaJNMcRuCQ1ygCXpEbZQhm2Qc8U6f1+ktY1R+CS1ChH4NIEON0Hszs27xhZHWqLI3BJapQBLkmNsoUySn4AKWmAHIFLUqMMcElqVF8tlCTfBX4C/Ap4vKqmk5wPfArYCnwXeF1V/Xg4ZTbAKwP2pffGEJdftHGMlUjtW84I/A+r6uKqOnl3+r3AwaraBhzsjiVJI7KaFsouYKbbnwGuWX05kqR+9TsLpYDPJyngH6pqP7CxqmYBqmo2yYaFXphkD7AHYMuWLQMoWWuF7ZTVW2wBkIt/1od+A/yyqnqkC+nbk3y73zfown4/wPT0dK2gRknSAvpqoVTVI932GPAZ4GXA0SSbALrtsWEVKUk61ZIBnuTsJL95ch/4I+A+4ACwu3vabuDWYRUpSTpVPy2UjcBnkpx8/j9W1eeSfAW4Kcn1wIPAtcMrU5I035IBXlUPAC9d4PyjwM5hFCVJWprXQlkr5l9nxcVEa4b3ANViXEovSY0ywCWpUQa4JDXKAJekRhngktQoA1ySGmWAS1KjDHBJapQBLkmNMsAlqVEupddTbqwgqR2OwCWpUQa4JDXKFspa1Xt1wp4rE55sl0z6PSi9Ap+0NEfgktQoA1ySGmULRVqDTteC2rF5x8jq0HA5ApekRvU9Ak9yBnAI+H5VvTbJ+cCngK3Ad4HXVdWPh1Gk1r7eueiT/gHravjhrAZpOSPwtwKHe473AgerahtwsDuWJI1IXwGe5ELgNcANPad3ATPd/gxwzWBLkySdTr8j8PcB7wB+3XNuY1XNAnTbDQu9MMmeJIeSHDp+/PiqipUkPWnJAE/yWuBYVd29kjeoqv1VNV1V01NTUyv5FpKkBfTzIeZlwNVJXg08EzgnySeAo0k2VdVskk3AsWEWKkl6qiVH4FX1rqq6sKq2AtcBd1TVG4ADwO7uabuBW4dWZWvuv+3JL0kaktXMA98HXJHkCHBFdyxJGpFlrcSsqjuBO7v9R4Gdgy9JktQPl9Jr4nzhW0e578QJAC7efO6Yq5Eml0vpJalRBrgkNcoAl6RGGeCS1CgDXJIaZYBLUqMMcElqlAEuSY1yIc968JRrslw6tjIkDZYjcElqlCPwdar3HpSS2uQIXJIaZYBLUqNsoawzFzxyxxP7P3zuq8ZYiaTVcgQuSY0ywCWpUbZQlqN3PvX2q8ZXhyThCFySmmWAS1KjlgzwJM9M8h9Jvpbkm0n+ujt/fpLbkxzptucNv1xJ0kn9jMD/F3hVVb0UuBi4MsnLgb3AwaraBhzsjiVJI7JkgNecn3aHT+++CtgFzHTnZ4BrhlKhJGlBfc1CSXIGcDfwO8AHq+quJBurahagqmaTbFjktXuAPQBbtmwZTNUait5FPr36WfBz34kvLfrYi899xYpr0mjd+dCdC57fsXnHSOtQf/r6ELOqflVVFwMXAi9L8uJ+36Cq9lfVdFVNT01NrbROSdI8y5qFUlUngDuBK4GjSTYBdNtjA69OkrSoJVsoSaaA/6uqE0nOAi4H/gY4AOwG9nXbW4dZ6MR5yk0SdDqna69o9BZrk6g9/fTANwEzXR/8acBNVfXZJF8CbkpyPfAgcO0Q65QkzbNkgFfV14FLFjj/KLBzGEVJkpbmSkxJapQBLkmNMsAlqVEGuCQ1ygCXpEZ5Q4c16t6HToy7BElD5ghckhplgEtSowxwSWqUAS5JjTLAJalRBrgkNcoAl6RGGeCS1CgDXJIaZYBLUqNcSq9l8fZo0uRwBC5JjTLAJalR/dyVfjPwceC3gF8D+6vq/UnOBz4FbAW+C7yuqn48vFK1HvVeVfHizeeOsZL17XR3st+xecfI6tBT9TMCfxx4e1W9CHg58OYkFwF7gYNVtQ042B1LkkZkyQCvqtmquqfb/wlwGHgesAuY6Z42A1wzrCIlSada1iyUJFuBS4C7gI1VNQtzIZ9kwyKv2QPsAdiyZctqah2P+28bdwWStKC+P8RM8mzg08Dbquqxfl9XVfurarqqpqemplZSoyRpAX0FeJKnMxfeN1bVLd3po0k2dY9vAo4Np0RJ0kL6mYUS4MPA4ap6b89DB4DdwL5ue+tQKhw1WyaSGtFPD/wy4I3AN5Lc2517N3PBfVOS64EHgWuHU6IkaSFLBnhV/TuQRR7eOdhyJEn9ciWmJDXKAJekRhngktQoLye7hvReN6QfFzxyx5AqGQ6viyI9lSNwSWqUI/B16p5fHFn0sUvP2jbCStQ6r1Q4Po7AJalRBrgkNcoAl6RGGeCS1CgDXJIa5SwUnWL+DJXHTpw1pkoknY4jcElqlAEuSY1avy2U3hs3bL9qfHU04JxHv/7E/mPPeckYK1mYS+y1XjkCl6RGGeCS1Kj120KRNHSLXSfFa6QMhiNwSWrUkgGe5CNJjiW5r+fc+UluT3Kk25433DIlSfP100L5GPAB4OM95/YCB6tqX5K93fE7B1/eiPTOSFETlnvzCmktWnIEXlVfBH407/QuYKbbnwGuGXBdkqQlrPRDzI1VNQtQVbNJNiz2xCR7gD0AW7ZsWeHbaTGORKX1a+gfYlbV/qqarqrpqampYb+dJK0bKw3wo0k2AXTbY4MrSZLUj5UG+AFgd7e/G7h1MOVIg3PvQyee+JLWon6mEX4S+BKwPcnDSa4H9gFXJDkCXNEdS5JGaMkPMavq9Ys8tHPAtUiSlmF9LaVfh/O959+cQZoEiy2xB5fZL4dL6SWpUQa4JDVqfbVQNHCTfrMHaS1zBC5JjTLAJalRtlDWgEmZaWI7RRotR+CS1CgDXJIatTZaKL0LdLZfNb469ITedgqMrqXidU/a5yKf/jkCl6RGrY0R+Bqy2Ajy4s3nLvmcUZg/sl7J6/yAUyvl6PypHIFLUqMMcElq1NpuoTRy9UE/eJsMvT+H3paV1K9Rt3gcgUtSowxwSWpUmy2U07VGRtw2ufPH31rw/I7zLhro+9hmWZ1+ZvdobVrLM1ccgUtSowxwSWrUqlooSa4E3g+cAdxQVcO7O/2Ezyh54PjPnth/wdTZi7ZWFnztL5587aVnbRtoXZNoscVAp1vgM6yFQIu1VpyR0p7TtUpW8poW2isrHoEnOQP4IHAVcBHw+iSDbfxKkha1mhbKy4DvVNUDVfVL4J+AXYMpS5K0lNW0UJ4HPNRz/DDw+/OflGQPsKc7/GmS+1fxnoNwAfDDMdfQD+scLOscLOscrKXq/O2FTq4mwLPAuTrlRNV+YP8q3megkhyqqulx17EU6xws6xws6xyslda5mhbKw8DmnuMLgUdW8f0kScuwmgD/CrAtyfOTPAO4DjgwmLIkSUtZcQulqh5P8qfAvzI3jfAjVfXNgVU2PBPTzlmCdQ6WdQ6WdQ7WiupM1Slta0lSA1yJKUmNMsAlqVHrOsCT/EWSSnLBuGtZSJK/TfLtJF9P8pkkE7WmO8mVSe5P8p0ke8ddz0KSbE7yb0kOJ/lmkreOu6bFJDkjyVeTfHbctSwmyblJbu7+XR5O8opx17SQJH/e/bzvS/LJJM8cd00AST6S5FiS+3rOnZ/k9iRHuu15/X6/dRvgSTYDVwAPjruW07gdeHFVvQT4T+BdY67nCQ1dSuFx4O1V9SLg5cCbJ7ROgLcCh8ddxBLeD3yuql4IvJQJrDfJ84A/A6ar6sXMTbK4brxVPeFjwJXzzu0FDlbVNuBgd9yXdRvgwN8D72CBxUeToqo+X1WPd4dfZm6u/aRo4lIKVTVbVfd0+z9hLnCeN96qTpXkQuA1wA3jrmUxSc4BXgl8GKCqfllVk3qh+jOBs5KcCTyLCVmjUlVfBH407/QuYKbbnwGu6ff7rcsAT3I18P2q+tq4a1mGPwYm6ZKMC11KYeKCsVeSrcAlwF3jrWRB72NuQPHrcRdyGi8AjgMf7Vo9NyQ5e9xFzVdV3wf+jrm/rmeB/6mqz4+3qtPaWFWzMDfgADb0+8I1G+BJvtD1v+Z/7QLeA/zluGuEJes8+Zz3MNcKuHF8lZ6ir0spTIokzwY+Dbytqh4bdz29krwWOFZVd4+7liWcCVwKfKiqLgF+xjL+3B+Vroe8C3g+8Fzg7CRvGG9Vw9HmLdX6UFWXL3Q+ye8y94P9WhKYa0vck+RlVfWDEZYILF7nSUl2A68FdtZkTdpv5lIKSZ7OXHjfWFW3jLueBVwGXJ3k1cAzgXOSfKKqJi10HgYerqqTf8HczAQGOHA58N9VdRwgyS3AHwCfGGtVizuaZFNVzSbZBBzr94VrdgS+mKr6RlVtqKqtVbWVuX+Ul44jvJfS3TDjncDVVfXzcdczTxOXUsjcb+kPA4er6r3jrmchVfWuqrqw+/d4HXDHBIY33f+Rh5Js707tBPq/c8noPAi8PMmzup//Tibww9YeB4Dd3f5u4NZ+X7hmR+BrxAeA3wBu7/5a+HJV/cl4S5rT0KUULgPeCHwjyb3duXdX1b+MsaaWvQW4sful/QDwpjHXc4qquivJzcA9zLUev8qELKlP8klgB3BBkoeBvwL2ATcluZ65Xz7X9v39JuuvcklSv9ZdC0WS1goDXJIaZYBLUqMMcElqlAEuSY0ywCWpUQa4JDXq/wF19eQfEU+/9QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "x1 = np.random.normal(0, 0.8, 1000)\n", "x2 = np.random.normal(-2, 1, 1000)\n", "x3 = np.random.normal(3, 2, 1000)\n", "\n", "kwargs = dict(bins=40, alpha = 0.3)\n", "\n", "plt.hist(x1, **kwargs)\n", "plt.hist(x2, **kwargs)\n", "plt.hist(x3, **kwargs);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Another excellent way to visualize distributions of data is via boxplots. The orange line gives the median, the vertical bounds of the box the quartiles, the \"whiskers\" show the range of the bulk of the distribution, and the individual points are alleged to be \"outliers\". To create multiple boxplots (almost always a good idea), pass a list of 1d arrays to the function call. " ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAATEklEQVR4nO3df2xdZ33H8c/Xjh0nDSWx6oauSRomIXRjbyrrFftBlczQTR1DYE2bhBtQtli4UpW7sG7KYPePlj8cqdnEhpxpVYLNiBZfNgELKEIbDbkIWRoMpw1QelmHCGkdCjFK1oIr107y3R+5sWrXjo9zHvvc5973Szqy7/H1c766N/eTx895nnPM3QUAiFdT1gUAANIhyAEgcgQ5AESOIAeAyBHkABC5NVkc9I477vDt27dncWgAiNaZM2d+7u4d8/dnEuTbt2/X2NhYFocGgGiZ2fmF9jO0AgCRSxzkZjZsZhfN7NnX7Ws3s6fM7H+rXzetTJkAgMUsp0f+z5IenLfvY5K+5u5vk/S16mMAwCpKHOTu/g1Jl+bt/oCkz1a//6yknkB1AQASSjtGvtndX5Kk6tc7F3uimfWb2ZiZjU1MTKQ8LADghlU72enuR9w97+75jo43zJ4BgFtWKBTU1tYmM1NbW5sKhULWJa2qtEH+MzO7S5KqXy+mLwkAkisUCnryySd18OBBTU5O6uDBg3ryyScbKszTBvmXJe2pfr9H0pdStgcAy3L06FE98cQTevTRR7V+/Xo9+uijeuKJJ3T06NGsS1s1lvR65GZWkvS7ku6Q9DNJj0k6IenfJG2T9IKkP3H3+SdE3yCfzzsLggCEYGaanJzU+vXrZ/e9+uqruu2221Rv91swszPunp+/P/HKTnfvXeRH77nlqgAgpbVr16q/v19nz55VpVJRLpfTvffeq7Vr12Zd2qphZSeAqO3atUvHjx/Xzp07denSJe3cuVPHjx/Xrl27si5t1RDkAKJ24cIF9fT0aHh4WBs3btTw8LB6enp04cKFrEtbNZlcNAsAQqlUKnrmmWfU0tIyu29mZkZtbW0ZVrW66JEDiFoul9Po6OicfaOjo8rlchlVtPoIcgBRKxaL6uvrU7lc1szMjMrlsvr6+lQsFrMubdUwtAIgar291yfUFQqF2VkrAwMDs/sbQeJ55CExjxwAlm+xeeQMrQBA5AhyAIgcQQ4AkSPIASByBDkARI4gB4DIEeQAEDmCHAAix8pOANEws9Rt1NvNJiSCHEBElgphM6vLoF4KQysAELkgQW5mf2Fm3zezZ82sZGaNcyFgAMhY6iA3s7sl/bmkvLt3SWqW9MG07QIAkgk1tLJG0jozWyNpvaSfBGoXALCE1EHu7hck/Z2kFyS9JOlld//q/OeZWb+ZjZnZ2MTERNrDAgCqQgytbJL0AUlvlfQrkm4zsw/Nf567H3H3vLvnOzo60h4WAFAVYmjlAUnn3H3C3WckfVHS7wRoFwCQQIggf0HSb5nZers+W/89kioB2gUAJBBijPxbkj4v6WlJ36u2eSRtuwCAZIKs7HT3xyQ9FqItAMDysLITACJHkANA5AhyAIgcQQ4AkSPIASByBDkARI4gB4DIEeQAEDmCHAAiR5ADQOQIcgCIHEEOAJEjyAEgcgQ5AESOIAeAyBHkABA5ghwAIkeQA0DkggS5mW00s8+b2Q/MrGJmvx2iXQDA0oLcs1PSpyT9h7v/sZm1SlofqF0AwBJSB7mZ3S5pp6Q/lSR3n5Y0nbZdAEAyIYZWflXShKTPmNkzZvZpM7tt/pPMrN/MxsxsbGJiIsBhAQBSmCBfI+k3JP2Tu79D0qSkj81/krsfcfe8u+c7OjoCHBYAIIUJ8nFJ4+7+rerjz+t6sAMAVkHqIHf3n0p60czeXt31HknPpW0XAJBMqFkrBUnHqzNWfiTpzwK1CwBYQpAgd/ezkvIh2gIALA8rOwEgcgQ5AESOIAeAyBHkABA5ghwAIkeQA0DkCHIAiBxBDgCRI8gBIHIEOQBEjiAHgMgR5AAQOYIcACJHkANA5AhyAIgcQQ4AkSPIAdSE9vZ2mVmqTVLqNtrb2zN+JZYv1K3eACCVy5cvy92zLmP2P4SYBOuRm1mzmT1jZidDtQkAWFrIoZX9kioB2wMAJBAkyM1si6Q/lPTpEO0BAJIL1SP/B0kHJF1b7Alm1m9mY2Y2NjExEeiwAIDUQW5m75N00d3P3Ox57n7E3fPunu/o6Eh7WABAVYge+bskvd/Mfizpc5LebWb/EqBdAEACqYPc3T/u7lvcfbukD0o67e4fSl0ZACARFgQBQOSCLghy969L+nrINgEAN0ePHAAiR5ADQOQIcgCIHEEOAJEjyAEgcgQ5AESOIAeAyBHkABA57hAEoCb4Y7dLj7856zKu1xEZgjygUqmkgYEBVSoV5XI5FYtF9fb2Zl0WEAX7xCs1c6s3fzzrKpaHIA+kVCqpWCxqaGhI999/v0ZHR9XX1ydJhDmAFcUYeSADAwMaGhpSd3e3Wlpa1N3draGhIQ0MDGRdGoA6R5AHUqlUdPDgQTU1NcnM1NTUpIMHD6pS4TamAFYWQR7IunXrdOrUKW3cuFFNTU3auHGjTp06pXXr1mVdGoA6R5AHMjk5KTNTa2ur3F2tra0yM01OTmZdGoA6R5AHtGHDhtke+Lp167Rhw4aMKwLQCAjygHbt2qVz587p2rVrOnfunHbt2pV1SQAaAEEe0MmTJ/XII4/o5Zdf1iOPPKKTJ09mXRKABmBpJ+Cb2VZJxyS9RdI1SUfc/VM3+518Pu9jY2Opjltrtm7dqosXL2p6enp2X2trq+688069+OKLGVYGxMHMamdBUA3UsRAzO+Pu+fn7QywIuiLpL939aTN7k6QzZvaUuz8XoO2aY2aJnzs9Pa3x8fEFf6dW/6EAiE/qoRV3f8ndn65+/wtJFUl3p223Vrn7otvIyIg6OzslSZ2dnRoZGVn0uQAQSuqhlTmNmW2X9A1JXe7+yryf9Uvql6Rt27bdd/78+WDHrTW1/KcZUKtq5XNTK3UsZLGhlWAnO81sg6QvSPro/BCXJHc/4u55d893dHSEOiwANLwgQW5mLboe4sfd/Ysh2gQAJJM6yO36mbwhSRV3/2T6kgAAyxGiR/4uSR+W9G4zO1vd3hugXQBAAqmnH7r7qKTkc/IAYBHLmd67UjZt2pR1CcvGyk5A128M0tXVpebmZnV1dalUKmVdUsO52dTepFuIdi5dupTxK7F83CEIDY+7OyF29MjR8Li7E2IXdEFQUvV4rZXXq+UFBXij5uZmTU1NqaWlZXbfzMyM2tradPXq1Qwrw3LV+2dvxRcEAbHK5XIaHR2ds290dFS5XC6jioDlIcjR8IrFovr6+lQulzUzM6Nyuay+vj4Vi8WsSwMS4WQnGt6NE5qFQkGVSkW5XE4DAwOc6EQ0GCNfAfU+TgfUqnr/7DFGDtwE88gRM4IcDa9UKmn//v2anJyUu2tyclL79+8nzBENghwN78CBA2pubtbw8LBee+01DQ8Pq7m5WQcOHMi6NCARghwNb3x8XMeOHZuzIOjYsWMaHx/PujQgEYIcACJHkM/T3t4uM0u1SUr1++3t7Rm/Co1ly5Yt2rNnz5x55Hv27NGWLVuyLg1IhCCf5/Lly0GuwpZmu3z5ctYvQ0M5dOiQrly5or1796qtrU179+7VlStXdOjQoaxLAxJhHvk8tTAPtRZqqFehrnfN+1Ob6v2zs9g8clZ2oqEs9SGv9yBAfWJoBQAiFyTIzexBM/sfM/uhmX0sRJsAgGRSB7mZNUv6R0l/IGmHpF4z25G2XQBAMiF65O+U9EN3/5G7T0v6nKQPBGgXAJBAiJOdd0t68XWPxyX9ZoB2M+GP3S49/ubsawCAhEIE+ULzud5w2t/M+iX1S9K2bdsCHHZl2CdeyXzWgpnJH8+0BAARCTG0Mi5p6+seb5H0k/lPcvcj7p5393xHR0eAwwIApDBB/m1JbzOzt5pZq6QPSvpygHYBAAmkHlpx9ytmtk/Sf0pqljTs7t9PXRkAIJEgKzvd/SuSvhKiLQDA8rCyEwAix7VWFhDqwkq3atOmTZkeH0BcCPJ5Qkw95MJLAFYTQysAEDmCHAAiR5ADQOQIcgCIHEEOAJEjyAEgcgQ5AESOIAeAyBHkABA5ghwAIkeQo660t7fLzG55k5Tq981M7e3tGb8KaDRcawV15fLly5lf5ybri66h8dAjB4DIEeQAEDmCHAAiR5ADQORSBbmZ/a2Z/cDMvmtm/25mG0MVBgBIJm2P/ClJXe7+65Kel/Tx9CUBAJYjVZC7+1fd/Ur14TclbUlfUrxKpZK6urokSV1dXSqVShlXBKARhJxHvlfSvy72QzPrl9QvSdu2bQt42NpQKpX08MMPa2pqSpL0/PPP6+GHH5Yk9fb2ZllaQ/HHbpcef3P2NQCryJZaPGFmpyS9ZYEfFd39S9XnFCXlJf2RJ1iNkc/nfWxs7BbKzV6oxR5ZL1qpV7Vw4+taqKFR1ftrb2Zn3D0/f/+SQyvu/oC7dy2w3QjxPZLeJ2l3khCPnbsvuEnSoUOH5uw7dOjQor8DYPlCXGKhHqWdtfKgpL+W9H53fzVMSfG6MT6+2GMA6SzWkRoZGVFnZ6eamprU2dmpkZGRm3a66k3aWSuHJb1J0lNmdtbMngxQU5TWrFmj3bt3q1wua2ZmRuVyWbt379aaNVzOBlhJpVJJxWJRg4ODmpqa0uDgoIrFYmNNNljsf62V3O677z6vN/v27fOmpibfvHmzS/LNmzd7U1OT79u3L+vSGsr1f9LU0Eg6Ozv99OnTc/adPn3aOzs7M6po5Uga8wUydcmTnSsh5pOdN1MoFHT06FG99tprWrt2rT7ykY9ocHAw67IaSi2c7KqFGhpJc3Ozpqam1NLSMrtvZmZGbW1tunr1aoaVhXfLJzuR3I0/7dx99k88ACsrl8tpdHR0zr7R0VHlcrmMKlp9BDmAqBWLRfX19c05P9XX16disZh1aauGM3GoO1lPMdu0aVOmx280NxbcFQoFVSoV5XI5DQwMNNRCPHrkAd1Yot/c3MwS/YwsdCJoOVuINi5dupTxq4BGQ488kBtToIaGhnT//fdrdHRUfX19kliiD6wkPnsJluivhHqctdLV1aXBwUF1d3fP7iuXyyoUCnr22WczrAzLwYyT+HR1damnp0cnTpyYHVq58bjePnuLzVohyANppClQ9Ywgj09TU5PuueceDQ8Pz/bI9+7dq/Pnz+vatWtZlxcU0w9XGFOggGy0traqUCiou7tbLS0t6u7uVqFQUGtra9alrRqCPBCmQAHZmJ6e1uHDh+d89g4fPqzp6emsS1s1nOwMhClQQDZ27Nihnp6eOZ+9hx56SCdOnMi6tFVDkAfU29tLcAOrrFgsLjhrZWBgIOvSVg1BDiBq/DXMrBVgDmatoJYxawUA6hRBDgCRI8gBIHIEOQBELkiQm9lfmZmb2R0h2gMAJJc6yM1sq6Tfk/RC+nIAAMsVokf+95IOSGLOFgBkIFWQm9n7JV1w9+8keG6/mY2Z2djExESawwIAXmfJlZ1mdkrSWxb4UVHS30j6/SQHcvcjko5I1xcELaNGAMBNLBnk7v7AQvvN7NckvVXSd6r3SNwi6Wkze6e7/zRolQCARd3ytVbc/XuS7rzx2Mx+LCnv7j8PUBcAICHmkQNA5IJd/dDdt4dqCwCQHD1yAIgcQQ5IKpVK6urqknT9ruylUinjioDkuB45Gkp1hlVqXLMcWeB65ICuB/D8rbOzU8ViUZ2dnWpqaprzeKHnE+KoNdzqDQ3vueee0+TkpIaHh2fv+bh3716dP38+69KAROiRo+G1traqUCiou7tbLS0t6u7uVqFQUGtra9alAYkQ5Gh409PTOnz4sMrlsmZmZlQul3X48GFNT09nXRqQCEMraHg7duxQT0/PnLuwP/TQQzpx4kTWpQGJ0CNHwysWixoZGdHg4KCmpqY0ODiokZERFYvFrEsDEqFHjobX29srSXN65AMDA7P7gVrHPHIAiATzyAGgThHkABA5ghwAIkeQA0DkCHIAiFwms1bMbEJSPV/I4g5J3PIuTrx3cav39+8ed++YvzOTIK93Zja20BQh1D7eu7g16vvH0AoARI4gB4DIEeQr40jWBeCW8d7FrSHfP8bIASBy9MgBIHIEOQBEjiAPyMyGzeyimT2bdS1YHjPbamZlM6uY2ffNbH/WNSEZM2szs/82s+9U37tPZF3TamOMPCAz2ynpl5KOuXtX1vUgOTO7S9Jd7v60mb1J0hlJPe7+XMalYQlmZpJuc/dfmlmLpFFJ+939mxmXtmrokQfk7t+QdCnrOrB87v6Suz9d/f4XkiqS7s62KiTh1/2y+rClujVUD5UgB+Yxs+2S3iHpW9lWgqTMrNnMzkq6KOkpd2+o944gB17HzDZI+oKkj7r7K1nXg2Tc/aq73ytpi6R3mllDDW0S5EBVdXz1C5KOu/sXs64Hy+fu/yfp65IezLiUVUWQA5o9YTYkqeLun8y6HiRnZh1mtrH6/TpJD0j6QbZVrS6CPCAzK0n6L0lvN7NxM+vLuiYk9i5JH5b0bjM7W93em3VRSOQuSWUz+66kb+v6GPnJjGtaVUw/BIDI0SMHgMgR5AAQOYIcACJHkANA5AhyAIgcQQ4AkSPIASBy/w/mnClya/c/UwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1)\n", "box = ax.boxplot([x1, x2, x3])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# 2d Histograms\n", "\n", "Histograms are also an excellent way to visualize relationships between variables. " ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "mean = [0, 0]\n", "cov = [[1, 1], [1, 2]]\n", "x, y = np.random.multivariate_normal(mean, cov, 10000).T" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3df2wc55kf8O+zq5G8VBIvXTOItZYsn+9i9xRFYs3EaoWitS6JfZGj8OIkOsc5BC1Qo4deEakOr3QjxLLhg4USObnAFSh81wAtrHPoHzkmrq6QE0goUPXkC2WSpyiR0ji25KxThDlpfTlxbS3Jp38sZ7U7nHd+7Mz+mJnvBzBgLpczw9Xy2Xee93mfV1QVRESUXLleXwAREUXDQE5ElHAM5ERECcdATkSUcAzkREQJt6YXJ73xxht18+bNvTg1EVFinT59+peqOuR8vCeBfPPmzZienu7FqYmIEktELrg9ztQKEVHCMZATESUcAzkRUcIxkBMRJRwDORFRwvWkaoWIKGumZsqYOHYeb1Wq2FAsYOye2zE6XIrl2AzkRJQqnQyYUa7pkW+dQbW2BAAoV6p45FtnACCWa2NqhYhSww6Y5UoVimsBc2qm3NPrmjh2vhHEbdXaEiaOnY/l+ByRE1FqeAVM08i3GyP4tyrVUI+HxRE5EaVG2IDZrRH8hmIh1ONhMZATUWqEDZidTnnYxu65HQUr3/JYwcpj7J7bYzk+AzkRpUbYgNnplIdtdLiEJz+zFaViAQKgVCzgyc9sZdUKUZb1Y2VGP7BfA/u1KQ5YUAX2T85i4tj5Va/ThmIBZZegHVfKw3ltnfo34oicKGH6tTKjX4wOl3ByfBcO792Od2rLqFRrxtep0ymPbmEgJ0qYbuV1ky7I69TplEe3MLVClDDdyusmXdDXqZMpj27hiJwoYTpdypYWWXqdGMiJEiYted1Oy9LrxNQKUcI4KzM6VbXSjcqYTp6jW69TPxBV7fpJR0ZGlHt2EvUvZ5MnoD6ajXMisBvnSBsROa2qI87HOSInyjDTiLidniVh+Z0j6mg9S7X2DOREKeYVzLxaq3pVfMQVIP3OEaXta6fbxvYbplaIEs4UWP1SFzsPHXdd1Vhaqepw+16xYOHdxeVY0iHtnL9ULODk+K5Ixw7y8/3KlFph1QpRgnmt8vRbEOM1IjZVfIjAeMypmTJ2HjqOW8ePYueh474rTb2qSqLWymet1p6BnCjBvIK1XzDzqrM2rXisLNRcf8b+AAnTNsBrVWXUGvAs1ZADzJETJZpXsPZrCDV2z+2uqRe7ztptxePEsfOuxwTMI3WvlEvzOey7iP2TsygOWLBygtrytdRvmBpwv98tbTgiJ0owr5Gn34KYdvqM3H3HUKjra/6g8Uq9OFNElxdqgNRz8u30QElLD5WgOCInSgi3SU2vkWeQBTFh+4ycODcf6prtDxq/KhK3FFFtSbF+3RrMPvqJUOe0paGHSlCsWiFKAK8KFCD66sWgJYW3jh9F0IhRsPK4/84STpybN6ZjigUL69etMX5fALx+aHfAM6ZfxxcEiUgewDSAsqreF9dxicg8qfnwc3NYVsWGYgGH925vK4Af/M5ZVKrXJjG9aq5NeXdbToBraW3F5PffRG3JHPor1VrLuZ3SOjkZtzhTK18G8CMA74vxmEQE86Tm0soddblSxf7JWUxfuIQnRrd6HssefXsFZNNEpVsqx2blBVBgeeWaqrVlz+vw4zU56XYHAWSjr4qbWCY7ReRmALsB/FkcxyOiVkFGpgrgmVMXPUv+micV/ZRXVlg2a55EBIC8CID6ZOL6tWtaqkyi8JqcdKudH3t+DmMvzGV216RYcuQi8gKAJwG8F8BX3FIrIvIQgIcAYNOmTXdeuHAh8nmJOqXf+nQ8+Kd/hZOvXQr03MEBCzNfc58gNK149LLzthtw5F/9Y9/nhcmfl4oFLFxdrFenuHzPa/VlmN8h6Ss5nTq2slNE7gPwC1U97fU8VX1aVUdUdWRoKFwJE1E39duemAemzgQO4gBcg6OtnZWNJ1+7hANTZ3yfFzSf/cUdmwDUr1Mc3wtS6x3md0jrSk6nOFIrOwHsEZE3AHwTwC4ReSaG4xL1RL/tifnsK2+G/hnTh067k4f2NUzNlLH9sZexefwoNo8fxfDjLzfOdfcdQ6sCs5vJv36zMaJWoPEzQWu9w/wOWZksjRzIVfURVb1ZVTcD+F0Ax1X1i5GvjKhHOtGnI2wfkmZLbaQ/H3vprOvjbouEgl7D1EwZY8/PtVSZXF6oYd/kLDaPH8Uzpy4GSq048+iKaymQIOkrt9/Bykl9srVJmldyOnFBEJGD39L2sKK2VM2LuAbz1lK/VpcXajgwdQYnzs275vn3Tc6G/j3a+ZmgwnxImhY6uT2WlaoVLggicgizc02QSdGwLVWdx9z8DwquOfIv7tiEo3/zc8+ceLN8TrA2L55lgVYOiFg12JZeTUr226S2H+4QRBRQ0L0eg460w6Rq3I556cpV7LztBpz66WUsqSIvggfu2tioF3/m1MVAv9fSsqLqUR4oACY+tx3TFy7h2VfebCul0w4rJ6HqxTu11VySN59gICdyEaRPR9Dt0IKkarwW6VRrS3jjb6t47clPtjx/+2Mve66KDEtx7fd+YnQrNo8fje3YXmrL2sjpN79unQ603djOrlsYyInaFHSk7ddS1S2V41SuVBuBdcDKobaksS2+sQ0OWC1fm3LznXB5obYqSHc60KZp8wkGcqI2BZ0U9UvVuAUsLwsdSmJfXqjhwT/9K/zw578KnHePU7W2hMdeOtt4XTodaOOe1O4l9iMnapNXv29nuSEAnBzfhcN7twMA9k/OYvjxl7H9sZdDr7TspJOvXYoliEuQgnIXlxdqjdLMtWvcw5MCoUs43fj1a08SVq1QasU5Uea1wbFbGZxbqqRg5bC4rJ7dANOgWLBwcM8WjL0w19bvWioWcPcdQ76TuO1u+twsLVUrDOSUSmFKCL2OYU9ACtCy2MXrWO30M0kL+3WZvnApcDWNkwDIBczPBy3h7PcAHVTHeq0Q9aOoy+ydXQKdIcV0rKmZcteCeJvZi46yN7o40mYQB+o56qCTrF4lnP3SK6cbGMgplaJOlAWZgHQeyw4gcQgSpPstQSMCHPzOWeybnG372uwcdT5gkt1tYtJrE460BnNWrVAqRa1ICBLwry9Y2HnoeOP2feHqYqjqEy/9FqSDUEXkuvZ1KxOcD9y1MVCO3G1i0msTDr869KSmZDgip1RqpyKhudIk5zMitHKCK1cXW27fe1GylzaVar2efOSWG1CwzOEpL2Kco/D6sPZKryU5JcNATqnUvJONwL9FqvOP2C9H+57r1qS++qRXqrUl7Juc9ewJ8/XPbzP+W/p1eDSN2PutfXEYTK1QagVZZm8LsyinVCwkcvVfWvhlz+1/84efm3P9QDaN2JO80pOBnBIp7lxm0D9WKydYuLqYyBx2Wijgm+seHS4ZSyDvvsN9h7Ikr/RkaoUSpxO5TNMfa7FgNdIzxYIFiPdWatQdQVIeJ87Nh3o8ySs9GcgpcTqRyzT9ER/cswUnx3fh9UO7sX4d8+L9xK9eP2yqJOy8Sj9haoUSpxO5zCA9yJOQK80Sv1rzdlIlYeZV+gkDOSVOmD/QMLl00x/x1EwZB79zlnnxPuNXWeTXPjhNGMgpcYL+gfptTDA1U8ZjL51t5LztZk/AtZH59QUr1s0bKD55EUzNlD0/mIFs7OPJplmUSAemzjS2I3NufWbz2itz7J7bXbvz5QDk88JceEIIgAd3bMITo1s7vi1cP3wgcM9OSo2pmTJePF1u3FovqeLF02WM3HIDgGsjMFMofqtSxcSx867BehnAMoN4V0XZiUhxrUHXi6fLHdkWbmqm3PKhX65UMfbCXCzHjgtH5JQ4ppH24ICFd2rLvgt7SoYcO3WXlRdMfHYbAPf+7WGYPgxMbW7DGH78ZdeS08EBCzNf+0SkY4fFNraUaM19UExB+PJCLVAwuHzl3bgvj0LKCVBb0kbJ6P13liK15TWN6OOoNDKtG+in9QRMrVDfC7I5cRid2vOSgrP3jS5Xqtg3OWt8Xl4Ey6rYsLJr0JFTF0NVDyVhVWYcGMip7wXpg2LlBevXrmGFScosq+L1Q7tbHgsazOMqNSwaKpeKBSvysePC1Ar1vSC3x+vXrsF9227qwtVQNzlH1E+MbsXhvduNi4FEEPuqzIN7tsDKtZ7PykmjVLUfRB6Ri8hGAP8dwAdQn/R/WlX/U9TjUvaYSrxMC4CavV2tGXtoUHK5jahHh0vYb0jHqAJP7d0eazVJEurRI1etiMhNAG5S1VdF5L0ATgMYVdUfmn6GVSvk5LVZMgCMPT+H2rL5vcpKlPTxqgrx2uA6jkqVftWxqhVV/bmqvrry/78C8CMA/fNRRV3TXFmy89DxUN0IvfZZ3D85i0WPIB52n0dKhssLNeP7yCv3ncWeOLHmyEVkM4BhAK/EeVzqf1Fby3rts6jw3sPy/jvrPVLaXVRC/cv0PhodLhknG7NSqdIstkAuIu8B8CKAfar6dy7ff0hEpkVken6eucy0idpaNsof35FTF3Fg6kxfVRFQfEzvo4N7tiS2f3jcYgnkImKhHsSPqOq33J6jqk+r6oiqjgwNue/QQckVtbWs3z6LXuxl2lcX46kzp/7j9j5Kcv/wuMVRtSIA/iuAH6nqH0e/JEqiqNtkOSsDwiZJFFzok2bO95GzwulwzJUqSRPHiHwngN8DsEtEZlf++2QMx6UESfI2WdRZYaag16/Nr6rZdr6POrHVX9JFHpGr6v9GuH8rSqF2am3tDRu4GjPdgt5dCYCzj9/r2zLWaz4mq6NyLtGn2ITZJmtqpuxbG07ZYqdP/N5HndjqL+m4RJ96YuLYeQZxahB414Y3M827ZLHs0MYROXWcfatcrlQjbSJA6WTv8hP0bs5tqz8AWLi66Ln1W5oxkFNgYbe7csuBM4iTk1fFidd7zvneurxQi21XoKRhIKdA/DYydj63eVNjIpNSseAZxL3ecxPHzq+aKM/qpCdz5BRI0JWb9h9fO0F83Zocy59SKufyDysA7r6jvjjQrU+P33vONLlZrlTb6veTZByRUyBBKwUee+ls6J18rLzgo5sHcfK1S21fH/U3t3ltRX3DZMB942TT+8h+z3m1N26uLwfSn2rhiJwCCVIpMDVTbmskXltSBvGMqtaW8Owrb7qOvE3s91yQtg5h+v0kGQN5xgVtPTt2z+2w8o5dUvLSUjKWhT8Yil+YCfDmVZ7OXismWagvZ2olw8JMYAJYvUTP8XUW/mCo+5o3YHZWSjUvHjJtNpGF+nIG8gwLs9TZbQFPbVkbo/CJY+dDN7oiCsJtA2Y3bvXlWen3w0CeYWGWOntVCOyfnGUQp46xR9R+6xiSsLdmpzCQZ1iY1rN+FQJEJuvX5vFHv1Pfe7U5yF55d9G3YZo9og6aBgzT7ydNGMgzLMytqGlZNJGT3Yah5DFqBtw33HYex94oYueh4+x46IGBPMPC3Io2P5e71ZObfE6QAxpzKeVKFWMvzOHgd87i7Wqt5f1lWvDTbFm18b5jx0NvDOQZF+ZW1H7ubY/8JXum0CpLywpnWK4taSN9YqdDpi9calkAZNKc4ou6A1XasY6cjEw15gzi1C7TAiAnZ4qPO1B544g8Y4J2MHRu/FCuVDH2/BwAYHDAYkMsaluQgYBzE+UsV6QEIdqD0dXIyIhOT093/bxZ5za5VLDyuP/OEk6cm2/5AzFtwVYsWADA7dmobX496UvFAk6O7+riFSWHiJxW1RHn40ytZIhpAdAzpy6u2sjWFKgr1RqDOLkqWHkULO+QUrDy2PFrg57fZ7okPKZWUsQvbRJ0hp8lhhTW4ICFRz+1BQBW3fUJ6msN7HJEU0+e5nJDCoeplQRrDtzFAQt//85iyzL6gpVv+cMw9aIgioMdqAHvXPat40eNi8gEYP7bA1MrKWPnu+2UyOWF2qpeKM4WnkHafhK1q3m15cnxXXj90G6cHN+1KiB7lQw2p/eysilEHBjIE8pvMYWtOZ3S3PbTC3fpIaC+wCcsr/7fdjlruVL1fY9lpY94XBjIEypovjsn0jKyGR0u4eT4Ljy1d7vr6HxwwMKDOzbFdp2UTIL6Ap92uL03m+8ggfrI2y+Yc9VmcAzkCRV0RduSquttqrMpf6lYwBd3bMLA2jV45tTFDlwxJUmUmTO396bbHaSiXs5qukPkqs3gWLWSUG5NrKy8YHFJV/0RNt+mOieh7HpdvwZGREGYygdNo+tKtYb7tt20ask+yxDDiWVELiL3ish5EfmJiIzHcUzy5jainvjsNuNIyp5ActaL2yP1djZNJmpWLFjG8kGv0fWJc/Or3sssQwwncvmhiOQB/BjAxwH8DMD3ATygqj80/QzLDzunnYZWXHJPcfBakTk1U8a+yVnX7wkQaAcg6mz54UcB/ERVf6qqVwF8E8CnYzgutaGdhlYM4mQTeE9C2i0a3HhNTo4OlzA44P6zzIVHF0cgLwF4s+nrn6081kJEHhKRaRGZnp+fj+G02WbqTOhXWkjkRWGe6LRygoN7trQdkB/91BZ2MOyQOCY73T7AV70XVPVpAE8D9dRKDOfNLLdtr8aen8NjL53F5YVaY0k0UVyKBQsH92xp5K3b2eSYHQw7J45A/jMAG5u+vhnAWzEclwzcSrlqy9pIkTCIU9zWr7sWKqIE5KzuqdlpcQTy7wP4DRG5FUAZwO8C+EIMx800rwZYXChB3ebc7NgZzO3yVgbp3ogcyFV1UUT+AMAxAHkA31DVs5GvLMP8dgz32tE+CCsnq/qyEPmp1pbw8HNzja+D7GpP3RFLHbmq/qWqflBVb1PVP4rjmFlm6htuj3qiNL8asHJspkJts1cKu607YH+U3uES/T7kt2O4czFQsWDByntH57wIntq7HYPr16G2xNE4ta9aWzKWrDLt1xtcot+HguwY7pw0mpop4+Hn5ox15MuqGB0uGRdlEMWBNeG9wRF5H2pnx/DR4RK+/vltxqzJhmKB/Z0pkCDta4sFa9V71MoJFq4urlrbQJ3HQN6H3PqoBOk9MTpcwoM7Nq0K5vaHAPOXFMR7161pvPcGByxYjsBesPI4uGfLqvQepL5KmJtDdB+3ekshU+mi1xZbRDZn7xO/vWAB8zaCXv1XKDxTrxXmyBMkyB8U4L7oYmqmjJxIW71YKFucee4gi3j8JuipsxjIE8JrWX5loYYNxQLuvmMIJ87N461KFdcXLIgAlYVaY2NmBvFsW782j4WrS77rENrpfRJkgp46hznyhDj4ndV1u/ayfDsn+cypi41+45VqrfE9t42ZKXuKA2sbGyKbmqsVC1ZbC3ramaCn+DCQJ8CBqTOoVNlqlqJpTnOYAu/BPVvaOna7E/QUD6ZW+tiBqTP481cugoNp8iMAHtyxCSfOzRvTJs51CEC8nQjZEKt3GMi7IOgkZbMDU2e4CTIF9k9uuwFPjG4F4L7/qluag4E3PRjIO8yvAZbzuXbA5yCcwnjjb6+Nwp2j7eKABVVg/+QsJo6dZw/wFGKOvMP8GmDZ7IBfZhCnNjjL/EaHSzg5vguH927HO7VlVKpcqJNmHJF3mKmOtlypYueh4410y5V3F7mLPbXNVObnNZDgqDw9OCLvMK86Wnv0Xa5UWZVCbfMq8+NCnWxgIO+wKL3Dm+VFsNanVS1l0/13mictTQMJLtRJFwbyDnOrrw2qYOXx1N7teOPQbrz25CfZR5xcnTg3b/weF+pkA5tm9cDm8aPG75WKhVVlinY1S5Tt3Sjd3N43tnbKX6k/sWlWHylYOVRry8bvH967HUB9WT43gqAg7A95uypl+sKlRt8dBu/0YyDvsqmZMhY9lmqWK1X8u8lZmMM8kbdqbQlHTl1slLG6rV3gKD1dmCPvsolj531z3QziFJXzHda8dsG5ZoG15cnHQN5lLPuiXrHfe0EXqVFyMLUSUdhbVL9e0ERh5KS+x2bzXZ5g9YgcuFZyyNry9OGIPIJ2blHjqisnAoD3XWdh4rPbWspbH9yxybPkkLXl6cMReQTtLH8eHS5h+sKllsko4NooqliwcOXqImvGM8zKCZYBLAXoX/x2tebaxXDklhuMd4pj99weqDsiJQcDeQTt3qKeODe/6tZXcW2j2qmZMh5+bo5bs2VEsWBh/bo1LUEXaO0VfuXdRdc2DqZRtFeL2k70IqfeihTIRWQCwKcAXAXwGoB/oaqVOC4sCUz57pwIbh0/avwD8fsAGB0uYT/rx1Oh5BGEgWu78pg20bYF7TEeFHuRp0vUHPl3AXxIVT8M4McAHol+Sclhyncvqa7KmU/NlLHz0HHcOn4UOXHvmbKhWGg8j2Px9Lhv202u75PBASvwdmjcSo28xLZEX0R+B8BnVfVBv+emaYl+c9VKTsQ1HZIT+G7XVrDy+Eebrsf/ee0Sg3jKFKw87r+zFNtKSy7myS7TEv04A/lLACZV9Rm/56YpkDe7dfxoqCDcXCa2fm0eV66yH3la2fMfQXgFalOKhaPzbDAFct/Uioh8T0R+4PLfp5ue81UAiwCOeBznIRGZFpHp+Xlzt7YkC1u+1Rz0GcTTLWiNtl9JKxfzkBvfQK6qH1PVD7n8920AEJEvAbgPwIPqMbxX1adVdURVR4aGhuL7DfoIa8TJJOiHvF+g5mIechNpslNE7gXw7wHsUdWFeC4puZonpIhsYapL/AI1F/OQm6h15H8CYB2A70q9EuOUqv7ryFfVQ0EmktyeA7TW5Q4OWLi8wO3bsigH4PoBC5WFWujJSFNJqx2ouZiH3EQK5Kr663FdSD9wTiSZ2n86nzP2/BwgaKzGZC+V5BIBos7/Xz9g4dFPra4Nbx4AFAcsqNZXZjYHe79AzcU85IYrO5sEWXLv9pxagKXUlAy/PrQe//cXVyId4/JCzXcA0Hy35jZg8ArUXMxDTgzkTYJMJHFSKb0EwE8iBnFbkAGA6fkM1BRWZgO5W577+oLlupRaAew8dBxj99zONrQpFua+av3aPBauLnn+TNgBAAcJ1K5MtrF1q9Ude34Ov3p30fgz9u3v3XcMrSoxtHICK+++7N4pb1ieT8ly5eoSrrPyeGrvdmOVUk6kUf8dpKqElSfUrkwGclOe269taLW2hBPn5ltKDPMiqC0r1uSCBehlVTy1dzsYzpPPTod49dyxF/P4rTFg5QlFkZhA3tx0aueh45H2F4xyC/tWpdqoLihY+UZvlWot2E6bG4oFjA6X2E8lJez3w5Of2ep6t9Wc+25uejU4YKFYsNgAi2KRiBx5kLLAMKLkue3bX7/JKzfNo64Sc+2pYL8fvFoPN7cnZrCmTkjEiDzu/hLtLqVvDsRBR/WmUReX8yefMx3CVZfUK4kYkcfdX8KtVnfh6qLrSsy8CJZVV9XzBhnV50Uw++gnWh5rrpa5vmDhOiuHykIN13OLt57JAcjnJdRrP+iy6IerLqlXEhHI/ZYtt8N5mxu2PajbH63TkuqqwN0crCvVGgpWHof3bsfEsfPGXWQoPoMDFoBrC3KKBQsH92wB0PrBfvcdQzhxbh7lSrVltadbALdx1SX1Smz9yMMI24+8Wz2Y/fqsOL/f/MfuZnDAwju1Zd9ceqlYwFsrpZDUOVZeMPHZbQyslFgd31gijHY2luj1riheHyYAXL+3bk0u0ChbEG0ClvzlBPjCXZvwxOjWxmO9fk8RhWUK5IlIrQC9n/F/7KWzxglXe+cXZ1AIuoGy/Xy/VI2taFiBmibt/o6WIde9rMCLp8sYueUGjA6XYq+EIuqlxATyXpqaKRtb0nqVlk0cO+87yrYnw5rzq34/s37dGry7uBS4dr0f5A37mbqxU00mpiCfl3oQN52rWlvCvslZTBw7jyvvLvo2SCNKikSUH3ZK0EVGXmWOXhOubiWGVk4wOOBekjg6XMLJ8V1449BuFCzzP025UsViwjouBg3iAjR62phUqjXXlbH2OfzOVa5UjaN99juhJMrsiDzMrbXXH7dXaZlzlG0v5x9Yu8ZY+WDzG23XlhRWDkjQoDwQBbBvctbzg8x+ni3MaN8Pa74piTI7Ig+zyMj0x10sWL634W7L+cuVKvZPzuLA1Jk2r76utgz8xvvXRzpGvwqaNioWLCzHFMRZ801JldlAHmaRkVuKpGDlG/XHftw+NBTAkVMXXdM5UzNlBG2SGFf/7KSqrOyw46ZYsBp1424GB6xG7xP2O6Eky2xqxVTu19x7vDl3DbS/0MP0oaGoV8M4a9NfPF0OvN1YsjLlnWFaUXnftpvw4mn3eY+ClfdNbxElRWYDuVe5n1u+PEr5o2nDCqC+wtCuiClXqjhy6mLignNO6uV9vTA4YBk/aE2NzfIiHH1TqmQ2kPuV+4UtRTswdQbPvvImlrRe/vbAXRvxxOhWTM2UceWqecMKJ694KD7f75VuBHErJ1gGWnrGW3nBo5+qp7fcPmhNdfzLqgzilCqZDeTAtT/+W8ePugbIoKVoB6bO4JlTFxtfL6k2vj5xbj6WRlhuQbyXI+E4Bfk9Jj63DdMXLrV8WO79yEbPgNyJHj1E/Sizk53NorYfffaVN42Px1GXbBqJtxPE+22nOSsvWLfG/204feESXjxdbqkVf/F02XODEdMkNStTKG0YyBH9D95Uw7y00v62Hc3VFHEOuqNU6sX9GTA4YAEarNTw2VfeDN2T3rkrDytTKK0ynVqxRa1KMS1IyYs0eq6EiZ+lYqHRvwWoV9F0o6GW/XuY7gDi/ECx9zw1tT5wMn1Y+t3x9LpHD1E3cES+wl4e//qh3Tg5vivUH/8Dd200Pj46XMKDOzatGs1aOcH6tat3CHK7Ewi7m1A7I+eClcfXP78NbxzajcN7t7dxhHDG7rk9VNrJbT9MgPluIoCBPBZPjG7FzttuaHls5203NFqmPjG6FYf3bm+5xZ/43DacffxePOV4/P47S5g4dr6l/0tzigBwD9R2nAubijH1fCkZAmRc6ZXR4VLgIFyw8njgro2h0l9xbtZN1O9iSa2IyFcATAAYUtVfxnHMJJmaKePVi2+3PPbqxbcbQRgw3+I3P+7V/6WZW6DecH09HbWTIhkAAAfUSURBVDM1U8bDz80F6j3iTOE0My2yuf/OEk6cm2+koNpJ+dirLd3OYad17DRPqSnNNXLLDYHSX2xRS1kTOZCLyEYAHwdw0e+5aeXVtyVM4DAd56t/cQbLCs9e5W9Vqo0AFiSIN49m3TZYAIDrrFzjnPaWaM7fZ/jxlwPnuW1//85iy4dc0LkJ+0PPvt79Ky1pnT8T178HUVLEMSI/DOAPAXw7hmMlUlybQ5uef+Wq/2YTG4oF40pGN3YqxW30OvbCHKBAram+sVKtYd/kLPZNzrYE9XaqYGrL2giqYScjg4y2496sm6jfRcqRi8geAGVVnQvw3IdEZFpEpufn56Octu9ErUNv9/nN7r5jKHCgKhULLaNhZ/CvLWlLEHeqVGsYe34OUzNlvN3mTkXtBtUgXSvj+vcgSgrfQC4i3xORH7j892kAXwXwtSAnUtWnVXVEVUeGhoaiXndfiWvhSZSFKifOzQcKVM7rajeg2qPqdoNjuz8XZLTNhUCUNb6BXFU/pqofcv4H4KcAbgUwJyJvALgZwKsi8oHOXnL/iWvhyehwCcWCue2ql7cq1dA7EgHRRqmmc/pVtjjz82GqS4KMtrkQiLKm7Ry5qp4B8H7765VgPpLFqhUgvoUnB/dsMVZyeNngSJcEXdjkVjli5WVVjjzMOe1WvG6/Q3MVSjvVJaZqGudomwuBKEu4srPPuAXGoBs42z8fJoCZgr/9mOncVk48zxmkVLCd6pKoq3CJ0kg0pm2ywhgZGdHp6emunzepvJbol3wCmVtpYdigNzVTxmMvnW2UGZpKEcMydZ0UAK8f2h3p2ERpJCKnVXXE+ThH5AlgSif45X3jWhjTqTQF28wSxYNL9BOg3cm7MBtMxynoBCarS4jiwRF5BHGkLYJyjortYOl17l4sjAlzF8B8N1E8GMjb1Mt+HkHP3YvURdgJTFaXEEWXqdRKnB3xepW2CHPuXqQuuDyeqPsyMyKPewTdy4AV9Ny9SF1wApOo+zIzIo97BN3Lfh5Bzm3ffdg7yR/euz30hhnt4AQmUfdlJpDHPYLuZcDyO7d991GuVKG4dvfRjc0VuDyeqPsyk1qJ+5a/lxUXfufudT9uTmASdVdmAnnQHh1h9DJgeZ2bE45E2ZKZQJ6WmuUgteuccCTKlswEciD5t/xBK286cfdBRP0rM5OdaRC08oYTjkTZkqkReS/FsZw/TO476XcfRBQcR+RdEFc5IPeiJCI3DORdENdiJC62ISI3TK10QVzlgGmpvCGieDGQd0Gc5YDMfRORE1MrXcCUCBF1EkfkXcCUCBF1EgN5lzAlQkSdwtQKEVHCMZATESUcUysd0s2NmYko2xjIO6CXGzMTUfYwtdIBvdyYmYiyJ3IgF5F/KyLnReSsiPzHOC4q6bixAxF1U6TUiojcDeDTAD6squ+KyPvjuaxk48YORNRNUUfkvw/gkKq+CwCq+ovol5R8XMlJRN0UNZB/EMA/FZFXROR/ichH4riopOPGDkTUTb6pFRH5HoAPuHzrqys/PwhgB4CPAHhORH5NVdXlOA8BeAgANm3aFOWaE4ErOYmoW3wDuap+zPQ9Efl9AN9aCdx/LSLLAG4EMO9ynKcBPA0AIyMjqwI9ERG1J2pqZQrALgAQkQ8CWAvgl1EvioiIgou6IOgbAL4hIj8AcBXAl9zSKkRE1DmRArmqXgXwxZiuhYiI2sCVnURECSe9yISIyDyACxEOcSOYiw+Kr1UwfJ2C4esUTKdep1tUdcj5YE8CeVQiMq2qI72+jiTgaxUMX6dg+DoF0+3XiakVIqKEYyAnIkq4pAbyp3t9AQnC1yoYvk7B8HUKpquvUyJz5EREdE1SR+RERLSCgZyIKOESH8hF5CsioiJyY6+vpR+JyISInBORvxGRvxCRYq+vqZ+IyL0rO1z9RETGe309/UpENorICRH50cpuYF/u9TX1MxHJi8iMiPyPbpwv0YFcRDYC+DiAi72+lj72XQAfUtUPA/gxgEd6fD19Q0TyAP4zgN8G8JsAHhCR3+ztVfWtRQAPq+o/RL1t9b/ha+XpywB+1K2TJTqQAzgM4A8BcMbWQFVfVtXFlS9PAbi5l9fTZz4K4Ceq+tOVvkHfRH3rQnJQ1Z+r6qsr//8r1IMUG+67EJGbAewG8GfdOmdiA7mI7AFQVtW5Xl9LgvxLAP+z1xfRR0oA3mz6+mdgcPIlIpsBDAN4pbdX0reeQn2AudytE0ZtY9tRPrsT/QcAn+juFfUnr9dJVb+98pyvon57fKSb19bnxOUx3t15EJH3AHgRwD5V/bteX0+/EZH7APxCVU+LyD/v1nn7OpCbdicSka0AbgUwJyJAPV3wqoh8VFX/XxcvsS947eIEACLyJQD3Afgt9otv8TMAG5u+vhnAWz26lr4nIhbqQfyIqn6r19fTp3YC2CMinwRwHYD3icgzqtrRdt+pWBAkIm8AGFFVdmVzEJF7AfwxgH+mqqu24MsyEVmD+gTwbwEoA/g+gC+o6tmeXlgfkvqI6b8BuKSq+3p9PUmwMiL/iqre1+lzJTZHToH9CYD3AviuiMyKyH/p9QX1i5VJ4D8AcAz1ybvnGMSNdgL4PQC7Vt5HsyujTuoDqRiRExFlGUfkREQJx0BORJRwDORERAnHQE5ElHAM5ERECcdATkSUcAzkREQJ9/8BEkY6GxG7J4EAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.scatter(x, y)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW0AAAD4CAYAAAAn3bdmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3da7BV9Znn8e/D4X4VOIDIRfDS6RCDURnNtDMdI2rZCaNpMzrdXaboaBWTqu6UzsSKGF9Mz1RNlVOZTmJ10jVDmdiktSdaiTYmlYuEiZ1Je0mAeEc7iaKiCBxE7ojIMy/2IpywnnXO2qy191577d/H2nXOXmdd/nvvw+M6z//ymLsjIiLdYUSnGyAiIvkpaIuIdBEFbRGRLqKgLSLSRRS0RUS6yMhOXLS/v99PP31BJy4tIhmicWRWwr5Frg/wy40bBtx9RpFz900+3f3IweHbcHDHj9z9yiLXarWOBO3TT1/APz+xvhOXFpEMR4+mw+aIEXEobmbfItcHmDBmxCuFTgz4kYOMed91w+536Mmv9Re9Vqt1JGiLiLSXgdUjG6ygLSL1Z8CIvk63ohQK2iLSG6xo1r0aFLRFhtCK3G1VRa8rK8/civegte+r0iMiIt1Fd9oiIl3C0J22iEj3MN1pi4h0FY0eEam/unY65lWf16+OSBGR7mEoPSIi0lV0py0i0i2UHhHpWc1MOMnaN+/xzch7rWau087JNS1lQJ86IkVEuody2iIi3ULpERGR7qI7bZH6K7pgVDP7vnvkaGpbX3B81jm7bxGnNtOdtohIlzBNYxcR6S6axi4i0i3UESki0l2UHhGpl3ZVGM86Z95Ox6jDMsuokfW4uyysxPW0zewU4C7gHMCBG4AXgfuABcBm4Dp331XKBU9Q2idqZn1m9ksz+15Z5xQRKUeSHhnukc+dwA/d/feBc4FNwEpgnbufDaxLnrdEmf8bvolG40VEqmdE3/CPYZjZZOAPga8DuPthd38buBpYney2GvhEi15FOUHbzOYCH6fxJ4OISPUcG/Y31AP6zWz9oMeKE85yBrADuDvJLNxlZhOAWe6+FSD5OrNVL6OsnPZXgM8Dk7J2SF78CoB58+eXdFmpiyosTPRe0IZ3gvzxuNHxHVn0Gva9cyS1bfK4UeHxefPXRfPU7VzEqjIs9+iRAXdfMsTPRwLnA5919yfM7E5amAqJFL7TNrNlwHZ33zDUfu6+yt2XuPuSGf0zil5WRKQ5+e60h7MF2OLuTyTPv00jiG8zs9mNy9hsYHtLXgPlpEcuBq4ys83At4BLzeyeEs4rIlIaMxv2MRx3fxN4zczel2xaCjwPPAQsT7YtB9a04jVACekRd78NuA3AzC4BbnH364ueV0SkLI1qY6Wlej4L3Gtmo4GXgE/TuAG+38xuBF4Fri3rYifSOG2phFYVAGjmvGGuOMgpNzNOOspfZx1/8N330m3qS7fpyOH0fpCdaz9RbfLUzTDDSnrd7v4kEOW9l5ZygWGUGrTd/RHgkTLPKSJShhLvtDtKd9oi0hMUtEVEuoiCtohIt7DkUQMK2tJ1ii7stO9QesILwJiCk1aORJ2hQadjuB8wcUyxf45FK980oxWLa7WSkW9IXzdQ0BaRnjBiRD1WPFTQFpGeoDttEZFuoZy2SD55c59FFzE6mDHhJNJM7jpaMCrKEwOMDLZH+etDwSQagHffCxanGpW/rmE0OSjK34/POQkH2lv5vdV0py0i0iXUESki0mXKmsbeaQraIlJ/pvSISC55c59RAQKI87R7Dr6b2pY1xjkqQpCRUg7bEC3YtDu4flYbjnr6nFk58RFBUIkWkYr2yxLl75vJRzezEFfVx24raIuIdBEFbRGRLqGOSBGRblOPmK2gLSI9wDSNXXpYKzqcsiqMh9cK/syNOhyzHM6oHLNrf9zBeKLZp4wNt0eTY8YGk2Oi/SCedDM6eF+yJgdFHZxR52oZVX5acXyrKT0iItJN6hGzFbRFpDfoTltEpEuYafSI9LCikzOiPGvWhJO39h9ObYvyxFl56mjfgb3pcwJMCCbHjBqZbteWtw6Gx0eiauzjRsc56WhxqSlB+6MJNxBP7onOmfVeR4rmqVuVPz8ZCtoiIl1Ea4+IiHSRsu60zWwzsBd4Dzji7kvMbBpwH7AA2Axc5+67SrngCeoxcFFEZCh2PK891KMJH3X3D7n7kuT5SmCdu58NrEuet4TutKVpUQHZLFH+NMqzRsUGIB6nHI2nnjg2Xth/2+53UttOmZDOMwNsD/ZtJv97ICjEcHr/+Fxtgvi1RotTTQny5FnGBQUPsvLMVV/wqQgDWpzSvhq4JPl+NfAIcGsrLqQ7bRHpAcPfZSd32v1mtn7QY0VwMgceNrMNg34+y923AiRfZ7bqlehOW0R6Qs6/GgYGpTyyXOzub5jZTGCtmb1QvHX56U5bROrPGumR4R55uPsbydftwIPAhcA2M5sNkHzd3poXoqAtIj3AaNxpD/cY9jxmE8xs0rHvgSuAZ4GHgOXJbsuBNa15JUqPSKJVkyCizrlmqqFvGUhPZJk1ZUxqW9aEmWgiz96D8WuNOkOjyuVZnZOTxqb/OUWdplGHI2R3pp4oa3GtqPL6xKBNzVQJKqpKHZkldUTOAh5M8t8jgX9w9x+a2S+A+83sRuBV4NpSrhYoHLTNbB7wTeBU4Ciwyt3vLHpeEZEylTFO291fAs4Ntu8Elha+QA5l3GkfAT7n7huTPxs2mNlad3++hHOLiBTXRM666goH7WR4y7GhLnvNbBMwB1DQFpFKMExFECJmtgA4D3gi+NkKYAXAvPnzy7ysDCHvhIms3GOUJ41y0tGEGYgX/I8WXIoWdoI4/xpNTnl+YHd4/A827Uxtu/6808J9RwX/qF/bk77WB2ZNCY+PqsQfDHL6U8ZnVI4/lN731CB/H10H4srx0effitx1N6jLnXZpn56ZTQS+A9zs7ntO/Lm7r3L3Je6+ZEb/jLIuKyKSS8nT2DumlDttMxtFI2Df6+4PlHFOEZHSKKd9nDX+9/R1YJO7f6l4k0REytVYe6QeUbuMO+2LgU8Bz5jZk8m2L7j790s4txSUd5xslHuFOH8dLWKUldPeezCdE49yr1FRW4Af/iY9sezyhen02vot+8LjtwzsT23beShesGn8yPQ/h+d2pM87eXS8YNO08aPD7SfKeq+mBgtZZS2kFYkWh4pUqTBBO9UkZpcyeuRn1KZkpojUVV3+p6QZkSJSf6b0iIhI12jDetpto6AtIj2ge4b0DUdBuws1U2EkbzX0rE6sqIMyWjDp4OG4w2zHvnSn35sH0pNrXtgRVzi/cM7k1Lb1b+YvvXftv0pPpNkeLOIE8PTr6fMunjMxtW3cyPi9iiriRBXS9wSdsxB/LtGCVc0suJXV6diLahKzFbRFpAeYOiJFRLqGxmmLiHQZBW3pmGb+zIsmwkTGHs03MQPiBZuyJuccOJLO374dLEI1LWMRpQeeTU+ueS2aMPPWgfD4+R9dmNr2q4FD4b5L5k9Kbftgf3pxqKzJMW8HufKRwWc1Y1K+STgQf35Zi2vlzXU30/9Rl5QCVCenbWZjgE8CCxgUg939v+U5XkFbRHpChe601wC7gQ1APD13CAraIlJ/1Vowaq67X3myB/fmwroi0lMaRRCKF/YtyaNm9sGTPVh32hXWzBjbrGKtUf4zKkyQdXxUmHZvkJPOWoQpWnDp3nUvpbZNmTI2PP6suemc8rz+CaltWQv73/foa6lti8+cHu4b5YSfCYorXHTatPD4CcGY7MjO/XER4ukT0rnuKePSY7+zgsu7weJS0fvSqwtGjajOrfa/Af7czF6mkR4xwN19cZ6DFbRFpCdUJ2bzR0UOVtAWkdqzCiwYZWaTk6pee4ucR0FbRHpCBbI//wAsozFqxPndJa0dOCPPSRS0RaQnlJmzN7M+YD3wursvM7NpwH00xl5vBq5z999ZzMbdlyVf05MHmqCgXWHN/JJldSSG5w3+TBzYG3eORRVlmlmw6MfB5Jg3Xklv2z89rnD+xuvpjsAzgo7ELVtStaQB6O8fn9r22NNbw32njk93+n14bnrBqsMZ1WSizyDqCM56/7Im7eS5DsQLeUXq3uEYMRojSEp0E7AJOPYLshJY5+53mNnK5Pmtme0xu4ZGh6QD/8/d/zHvhTXkT0R6wggb/pGHmc0FPg7cNWjz1cDq5PvVwCeGOP5vgc8AzwDPAp8xs6/lfR260xaR+rPc62n3m9n6Qc9XufuqE/b5CvB5YPC6B7PcfSuAu281s5lDXOMjwDnu7o2m2WoaATwXBW0R6Qk5B48MuPuS7HPYMmC7u28ws0tOsikvAvOBV5Ln84Cn8x6soF0R0cSIrBxlVKE7a2GoaBH+aMGnLI9u2Zlrv8deTueeAS46uz+17fHv/jS9Y0ZOe/vGnwdbL0xtmTBxTHj8/mAiy8J5p4T7RqaPTZ83b+4Z4Kjn3zfKdTeTf1bBg2xGaZNrLgauMrOPAWOByWZ2D7DNzGYnd9mzgVTHjZl9l0YOewqwycx+njy/CHg0bwMUtEWkJ5TRAevutwG3ASR32re4+/Vm9kVgOXBH8nVNcPj/LNwAFLRFpAdY6xeMugO438xuBF4Frj1xB3f/pzIupKAtIj2h7LVH3P0R4JHk+53A0lIvkEFBW0R6Ql1Gpytot1DUuQjxymtZq9RFoskVUYUUiDsdo1X+sirPnDYxvfreX6/9dWrbzp1x5ZhvP/ZEatuSZR9NbXvxuTfC42een+50HNGXfq8u+MCs8PjIv14Yd3qeNzPdQbntQLrKzeLZcUdm9LmMCtqa1cEcHR/lYdXheHI6vfZIWRS0RaT2GqNHOt2KBjO7GPgr4HQaMfjY0qxae0REBABra5GD4Xwd+E80Fo6K/8QdgoK2iPSECqVHdrv7D072YAXtFsrKU+edSHMgI88cLeKUVaF74L305JJowaM9h+PJOfdsSOeaF50+NbXtpYyqLQc/dH5q2/rv/SS17dRgP4D9u9OV1889f35q24xJ8eSatw+kX9eug+nKOwBTg8ox0ecybnT8uUb563Gjg8pBGX0dRVXoTrJyqpQeAX6SjOt+gEGFfd19Y56DSwnaZnYlcCfQB9zl7neUcV4RkbJU6E77ouTr4OnyDlya5+DCQTtZV/ZrwOXAFuAXZvaQuz9f9NwiImWpSsh29/TwqSaUcad9IfBrd38JwMy+RWOZQgVtEakEs/zrjbeuDXa9u99jZv85+rm7fynPecoI2nOAwSWvt3D89v+3zGwFsAJg3vx0TrKOmsldRmN0o9x1li1vHQy3R5XTJ41Nf+wPvTiQ+1obX9yR2hYtzASwa/uu1LYxp6ULd+x5Ky5icODpf05t2/e+61LbfrMtXfUd4K096XHW/2HxqeG+u4LXMHtqepx6VhGCqFshGv8e5bmbodz1yalAemRC8nXSkHsNo4ygHb0Tqd/qZE3aVQAXXLBEswNEpK06HbPd/X8nX/9rkfOUEbS30FgP9pi5QDy9TUSkAwwrfe2RTimj3NgvgLPNbKGZjQb+BHiohPOKiJTDjq/0N9SjGxS+03b3I2b2l8CPaAz5+4a7P1e4ZSIiJapATrsUpYzTdvfvA98v41zdIFqwJ+ocKjq5ZnTG8dEkjn2H4k7LzW+nJ6f89LV0NZoxffEv9CtvpDsIn9v4UrhvxPe/nd44dmJq06kL54THL/7zP0tf/6ktqW1//KkLwuOj97B/fDwRJ5pcE8mayJS3Sn1WB3X0O6BOx3IY0FeRoG1mNwF3A3tpFAc+D1jp7g/nOV7V2EWkJ5RVjb0EN7j7HuAKYAbwaRpFFHLRNHYR6QkV+qPlWEs+Btzt7k9ZE7kbBW0Rqb1GR2NlovYGM3sYWAjcZmaTgNyTOhS0T0LRPGOU644mYWQtGPXukfSEmazZXm/sS08uWbVmU2rbzFnxeP9XfpMqKs25F56Vvs7rcTX2M878/dS2yePTueP3z5kcHh8tbrVofnrBqnNnxcePGpF+r08LJsxAXOV+fDARppkK60Xz1Hn7T2R4FXrbbgQ+BLzk7gfMbDqNFEkuymmLSE+o0JC/te6+0d3fht/Wl/xy3oN1py0itWfAyA6nR8xsLDAe6DezqRzPbU8GTst7HgVtEekJFUhp/0fgZhoBegPHg/YeGiul5qKgXZJmFoc6EuQpo2K724OivBAvOLRjX7zvv53fn9r2f89Lj4k+8E5cGOC0memydS+/lh57vXtnnNN+4VC8kNSJzpz1gXD7/iCvf9lZ6Zx2liinnbXgU2RiRnGHSCtyzcpfl8Os89PY3f1O4E4z+6y7/83JnkdBW0R6QhkxO0lx/BQYQyN+ftvd/4uZTQPuAxYAm4Hr3D29xCXg7n9jZn+Q7Dty0PZv5mmDgraI9ISS/mh5B7jU3feZ2SjgZ2b2A+AaYJ2732FmK4GVwK3RCczs74EzgSc5XtjXAQVtERFIprGXELXd3YFji7ePSh5Oo/DLJcn21cAjZARtGmXGFiXnapqG/IlI/eWYwp7E9H4zWz/osSJ1KrM+M3sS2E5j+N4TwCx33wqQfJ05RGueBeJKHDnoTjtRdBJD3okZEHdERgs+ZS1gFE04eXlPemEoiKusjwxeV1ZH6qagSs3iRbNS2xadOT2+/sGgGvru9ISfs/rjCS+nBFV2po9NL/g0Zeyo8PhZU9L7TsjoXIzel0jW74UmwlSb5asSOeDuS4bawd3fAz5kZqcAD5rZOU02pR943sx+zu9WY78qz8EK2iJSe0b5MyLd/W0zewS4EthmZrPdfauZzaZxF57lr4pcV0FbRHpCGUHbzGYA7yYBexxwGfA/aBR+WU5jtb7lwJqsc7j7PxVpg4K2iPSEkhaMmg2sNrM+Gn2C97v798zsMeB+M7sReBW4doh27OV4Hd3RNDoz97t7vIDOCRS0E1HusZnF6qP8dZTnBtgfTGQJ5oCweeeB8PgoT53l8S3pIgafPDedk77vl2+Gx0+dOi7XdaZk5N+3vZ2uEn/W3CmpbfMnx9eZOT6d644mF2UVJpgWtCtrIa5o0s3EIKeeRfnr6jKDoHZI09z9aRpFC07cvhNYmvMcv7M6m5l9Argwbxs0ekREesKIZFbkUI9OcPd/BC7Nu7/utEWk9lrREXmyzOyaQU9H0Bi3nXvMtoK2iPSECiwYdcy/G/T9ERrT3q/Oe7CC9kmIcp/R2OushYmiccLbgsWhZkyMC9A+/XI6Tz11XPxRXnFGevz0+q3pBZ/mTx8fHh+NXR4ftH/+1Lit1yw6O7Vt8uj0mOrJ4+Jx1lH+P3r/oj4BaG5xqKivIhp7nUU57SozRuQbp91y7p674EFEOW0RqT2jOkUQzGyumT1oZtvNbJuZfcfM5uY9XkFbROrPGn81Dvdok7tpjOs+DZgDfDfZlouCtojUXpXutIEZ7n63ux9JHn8HzMh7sIK2iPSECg35GzCz65OFp/rM7HpgZ96De7IjMu/CPlHV9EzBRJpd+/NPgsnb4QbwB3OmpbatfTm9sBPAroPp8765N92ud96LO9zOn5+epDVrUnrCyviR8eSWqWPT+44O3teocg/EHZSjR6Y/qykZHZmRrIW8IqqcXh8VGj1yA/BVGsV8HXg02ZZLTwZtEektRnXSCu7+KpBrRb9IVV6HiEjrWHXSI2a2OlnW9djzqWb2jbzH605bRGqvMSOyMvmRxe7+28kS7r7LzFLrmWRR0B5CMxMrDr6bXoQoyr0CHDyczt9Gi/XvPhBXSI9y3XOD4yGuZv6HC05JbdudsQjVuTPT+z4/kK68fva0SaltAAeD60ftz5oEM2dqXBzhRE31P7SI8tfVVqFPZ4SZTT1W+DcpCpw7FhcK2mb2RRpTMg8DvwE+Pfj/ICIiVVGdG23+GnjUzL5NoyPyOuC/5z246O3JWuAcd18M/AtwW8HziYi0gGE2/KMd3P2bwCeBbcAO4Bp3//u8xxe603b3hwc9fRz490XOJyLSClUaPQLg7s8Dz5/MsWXmtG8A7sv6YVLVeAXAvPnzS7xs86Lc475D6fxx1njeqLjBqGCF9ahYL8QLEw3sPZzalrWw/xOvp8fh/94pGTnlI+k2REUUTg2KDUBcRHhRf7qIQZSThzh/Hy3uNGV0/KsYLcQVTTfO6n+IPmstAtWbKtQRWciwQdvMfkxc7v12d1+T7HM7jSUG7806j7uvAlYBXHDBkvz/akREirLSyo113LBB290vG+rnZrYcWAYsdXcFYxGpnKqlR4ooOnrkSuBW4CPuHhc0FBGpgJ650x7GV4ExwNrkDXnc3T9TuFUiIiWrR8guPnrkrLIaUlQznUuRqOp2NDEE4FAwkWZiMGFk4ti4I/GtfemOwKhzctrEeBGkJbPTC0YdzchMTe1LL9i0c1u6Ss6ZkyeGx0cmZVTJiRRd3Cl6XyJZHYbvRp3GFZiII+1lQJ/utEVEukdNYnZtcvMiIkOwXP8NexazeWb2EzPbZGbPmdlNyfZpZrbWzH6VfJ3aqleioC0iPaGkyjVHgM+5+/uBDwN/YWaLgJXAOnc/G1iXPG+JrkyPRPnraMLLuIzJMdHx0bYxGbnPaHGjvBNuIM5VHwjy59HEFojzvKdNiSfHvPDG3tS2aHGnaHIRwOxgwaYof5+VU4/2jd6rvLnrZvdV/lrg2JC/4vkRd98KbE2+32tmm2jUebwauCTZbTXwCI2RdaXryqAtItKU/HfS/Wa2ftDzVcnEwPQpzRYA5wFPALOSgI67bzWzmYXaOwQFbRHpCTmnsQ+4+5LhdjKzicB3gJvdfU87x4Drb0cRqb1GEYThH7nOZTaKRsC+190fSDZvM7PZyc9nA9tb8DKADt1pO+kccjML80T7ZuWvI83kvyNRTjUsTBunpNm2Oz1OOlpwKSun3R8U1s0ydUJ63/nTx6W2RTn1LPuCIsTRIk6QNX5df+BJ++UZHTLsORq31F8HNrn7lwb96CFgOXBH8nVN4Ytl0L8eEekJJWUwLgY+BTxjZk8m275AI1jfb2Y3Aq8C15ZytYCCtoj0hDLutN39Z2TPiF9a+AI5KGiLSO0dy2nXgYK2iNSfWe8UQWgFoz0VQbIWkYo6HZtZcCo6fszR9ECct/anq9FAvOBS9As1LehEzDpv1oSTqNMy6nTMmkgUnTeq6JP1eeadyKQKMdJqdfkN0522iNReIz1Sj7CtoC0iPaEeIVtBW0R6RU2idlcG7XblRJtZWD+SlZPevjc9uSbKHUcLUwGMDvLPzfzpF00EamZxrKKUv5ZOUHpERKSL1CNkK2iLSK+oSdRW0BaR2jPKmRFZBbUO2kVzp1m56yj/HF0ra+z39CDXHS1ilbUIU1RcISsnHS3uFBXWzcpdFy0iULTgskgp8q+nXXm1DtoiIsfUJGYraItILzDaWaiglRS0RaQn1CRmK2iLSP0ZSo90VKs6GE+UtQhT3k7HZiamRB2JWZ2AfTmrwUM8oaCZdkXvVd6O2KG2i7RdTX4VuzJoi4g0S0P+RES6iHLaIiLdQuO0u1szOdlI3vx1Vk46yhM3M4klamvW5Joo11309TdDBQ+kKuqSHik23S1hZreYmZtZfxnnExEpk9G40x7u0Q0K32mb2Tzgchpl40VEKqlLYvKwyrjT/jLweUCLTIhIdVmORxcoFLTN7CrgdXd/Kse+K8xsvZmt3zGwo8hlRUSaNiKpyD7UIw8z+4aZbTezZwdtm2Zma83sV8nXqa16HcOmR8zsx8CpwY9uB74AXJHnQu6+ClgFcMEFSzp6V97pjrCsSTsnylohL+pczOqIjLYXff3NdC52+r0WOabE38S/A74KfHPQtpXAOne/w8xWJs9vLe+Sxw0btN39smi7mX0QWAg8lSzEMhfYaGYXuvubpbZSRKSokqK2u//UzBacsPlq4JLk+9XAI3QqaGdx92eAmceem9lmYIm7D5TQLhGR0jRRBKHfzNYPer4qyRIMZ5a7bwVw961mNnO4A05WT47TFpEek39I34C7L2lxawopLWi7+4KyzlV1UZ62ndXgx43uK/1aWfnzdr1WkVZr8W/tNjObndxlzwa2t+pCpUyuERGptkYRhOEeBTwELE++Xw6sKdzkDAraItITypoRaWb/B3gMeJ+ZbTGzG4E7gMvN7Fc0Jhve0arXoZy2iNRemXNn3P1PM360tKRLDKkng3a7FjFqJk/cKnlfq/LUUns1+RXvyaAtIr2nLqv8KWiLSE/ollX8hqOgLSL1Z1CXDKCCtoj0iHpE7UoH7VZ15LWr0y3rOu2s5qIORpHjRRDqoNJBW0SkLDWJ2QraItIbdKctItJFCk5Tr4xKB+1uyseqMIBItdXlX12lg7aISBm6qdr6cBS0RaQnaEakiEg3qUfMVtA+pujYaeWpRaqtLv9CFbRFpAcYI2qS1FbQFpHaq9OMSFWuERHpIrrTFpGeUJc7bQXthDoSRepNQ/5ERLqFJteIiHSPOnVEKmiLSE9QeqRm2lmYQETary532hryJyI9wXI8cp3H7Eoze9HMfm1mK1vS2CEoaItIbyghaptZH/A14I+ARcCfmtmi1jQ4pqAtIrVnwAizYR85XAj82t1fcvfDwLeAq1vZ9hN1JKe9ceOGgXGjbD8w0InrV0g/eg9A7wPoPYDs9+D0oifeuHHDj8aNsv4cu441s/WDnq9y91WDns8BXhv0fAtwUdH2NaMjQdvdZ5jZendf0onrV4Xegwa9D3oPoLXvgbtfWdKpotvx9CiGFlJ6REQkvy3AvEHP5wJvtLMBCtoiIvn9AjjbzBaa2WjgT4CH2tmATo7TXjX8LrWn96BB74PeA+iC98Ddj5jZXwI/AvqAb7j7c+1sg7m3NR0jIiIFKD0iItJFFLRFRLpIJYK2md1iZm6WaxxlrZjZF83sBTN72sweNLNTOt2mdun0dOAqMLN5ZvYTM9tkZs+Z2U2dblOnmFmfmf3SzL7X6bZUWceDtpnNAy4HXu10WzpkLXCOuy8G/gW4rcPtaYsqTAeuiCPA59z9/cCHgb/o0fcB4CZgU6cbUXUdD9rAl4HP0+YB6lXh7g+7+5Hk6eM0xn32go5PB64Cd9/q7huT7/fSCFpzOtuq9jOzucDHgbs63U6dSMMAAAF5SURBVJaq62jQNrOrgNfd/alOtqNCbgB+0OlGtEk0HbjngtVgZrYAOA94orMt6Yiv0Lh5O9rphlRdy8dpm9mPgVODH90OfAG4otVt6LSh3gN3X5PsczuNP5XvbWfbOqjj04GrxMwmAt8Bbnb3PZ1uTzuZ2TJgu7tvMLNLOt2eqmt50Hb3y6LtZvZBYCHwlDVW15oLbDSzC939zVa3q52y3oNjzGw5sAxY6r0zcL7j04GrwsxG0QjY97r7A51uTwdcDFxlZh8DxgKTzewed7++w+2qpMpMrjGzzcASd++plc7M7ErgS8BH3H1Hp9vTLmY2kkbH61LgdRrTg/+s3bPLOs0adyyrgbfc/eZOt6fTkjvtW9x9WafbUlVV6IjsdV8FJgFrzexJM/tfnW5QOySdr8emA28C7u+1gJ24GPgUcGny+T+Z3HGKhCpzpy0iIsPTnbaISBdR0BYR6SIK2iIiXURBW0Skiyhoi4h0EQVtEZEuoqAtItJF/j+JS4xEDGG0WwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1)\n", "h = ax.hist2d(x, y, bins=50, cmap = \"Blues\")\n", "plt.colorbar(h[3], label = \"counts in bin\") # h[3]: annoying syntax thing to memorize or look up" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "While there's often no problem with square grids, hexagonal grids often look a bit nicer. Note that the relevant keyword for hexagonal grids is `gridsize` rather than `bins`. " ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW0AAAD4CAYAAAAn3bdmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOy9aZBlW3bX91v7DPfmzbHGNw89qSXRiBY8N8JyYEAGyyAjBQ4RDlsOWVZE2xGYELYJaPkT+JMctgk6CBzhDinkxoCNAoMlCAmhaCxsjJDUDd1IPUit7n79pnpV9aoqxzuds/fyh33mczLzVuV0M+v8K25l3p377L3PtM46a/3XWqKq9OjRo0ePywFz0Qvo0aNHjx6LoxfaPXr06HGJ0AvtHj169LhE6IV2jx49elwi9EK7R48ePS4RwouY9ObNm/rqq69exNQ9evS4ZPjc5z73nqreOskYwcYrqunk2H46uf+Lqvq9J5nrrHEhQvvVV1/ls5/97EVM3aNHj0sGEfnmScfQdMLgw3/62H7Tz//1myed66xxIUK7R48ePc4XAnI1rMG90O7Ro8fVhwAmuOhVnAquxqOnR48ePY6DyPGfhYaRLRH5uyLyFRH5soj8ARG5LiK/JCJfzX5eO6vd6IV2jx49ngJk5pHjPovhk8A/UtVvBX4P8GXgE8BnVPVDwGey72eCXmj36NHj6cApaNoisgH8QeCnAFR1rqrbwPcDn866fRr4gTPai+W2aauCyz6h8cdUFazzfw8qbakDI/5zWJvLtg2Mb4N2W3XOrrbmOpR6W+q8+eys11ad8zTXls951mvrz9/Fnj/rlNRBFIAROafzpyTWjxXkCzkvCItq0jdFpEpt+5Sqfqry/f3AfeCnReT3AJ8Dfgx4RlXvAKjqHRG5fToLb2NphXbq/CeHtf64V3MSpo02m/1ybJtdvK0556Jt57G2Ji5qbV1zLrK2p/38XcjaVHFatqUOAlHf0f936ucvzyTqtJxTUIYhmHMT3rKozfo9VX3tiL+HwO8F/qyq/qqIfJIzNIV0YWnNI1WBnUPbTWfedhFzPk7bMqxjkXWdxzq62pbpXC3D2qweto7jBdpJ1uG03SfpuMfPFCY4/nM83gLeUtVfzb7/XbwQvysizwFkP++dyT6wxEK7R48eF4erl2dfTsURqarvAm+KyIezpu8BvgT8HPDDWdsPAz97FnsBS2weab6K9ejR4+TIrRGaadylfu2/aWY+yTvl/WVBOtxR82o2b45ztWoLi5pHFsGfBf6WiMTA14EfwSvAPyMiPwq8AfzgaU3WxNIK7Tjwr1Tn/grVo8cVRil81csx3+hbqgI7g9NS0J90TlFFs7Hi4CKckadjWFDVzwNddu/vOZUJjsHSCm0RCHIvdK9y9+hxqvCCtK5rd/Y73UkRlEEgp6j0Ljz5qQnti8bSCm0oX+F69Dg/HC/IepwE4umY532IBQiuRhj70gpt63rTSI9zhGpFQcj4FBUa3NWBUvcx6pFa71k8wpKMzus54qc8+FE4f/X+TLC0QrsX2D3OEw5oEiYMV+Y+L+Cy+6q6q4IPsClNklqzZZ/UCdkFJQsCOjfltzeP9OhxpXDlGG6HoGs3myJZRM7ffHEeuCJP4KUV2j3lr8d5oklJO637W1VrvpmTaq9NhocPEV98rGqIPbQpf1cavaZ9togD//rUM0d6nAckYzZo8R1OKsROm0LnnHZHM6ouLLjzfqYy0lmYP45DZC7Ant1r2mcLEQgDwPaCu8d5QRBZXo3zMNPGkwjdRSl/Z4HYgLkIpfeKFEFYWqENPeWvx0VgUSHWJfCWReAvyzq6cTH3dO+IPHM4hcT2QrvHEkF9NrwyBLxMtaSUbcdZVkqTRt6ymIA10k685NSbOso4x651VMc/ivJ3PoI+cT5xVWTO2WLRm0fOFvOO9JA9elwkckdlLvRyfVal3ZY7CKsUumKcTOg/LqVQRApbdG4bz00jmvHMu9ZRhcukelVuXwS10WnJ1T4XCL2m3aPH04gmNVCL/5oopaCItDfs6LcIciHdRclbhLZ4eJ+roYUejt48cuboKX89LgpVzThPDlewLirmkfy736Zs813rduXjqHZd1MDHdTA219a1+VNN+esdkWeLOPCvT7aX3D3OEU1anVb+a2bIgyo1sNrWFn45W8O02o6gBvL4VL7utbX7tSl/53ujnTvlD3qb9llDxNu7pKf89ThHdNLqGlpvF11uMQqddGrh3T2fTO9dOHtfq9/5CbQLofxJbx45N/TyusdyokvInbXgW5RmuOg6robmuTB6Tfts0VP+epw3VLXTl/K4popDRm+wNsq5upyDOXUvj3Ysy3/VaYbVsY7lGi4J5g7MBVD+LiLy8yywtEK7F9g9zhNOm/xlj9PKdNdFtcstFIdd5y6X6FJfnxGfla+ZqS93ml4GnDfl70kjR5cRpya0RSQAPgu8rarfd9LxeoHd4zxxmMA+rRu9U5tehKKHd45W4VqO0aeC+3EyiCDn7vk8G5ympv1jwJeBjdMYrKf89TgPqCrWefZGVVPN20IjhIF0Cu8q6yOXB1ZhlvjIsEEUEJoKVZBKFsFsrjJgp5vyJwKRkWI94GsrHkfva65N1feptl0VzXNRXJX9PRV3qoi8CPwJ4CdPYzzICn9ejWPcY0nhVElsKdwUSG29LXX+uzbU4qoQ9WPBeG6ZzC1O/ffJ3BZasohkZo1cYJZtgh+rSTU0AoH4B0ZghCjwn6pAr48ph65Nod12gUnEQ4HwnMkckh3Loz6XAaelaf9V4C8A64d1EJGPAx8HePnll48dsKf89ThrHOYAbCIwi93QTa41NDXgLspfya+u96wP8OQ0w+XDRWX5uyxC+Tic+NCJyPcB91T1c0f1U9VPqeprqvrarVu3Fh6/l9dPM05y9p+WK+dqCKIzhyz4uQQ4DU37u4E/KSJ/HBgCGyLyN1X1h04yaO5d7tJeelxllGHd+ffFs9B10OqO2LZL8erypVinLRZJV7/A+MIdtW1VCTr7V6l8/i9aKWZQJH7qoPzJYYu/JJg7CNSbR85rN4TLY/44DicW2qr648CPA4jIHwL+/EkFNvSUv6cVmjEjqg9rIXslPOaeK2zTlW2PymBnRIgCv53LBHMU+ZfPxDpfeFa84Eys97EYk9mijbTCzwdhgFNlnnrJHRjBOVAUY0p7M9majPh508oggZROytw5CnUnpVcKD+eNB4215fburraLglXgPLP8AeZCKi+cPpaWp90L7KcXzbcrVVBZQM8+9KI5fMs8fWrQoINFgaF5Ffogm/q20shPbUSIQ1NrU0phXV3rvKmW0xamXcE3i1ixu4rzXtmCvQui17Q7oKq/DPzyaYzVU/6eXjST/UtDYB+mMTZpdZUtyE0uuVYtcriDMe/XpNXl7JLAlHS9JHWoQhAIpkibKqiUyaceZ72JVYxo8RApCy6UjJPmiruohycpHAxXkBp4ijZrEXkd2AMskKrqayJyHfg7wKvA68CfVtVHpzNjHUv7vhAFF5AFrMeFI6deValsJpdYHEJn07rJoYsGpxm9z2aC1KkXwE3am6uMn9PqqlA8DXCWOGaJK+h0+Vj5vF57r4/XXm9O+avDp3BQHwVZmTenAVb3q3N8nozOVzzUzoEaGFx+yt8fVtWPqupr2fdPAJ9R1Q8Bn8m+nwmWVmgb8VztsBfcTyVygelvpEUvAt+3vq1v75I7Xfdos9thN/KR1MDsI41Q88PWfJhDtAkjlfGPWIccse7HxVncfrHJ6LzneG8LxwvsEx6z7wc+nf3+aeAHTrzoQ7C0QrtHj0vDwepxKSBGjv0AN0Xks5XPxzuGUuAfi8jnKn9/RlXvAGQ/b5/VfiyvI1J9AdCe8tdjEWilUG6pVbdpdZUtcApBZuuu0urKsbrbDoNz6tkl1E0UVSrfwvtDcy7NGC5ezzqMeuhrRZatx615kXWc9qPzIih/LG6bf69i8jgM362q74jIbeCXROQrJ1/g4lhaoT3vKX89OtCks0HFaUbbsSeZXbhK1KgWwbXOERghDBr24UxIV9ty+qARMBkNI7Va2JpTpxgt6XrjuS043oPQlDbsivAwAkHgBXHqFGs1W09pw7eupPylzhEFkjlD69RD5xSrWqSYjULzWEV788LB1WN4Vo7Ii6D8nV7yL30n+3lPRP4+8DHgrog8p6p3ROQ54N6pTNaBpTWP9AK7x2HochBC7jSrt6WZ0Duqn3VKkrpWxfTUaqutzPuR5QQJpLBl5zbsSeLYm6bFvLlzsMlWCQMhDEyxbRQYhnFAFJZtgZFWXpK0UYMvHzN1Jf2wi2a4CPyDxc+7aPj+ZcFp2LRFZFVE1vPfgT8G/Cbwc8APZ91+GPjZM9qN5dW0e8rf1UGVRlalrnX1yzXKKnOoK0ikxV0+hFYn2YVU7R4awXVo61Kpr5i6kvZXdexVCyKoes1YG3OGRjChYWZdsc556rDOEYcBRvz8k7kjNJ7XDV7IzlNHGEjGE+8WvE1OuaqSdnC+u/otc8DNWeIUIyKfAf5+NlYI/G1V/Uci8uvAz4jIjwJvAD94GpN1YWmFdhT0YexXAU1TRv7q3awE0yyo23Xei6IAR4yf98vHNxn5pGYyEUEUjCppzWySUQqbY2UUvny5TsFZ11pvLvQDQ6atGsZzV1AMnYXEWuLKYHOrzK2trd+myiy1RVrXHEYoBHx+7BLrWpp3WDGv5GvuettoHkvXcV7OChdB+TsN47yqfh34PR3tD4DvOfkMx2NphXZO+UtsX5H9MqOTktah0Z3dKa4Lr3INUhQTOG4dQcer83FUu8Ip2tFnEede19+7zBVdmvhh+b+XBRdT2LcPYz83LO+l16PHWeMsuBtPgmVZx8mwzA+yx8HSPnpUvZbd59K+3Oi6TfJKLbWGI1FS5/JP2dbuW6XaNber/61ukmlul/9uG9se9V7Q3La5/49zOTf3LY9OrK6zK2rYs03KY/TkyPe7YzEnwNxBak91yMUgC3wuAZZW0+4pf1cDBY2s4gysOsHg+Js3z/xXG4OsrdLPqecyzzImSBQIwyggtY7J3GIVhpFhJQqwqgVdL3eO2qxKjV+jt+3OE8cs9bTAtWHAoJEMitr82T5nlMLEetqgUHemOvV9jPF5RHLBm7NE8raqZijkofhZuTFyVoxgpF6KLB/L0/3yPvU15t+r8rgZ9l/9m5dpx3PVF0WqoJeU8nfRWFqh3Qvsq4PDKHqwuKO5y9nYxCxxtTSniVWmSVLrM01csaYcmm1bReqU6bx0EFqnHMy8c/C4mz+p0Azzvs3sejnjxFTGOiwDX1OAW6dUXZd59r5WhkAtbeuHreMwGdzKtFh0u5yC7xTC1JcGSyu0e8rf5UDO3jiKyncUpMnsyP5rarPVfodR1/LAF1fRDqPAOxwLznQmfOPAEGWJbeapYzyzxKEpWBdJ6kisy4rollQ+qxBUKH/zLMtfHJoiZ3ZeeDepeNCb1EOnynRe0vskowH6YJx2EE7tWBxyjJ6E/XHY+euiSl529EL7jBEFfSGEZcaiVL6jICLZw7l+lgV8atNc+Er5EK8GkECe9S8zc4iPbKzT+/z3wMDB1JYFe61FZ/Vow8TaVsXG1CqDUFgfhcXfrEKS2LpWP7eM4oAokGxOr/0mrhICL4Jzjnla7nFO74sr/DeX9R+GuQAtj2eZepWCSlhFLrgD007hWsVx5y9PG1sd4zSFnrkAyp9ckbShSyu0fehvT/l7GiC5WlctWgsV6Vtqf11pQpsiv03vA2e7+dxdZoAmVuKgZsaAbqpdVKHald2bDzFpP6S8wbrGSwykLbCbbzLFdg2cTiSj1PbhNM0ifWHfk2FphXaOq3GYn3YsB2XsiU1u57H0iz88R2CpF7cYnsB0t6zoKX89zhAVytiR4lIrNMA2Xa+77YjRWhS90pzQnHcR+Oo0Wmi1ubmjRQ10nr1SbWuvt2NObSvMrmtbbc/ZHv/sChecFubOJ/A6zyWWNvujP5cBS6tp95S/5cZRVD5YjDJWCJ1qP/U0tdT6grcCGTVOmSaWudUinDufLyDP5OdNILPEYZ0SBT6vh1NfaSYvH+YKadGky/nakAqkqSvMcgczyyxxrA4DosD4qjWpQ8RvAzBLLDvjlFEccG0twogwT/0+5NVmVJVZ6oW739bbjr0jsl0xZprkjsrS5FKlSlar+OS2aMmOq0OytvJ4P875O2skzjNZzo/y17NHzhy9wF5+PA6VrzTZypH9rCpp5fVKgVnqigrn+TapVeKKo06ASWqZV7ZNrJJU8nrk662On9PlavZofFpTW6EBpk7ZHacZ4yR/6MB0bmsskfHcovvKaFDeWv5BYuuOv4xxUk3qFEjdVq7ZPgzCtjBNGo4eJc+PUvZrmMlbOOr8XUWYK+KIXFrzSI/LgapZoPq6nkflHfaanvdJrauZAlKrJJU2KKl7Utl2llgOZrZlRghMfY7UukqEYIXJUYF1ymRuSa0r1pxabW1XBOm4sm0yT5kmabEOEe+4rCZ76kpqBZBaV8wJno6YWlc7llXKYRVNdoinCdbn9Me3vg+uca6eFPlYaeNcLS0WMI1clgfY0mracU/5W2p0UcZyYZA3538PTWlTJOs3T8ssec4qqo7Eaq0tyKqSBwaMCqERDmYp00SLWSdzx0pc5p8OAx/NOJ2X5g3UB6MExjtETcYhtK4szgue8pcHvFT3wVpbK6IwT1NUHXNbmnfmNuWZ9Zjr64NCmIbGm1bSRtZUkwkIn/nPR2FWnaRO/X5vDoNCE8+JIrlLN5CyEENVWxcaebWzOfIIz/L8PXlGPx+ZWZ58Z9XTDB8jUZWR+gP2rJGb2a4CllZo51n+UtdT/i4TmqeqmY8aKCqiV5F2tBnTptDNkuMuBt+x+5ppmGY6Uv+Wbs96vybmtr3e62uDurYLLYEN3UFIzbFy00ktfW1h+C8FeZcd+qydkF1vDfIYNMPInK/AznFZNOnjsLRCG8qgih49LhRPzBU8fLinGRel8F4VR+TS2rRVfSawnvK3nFhUmWtR1zpsytAtyFzDJq5ZkqUFVtduOYQGeFi/o4drU+0EmFfs01mH7iyHRy8eoGY3z8cqv9edlc31HnaEzlIDP8p/0cTMnj/lr7dpnwN6yt9y4qgbM6eRVXnCTsE5T10LTFmfsZojxGTOO5uxLHLziU29fTrOKHFJxhixDpIsFDyqZEDK7ezjmWeRmIxWlzpld5qQOsfGIGJ1EGb2b8NIhHmWZ2SSWh5N5xgRbqzEDMPAO/OK5E9eIB4klnGSEhhhlJUPWx0E7E1SkihgNPA8tvHMFbTFwp6vXmAFmYmgS/vzVV2EWVI/brkgzzMcevqfv08M2QNCfNUaVWph9nl+kXy7kwgob7ahSAGQH+fH0WTPm/InSF8EIYeIvAT8DeBZwAGfUtVPnnTcXmAvJ47LypdT6OYN2l5ic36y1PrVait6mgj707Q237RiGJZMKIUdN/vB1DJrUAMfHcxqwmt7mrASB0WYuDHCMA742t39ykjK3YMpt1eGNUE0Tx2786TMS+IUFbi1Pij6zVJXW0O+/81nnXXaquEIsDYISq57x3HL96u5ad6W+w9EIKBNIcyL9p4URqRwhF4WXBZN+jicxqMnBf4bVf024LuAPyMi334K4/Y4A+SvsVU622H9mvSwo9raBQbqcBlNrypE5qljZ5yQZNQMVWWeeE521axyMEvZPpjXTAZJ6spIRbxw2pnOmSRp0TZLHd/cmfBwMi/ajPh5q1Q765SNOCKq0fRgb+618xxRYNhaiQqBDxAH5thjCXndyHpbU4gUGusR4+TnbpbYWvQltCMhc0dlbU4/ypFrvao4jWrsy4ATa9qqege4k/2+JyJfBl4AvnSScSPjX6F6nB6aND2XqYBNBkKzyK6njblWW1V9VK3YjRtz5kUJwGfD0yxCMS9CsDe1WaCMFHOnWWDMziQthNP2OGFrFBFX1OzEKrvTOQ8n82xC2J3NubM/5529GeAF1SAQXnv+GqMowDpP95MsItAqbAwi1uOQcZqyPfXa9EFqOUgtm4OQ22srhQDcWomZzFOiwGSV1TOKolEv0CvH0ptB8sx/VNorJgqFlUFAHFQz/ZUCuKgE7/yxzIdJ5r5IcByaso+CaJmm1keM+odmF9/7oiCcM4PkEtmsj8Op2rRF5FXgO4FfPelYuWbSU/5OD12mjU7KWMe2i56CZr9qqHUOb0OuN6a27aDcnaStbHph405XVR6M57W2aep4e3dWW8sojhiGQW1frdYFqYh0Uvk2hnGLCz2M2reO6fBmVW29+Z+axzwKJbPbV+dQmi7FxLnW2rq40VVKYf6XoLqAC8ZFUP6atNPLjFM7dCKyBvyfwJ9T1d2Ov39cRD4rIp+9f//+gmP29KgeT4Yus21bHC+G074GT1N2LDzUEt1IF0f5uxrskVMR2iIS4QX231LVv9fVR1U/paqvqeprt27dOnZMzbzsPeXv9HAo/ewJuFddFLoTraNhj4Xu1/mm/Vi1o7ittt8q5qlDkMa2bft7l29tltqWxp/P/URobJYPXY6nnTTDTlrkYSbEBvWwac66SMzs4mXmThPGyLGfy4ATC23x7xw/BXxZVf/KyZfkMbe9Tfu0kds5a235HyowjTccwYdk55nq8jDmSVaTsWp7LQrG5rkpGqlKVZUgMAwiUxsvzfrm3OzJ3NYclKpK4hxv74x5MJ4V+U4eTZNCAFinTBLLr7+1w5fu7PFgf1YI9bla/uW7D7k/mWV2dsu7BxPe2Z8wzpyXkyRlL0lJnavN+TuP9vnaoz2maW6PtxzMUmZJ6cwcRIZhbAoKYr7/k3mZ7MqbVQzDyNTMA9Ype9PUh/FrPfNeNRdK0siHEgVClGU77DrZ+bZ5sYfSTn7xwntufRzGuUF6R2QV3w38J8BviMjns7b/VlV//iSDXvxldfVQLe7qEyd1X6Q5ZazuCBMMPsdGM5Oeyeh7ed9AYH9q2+HgDUeciQz70/LOVYVpYtmfprVtbSY887a9WcreLKWJX3n9Eb/93rjo99b2jBtrEaMsN8nMKb/9cJ+d0by273fH01bQ49w5ppWEI4+mCdN0j5fWRsUbQKpKrLA5CgoOcGDAOltTOFKrjGIYRqVNfWB8RsDqsRnPbY3yl2N/1pZuG8Ogxjs21J2W0KYa5r8ezU85P5znPe5t2uc44RniNNgj/4ylspj1WATHaRW5hqaqRRCIZhF+gZTOYVWf/CkMpKhUnrMcfLKnTMBZx940ZSX2OanBC+idacJqHBBmAii1jlQVQznn7izhILVcG0aFwKwQLwDPjpgrDKKAaVIKub2pRURYzQS3qrIzSxlFJVtjklju7ifcXo0YxSUzJRQhrUi92Pjc3FWnlhFvxpNKQibvLC2drf7NwaEqRXIrgDgQX1k9myJ/cDTPTBRILQlUVyDLUUWAewUIOOV82iISAJ8F3lbV7xOR68DfAV4FXgf+tKo+OrUJK1jaiMhQenv2RcE6z4MuqGW2DAQJjMkKB2TpTB2AMks948E5mNTsWsp4bplkWuXOOCEKhHFqC01zewIrkU9nmptTHHAwS/jmzoR5Zrh9d3/GSxtDrg/jMgBFlW9sT/jqgwkWuLEWY51/QATG8HBseTSxDALhxWs+oVNqLXszixHhvf0574291v61BxOeX4/58G2vTceBEGVUuWvDmMhkBRJUiY2wsRJiRLK83cogMlnhBSEOA1Rh+2DOLFWmiWNv4qMon7s2LAo0QG660Fo61xwiwiA0DPBvHJEpK8YfBwNoZupSrWqaXY+G88W5U/44dU37x4AvAxvZ908An1HVnxCRT2Tf/+KpzphhaeM6wwAGwRIv8AqjmjY1Rz1nhyBIK4Nd2hLYvvrKeG6L7HkK7M5SDhpt1rma/Rvgjd0JU+tqNtnNQdQyAXzp/pgkoxZ6u7ohDLyZQcm4y5WglXzOB+OE+wdpbfzVOGhVIL82iBkE9eK+68OQwNRpenkhhcI+KjBrRIYa4yM6q/0CI0TmcLtq3j4ITG3bY1GZwxgq212swA7FZ/A8V7+fnJ4jUkReBP4E8JOV5u8HPp39/mngB051/RUsraYNmePMHOEh73EpcLqv6AveWJ1zNvjomfZZI1VI29TQJR8X0toO2+kuRfcxhPCTYXksmN7cdr5znjJP+68CfwFYr7Q9kwUaoqp3ROT2aU3WxNIqsgXlrxfYT4RF6HjdtL1ualnH1q2WThJDR89O6loH/a7LXJAzO4pVqPoyW9XxpS0Xu1gToZEWC64aQl/O2d7XrtD15lhdcQbq8iCiJ3uMHU8zzNbVOKfLwhoBzxy5CMrfguyRm3k8Sfb5eGOM7wPuqernzn8PPJZW0046EtT3OB7V7HpwCKsg65NmFWQCk4VC42/rMHOQVbPw5dVc8m2TrMBuVdjm2zqnRQY4VZ+fI3WuyJi3k+X0GIYBBu/s256lxIFha+CdjY8mCdtTT/sbhKWG+ZUH+9weDbi9OgCUrz2ckFSyzeU332ockFhfXSY0cG0UFtnvqvu2MQwYJ95puBIZMD6p1ObAmz+GgWEY+gRTNnvAWVXu7E5YG0RcG8UYkcKeDaWTNjTC89eG7E6SzCEKo2HALHVZCHwl416Dmh1kVgznygdQNdLxsPOeM0b85t1tfoiL07wVL7hDoTPx11lhQUX7PVV97Yi/fzfwJ0XkjwNDYENE/iZwV0Sey7Ts54B7h69DBsB/gHdaFjJYVf+7RRa4tEK7F9hPhuZx68oIpwoHVbqZVayzDMJq2a4yEKUqKMbzsnhuTu+bpyWzIc/+tj+dV4S+EAcBX9vbraUk2JsnPpQ8+56mlkfThPcOkqKUF4CmykqWw9MpvHsw47cfHLA9sTVndVypI2mMMDABL18LGEZ12/PDg7RYRxQYNgPDzbWS1TJ3/iHykZubRBVanVVlaksq4+4sIQqFV66NCCpetaBx3G6sD9hYcTXTS+qUQaWifL5wKVKn5o7fo+mZVVTPfZPuR/5dnk7KH5yOeURVfxz48Wy8PwT8eVX9IRH5H4AfBn4i+/mzRwzzs8AO8Dlg9rhrWFqh3ePJsIj92GaacDUHhg+Y8RpgTo0rb3ot2sWDUGsAACAASURBVKpmgbxtPE8JxJRsCKeMk5TQmILK51QZBEHmWPTbz1PHo5ljc1gWwu0qRWYzCmGcOeHy4J6aAUCVvYlP4rQSB0W/vYwvnlP5nLZzoUTG588+SCw2W5sRmKaWICoL5ho87W9WcbIYEWbWsWLqawOfK7qgSnaYepzTrEwXxdqS1DGIgoL3nh/zau7tnI5ZP3/d57oJH7jTfhhfeciZ29F/AvgZEflR4A3gB4/o+6Kqfu+TTrS0Qrun/D0ZcqUtP3T5azV44XcwSwuhYqFW8DV1PlIw7khCZJ1jmtiawNufzNkeJxX+sKfq7UyTIqovDjzrwSpsDmI2VNmZJXz5/gEPJ95k8M4O3BiFrMZC4mAYGgaBj7icpV5AzVPHAT4yc2oV67yWHKryaDzn4V6ZfjUKDNfXYqzCjgC7MIoNN0cRlRcMQoH3X1/hmbVBcYz25ml2DAzbsznbszmbccTWMMYEhigwjPDJm9YGIVFg2JtYDqaWYWRIbN2stDYMD2UlJE5JnCJk6VazC35/5t96hpGpFB3OBDdl0eHm+auiK9928/cnLex7Wjjfwr6nH6auqr8M/HL2+wPgexbc9J+LyO9W1d94knmXV2gHEGhfweZxkd+EUhPbHj40/MmOZjV6D7xD8OFBUmubpo6DpN7POrCVm1NEuHeQ8mDi++UC30c9Sq3fPC0FYL7q/Y7Cvvd36m+YhTCU0tE4T5VpWhdSt9Zinl2vF+ONTWlozWcahO0oxa2VuCYEnHp6Y7Wf08W0u9TVi0b4bRcrftCFrsyNy5J3BPzD8iIYJEuUmvbfAv5TEfkG3jwigKrqdyyy8dIKbfAnNTA9g+TJ0L5AT/W21cOZIcfyGzoESJcwkkUG68BhLJYWle+QvousbZENl0ZELBkuQmDDUoWx/3sn2XiphbbTXmC3oN4ieRwDoMuJaDokYJMal7dxjPOrW2BLq00zull1rLw6S7VvYpXIdBTu1ePV1fyvpQm+vZ+OdiGANGORHBdg6Cl/9SOlOR3jqHNwWPuCjsVF4FfVPounOcdpY259cM15Lu8o5s35rUE2srTVeycZZ2mF9kVxOZcVBaUrsxWXxV2hJUwoHXqZbxAFotCgOM/2wGuQcWhQpcjWpwrj1BGFhjwNh+Kz2KXWV5sBL2TX4pBpagsec2SEjThkklrm1uHwiZYAbo5iIjGME8s3H47ZHs9ZG4YIwnia8K9+e5vrqzHf+YHrxKHh0XjOl97aZXUQ8eJNH1Y+SxyPJgnDKGAlW9wssWyMIqZzXx9S1dvtdw/m3N5aYWUYogp3dia8857yvttrfl4R7h0kWHfAB66vsBKZjB3iSJwPU8+jEIOMz1295ydzV4Si52ybPJdIrmjkOUO8LTo/P/n5K4XqMBJWopCDmT+WRnxyqSiQorqPUPLWq/eFVP7PUZhltLw+qoWULxqKT8963pS/Jci8+reB78OzRppPWgXev8ggSyu0l+UCWxYcTuWrX4lNm3X9u8+JEQUZr7l6FaeW8bzsO08d1klRBsyI33YyT5jMbBGSPYpDpvO0CCEPgDgI+K2Hu7UseW/sTvjCW3vsVLL67U+nfOOtbcZZId/721O+dneP29dXCxv6g/2Eh5kAznclsSn70xQRinlXBiH74zF747Lw7hv39hnGQY2R8vnXH/H7P3SjKFn2YJLw4O2Eb7u9QvWlLlF4dX3EIDhcqiRWCQNlfRgQ5Fn+AojVvzFUmR3Ny9mHzJsiyRb4ZFdNVk9gSqGdtznniv0+DPnLRt4lp2d2vYFdFM77Fr/ofNmq+n3Zz/edZJyljYjs8WRY9F6Uyg2ctbQZB87TAEst3+d1rkYMakb+bQbvhMbUCuCm1nF/f14rRpsklu3tSW0dzinb+/Na22Sacu/huNY2Sy3TrN5kPud4kjBvOkK1TQ189+GY6TyttT0YJ63Ix2qu8LxfV6Spc/VjaTrCzU0H5azpGMuzKDbRtvY8eTj7wnlLrhi86nH8v3Nbj8ifEpG/IiL/k4g8Vp6SpdW0q+k/e9Rf7Y7K2FblOze185zfW4NzRYBLzrOeZhXRVYEEjDiMCPuztDCPiHrhnVb40qrKNLWMU8vmIGIjDtmfp/zaG4/4rbvjgiMdG7j77h7v3PWmPd2bMlqNWRkNfJbAyZzJxNu+xwdzJpMEY4QoDHjl5S3CKPCFBQTPA7fKnXv7zBIvxOPIsL4xZGUYFgdKVTk4mDObpfzm3gzVR7x0e8SHX7oGxvD2ToLuJtwehbyyNfSUv+mcbWAjDhlFUe2wBSjD2CeRmiYOSWElNgyiRsIpvFkqb8sz+nmzlOIQNON9180eWpxdOJreV4nNyfZ1qZxuh+K8Fd8lMI8AICL/M/BB4H/Pmv4LEfmjqvpnFtl+aYV2FEDYU/4KlDdrpo0dcldW+zXRZXKqMs3ybQuBXdlud1Kn97lKEEk5lnJQKUciIryzM+O37o5rWuxbd/a4d3e3NocxpiL8/c9HDyfYzMRirWJtytwqVlyxi/PU8tZb9ZKkqVWGFYENkCSO2Sz1YeyF3R9UvH6VZZgt8mznGQIBjLRfSFcGYa3Yr6o3bzS15ybnPRAlEtMSvq0HbGvGbqdxUWWoMp7k1J4l5a8EeFv2uT5YlusN498GPqLZ65mIfBpYmLO91OaRnPLXo4qO9+yuXhdxgWpbTDjXruFYj7b0CDrUoK5dWGS3uve9o+aktF+IDzNttObomvf4pS3c67HQsV/LKrDhAgR2BpHjP+eE3wJernx/CfjXi268tJo2eO3D9pS/GrocSUe1dYxA64ZuvE8LC3CtD9EFm62BkVZRXGPycmZlW2odYWNezdgP1X7WKSYo9yFfdXXNbeKhF8ausQ512gpYye3Y1WPpNfMG5a81Q77O46h2XfryY6DT9nHCMc8Zic1D/M9vTuHig2tE5B/gT9Ym8GUR+bXs++8H/vmi4yyt0E5sb9OuopqZzzklDMqCrrlzsGo7zXUtzbbVbFtUi7cXpzCd+3JccXYluIyi5mw5r8u2cc6bRRR4OEl4MJlza2XAKPJVWt7am3F/POW5NW8Tts6HmAdGcDYv2guDYcTa+pCD/Vk2n/g0vNZ5D79CmqaM98aIEaLBwK8tSfjyv/o6L7x8k2u3N33U5DQlCAStluMKDOPxnNFKhGSsjpVBQHh9xO7elCT1NnoV2JkkbK5EvhBB4D9VEWgE5tYRiBBlRQ88BbAtKHfGCaNBUNAAs13Jw9388cQfw6DQ7ry272tLlmN5tki32aQlsi+JDTuHI6P8mZKSeh64aPYI8D+exiBLK7R7gV2Hz/Nc/97UQmeJK6qnkAkN55RJw0adJI7EuoqQUOYpFS6vEBhIUlepYuOrn7y9M+HeeEqSLWZ3njCZO/bntmh7MNnnrfcmvHcwz2zZfryHDyYkGYl5tDZkOIoZH8wxmQ3ZOZhO5uw93GU6LW3ounuAUUuSMT6++sU32Lp3jc2bW4VjFGA0ioiiAJPxqg/GKdc3h4SZEI1jGK1EDALDcOBzgsxTx/29GX/swzdYj+vh6qthUNRjTNXnCNkahFkNyDac+oLGa0Nf2aYZgl5FmjkjjVA8WIwcTu+j0galUF8iO+1j4zwj68/Z/NEJVf2npzFObzG+YOSa7HFoMgSg+4XYtShpGcOjSqtTzdKp1ql2eXBKtc021pY4V2tTVe4fzBknttE2I6mojqrK9v1tnK07KocrcU3wqNNCsNfW0QyNlfZNv7kWMxrW9ZAoNDUNS0TYXItrNvRAhBsrUUsANo09XiOWVltTgeu0xbebWiYsr53Xbeo5Re8yC+dlgcmO71Gfy4Cl1bSvOuVPVWsVtsEXdj3sFS4KDAQllS+s5I5OrSui8KwDiyI4ZoljPC/rPQbG951Wki4ZwRcLyGgkeTReYstq4iLK3DoeTGYkzrERh1hVvvlowm/fmzBJLKqwPvTsiTvbU+bWa/eiljvfuMM3vvJNnPVrufniszzz6vPEg8gzNxT29qbsPDxgMpkTD2NQZTaZ42wKCBpEmECJIsOzL95kuLqSk80ZRiEfef911kZe8O7sz3j74YSt9bjQilPrCAPDzfVBIXjnieV9WyP+jRc3i1za+/OkYMCMU4fgGIUB10cxo7jUnkV8xGM+lqp/uKwOQ6JAvHacmUGM8TnGIffTlAUknFWMaGeF9eo5yrctr5bMLJLZ/S+jUD9va8XlO0LdWFqhfdUpf13Z9vx91+1UqibFD6CljTWP0iRxTOZ17XQyty3H7mReT7fqtJ3Rz9urx7X50lT5/Nv7NW330UHC/jStreTe2+/xtS9+o6DuAQSBMBiUlDwjMN6fMp3OawfDWotWiyEgvPTBFwijsNJN+APf8WxNS9pYi5k3Du/qIGJrta5Nf+vtdX7vsxudzJVyTthqCGzAFzCoafCwsRIWa8phGsLYoC1pddT1XRPGHW9ky/Da/7gwnL8jEi7ng60LSyu0oaT8PV1Jo46+sLqDKzr4Hmf8pPPJl+pvQ0qh/BawzjsXK1aRwuZehXTQAI1ASqPNtC16hznmGgGfC0Updh22oMs8cchpuhjBcLmE0YUIbJYquOa7gb8EvIKXwZmvWi937hF4Sil/T0AF6KT3nfEFakQ6g3WabYExLcpfUXi2asumTTXM6Xi1NwHnCIJ6vsLmu0k1j3a1U5PK17X+rsPWRfk7nGV31vS7Tv7IGc95urgIyh9y+kUQToCfAv4rfOIoe0zfFpbWEZnYq2saATKWR70tT2y0aOXsPA9GMxeIz/2RJ3sqsToIWIlNrW9oypwj+Vhza7FZ1XOnyiRNMVCkWbVO2Z2lbA6jWlmsNLWZGaTM0bGyuc7Nl58niEJMYDCBwTmy6ETFWUs6nTJ9+xsk999Gna9mE0UBz798k63r695BZ4TVtQFGlDDjiakq81nCZ37lG7z38ADnHALcXA35jufXWB8EBd1uMk95b3eGtQ4jPtx/fRAWYf/5vgpS3BSC7ydQy0thpOSaU+k7Sz0jp/oQzcuzVbn0LcflAue4PNf18XPb9mVCTvk77zfoqlP3sM85YUdVf0FV76nqg/yz6MZLq2lfZSck+AsoCiQTjhWKl+ZMkeMvoCYNMN/es/aEQRQQh3mNQX9RDoFBZHm4lxT9wsBT35IsnaofRzlIEvaThFnldef1h1Pe3pkxTnzbIA55b2fCeJYWa7FO2d+fFWHjG7dvsn7zBvPxmJX1NUxgSFJlOj7g3X/56+zdvVdInuHBNi987LsYrg490+PaGrdfuI6gjNaGxToOHuyz8/CgSBD1C/f2+I//2If4Xe+7wSAT6i9sxvyT337Io7HvM0sd+7OUH/nYy7y46bnk4IXgO7uTgrIYGEMs8OzGCuuDivMRz48PK6HpOQc+t6mnzj8II9N1Tstzb7K0ul0RmE3UHdb+/IbSU/4eB8tkHgH+b/GFgP8elcK+qvovF9n4VIS2iHwv8Em8j+wnVfUnTmPcpwFdN52q12rrr/Ja4++q+iAbbbT5ElVtZ1gTrqL5+Z+07BOz1HIwT2u85O1xwu603mZy43AeQKLK/t7Ua9a5DVogHo1qBW5VlYO9Sf3uVW3ZrUUtxiVAKbTjMCCOTC2r32RmfWHcTGiLCJNZirWuVi19fVAWEj4MRoS1Qdg6Py3Ttkj2FtJuPwp5UM1C6DLjXHKBfRFYouP1+7Ofr1XaFPgji2x8YqEtIgHw14E/CrwF/LqI/Jyqfukk4151yl8OL4j9763iq5ntN6+wAnmpJqlr2c6V3Gvyfu3E94LjYGbZnXgNWIBcPhoRhmGAU2VnlvL69ph7B57NMQx9xOEX3trlrUcTVGF1GHJtbcAwDri9ucKtTdg9mPH1Nx/x9hsPmGY5stc3RwxGA+bTFM1qHg5HIWmSsrs9Zvj+3wVpwuzum8SjVQbPvMSDBweE4ZSVgbD/xtd4+I3XEWD1xjVe/dhrvPr+51l5bgPV2zx8NOarrz9gtDbgC3cO+M27E77lmVWe3xrwa7/zkLs7M1SVrfUBH3nlGh99aZ139sa8ezDm1mjIjeHA0xBHAx8cM08YRAE3VweFUK1mR0znLqNFGoZxUCtMkFpHYKTGRpHsXOTnuGrOKH7Xeg7tJsLMBmWLN4GlET4nwkU4I5cBqvqHT7L9aWjaHwN+R1W/DiAi/wfw/cCJhHaUFfZNrrBdG6o3qrZ2VLUjk55tB7xY7SgM6zTzWJSX6vY4ZTyrBMFQd875RPnCb9zbZZKU/O5Jqnzmi3drAmcytzwbBYAUr/6aOr7+1bu1HB+T8Qylbkff3R6TzLMq7mIgGrD66ocJwwDFOxGTxPLwC5/DjfdQ59eyd+893vfyDUYreZpU4eaNVR5UCgxbp3zlnT1+/at1E+F6bPiu920Wbwi+EK/FDbSwZwYCt9eGDCLTeMupnxenFAK7oGKKYsRkgqjcNix8F4Xx/7Gu51qGQNNuu4y4CMqfZ6Jd7HETkR9S1b8pIv91199V9a8sMs5pCO0XgDcr39+iVP8LiMjHgY8DvPzyy80/d8LI00T5a3InTmvMEl0sk062hGuvpCm48tVWt09TH8Ayd9XoyC4GSJveFwaGZgkAcRbn6ic/bOS17oLVdqKp0JjWnIfmpl5AmrRt0W3HcvVvp4HLLqxzXATlD5bi+K1mP9dPMshpCO2uI9GSDqr6KeBTAK+99tpC0kmvMOUvz2B3DHdg4bEO3/7oC9VBFlVXeaXv2KQpBLtG9nzs+gkT6RD40i7s67P3NdbWsRCXsVOO2i+hTfkrj3ettbWtFv8ffdw8bZHGwbpc1LuLQmovJj3rRctsVf1fsp9/+STjnAbl7y18PtgcLwLvnHTQxHpa0FUyjdTpW7nJul3CqvzepocBNZqdqs/VMZ3bWt6R0iRejh8G5XjW+dD0r7y3z539GTbLMzK3jufXY4YZXdA5ZZpYnt0aEgWlw3M6TfjKNx4ynaVY53ze7NBw/dYGQWAKNkUQVDIPZr+srQ9ZXV8pGRWqzKZzbJIUwlVE2PzW72B06xkky64nQcj/8wv/ggf3tn20ZOaMvbY+KPYtELi9NeSjH7hehJSHgbAxjAjFFNkP/UdaeV+k8n+OQSi1vO4iFIWLq8c3p/tVz6X3R+RtWmz/NMPi7+3zVMiE4/OOPE25R34d+JCIvA94G/gPgf/opINeVSdkLX/EIftYcx5m9LDie15AVpX7e0kmrP3fpoljazWqObVUlf1pWti8w8CwN0v5zbu7vHswK0a+vhLywsaQ/cQ7EG+sRnzz4YQ3Hk7Zm1lAuLY+5K27ezx8NCHJKH937h/w7M0Rs8TXbByuDnluNGD74QHGBJgKrW44DAmzLHwrqwNW1ga88837aLYDaZISxiGb1zcIwgCRIStbrxHdfZe9t9/GBTEPd2f80j/8Vb7tOz/Eyx980VcvDwwbawOeXR9wa2PAxsibUL71xU0i63jp+ogba4NiHYMwYDUKfamyDGvDwCeXqhy33KkoIsSQOVIpUq/mynbz3FWdy6reDxEYySIrc9NS5ZxeEmFx2nDq6WbnArk6D8sTa9qqmgL/JfCLwJeBn1HVL5503MuOw4rAPgkOC113HXUgmywEzVglzbVtz5KawJmkjt1ZvaRYYh17s3oguU1dIbBzbO9NmVT6iQiDQVwI7LK9fuOISCGwq4gbeT6i0Yjh1lZt49kswVbs3SLCjfWYzVFUa3v11irXVuNa28YgqglsyIKdamttv+UYI4XAfhJUr4WnVVBfJE4juEZEhiLyayLyBRH5ooj85az9uoj8koh8Nft57az241R42qr688DPn8ZYOZrOq8uCwvRRaatqytWMbdWAC4VOIdwqxpvxtRW4uR6jCruTlHnqWF8JC3qYz+bnamwRp8qjyYx74xkvbsbMrXLvICEQz11OnGJUmaeOL769x2+9u49VWMkK2L73YMyDRxOgdERGocEnxVOCIOOKW0UCQxj4lKrJPGU+nbG/4/nSaxsrOOfYeeAL++a2/a3ra9x8dtPXi3TKZJKyvh6z8eoW6MvsvrfNW197mw988Hle/dALRT5stcq3PjvimY0BIjBNlcncEQbCm3szzP6Mm6OYF9dXuDGKi4x+OTNnEBqSbB/CoHzw+drBvgBEGEjGFjn83Au+X1ChAea8+byIb/X8Fdvp5c3UdxKc594KPofMKWAG/BFV3ReRCPhnIvILwJ8CPqOqPyEinwA+AfzFzrWI/Bjw08Ae8JPAdwKfUNV/vMgCljYiMg78RT9/7Mj8i0XXK7OSC+tSRBtKYeX/V9LjBHZlPMg1B59dzjacifNU2Z/WD97d/SkPp7MiIGcQCjdHYc0U5RR+8TfvczBLi/ZpYnnzzZ2asBFpB+1Y60hrOyGos0wPJsV21joevbeLTW3JV1Z49oVrbF1fK7RpCYRnn10jisoX6K1nrvMt3/K850Zn/QZRwMc+sF6jcw1CSF353WUmittrw1a/pobVzE+i+FD2Kr3vMAyiMqgHwOTUmWNKlT1t7suLyvJ3Gow/9a9L+9nXKPsonub8h7L2TwO/zCFCG/jPVPWTIvLvAreAH8EL8YWE9tLmHgGyHBEXvYrHQ9fLQdu5lQsAafV6InRoaV3rsOo625tIrKtn78ucps3gH5F6NcZCRlXndG1+bBflL47D1l0cBO1jYioCG+gOT9b20QyMtIoalOeh3tZcx+NowbV+0vqlBxdH+TNy/GcRiEggIp8H7gG/pKq/CjyjqncAsp+3jxoi+/nHgZ9W1S/wGBfJUovEq0L508r/ZzZBm+O28KYtnOCG6uJzN4OBugSgde0T3bkLi+yWHKbNPqEteuGJeywCz7I53znzB+8CNu2bIvLZyufjzbFU1arqR/FMuY+JyEceczmfE5F/jBfavygi68DCkm5pzSOppWUuuAzotMV3UXo7EBrJaip6iAihUGurDZtd+bOsluPKwJsSpOxQy03y7NoKADuzpKD4vbM9YxgZ1ochzimJdZhZynzuiKLA5wpRZWtryP6eL8SLCMM4YGM9ZnvXt1mrpPOE8XjGytoKgfGW/NXVAcbA/u64oNgNBjE2tMym8yyoBjZXQjbWYnYP5p7BIXBzNcIi7M/SImJxFAkqQmIVI8IgENQJYVj6CUZhwK2VAW/vTYtCxKrKLLWsRPU6kE1ziLdJ15UFH54eIJV8MIGQUf7KfklWHSdnj3g+evutovtcAnl4Uc2UdvWQqr+/4+B8kzgtONd7qvra8d1AVbdF5JeB7wXuishzqnpHRJ7Da+GH4UeBjwJfV9WxiNzAm0gWwvIK7UsosIFCuOTZ34TFX62NESIpBUF+44eqzJK2aWN7nDBNXFHcdm+SEkeGWeJqOaxjkzMfQt5/fZ1396b8/Ffu8tb2tDCDuFnC/u6EL31zpwhD37o+YnNjiHVKFAVsXVthEBpWVyLWVn1tx+efUT77uW/y4ME+ScY+2Xm4yysffJ54EGKMYTiK2dxa5d67OwRBUCSNWl0f8eFXr/HKS9e8eQRIkhSXOp67vlKYVQ6mKdYpN9bijGqnRMawOQi4XWlbCQPWBhEbA88g+cD1Ne4dzFiNQtYH2fhWGYTSMnmsDoJaNRrrfLm1nPLn1AvqvMRYvm1iXSVsXzy9TygSVBmpm2WqVMyu8mFGnvyN4LLBOloBVWeJ0zDJiMgtIMkE9grw7wD/PfBzwA8DP5H9/NkjhvklVf2e/IuqPhCRnwG+54htCiyt0L7syKP+qqjmVD6szf/sblOtb9ssH6ZQY4vk46cKAyktYaEx7EzSmt16Z5Ly+p39Wt6QZG5r30WE0TBidRTX1ibOkSYVaqC2H1YmMERxWNM6jRGef36zENgAK4OQG1th7V1xdRAQGlObc3MYcrtSPkxEWI1D1ipjGRGeWR3QRJeNOhfk1W3DLpt6w3FpRNBGt5bl/BDK5rnbCJ5iCBCejiH9OeDT4hPlGTzF+R+KyK8APyMiPwq8Afxgaw0iQ2CEN8Fco3yd2gCeX3QBSyu0Lyvlrwt55J4tNGivjbmKc8+Iz/SWOi0EcWC8Ay1JlZwaLapY67zQbVjBjFAUkE2sY3eW8u7BhKl1rEYBL22sshoFrEYBH//YK7y5M+GffO097h/MuXV9xM1rIx49GvP1Nx8RDyM2NsuoxSR1OOt498EYeTjm+saQzXVv+viOj77Ct1vL73z1Lm+9+ZB4EHH3ziPC0LBxbY3R6gAUNrZWUKccHMyxqX8b+Cf/3+tc3xzybR+8yYdf3ODVrQFRIEwSx7v7cx4cJOxOs8LBg4CbaxHXVkIGEewmKbERBoGnJD6YzHk4mbM19Nq2No7vIDSMsihJVa91qypRaJimiqDEoSE0EIaGGL99Yp0PsAm84zUvzJtzuU1WyDfPW+6V9eM5IWXxiSrDaLFtrwLOew9PQ2ar6r/GU/Sa7Q84XlP+z4E/hxfQn6M8BLv4TKkLYWmFdmSyG+YKOCJbSey1bf5x2i72m1pl1uiYWMf93TlNNINqducpb+weFBrrQWLZmc6JZOCFjRFe3lphayVie5oW9uAbN1aZqbfjli/9MJvZmpDZ2Z+xvhaDE8QIkQl56aUb3H9wUNjw09QxOZgyGEalsycQ4jhgkrrimDzYmXJrYPjA9ZKSN4oDHh6k7E7LiM+9meXbnhmxEpfv1IlThpVXbAV2ZwmjsB6gExhhs6aZQ5zZjguTRdY3DEqtPs/e57dp01SKfmjxwHx8pokWdvCnRVhH52zP9g/Yiz22qvpJ4JMi8mdV9a896ThLK7RFvA1RuQJZ/k7xjUG1DGypoikoUuda7ugm11jEa7Tt1KN1cpwW7e22ar/UeUdcUj1h4sO3q2txHa9QG6tRixqYuLbdN2qYLLopltLSVY2027rNFs32BbP3dZjDFsXTFlQTX1iWv/Ofswuq+tdE5N8EXqUig1X1byyy/dIKbci4wUstsKu66PnP+tjbaZstcZqLzUgc8gAAIABJREFUEKRVxBfaXKZOyl9mqjhubYvue1snXvwsndox6tEJz92/nME1pwER+d+ADwCfpyzsq8DlFtrLTPlrFlb1VLC8RVp9qt+PEgZVO/7h2/uw61lazbtR76fAShgQB8LcllmqD5KUjTgiMLk9XXlxc8jd/VmR6Ajguesr3N+ekuQmEoXRSsh0ZgsaoaLsHcxZXy2dkpsbQ156YYs3397GZG9JW6OIW1sDHuzMAa85j1YjDDCZJoUG/IU3dnjxmTU2hmGhcX/kmRFfujvmYG59zAs+ZP/GWlxm6hMhNgZb0Bu9KSQMhUpa704WT26SqD5nUuuIAskiVg8/V77GRHlOngazxmkhVbD2fE0kwsUXQajgNeDb9QkTEy2v0F5SgQ3tUHWF1k2uZFVmGpqnoa3F5c4+MYa8Wo1rzJGkmjkpveMsDIRZ4kDKorKqyr2DGftzS+IcK0HIwCixMazFPqvdzDr2Jym/dX/Ml9/bZ5b6fCAB3nwShQaRmGvrQ959NGF3PM8qjwurq8p4PPeaUiDMUsdsZ8qtzaFPwSrChz70DK++epPJ7phnbq5z49oIgNnc8i++/C5B4LP8jUYDrHWsxQHP3l5jEId89q19nluP+ejza0TGMFoLub064O7+nP15yjPrsY9sVGUUeqbIehQVlD8JYBSHbAzDos05ZRAFhR1c1R/HILPrg394lU5EwTqwtKmBeU3OHEVGP+hl9mPCO3TPkfIny6NpA78JPAvceZKNl1ZoLzO6niddr9RdpoIuVLVAL2xsESaeI22Elud5R2rjiLA9nVOtAGNECoGdI7GOL97bK6qPg6ffjeKgeFiKCCtxwGRuCk1UxDsRk+YTVep5SKIo4KX33ySozDmIAwZxWNNqoyjglRc2a68K48RlZbvK43FzNWJrJSjMLCLCIDBsxBG5tMxpgCtxUNt2fSWorUNEWnZxI4KYtmbdtG/3LL3LjSXiv98EviQiv0a9GvufXGTjpRXaXc62ZUE1U181SEKcw2SUL1XvNHPqNbT8tR3aUed5EV9BmSaOSUaZCY2PzpunWhPYqXVsTxIO5imBCOvDCFXlzv6Eg9T6hDzGkDhld57y5t6E9Sjk2dUB48Tx5fsH3FiNmKWOvZnlla0hr72wyUoY8I3tMf/67j7392bsTpKMl5zR2QTCQcTKwGvOqXVc3xiwMfKpT+cZjS8KhAf7CaERVochqsrD/TlbG0OcU/bHCesrId/6wjrrK34dD8cpN0YRH3lmjapMHYUBG4MhRoSd2Zzt2ZytYczz6yOMwCwpo0GjwGBdRqk0sBIFRUa/rkRevr1ybDULjDFSbNd17vPzVzA+lkYWXC6c52ETlkrT/ksn2XhphXae5W8ZKX+F9qWuxW7wBXVLTc9TxupUMEE7izxsZ1XScyTWMa2nuCa1jrd3JsV3q8q9/QkPp/OS5YGn/O0npVF3L0l5590Z47krtPhhFPBdL21xezUuNPEPXBvx/37tIbsH8yKYx2TmG2PKMl6jQcj1jUFNE41DU0Rngqc6PtibcZBVZhcRgkD4wHNrvHRjVDzEhlHA73thhZujuGZ3vDGMGYal5rw1jHl+Y4RUTELD2LC+EtbPC7CZFf8tz1U7XZTryKJ4VM7s6vnTRluPx8N5h7DD8ghtVf2nJ9l+aYV2Qfnr4DQvDzreB1oXRhdlrPs9YpFXb9uRj9l20ACbZbTytqaQWo2CmulERGppWfN1ebNNYxea66djHdlbRnXeOAxaTqFhaFpt1ShIP6W0+OhV80iOfG9q/RZ8dWsX7O2CIPI4fJQeVQye3sK+AIjIHuXVGOPTux6o6sYi2y+t0IbDX2mXH80b+uxv8EWOU9dDodNkcMrruGx24LZ/ouv8LYcAuIzI87icJ7x58nznPAyqWqvGLiI/AHxs0e2XZDfaSF1W/HOJb/hGamcg50JTKzdmHUVx12oJMm20xR1XcvVCU1XCzAlXRWyEuNKWa9nVvNUC3BiFrEamMHcIcGdvinWuWINT5WOvbBbFbLtFlZJa5VGW4S9vzxMq5d+N+AINN9aiIl9xIOCcIzTlGoz4Qs5hFqCSf1J1LTrjpFHAuLq/5fcyRP2oJ0aXOTp1rsgKWB4TKudtiS/IS4LEwSw9/zQVy1rYV1X/L+CPLNp/aTXtyxAFmTMRVL0Qq7INNKOR1Yv0eiFRL+7rv1j1meGGRotUrLl5QFXZPkiYpT6j3zAKiUNlZzov+MnrccTIKa9vjzlIHPPsaWdEeXYtYjUOiYzh+XXYm6XszlJurcbEgeHRzFdBn1vlwXROPDD8wQ9f54tv7/PO9owqxpOEySxlktmp37q3z+96/3VWBmGZEc8pa7GwtRKyPgwyk4vl3s6M57YGhb3ZOmUjinhxY4XV2DssE+dQlGujAVHgKZDTxFebV8BZJZmkrA4C4rCecS80PgQ9zHKEzK0SmdKkkmdfrLJh8jYRf7xzu33T1JP/LCsO9TgJzpvyt0yOSBH5U5WvBs/bXvgRtrRC+zLB5/Jot7frPbbPTG5nrtLU8mo91bbUupqTLw9KqY4XGGF3ZlvlwzZinyI1x/ogZC2u55WepDZzZkox/uZKyN3dWS0x1TyxhcCGzPGqSljR9I0Rntsc1Kh1q4OADz+7WndcBoYPXl8tKHkiwjAMWB0ExXpFpDMowik1gQ1eYEfHlDrqyq6XP3ybNs/DTD1N7b8Yt8dSY4lO0b9f+T0FXseXK1sISyu0l5nyV8VhdRyhpPzlfO0i4RJe2FmnNY04pwrmQjcQ72EczyypKyMmBYgjw6srq8yt472DGdPEMgwC/uDLN9iZJnz10ZhJatkcBuyllsg4BoFXaxJrcXhzRCDCnd05n3tzh71ZykvXhrx4fYiIsDWK+X2vRLz5cMK9vRm31ge89uoW1jp+45vbvPtowre/uMF3vLKBdcq9/ZSZVV7aiHlmNcKqcpA4QhFe2ljhxnDA3jzl7f0JiXWsxyF3DqYZrS8iCrwzcpYqgbG+FmT2EBkNAuapz3Edh4ZB5PnjeVHcldgLe6/BaUWot89LlbIJ3ZFy1XNVS/aVjZ/T/QqtXZ/O4ryXB4JZEj+Eqv7ISbZfWqEdB/71aXmZI0cL7GqWOJCac0tQpqlracTNsVKnPNovOX8ivpLNaBAUodODMODm6oDdcan9XluJ+BYz4t2DaWG0TZxiXVoTUKkqv/L6Nm9vzwqTzFuPpqzEAetZZr4wEF69OeLVm6Nyv0LD7/vgdTay9HrGeG34pa2IrWFU2AZDEZ4Zhby4vpqZhoSNQYiRFbZnc/LFjVPLIDS1qjK5HyBfrogwiAJWB/UoRQVurMWFCUSAPHW4P0SHU/cCNIv2LHpXqs74b4fRM7uEuUBbFe/Ricicc5Y/lue0iMiLwF8Dvht/6fwz4MdU9a1Ftl9aR6QIhAGES3Kgu7DY86RNSQNZyAmjlU+xZTGM1DrWD5P4uoyNY9d1k+xPba2cmVVf4bw+Z62GHlDa26uRkIFp85vjIMgcfnl7trbG4gZBe9uuU9+lzTZpgMVaj7lLy6jSSr9irHK9j4VlkQxLjEFwAQmjxAerHfc5J/w0vtLN88ALwD/I2hbC0gptuGjK3xKr+B1YiPK36HYL73qXCvpkx+1yHe0eJ8FFFDepvkEd9Tkn3FLVn1bVNPv8r8CtRTdeWqFtL4TyV6XgUaF4NReR08kWG08q3/O2ZimrnG5XpQOiStQ4QzYrvlusLRurqiWoKqMgIDKlfquqzKzLCjKUc3zo9gpx4J2fOSXPWkcgeb3COnIvvAgM8yK2QJCZJoZhUGxjBBJ1NZOMqjIMgxpPV4BpB12oMHdU+kWB1L4LPnL0SVCSFSttxXkvsagC1rVtjzZ6yh/vicgPiUiQfX4IeLDoxktr076I8PU8d3f1WupK0ekOof/mr+m50LeNfrnd2jrNGAueT1zNHudpbQ6rSpIqYgyR+MK+OdMksRYjltWhP32qsDoMsU65vz9lbh2JUzbjmJm13NmfMU5cURlnJTKsD3zGv9Ew5A98cIudcYpa+N3PrbM2CHGqfP3RhO0KUwRgPQ7YHIZFyHnqHJPUsR5H3F4dYESYWcvePOHaKOb6ik/dmljH3Z2pz1Soymrk5wgDw9aKD1fPBZ4xsBKHRUY/57zteWsUEQYG9/+3d7ahsl1nHf//136ZOeee+5amaVOTNBUtKLUgXKtQxGqqRC2tiIqKEq0QFBULFm2bD34qFEWrUEEuWoxY1KLVFqFoWhQRqTaJrW2N72htjU3Tpia5956Z2Xs9flhr7b32y5zZ58ye2XvmrF84OXfW2bP2s1/mmbWf9V/PI4Lnbi2K41beOQfKEIprc/lgXJspF8bKdi77wLLl6y4Nq79d2z4D3di25A8YVeTqjQDeBeCdMKfib2xbJ9Zy2iR/CUa+Mgfw7wB+VES+tE6fQ9L1i/8khw243wJp+eLxM/+ZvB6CvLZNlmtUi7+YfnXtCyDPpRZTJm5leWW7RCn833F1D8eZrsQUFYm7rkzxkqNpMdpQNJOGz86ySn9XpjHuOCqL5cZK4aWXp0ij8tM3iSK88LZJNfNfpDD3S5iRZqLyaFo5b6qYaC3bDlKFo2lZPkyRuHiQVK7DsqyHUS17n5tcbWy34gNdTiyfvM/AOCHGE1YQkU8D6JTRr411j+MRAK8QkVcC+BcAb12zv4IhPgtu1Z6PGxm7sEmWm0x2ua6vyiu30yJFxjvXlmvBzLZpLdDaZKd77laGm7PcFM7V5u8CVkITisDRNMalg7hwONPEyN5iGwLJteCZWzNE1lEBwEJrfOHWDEcTYhq78IXCvZcPcO/lQ9w2MSsVE0UcxArPzzMcZ7mxLc+hlOAll1IcpSYMcvflA7zixZdx58UDHCSlk37m1hxfuHmMWZab0I8iFrk5By7ks8g07rx8gNvsCD1SxFX77/LLA5jU4kEu494iN+fenaO5ze5Xvw4+1bBU4FzD8YRHSD5M8or3+irJd3d9/1ojbRH5c+/lRwB8zzr9+aSRWRW5zZh2MbqTZka4+upGoKxe4o/kFjUpHwDMsmY2wBvHeaWQ7yIXm1a0tIWsjuAViIvTukoCmOc5/udZm/mPZrsbWYYnbxwX2yUx8IKDBLcdpMXNeWWa4iBWyKwcUQDMc40bXmwqVsQLDhN840svI43KpE5XD1IcL24VEf+FFjw/X+Allw+9ogyo5N5WNDm7ncMvZXVGKZTUsutNbbEH/0nDnDNP7aLLgro+vlQwMD6GkPwNXdjX45V+REJEniHZqPC+jD6fGN4I4IN9dUaackRDSP66Xls3I+3TqtBoaexaIKGxzxb5XWYz6fk9LlwODa9tGqvGjSu1WOwymWFbFr46qtOn0Le/GtppSPla2k5DcNjjJFXD1Ihkh58toUheLewib8MpBtArNyT5IZjSOHUeEpH3220eglmO+Z4T+nkQwIMAcM8993S1b/wPtk3J8SB0ES2sJ+/r9a0jon4BR3JB95ih7psRfYf/MoC/IfmHMKfj+wC8veubVzptEXntSX8n+QCA1wG476RClSJyHcB1ALh27drK62YehcftGIoRqQsv2KXMpdFSVDipH8c0Ubg110U/QPtN5ffv+nOL7lwejUmsbNxXFy7nYhrjuXmGudaF4mGe61ICaH9PIoXjvGpHQiJzskDb9sVbM9x2MDGLT2AK3x7EEW5lZpKTsKEhLa2lu06CMIUb4tp5yvJq6TG3bVOAeVo3K7UvufL6lT2W4ZtAvyxs2DPZ6mh7vae2PhGR3yH5KExmPwL4bhH5x67vX1c9cj+AnwfwTSJyc52+6gzpsFdlhKurOAAAhc7avZRyctL1C/NYmMYRLk6N5O1LN+aYLTRiG7oQEcy9OLCTg4vWEEGxepE08WazD+CFF6aYZTmemy0wiSKkkcKdRwd45niOz92Y4UUXpriUmst9nOWY5RrTOEJiiwk/O5vjWGscxqamohbBjUUGLYKjJME8Ezz1/DEuT9Niv5enKS5ojZuLDLFSNlOfKWgcR+2Oux5BSWOFC5MYk0TZYy/zi7TmBFGlcyXKOLg73y5kddIHVDuJu3/NG448ZPTbJFrMnFWyzSx/29lVJ6yT7uyofdbVab8LwATAI/ZD8hER+fE1+wQwjhH2soxwbbGIthh1VmtTNKW1fOna4SSGL/pzMsB6d/WJUCkm5UomceQtFzdcnaa4PEkrbdM4qkj0SOIoTZDqUpKnSFxMk8r7tACzLIdKytsmVgpXpmkjLt5GY7k5gKsXkoq8L40V4hUfZCcN9OPgXWR7hX0tBoY1MftPHxORJO8G8DswIWMN4LqI/JqNS/8BgHthsvZ9n4g8s/YOW1hXPfIVfRlSp+0xuA/cKM05QJdZj6y2mW3b+/CLuwJu0rSa0Q8A0ohG9mfDE2ktbWiuxY5KWWT9E2nm286sxA0AoshIkyaJwkEaIcs1bs1zozyBkcxp2y9gRqJJrDDPNOaLMnzijtsoPzSenS2w0BqTSCGNFGKlbLEFYp7nyERwmJrkVCRwPDerK52GOteCGzaPSRorTBMzWl/kZidpojCxNSTnmSlucGFS9c4uL7mzkVam5WSB9S8up94x14CNvrQ997QOvtSgN6+fu9/8tsAeseLp6xRkAH5WRB4neRHAYyQfAfAjAD4sIu8g+RYAb4GJQvTOaFdEbkry1ybbA5orF0/CLK/2i7sCABFBINb5uxskjoCUZU5o135rXiZqUiSoYFc9lvsREdyc51XbtOAFl9IiBBApE0Z43ltAoyIiVlJRX0wTBYrgODNezDnFLx3PKgWAZ7nGhSSuFNSdMsLVoxRxVDq+C1Papwa7YEUJkgssRv9mFSERR1JUNieJlMA0jTGJyzbAxK/rssjIJvGhPcmRAiioZEx08XofcV8WxWsTskkiT32DatgKS65pYDPEarslx/oKj4jIkwCetP9+juQTMEmf3gDgNXazhwH8Jc6b03aSP+ZjTc/ashquOgtZbNdWLLYeTiHbM/814qx21U3ZH1uL+NZlgQBbvwDntbwdApuZzw9jsFkkwGnEi7YiTFGzl7VMemSrvK/tEkds7hNSr9/YjeZbnF3+NGYo2LsNUoXWoiGbpu+JSJL3AvhaAH8L4EXWoUNEniR5R6878xit03as56/38wO4+SPaz/PWTuObdxArApun45W93So7HNet8q3aF3kE4I8AvElEnt2mMmW0TvvMkj+rsSsHn6XSoBR6yEYef6sZ/cq26kpHY1sSsaha40giNiYX44iVMmN5bmLffk6N2CZEWhXWccvLfQ6TGIuZKbTgXPVMaxzWhkKzTGNSW7FYXxHadjadUsP/Y64FsarGohWroTCXMiChC6Ms24Mviyzb6s88YaJxPMw1EIkJkWzL1xHmya0DT4vItRP7IhMYh/0eEXmfbf4cyTvtKPtOAE+tZfAJjNZpn1Xy5yaT6jpcgJVsekbOdfYRpYsJO7Jciix8gHE6zpZcBAttisxGyiwXjyOFSNnsfbaPaRoh1WaiDjDSwEvKZMh77uYCJHEwiUyWwVyQxG4JOHHlQoJFLrg5yxArFrUSsyIHiCkcfEER84XJAggAR2mCwyTGzUUGDcEdF6c4TGJoLbg5z6C1mUDV2kw+pkkpxXMTimmkivh5GptwjdaCiGVyJkH5paIFuLUwKVvTqIzNK5qJyNzJJXPBLM+tDp3FhGJddukmVJ1EzyzbN3a40mBtIZnAcOQCYIuSP6CfLwiam+i3ADwhIr/i/ekDAB4A8A77+/3r762d0TrtdQZG9VGVUWRIo+2sF7EthWemmyn92sqH+XFgl3HOH3ErRUzT6p2cRAqXDpPGRB1RXTqexkT9ksYRcFxWLLMOTWEx9yYuSVw9SKuZ9BRxkMaV2LvATIRW82M3r5WyTrMu76tfg1xLZTbKZeVbZM3tktp2bbpqc07K7RTZmpckcB5hQw57Rl4N4IcBfILkx2zb22Cc9XtJ/hiATwP43j521sZonXY9nOGcgxvhnjRq8uVs5jXt43e5QtHJyHxH6EvB2qSBpH/Zq6P0SaKQ51LRZke2f1+RYBaAeGEUEoky9RqN9lrj5sw41MNJVCgvlCIo5QQmYb4ElEgxiiRMjpDcU09MEyPJu3Gc4/njzEw0JgrTVOHmLMfxQiOJiDsuTXCQRrgxyzHLtNGUW0nezXmOLBekMXHpILGSv7Lqu5HiSVGIIdeCY21ep3bRkKKxLdPlOUoj++ThXytlkkrNM118ySUtjtfdB/71K8f0wVEHmvQx0haRv8byG+y+9fewmtE6bSf5y+qLSmCcaj27nqNsKx1jIemqf8BZlQG5EZyWZjFXLSiqyLSv9CMYmf1mutwmop8o37oV35EDoCISAb7w/KIIjQDAszczXDmMEUU20RObsVoXKpkmpUqDNmbuS+0uHRKTmLgxzws7Lh4Qd1yKcDCxdRxJXIpM+CT34v6XolI94o78wiTCjbkurBFBEYZxLHKBFo3LB3FhR0LjrN0TglueX1wru89poqq6jto5dy/r128cy7ICq4hoYtrbgsBoqrGvy2idtpP8AcA8b/5tVXyy/ABzRRtqbb4UDLXtTt4vyRaX0S1RPsmKwwaa4Q9nR0MGSP8dpq96alJFFilY/X2mSTXznyKhgcZ2ddliacLJBxdXnHOLVBDLrxVPHDW7/pw1la/owIgZRPLH7U16bprROm3HeuOmtqvU95ULj+P9sM61Cuc/sJoR5dNeizHlUKmgxYyw22q2+jHns2O0aKuL+Jb7rL6uv9dwlhvD9VVf5i5Aa2WW+h5cUio/S52b9POL+Cb1UTvQGN2LSGty+vrxE+0r2upNWcsFdOcsFME9n8y1VYdt8fITZWWqk352gdGOtJ3kzyyHLuVzXSYiu1Bkz/MnIr2+k8jK9exkmC8Z07oq7zPJi8x7nZRvkWtoXRaVBVDk1fBvjoXVXbvJxyQicjFyuWkaFSslY5pHyiQy8edcC24tXDkv4MYsR6xos+WZY8u0qaoeKRN6meelQ1YkLkyi4ovC5fVY1JaTK5aSvqKQLY2iZZJExcSpSYYVgzTHdGueI4nKY6jjMvJhydxEYL8ZRPK3J09ko3Xa/pdwMQHW4zmvO2zXRi8uTLZXY6ln79MC2AJaxXuTSEHXnmPiqCkDzLWpE1kqTEy9xqQ26o69iUX3Wmp5STItiHXVCWbWuTtcvPvKYVJzlqaGZWWfqrp8nSSUqj5NJJHC0bR6jtKYjaeGNkJgKbBN9mVsMNrwSF/n1xXVrYcZ3Ii9um17OKJO7I2eAX/kXIYi3D51LQxQ/w6YxApXj5JC1tb2ReGK/NZtO5rGOEij4lwtKwdWt5c0Mr5FrotseDdmGWaLvNgHW/oTEeS5FO8DUCxg6XLe6uzJZyiwI7DDf7vAaEfaiZX8nbGUIgBUVkACVY20VdBVdNT+dsskhUCpmfaf7Ny2udaVOLy/nN7J2ZS3xygmYiHSiyluHGdWm1zudxIrqxEvbXPSw0gpxEowiYljO5puSuPsE0OLXvx4oZFl5epIwMS4D9MI06Saha+S4EqaBXa7nDefUHj3fDOI5G9PbrXROm0TSzWx7bOmZ217m3HYvrSsubKuC74e3B8ztvXVLnGr/R2wuZ9bRtm1lYV+HyTt6s4VdyTbxxFZi8H1HCN908ecRGB3GUbyx71Rj4zWaTv24zQHAoGh2RdfMlqnLbJ+EQSiOdoW2/m6X7rV+K0NCSwZspvmMjxSia17bbFiY5KzniOlEPVJWey3bUlPVyKyMdpe5BopVcNev/jAOhh1S1CNnFeGyvIXRtobZt5DYV8nU1uVEc7FmOt5RpbhZHYuzuskfy5MoeAVBPb247TgjWXsdvuDNDJ5OxZO/qdaJyV9uZzYXtIIRWmzNpz0D7AZCe2+p1YiOMvM5OLhJC4UI9qmS3Q5W1w/hCl7RpSTt2C5+rFNmePoS7IZ2G2GkfztB6N12n3p7utOtOjfG8GWeUK69ZnVPKMW8z8/zrysrzY7/KZImTJezeOojhTqjtFlDMyz5mKWtFZ1pi49JE1yqriWvtT/YvK39aWHbWoXE59vXsFVX4aBwEbZk1tvtE67LbSxVn+sO7nV76nnYwaWjyCrypSyBmPdSZ8QRSlw+/Iz+kV2trJxDDWnHytWFv4ARukRqarTbCua0ChqYEfnlfSsAmSZNiNtlk8y9RF0m22BwJCE8MiGSaKzF0Ko49QZUmtbhl8VHPZ9ukNw3a0WdKPsRqjEs6NNyuj+7hyeisqUSTa83DwGAlqX7ZEilJgvHDfoFpS5vBOv4o0vPfQlhbQxf2X7V3ZVpgvvaBA6FyiWDt6N/J3Td2lu/dMcRtkBx7Ylf8DeDLTH67QVgUm8nuSvxBv9dVyHd9Zd1rPhVSco/VFocw9ulO5eFM7a36bF9Lps0ckA60fRWii3haZEETh9QV0/u2G3cx44HwxV2HdfbsHROm1H/+d5T65cIBDojBkA7cdnf7ROuw/J3xC4pEr1BTFnOYyqNPB0N5wCUEtDDi0CVbGtWaDBtOIMI+s2vP4h7Y8JgXPHXAOxlEnWtgL35/YbrdPuQ/J3VlyWP20ldG5Cjiiz4SlWJwz9tsaqxhbpYb1AbVksoJlxzy+A64QavuokUkRUsyNJFFIA80WOhTZtTh1S6q7LPlwsuoyrl8egCJii6EZH7iR/LurhHLM/0ekXNjbvtBO1LjgfONdkAkiQ/J2J0TrtoQfYJBFFRP2eilu0fG1tbf3VpYckTcrVYvbRkGld1HgEbKKoiIVjB9pHDXU7CGCSRFAtjysN6SGsBrv2hKCss3btaSXzVO33iv65bONAYOM0B1O7Si/TASTfTFJI3t5Hf8A4PtouW58+ZQY79z4zUpairUvyK7fqMKo5aFXTCrbZpsXl8XYFGlyoptp/26GcNZnOsnPkRuyBwFggV//sAmuPtEneDeBbYcp5v/cKAAALEElEQVTG90afkr/T0lhFiaqc7SS0lQs6stwsM29bgFI4NrptdaF5jpVReyjljbDZbpsrylDs0/7bOX5FQiB2u6q9ZEsdR6x2uqvOUbvk74QOA+cKtWXJX5sSa1fp47S9E8DPoWf/6iR/XVcpjgXdMpxul+lZJ8nydqq+1TzOuTg0ebrbzkaWvf21ZyhxMsC6xK9u2+kp+y0lhDt2MQMbIVWwczRb3jE7/OwAazltkq8H8FkR+XiHbR8k+SjJRz//+c9338c6BgZGQriKgeE5N0UQSH4IwItb/vQQgLcB+LYuOxKR6wCuA8C1a9dWjsp3VfLXdtnbJHT9yeraceGKVd0vO719ZfQLBOrMNZDAqY22t999uZ1XOm0ReW1bO8mvAfAyAB+3H+67ADxO8lUi8r/rGja05E/VsvF1TXYURQqqtgy+LqED/GyDpQwuiWi/rMw745bZwTbb6pK/Zfamscno56SEsWKhOPFDM0UfJ0j01jlHgcBCm9Dn1iR/OzTRuIozT0SKyCcA3OFek/xPANdE5Oke7BqF5O+sF9lk3OtWKk2knIx0+0xXSDmW2RZHbGTlA6rOlAQmqvlJWZUJ8TR2BAJjZFfCH6sYbWHffaDu0NjSVscV2j2tzHDZPrf13kBgzBSKqB4kfyTfTfIpkp/02m4j+QjJf7W/r27oUPpz2iJyb1+jbMBmo+urs4FQMKNc90OWbfTa3IE6SZ7JpQcvs94p9snmT5ebkVhiWyCwAQizjH3b++xJPPLbAO6vtb0FwIdF5CsBfNi+3gijHWlHVha0a5I/h5PM+ZI3v035bQ2B3lp7PpvUboVtgUBfJMrIebdeHb0nry0ifwXgi7XmNwB42P77YQDftba9SxjtMnZg9QKPzdJnOtGhjqLLftuOMzjqwObYurMu9ttpx7eTfNR7fd0q31bxIhF5EgBE5EmSd6x6w1kZrdM2pa5MYpnt79vsdHOZ6dxS8qFldWOxI3CemOVmtL11yV+3zZ4WkWubtWQ9Rhsemefbd9h+vg4tgEZZraXP/rXrX6r9q9qTxabu56V2nGHiMxA4CwsNtJQz3Sw9BrVb+BzJOwHA/n5qrd5OYLROewj3UTgwZ0NhRH/uU9d8o/+dQJtrpJhEbEnzuik7/OMOBPYN45M3uiLyAwAesP9+AMD717V5GaMNj5wFP090OQnXHbf1Jp0XV/RfzcK3ubDFKjsCgb2ix1AMyd8D8BqY+PdnAPwCgHcAeC/JH4NJnve9/eytyWiddqxO9/jUyDpn13GfZpWeUuViFzc9V+o6+rnibnK10qMn+fMdqQhAaWYI3JQdIawd2BZDSf76QER+YMmf7utpFycyaqcdsXv+kbbVh6dfVl2u8GOv6hHbp59Jr6X/tsPchCNdZUcgsEkStX2H7aSw+8BonTYwtORvPy7was7LcQbGwlCSvz3x2SOeiOxB8ucy6Y0Pm2hplLYFAptllnfLy9MnXYQju+LTRzvSPm2WP/ftXanBCIzq67WYKIWXKMrLpKdYLbY7HssDgX6Z50BMIN5iYd99+UCN1mmf9ovYxauUfedY41f1El1oyfJXLO4Z6TEEAn2w7efMfcnyN1qnfVY27ehEpFSWnGFfXaR2wVkHAv2zLx+r0TrtmMMsYV9GWyHb00oKAU9q56+Q35ObKRDoytYlfxxuArRvxuu0IyASU5F926tdu3Lae6BaPDdI7QLnk5jGYW9/5Lsfn7fROm3AXFSlAD1Wr411brz9uIECgdMyhMM24czt7nNTjF/yN2KHDYxUURgIjJj5AJI/IEj+Ns5CD3Nhl+EK2VZs2pWrHAiMCMEwkr99GWmP1mmPyWE7SCIKkrxAoBe2Lvnbk8/raJ32mNmXix8InCf25VM7Wqc9NslfIBDol21K/k5TbX3sjNdpW8nfaZezBwKBcTOU5C+siNwCtBd37AqSQCDQnWE02tib+MhoJX+AmYwcjcO2dRUDgcB6LPJhpLJB8rdhhtJy1vFLmAn8hFTA7lzmQGA8aJj0rNuV/BFqT4Lao3XaY3DYQNMOLS6HwX7cAIHAUGzzIx5WRHqQ/GmS/0zyUyR/sQ+jAoFAINDOWiNtkt8M4A0AXikiM5J39GOWqQ/ZpTbkplmeSjUkfAoE1mHbWff2ZaS9bnjkJwC8Q0RmACAiT61vkiGJgNhORApMoV9X2WWhjbuM7Sy0m7A8qQ0o2yKWGlHX5voHTM6TXGxbZNz2wsbYkwiIFCEttnW1Q7DcNt+Otra289GlbaPn6JS2LezkcrKN63cK27rYscq201w/50N27fq5tnWvX5D8nY11nfbLAXwjybcDOAbwZhH56PpmGUjjJOttaa1NdWzzb+ST2uKofmKItHam2mzragdxdtvazkeXts2fo+62bf/69XuO1rGtvs++bdvG9evTtq3BczTSJvkhAC9u+dND9v1XAXwDgK8D8F6SXy4t2jiSDwJ4EADuueeedWwOBAKBU7FPE5ErnbaIvHbZ30j+BID3WSf9dyQ1gNsBfL6ln+sArgPAtWvXRhCtDgQC54l9CY+s+7DyJwC+BQBIvhxACuDpdY0KBAKBvnH5R0762QXWddrvBvDlJD8J4PcBPNAWGgkEAoGhYYefTv2Q91uZ87+RfMtGjD2BtSYiRWQO4Id6siUQCAQ2Rw8jaZIRgF8H8K0APgPgoyQ/ICL/uH7v3Rh17pFAIBDoAwJQ5MqfDrwKwL+JyH/YQevvw6xV2RqDLGN/7LHHnib5X6d82+04X/Hy83S84Vj3lz6O96XrGvH444/92UHC2ztsOiX5qPf6uhVROL4MwH97rz8D4OvXte80DOK0ReSFp30PyUdF5Nom7Bkj5+l4w7HuL2M5XhG5v6eu2objW53HC+GRQCAQ6M5nANztvb4LwP9s04DgtAOBQKA7HwXwlSRfRjIF8P0APrBNA0abmrWF66s32SvO0/GGY91f9up4RSQj+VMA/gxABODdIvKpbdrAIKsOBAKB3SGERwKBQGCHCE47EAgEdoiddNok30xSyE66y52E5C+R/CeS/0Dyj0leGdqmvhl6OfA2IXk3yb8g+YSt8vQzQ9u0aUhGJP+e5J8Obcs+sXNOm+TdMEtIPz20LRvmEQCvEJFXAvgXAG8d2J5e8ZYDfzuArwbwAyS/elirNkoG4GdF5KtgUhn/5J4fLwD8DIAnhjZi39g5pw3gnQB+DlsWtG8bEflzEcnsy4/A6EH3icGXA28TEXlSRB63/34Oxpl92bBWbQ6SdwH4TgC/ObQt+8ZOOW2SrwfwWRH5+NC2bJk3Avjg0Eb0TNty4L11Yj4k7wXwtQD+dlhLNsqvwgyu9NCG7Buj02mvqJTzNgDftl2LNsdJxyoi77fbPATzaP2ebdq2BQZfDjwEJI8A/BGAN4nIs0PbswlIvg7AUyLyGMnXDG3PvjE6p72sUg7JrwHwMgAfp8nGdReAx0m+SkT+d4sm9sZJVYEAgOQDAF4H4L49zFM++HLgbUMygXHY7xGR9w1tzwZ5NYDXk/wOAFMAl0j+roiENM49sLOLa0j+J4BrIrKXGdNI3g/gVwB8k4g0yrftOiRjmAnW+wB8FmZ58A9ue3XZtqAZaTwM4Isi8qah7dkWdqT9ZhF53dC27As7FdM+Z7wLwEUAj5D8GMnfGNqgPrGTrG458BMA3ruvDtvyagA/DOBb7PX8mB2JBgKnYmdH2oFAIHAeCSPtQCAQ2CGC0w4EAoEdIjjtQCAQ2CGC0w4EAoEdIjjtQCAQ2CGC0w4EAoEdIjjtQCAQ2CH+H+SYsfqZD65rAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1)\n", "h = ax.hexbin(x, y, gridsize = 50, cmap = \"Blues\")\n", "plt.colorbar(h, label = \"counts in bin\") " ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }