{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Gaussian Plots" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "data2 = [2*x for x in data]\n", "data2" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As you may have guessed, x acts as a variable that takes on each value in the list data. The first part of the command computes 2*x for each x and assigns the resulting value to the same location in data2 that x was drawn from in data." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "import math\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def gaussian(N,s):\n", " return math.sqrt(2/(math.pi*N))*math.exp(-2*s**2/N)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.252313252202016" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "gaussian(10,0)" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "ins=np.linspace(-10,10,100)\n", "data10=[gaussian(10,s) for s in ins]\n", "data20=[gaussian(20,s) for s in ins]\n", "data100=[gaussian(100,s) for s in ins]" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEWCAYAAACXGLsWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl4HGeV6P/v6dZmW7IkW5I3yfsSb7Idy85G9gUDWWAmEGfITLiEyQ0hszwM8yNzYQKEy8wAM3DDEG4IJJDLlgQIYIiTkM0JIXFiOZYt75ZXyfKu1dauPr8/qtrpdFru1lJd6tb5PE8/qq56q+uo1KpT9b5V7yuqijHGGHMuAb8DMMYYM/xZsjDGGBOXJQtjjDFxWbIwxhgTlyULY4wxcVmyMMYYE5clC2MAEXlGRG73cfv/S0R+eI7lnxCR15IZkzGRLFkYX4jIahF5U0TOiMhxd/puERE/4lHVD6jqY35s293+v6nqpwBEZLqIqIhkDPTzROSAiBwTkTER8z4lIusG8FlfFZFqEekRkS/HWP5XInLQ/Vv+VkTGDTRuM3xZsjBJJyL/BDwAfBOYCEwA7gIuAbJ8DC3dZAD/MASfUwP8f8DT0QtEZCHwfeCvcf6ObcD3hmCbZpixZGGSSkTygfuBu1X1V6raqo5NqvpxVe10y31IRDaJSIuI1Eae0YrIFSJSF/W5B0TkGnd6pYhUuuseE5FvufNzROSnInJKRJpEZIOITHCXrROR8Jn9LBF5yS13UkR+JiIFUdv6nIhsEZFmEXlCRHL6+H0Pishyd/o294phgfv+UyLyW3f6yyLyU3e1V92fTSJyWkQuivi8/xSRRhHZLyIfiLO7vwl8LjL2gVDVx1T1GaA1xuKPA79X1VdV9TTwr8BfiEjeYLZphh9LFibZLgKygd/FKXcG+BugAPgQ8GkR+XCC23gAeEBVxwKzgCfd+bcD+UAZMB7naqY9xvoC/DswGZjvlv9yVJmPAauAGUA58Ik+YnkFuMKdvgzYB1we8f6VGOtc5v4sUNVcVX3DfX8BsAsoAr4BPBKn2q4SWAd8LtZCN9k19fFK9OpgIbA5/EZV9wJdwNwE1zcpwpKFSbYi4KSq9oRniMjr7gGqXUQuA1DVdaparaohVd0C/IJ3DrLxdAOzRaRIVU+r6vqI+eOB2araq6obVbUlemVVrVHV51W1U1VPAN+Kse3vqGq9qjYAvweW9hHLKxHrXoqThMLvLyd2sujLQVX9gar2Ao8Bk3Cqfs7lPuDvRKQ4eoGqlqtqQR+vuxOMKRdojprXDNiVRZqxZGGS7RRQFNl4q6oXq2qBuywAICIXiMjLInJCRJpxrgKKEtzGHThntjvdqqbr3fk/AZ4DHheRehH5hohkRq8sIiUi8riIHBaRFuCnMbZ9NGK6DeegGcsrwKUiMhEIAk8Al4jIdJyrnKoEf6d3bVNV29zJvrYbLrcV+ANwbz+20x+ngbFR88YSu8rKpDBLFibZ3gA6gZvilPs5sAYoU9V84CGc6iFwqqhGhwuKSBA4e+asqntU9VagBPg68CsRGaOq3ar6FVVdAFwMXI9T1RXt3wEFyt2qrNsitt0vqlqDk0z+HnhVVVtxDvp3Aq+paijWagPZ1jl8CfhbYErkTBHZ5raJxHo9lOBnbwOWRHzmTJxqxt1DFbwZHixZmKRS1SbgK8D3RORmEckVkYCILAXGRBTNAxpUtUNEVgJ/FbFsN5DjNoJnAl/EOUABZxuSi90DcZM7u1dErhSRxW5yacGpluqNEWYezhlzk4hMAf55kL/2K8A9vFPltC7qfbQTQAiYOcjtAmcT1hM4CSty/kK3TSTW665wORHJdBvwA0CGe6NA0F38M+AGEblUnNt07weecpOiSSOWLEzSqeo3gM/i3I55HDiGc/vl54HX3WJ3A/eLSCtOvfuTEes3u8t/CBzGudKIvDtqFbBNRE7jNHavVtUOnNt0f4WTKHbgHKx/ynt9BTgfp+79aeCpQf7Kr+AkoFf7eP8ubhXT14A/u205Fw5y++AcxMfELRXbD3BuBLgV+II7/ddurNtwqgh/hvO3zMP525g0Izb4kTHGmHjsysIYY0xcliyMMcbEZcnCGGNMXJYsjDHGxDXgXi2Hm6KiIp0+fbrfYRhjTErZuHHjSVV9zxP+0dImWUyfPp3Kykq/wzDGmJQiIgcTKWfVUMYYY+KyZGGMMSYuSxbGGGPismRhjDEmLksWxhhj4rJkYYwxJi5Pk4WIrBKRXSJSIyLvGXxFRD4rItvd4R1fFJFpEct6RaTKfa3xMk5jjDHn5tlzFm5/9w8C1+J0H71BRNao6vaIYpuAClVtE5FP44wrfIu7rF1V+xqq0phhr76pnV9trKOn1xnfqHBMFh+/YBpZGXZBb1KPlw/lrQRqVHUfgIg8jjM62tlkoaovR5RfjzMimTEp70hzOx/7/hvUNbYj7hh7qlB5oJEHVi8lI2gJw6QWL5PFFKA24n0dcME5yt8BPBPxPkdEKoEe4D9U9bdDH6IxQ+9Eaycf/8GbNLd1s+aeSygvLQDgB6/u42trd5CdGeA/b15CIDCgkVqN8YWXySLWf0LMkZZE5DagArg8YvZUVa13x/R9SUSqVXVv1Hp34oxlzNSpU4cmamMGofFMF7f98E2ONHfwkztWnk0UAH972Uzau3v51vO7yckM8rUPL0LEEoZJDV5eC9cBZRHvS4H66EIicg3OUI03qmpneL6q1rs/9+GMWbwsel1VfVhVK1S1org4bj9Yxnjufz+9g/2nzvDI7RVUTB/3nuV/d9Vs7rp8Fj9/8xDPbj3qQ4TGDIyXyWIDMEdEZohIFrAaeNddTSKyDGfs5RtV9XjE/EIRyXani4BLiGjrMGY4OnjqDL+tOszfXDiNi2cXxSwjIvzz++cxs2gM33mpBhvW2KQKz5KFqvYA9wDPATuAJ1V1m4jcLyI3usW+CeQCv4y6RXY+UCkim4GXcdosLFmYYe17L+8lGBDuvGzmOcsFA8I9V81mx5EWXthx/JxljRkuPO2iXFXXAmuj5t0XMX1NH+u9Diz2MjZjhlJtQxu/fruO2y6cRsnYnLjlb1wymQde3MN3XtzDNfNLrO3CDHt2/54xQ+D/vrKXgAj/8/JzX1WEZQQDfObK2VQfbmbdrhMeR2fM4FmyMGaQ6pva+WVlLR9bUcqk/FEJr/eRZVMoLRzFAy/usbYLM+xZsjBmkB55bT+qcNfls/q1XqZ7dVFV28T6fQ0eRWfM0LBkYcwg9IaUNZvruXp+CaWFo/u9/keWTWFMVpDfVR32IDpjho4lC2MG4a39DZxo7eT68skDWj8nM8i1Cybw7LajdPWEhjg6Y4aOJQtjBuH3W+oZlRnk6vklA/6MG5ZMpqmtmz/XnBzCyIwZWpYsjBmg7t4Qz249ytXzSxidNfC70C+dU8zYnAx+v+U9HRwYM2xYsjBmgF7fe4qGM13csGRgVVBhWRkBVi2ayB+3HaOju3eIojNmaFmyMGaA/rC5nrzsDC6fO/h+ya4vn8zpzh5e2W3PXJjhyZKFMQPQ2dPLs9uOcu3CCeRkBgf9eRfPGs+4MVn8frNVRZnhyZKFMQPwp90nae3o4YYB3gUVLSMY4AOLJvLijuO0dfUMyWcaM5QsWRgzAGurj5A/KpNL+uhddiBuWDKZ9u5e6/7DDEuWLIzpp1BIeXXPSa6YVzyk42lXTCskLyeDV63dwgxDliyM6aedR1s5ebqT9w3hVQU4VVGXzCriT3tOWl9RZtixZGFMP71W45z5Xzpn6EdnfN+cIg43tbP/5Jkh/2xjBsOShTH99Kc9J5lTksvE/PjjVvTXZW4C+tMee5rbDC+WLIzph47uXt7a3+DJVQXA1PGjmTputCULM+xYsjCmHzYcaKCzJ8Slc4a2vSLSpXOKWL/vFN291rGgGT4sWRjTD6/tOUlmULhg5jjPtnHpnCJOd/ZQVdvk2TaM6S9LFsb0w5/2nGT5tMJBdRwYz0WzigiItVuY4cWShTEJOtHayfYjLZ61V4Tlj8pkSVkBf9pjz1uY4cOShTEJen2vc6bvZXtF2KWzi9hc20Rze7fn2zImEZYsjEnQq7tPUjg6k4WT8z3f1qVziwkpvLHXqqLM8GDJwpgEvbn/FBfMGE8wIJ5va2lZAaMyg6zf1+D5toxJhCULYxJwtLmDusZ2Vszw7i6oSJnBAEvLCth4sDEp2zMmHksWxiSg8qBzhr9iemHStrlieiHb6ps53Wldlhv/WbIwJgGVBxoZnRVkwaSxSdtmxfRxhBSqDtnzFsZ/liyMScCGAw0sm1pARjB5/zLLphYQEGfbxvjNkoUxcbR2dLPjSAvLpyWnvSIsLyeT8yaOPVsFZoyfLFkYE8emQ02ENLntFWErphey6VATPdZPlPGZJQtj4qg82EhAYNnU5CeLiunjaOvqZceR1qRv25hIniYLEVklIrtEpEZE7o2x/LMisl1EtojIiyIyLWLZ7SKyx33d7mWcxpxL5YEGFkweS262d/1B9aXCvZqxdgvjN8+ShYgEgQeBDwALgFtFZEFUsU1AhaqWA78CvuGuOw74EnABsBL4kogk/7TOjHjdvSE2HWqiIsntFWGT8kdRWjjK2i2M77y8slgJ1KjqPlXtAh4HboosoKovq2qb+3Y9UOpOvx94XlUbVLUReB5Y5WGsxsS0vb6F9u5eVkz3J1kArJg+jg0HGm1cbuMrL5PFFKA24n2dO68vdwDP9GddEblTRCpFpPLECeuh0wy9cPVPhQ+N22HLpxVyorWTQw1t8Qsb4xEvk0WsDnRinhqJyG1ABfDN/qyrqg+raoWqVhQXe9tttBmZKg80UjZuFBPGDv1424kKX9VsOGBdfxj/eJks6oCyiPelQH10IRG5BvgCcKOqdvZnXWO8pKq8faiR5T7cBRVpTkkueTkZ1k+U8ZWXyWIDMEdEZohIFrAaWBNZQESWAd/HSRTHIxY9B1wnIoVuw/Z17jxjkuZIcwfHWztZWlbgaxyBgLCktIDNNsyq8ZFnyUJVe4B7cA7yO4AnVXWbiNwvIje6xb4J5AK/FJEqEVnjrtsAfBUn4WwA7nfnGZM04YPzUp+vLMDpsnzXsVbau3r9DsWMUJ7eOK6qa4G1UfPui5i+5hzrPgo86l10xpxbVW0TWcEA8yfl+R0KS8oK6A0pW+ubfb0zy4xc9gS3MX2oqm1i/uSxZGcE/Q6FJWXO6HxWFWX8YsnCmBh6Q0r14WaWlno/hGoiSvJymFIwik2WLIxPLFkYE8Oe4620dfWydKq/jduRlpZZI7fxjyULY2IIH5SXlA6fZLGkLJ+6xnZOnu6MX9iYIWbJwpgYqmqbGJuTwYyiMX6HctbSMueuLLu6MH6wZGFMDFW1zSwpK0AkVmcC/lg0ZSzBgFiyML6wZGFMlLauHnYdbWGZzw/jRRudlcHcCXnWyG18YcnCmChbD7cQUufZhuFmaVk+m2ubrAdak3SWLIyJUlXr9ME0PJNFAS0dPew/ecbvUMwIY8nCmCiba5spLRxFUW6236G8RziBba6zqiiTXJYsjIlSVds0LK8qAOaU5DE6K0jVIUsWJrksWRgT4dTpTg43tbNkmDy5HS0YEBZNzmfL4Wa/QzEjjCULYyJUuwfhxVOG55UFwOLSfLbXt9DdG/I7FDOCWLIwJkJ1nZMsFk0Z63MkfSsvzaezJ8SeY6f9DsWMIJYsjImw5XAzM4vHkJeT6XcofVo8xakiqz5s7RYmeSxZGBOhuq6Z8inDs70ibPr4MeRlZ5ytMjMmGSxZGOM63tLB0ZYOFg+jzgNjCQSERVPyz1aZGZMMliyMcYXP1MuH6Z1QkcpL89lxpJWuHmvkNslhycIY15a6ZgICCyYN38btsMWl+XT1hth9rNXvUMwIYcnCGFf14WZml+QyJtvToemHRLl7a+8Wq4oySWLJwhhAVdlS1zysn6+IVDZuFPmjMu2OKJM0liyMAY62dHDydGdKtFcAiAjlpfl2ZWGSxpKFMbxTnbM4RZIFOM9b7DraSkd3r9+hmBHAkoUxOM9XBAOSEo3bYeWl+fSElJ1HrZHbeM+ShTE4T27PnZBHTmbQ71ASFn4epNq6KzdJYMnCjHiqytbDzSwexv1BxTI5P4fxY7LYbO0WJgksWZgR73BTOw1nuob9k9vRRJwnubdatx8mCSxZmBEv3G3GcO8TKpby0nx2H2ulvcsauY23LFmYEa/6cDMZAWHexDy/Q+m3RVPyCSlsP9LidygmzVmyMCNedQo2boeFuyu3qijjNU+ThYisEpFdIlIjIvfGWH6ZiLwtIj0icnPUsl4RqXJfa7yM04xcqkr14eaUeRgv2qT8HIpys+zhPOM5zzrBEZEg8CBwLVAHbBCRNaq6PaLYIeATwOdifES7qi71Kj5jAOoa22lq62ZRCrZXgDVym+Tx8spiJVCjqvtUtQt4HLgpsoCqHlDVLYD1s2x8sfXsmNupmSzAiX3PcWvkNt7yMllMAWoj3te58xKVIyKVIrJeRD4cq4CI3OmWqTxx4sRgYjUj1JbDzWQGhfMmpV7jdtjis43cdnVhvONlspAY87Qf609V1Qrgr4D/IyKz3vNhqg+raoWqVhQXFw80TjOCbXUbt7MzUq9xOyzcn5WNnGe85GWyqAPKIt6XAvWJrqyq9e7PfcA6YNlQBmdMuHE7laugACaOdRq5qw/b7bPGO14miw3AHBGZISJZwGogobuaRKRQRLLd6SLgEmD7udcypn/Cjdup1NNsLCLC4in5NraF8ZRnyUJVe4B7gOeAHcCTqrpNRO4XkRsBRGSFiNQBHwW+LyLb3NXnA5Uishl4GfiPqLuojBm06jRo3A5bPCWfmuOnaevq8TsUk6Y8HT9SVdcCa6Pm3RcxvQGneip6vdeBxV7GZsyWOqdxOxWf3I62uLTAaeSub6Fi+ji/wzFpyJ7gNiPW1sPNzJuY2o3bYeGro2p73sJ4xJKFGZHSpXE7bMLYbIpys+2OKOOZhKqhRKQUp4H6UmAy0A5sBZ4GnlFVe6jOpJTahnaa21P3ye1o4TG57crCeCXulYWI/Ah4FOgCvg7cCtwNvACsAl4Tkcu8DNKYobbZHV1uSYqNYXEui6fkU3PiNKc7rZHbDL1Eriz+S1W3xpi/FXjKvS126tCGZYy3qg83k5URYO6E1G/cDltSlo8qbDvczAUzx/sdjkkzca8sYiUK9zmIcnd5l6rWeBGcMV7ZXNvE/EljycpIn2a7xVPcMbmtKsp4IOH/FBFZJyJjRWQcsBn4kYh8y7vQjPFGKOSMub0kxR/Gi1acl83k/Bwbk9t4oj+nVfmq2gL8BfAjVV0OXONNWMZ4Z9/JM5zp6k2bO6EiLS7Np7rOnuQ2Q68/ySJDRCYBHwP+4FE8xnhuS7hxuyx9GrfDyksLOHCqjea2br9DMWmmP8nifpyuO2pUdYOIzAT2eBOWMd7ZUtfMqMwgs4pz/Q5lyIVH/LN2CzPUEk4WqvpLVS1X1bvd9/tU9S+9C80Yb2ypa2LRlLEEA7F60U9t5W4j9xbrVNAMsQHfCiIid4vILSLiaf9Sxgylnt4Q2+pbKE+j5ysi5Y/OZNr40WyptSsLM7QGc9+gAO8DnhqiWIzx3O5jp+nsCZ2trklH5aUFVg1lhtyArwpU9cGhDMSYZAiP+ZCuVxYA5VPy+f3mek6e7qQoN9vvcEyaiJssRORvzrFYVfUnQxiPMZ7aXNdMXk4G08aN9jsUz5RHDLN65XklPkdj0kUiVxYrYswT4AZgCmDJwqSM6rpmykvzCaRh43bYwin5iDh3fVmyMEMlbrJQ1b8LT4uIAB8HPg+sB77mXWjGDK3Onl52Hm3hjvfN9DsUT+VmZzC7OPfs8yTGDIVEuyjPAD4B/BPwJnCzqu7yMC5jhtz2+ha6ezXtuvmIpby0gFd2H0dVcc7xjBmcRLoo/wywHVgOrFLVT1iiMKmoqtY50146NX0bt8OWTi3g5Oku6hrb/Q7FpIlEriz+GziOc5vs7yPOUgSngbvco9iMGVJVtU1MGJvNpPxRfofiuWVuVyZVtU2UpXFjvkmeRJLFDM+jMCYJqmqbWJqG/UHF4owtHqCqtokblkz2OxyTBhJJFodUVc9VQEQkXhlj/NRwpouDp9pYvWJkjNOVGQyweEr+2ao3YwYrkSe4XxaRvxORd/2XiUiWiFwlIo8Bt3sTnjFDY7N70Fw2AtorwpaWFbD1cDNdPSG/QzFpIJFksQroBX4hIkdEZLuI7MPpcfZW4Nuq+mMPYzRm0DYdaiQgpOUYFn1ZOrWAzp4QO4+2+B2KSQOJPGfRAXwP+J6IZAJFQLuq2vWtSRmbapuYOyGPMdkjp9/LpRGN3OncvYlJjkRunc0RkX8Uke8C/wM4YYnCpJJQSNlc2zSiqqAAphSMoig3m6pD9u9qBi+RaqjHgAqgGvgg8F+eRmTMENt/6gwtHT0j5k6oMBFh2dQCa+Q2QyKRa/IFqroYQEQeAd7yNiRjhlb4zHppWaHPkSTf0rICnt9+jOa2bvJHZ/odjklhiVxZnB3MV1V7PIzFGE9U1TYxJivI7JL0G0Y1nrMP51k/UWaQEkkWS0SkxX21AuXhaRGx2yzMsFdV28SSsoK0HEY1nsWlTg+01m5hBituslDVoKqOdV95qpoRMT32XOuKyCoR2SUiNSJyb4zll4nI2yLSIyI3Ry27XUT2uC97jsMMSEd3LzuOtIy49oqwvJxM5pTksqm20e9QTIobzLCq5yQiQeBB4APAAuBWEVkQVewQTm+2P49adxzwJeACYCXwJREZeRXOZtCqDzfTE9IRmywAlpUVUlXbRChknSyYgfMsWeAc5GtUdZ+qdgGPAzdFFlDVA6q6BYh+xPT9wPOq2qCqjcDzOA8HGtMvlQecM+rl00buucbyaYU0tXWz7+Rpv0MxKczLZDEFqI14X+fOG7J1ReROEakUkcoTJ04MOFCTvioPNDCzeAzjR/BY1BXTnUS54YBVRZmB8zJZxGpNTPQ6OKF1VfVhVa1Q1Yri4uJ+BWfSXyikbDzUSMUIvqoAmFE0hvFjss5eZRkzEF4mizqgLOJ9KVCfhHWNAWDvidM0tXVTMX2c36H4SkRYPq2QyoMNfodiUpiXyWIDMEdEZohIFrAaWJPgus8B14lIoduwfZ07z5iEVR50zqRH+pUFOFVRB0+1caK10+9QTIryLFm4D/Ddg3OQ3wE8qarbROR+EbkRQERWiEgd8FHg+yKyzV23AfgqTsLZANzvzjMmYZUHGhk/JosZRWP8DsV34aurjXZ1YQbI0y44VXUtsDZq3n0R0xtwqphirfso8KiX8Zn0VnmwgeXTCokYCnjEWjQ5n+yMABsONLJq0SS/wzEpyMtqKGN8c7y1g4On2lgxwtsrwrIyAiwpKzhbNWdMf1myMGlpY/j5iunWXhFWMa2QbYebae/q9TsUk4IsWZi0VHmwkeyMAIsmj5yR8eKpmF5IT0ity3IzIJYsTFqqPNDAkrICsjLsKx62fKpTJVd5wBq5Tf/Zf5JJO+1dvWyrb7FbZqPkj85k7oRca7cwA2LJwqSdTbWN9IT0bDcX5h0V08fx9sFGeq1TQdNPlixM2lm/9xTBgNidUDFcOHM8rZ09bKtv9jsUk2IsWZi08/reUyyekk9ejg0jGu2imeMBZx8Z0x+WLExaOdPZQ1VtExfNGu93KMNScV42c0pyLVmYfrNkYdLKhgMN9ISUiy1Z9OniWePZsL+Brp7oYWSM6ZslC5NW3th7isygUDHN2iv6ctGsItq7e9lcZ89bmMRZsjBp5Y19p1hWVsiorKDfoQxbF84ch4iTWI1JlCULkzaa27rZerjZ2iviKBidxYJJY3l970m/QzEpxJKFSRtv7j9FSLH2igRcPGs8bx9soqPb+okyibFkYdLG63tPkZ0RYOnUAr9DGfYumjWert4QG+1pbpMgSxYmbazfd4oV08eRnWHtFfGsmD6OYECs3cIkzJKFSQsnT3ey82irtVckKC8nk/LSfGu3MAmzZGHSQvghM0sWibt41ng21zXT0tHtdygmBViyMGlh3a7jFIzOZEmptVck6rI5xfSGlD/vsasLE58lC5PyQiHllV0nuHxuMcGAjbedqOXTCsnLyeClncf9DsWkAEsWJuVtOdzMqTNdXDmvxO9QUkpGMMBlc4tZt/sEIeuy3MRhycKkvJd3HkcELp9b7HcoKeeqeSWcaO1k+5EWv0Mxw5wlC5PyXt51nGVlBRSOyfI7lJRz+TwnwVpVlInHkoVJaSdaO9lS12xVUANUlJvNktJ8Xt5lycKcmyULk9Je2X0CgCvPs2QxUFeeV0JVbRMNZ7r8DsUMY5YsTEp7eedxSvKyWTh5rN+hpKwr55WgCq/stqsL0zdLFiZldfeGeHXPCa6cV4KI3TI7UIun5FOUm8XLO0/4HYoZxixZmJS18WAjrR09XHme3QU1GIGAcPncEl7ZfYKeXhs9z8RmycKkrD9uO0ZWMMAls4v8DiXlXT2/hOb2bt460OB3KGaYsmRhUlIopKytPsLl84rJy8n0O5yUd+W8EkZlBnl6yxG/QzHDlCULk5I2HmrkaEsH15dP8juUtDAqK8jV80t4dutRq4oyMXmaLERklYjsEpEaEbk3xvJsEXnCXf6miEx3508XkXYRqXJfD3kZp0k9T285QlZGgKvnT/A7lLRxffkkTp3pYv0+q4oy7+VZshCRIPAg8AFgAXCriCyIKnYH0Kiqs4FvA1+PWLZXVZe6r7u8itOknl63CurKecXkZmf4HU7auGJeCaOzgjxdXe93KGYY8vLKYiVQo6r7VLULeBy4KarMTcBj7vSvgKvF7oE0cVQeaOB4aycfKp/sdyhpJSczyDXzJ/Ds1qN0W1WUieJlspgC1Ea8r3PnxSyjqj1AMxAevWaGiGwSkVdE5NJYGxCRO0WkUkQqT5ywe8RHiqerj5CdEeBqe2p7yH2ofBKNbd023Kp5Dy+TRawrhOh+kPsqcwSYqqrLgM8CPxeR9zyiq6oPq2qFqlYUF9u99iOBUwV1lKvOK2GMVUENucvnFjMmy+6KMu/lZbKoA8oi3pcC0ZWhZ8uISAaQDzSoaqdrSyF6AAASs0lEQVSqngJQ1Y3AXmCuh7GaFPHW/gZOnu7kQ3YXlCdyMoNcu2ACz26zqijzbl4miw3AHBGZISJZwGpgTVSZNcDt7vTNwEuqqiJS7DaQIyIzgTnAPg9jNSniN5vqGJ0V5CqrgvLM9eWTaW7vtm7Lzbt4lizcNoh7gOeAHcCTqrpNRO4XkRvdYo8A40WkBqe6KXx77WXAFhHZjNPwfZeq2v18I1xLRze/33yEG5dMZnSWVUF55Yp5xZTkZfP4W4f8DsUMI57+x6nqWmBt1Lz7IqY7gI/GWO/XwK+9jM2knt9V1dPe3cutK6f6HUpaywgGuGVFGd99uYbDTe1MKRjld0hmGLAnuE1KUFV+/uYhFkwaS3lpvt/hpL2PVTjNjU9sqI1T0owUlixMSthS18yOIy3cesFU6448CcrGjeayOcU8uaHWuv8wgCULkyIe33CIUZlBblpqD+Ily60rp3K0pePsaIRmZLNkYYa90509/K6qnhuWTGKs9TCbNFfPL6E4L5tfWEO3wZKFSQG/qzpMW1cvq61hO6kygwE+uryUl3Yep76p3e9wjM8sWZhhrTek/PBP+1k4eSzLygr8DmfEuXWl00b0yGv7/Q7F+MyShRnW1lYfYf/JM3zmytnWsO2DsnGjuWnJZH7+5iEaznT5HY7xkSULM2ypKg++XMOs4jGsWjjR73BGrLuvnEVHTy8/+rNdXYxklizMsPXSzuPsPNrK3VfMJhCwqwq/zC7JY9XCifz49QO0dHT7HY7xiSULMyypKt99uYbSwlHcaLfL+u7uK2bT2tHDT9446HcoxieWLMyw9MbeU2w61MT/vHwWmUH7mvptcWk+l88t5tHX9tPe1et3OMYH9l9ohp1QSPnmH3dRkpfNR5eX+h2Ocd1z1WxOneniUWu7GJEsWZhh56lNh9l0qIl/fv88cjKDfodjXCumj+P9Cyfw3Zdq7LmLEciShRlWWjq6+Y9ndrJsagF/eb5dVQw3X/zQAkKq/NvaHX6HYpLMkoUZVh54YQ+nznTylRsX2h1Qw1DZuNF8+opZ/GHLEV7fe9LvcEwSWbIww8aeY6089voBVq8oo7zUntYeru66fBalhaP4yprt1iPtCGLJwgwLPb0h/uWpakZnBfncdfP8DsecQ05mkH+9fgG7jrXy/VdttOORwpKFGRa+/cJuKg82cv9Nixifm+13OCaO6xZM4PrySXzr+d1sOGAjHo8EliyM717dfYLvrdvLxypK+fCyKX6HYxIgIvz7XyymtHAUf/+LTTRav1Fpz5KF8dXxlg4++2QVs4tz+cqNi/wOx/RDXk4m3731fE6e7uRzv9yMqvodkvGQJQvjm/auXj7z87c53dnDgx8/n1FZ9kxFqllcms//+uB8Xtx5nO+8WON3OMZDGX4HYEamzp5e7vrpRioPNvKd1cuYOyHP75DMAH3i4ulUH27m2y/sJjcngzveN8PvkIwHLFmYpOvpDfEPv6jild0n+PpfLuaGJdZRYCoTEb7xl+W0dfby1T9sJzc7yC0rbFTDdGPVUCapOrp7+eyTm3l221Huu36BHVTSREYwwHduXcblc4u596lqnthg43anG0sWJmmONndwy8PrWbO5ns+vOo9PWnVFWsnKCPDQbct53+wiPv/rar68Zhvd9tBe2rBkYZKi8kAD1//3a9Qca+Wh25bz6Stm+R2S8cCorCA/+sQKPvW+Gfz49QP89SNvcvJ0p99hmSFgycJ4qrm9my/9bisf+/4b5GYH+c1nLmHVIhsiNZ1lBAN88foFfPuWJWw61MS133qFx986RChkt9amMmvgNp7o7g3x202H+fqzO2k408VtF07jn66bR/6oTL9DM0nykWWlLJyczxd/s9Vpx6is5V+vX8D5Uwv9Ds0MgKTLgzQVFRVaWVnpdxgjXnN7N4+/dYgfv36AI80dLCkr4GsfXsSiKfl+h2Z8oqr8ZtNh/m3tDk6e7qJiWiGfunQG1y6YSNB6FvadiGxU1Yq45SxZmMFq7ejmpZ3Heab6KOt2H6ejO8RFM8fzqUtncOW8Eutq3ABwurOHJzfU8uif91PX2M6Esdm8f+FEVi2ayMrp48iw4XN9MSyShYisAh4AgsAPVfU/opZnA/8PWA6cAm5R1QPusn8B7gB6gb9X1efOtS1LFsnR3Rvi4Kk2dh1tZePBRjYebGBbfQs9IaUkL5vrFk5g9YqpdiVh+tQbUv647Si/rTrMul0n6OwJkZedwbJphVRMK2RpWQFzJuQycWwOInai4TXfk4WIBIHdwLVAHbABuFVVt0eUuRsoV9W7RGQ18BFVvUVEFgC/AFYCk4EXgLmq2udI8ZYs+q83pHT29NLZHaK9u5e2rh5Od/bS2tFNY1s3TW1dnDrdxbGWDuqbOzjc2MbBU230uA2V2RkBlpQVsHxaIVefV8L5UwvtKsL0S1tXD+t2neC1mpNsPNDI7uOthA9JudkZzCgaw6T8HCbl5zAhP4fC0VkUjs4kf1QWeTkZjM4KkpudQXZmkOyMANkZAUsw/ZRosvCygXslUKOq+9yAHgduArZHlLkJ+LI7/Svgu+L8pW8CHlfVTmC/iNS4n/fGUAfZ1NbFRx8a8o89p/6m58iE/q511Xmvqu5PUNT5qRBSdV9OYgi/untDdPeGSOTmFBEoys1mUn4Os0tyuW7hROaU5DKnJI95E/PIyrCqAzNwo7My+ODiSXxw8STAafPaVt/M3uOnqTl+mv2nnBOU9ftO0dLRk9BnZgSEjKCQGQgQDApBEQIB96c4T5yLON/tgAiCOw8gIs9Ep5z+JqFkpqzzJo3lv29d5uk2vEwWU4DaiPd1wAV9lVHVHhFpBsa789dHrfuevqtF5E7gToCpUwf2JHAgIMyZkDugdQdD+vtV6uNLLGe/7LzrSx8IuP8YCIEABN1/lmAgQGZQyAwGyAgKOZlBcjIC5GQGGZOdwZjsILnZmRSOzqRgdBYFozPJtLpkkyT5ozK5eFYRF88qes+y9q5eGtu6aGzrormtmzNdvZzp7OFMVw8d3SE6e3rp6A7R0xuixz0pCp8ghVQJhZwTqF7nrArFea8RJ11h7zmP6ucZnvb7lHBwygpHeb4NL5NFrKNh9B7sq0wi66KqDwMPg1MN1d8AAcbmZPK9jy8fyKrGmCQalRVkVNYoJhd4f2A07+XlKWMdUBbxvhSo76uMiGQA+UBDgusaY4xJEi+TxQZgjojMEJEsYDWwJqrMGuB2d/pm4CV1rgXXAKtFJFtEZgBzgLc8jNUYY8w5eFYN5bZB3AM8h3Pr7KOquk1E7gcqVXUN8AjwE7cBuwEnoeCWexKnMbwH+My57oQyxhjjLXsozxhjRrBEb52121yMMcbEZcnCGGNMXJYsjDHGxGXJwhhjTFxp08AtIieAg4P4iCLg5BCFM5Qsrv6xuPrH4uqfdIxrmqoWxyuUNslisESkMpE7ApLN4uofi6t/LK7+GclxWTWUMcaYuCxZGGOMicuSxTse9juAPlhc/WNx9Y/F1T8jNi5rszDGGBOXXVkYY4yJy5KFMcaYuEZUshCRj4rINhEJiUhF1LJ/EZEaEdklIu/vY/0ZIvKmiOwRkSfcrteHOsYnRKTKfR0Qkao+yh0QkWq3nOc9KIrIl0XkcERsH+yj3Cp3H9aIyL1JiOubIrJTRLaIyG9EpKCPcknZX/F+f7fb/Sfc5W+KyHSvYonYZpmIvCwiO9zv/z/EKHOFiDRH/H3v8zoud7vn/LuI4zvu/toiIucnIaZ5EfuhSkRaROQfo8okZX+JyKMiclxEtkbMGyciz7vHoedFpLCPdW93y+wRkdtjlekXVR0xL2A+MA9YB1REzF8AbAaygRnAXiAYY/0ngdXu9EPApz2O97+A+/pYdgAoSuK++zLwuThlgu6+mwlkuft0gcdxXQdkuNNfB77u1/5K5PcH7gYecqdXA08k4W83CTjfnc4DdseI6wrgD8n6PiX6dwE+CDyDM3rmhcCbSY4vCBzFeXAt6fsLuAw4H9gaMe8bwL3u9L2xvvPAOGCf+7PQnS4cTCwj6spCVXeo6q4Yi24CHlfVTlXdD9QAKyMLiDNa+1XAr9xZjwEf9ipWd3sfA37h1TY8sBKoUdV9qtoFPI6zbz2jqn9U1R737XqcURX9ksjvfxPOdwec79LV7t/aM6p6RFXfdqdbgR3EGNN+mLoJ+H/qWA8UiMikJG7/amCvqg6md4gBU9VXccb6iRT5HerrOPR+4HlVbVDVRuB5YNVgYhlRyeIcpgC1Ee/reO8/03igKeLAFKvMULoUOKaqe/pYrsAfRWSjiNzpYRyR7nGrAh7t49I3kf3opU/inIXGkoz9lcjvf7aM+11qxvluJYVb7bUMeDPG4otEZLOIPCMiC5MUUry/i9/fqdX0fcLmx/4CmKCqR8A5EQBKYpQZ8v3m2Uh5fhGRF4CJMRZ9QVV/19dqMeZF31OcSJmEJBjjrZz7quISVa0XkRLgeRHZ6Z6FDNi54gL+L/BVnN/5qzhVZJ+M/ogY6w763uxE9peIfAFnVMWf9fExQ76/YoUaY55n36P+EpFc4NfAP6pqS9Tit3GqWk677VG/xRnO2Gvx/i5+7q8s4EbgX2Is9mt/JWrI91vaJQtVvWYAq9UBZRHvS4H6qDIncS6BM9wzwlhlhiRGEckA/gJYfo7PqHd/HheR3+BUgQzq4JfovhORHwB/iLEokf045HG5jXfXA1erW2Eb4zOGfH/FkMjvHy5T5/6d83lvNcOQE5FMnETxM1V9Knp5ZPJQ1bUi8j0RKVJVTzvNS+Dv4sl3KkEfAN5W1WPRC/zaX65jIjJJVY+4VXLHY5Spw2lXCSvFaasdMKuGcqwBVrt3qszAOUN4K7KAexB6GbjZnXU70NeVymBdA+xU1bpYC0VkjIjkhadxGnm3xio7VKLqiT/Sx/Y2AHPEuWssC+cSfo3Hca0CPg/cqKptfZRJ1v5K5Pdfg/PdAee79FJfCW6ouG0ijwA7VPVbfZSZGG47EZGVOMeGUx7HlcjfZQ3wN+5dURcCzeEqmCTo8+rej/0VIfI71Ndx6DngOhEpdKuMr3PnDZzXrfnD6YVzkKsDOoFjwHMRy76AcyfLLuADEfPXApPd6Zk4SaQG+CWQ7VGcPwbuipo3GVgbEcdm97UNpzrG6333E6Aa2OJ+WSdFx+W+/yDO3TZ7kxRXDU7dbJX7eig6rmTur1i/P3A/TjIDyHG/OzXud2lmEvbR+3CqILZE7KcPAneFv2fAPe6+2Yxzo8DFSYgr5t8lKi4BHnT3ZzURdzF6HNtonIN/fsS8pO8vnGR1BOh2j1134LRxvQjscX+Oc8tWAD+MWPeT7vesBvgfg43FuvswxhgTl1VDGWOMicuShTHGmLgsWRhjjInLkoUxxpi4LFkYY4yJy5KFMYCI/B8Rucyd/qGILOjHup8Qke/2c3sHRKSov3EOhIi80FfPpMYkypKFGfFEZBxwobrdTKjqp1R1u89hDaWf4PR2a8yAWbIwac19Svhpt8O3rSJyS4xiNwPPRqyzTtzxTkTktIh83e3o7gURWeku3yciN0Z8RpmIPCvOWBZfivis37rrbuurA8O+yrjb/pob+3oRmeDOnyDO2B2b3dfF7vzbROQtccZX+L6IBN2PWoPzNLIxA2bJwqS7VUC9qi5R1UVEJIUIlwAb+1h/DLBOVZcDrcD/Bq7F6Q3g/ohyK4GPA0uBj8o7g2t90l23Avh7EYnVw2xfZcYA61V1CU5/SX/rzv8O8Io7/3xgm4jMB27B6ZhvKdDrxoM6XVRn97FtYxJiycKku2rgGvfq4FJVbY5RZhJwoo/1u3gnwVTjHKS73enpEeWeV9VTqtoOPIXTxQY4B/9wlxBlxO6ZtK8yXbzTYePGiO1dhdMLMKra6/5OV+N0PLlBnNEVr8bpTiPsOE4XKMYMSNr1OmtMJFXdLSLLcfpC+ncR+aOq3h9VrB2nz6ZYuvWdPnFCOP2Koaoht9fYs5uK3rSIXIHTKeRFqtomIuuitxOnTOS2ezn3/6sAj6lqrO60cT+z/RzrG3NOdmVh0pqITAbaVPWnwH/iVNtE2wHMHuSmrhVnbORROCOX/Rmn+/FGNwmchzMsaLREykR7Efg0gIgERWSsO+9mccaFCI/TPM2dFpzxQA4M6jc0I5olC5PuFgNvuVUzX8Bpc4j2NO/u+38gXsO566gK+LWqVuJUX2WIyBacAaPWx1gvkTLR/gG4UkSqcaqnFrp3b30RZ9S5LTjDaIa7lV+O0/bRE/PTjEmA9TprDCAirwHXq2qT37EMNRF5AFijqi/6HYtJXXZlYYzjn4Cpfgfhka2WKMxg2ZWFMcaYuOzKwhhjTFyWLIwxxsRlycIYY0xcliyMMcbEZcnCGGNMXP8/tR+KRWjkGsgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(ins,data10)\n", "plt.ylabel('P(N,s)')\n", "plt.xlabel('s (imbalance)')\n", "plt.title('Gaussian with N=10')\n", "plt.savefig('Gauss10.png')" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEWCAYAAACjYXoKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XecFeXZ//HPtRVYYGkL0nsRFBFWwIYYS2yBJFawYKIxakzPLzHxSTNPqtHEJPqoMRpjwRoVFWOL2CmL0hFY+sqC1KVvvX5/zCwe111Y9pyzs+X7fr0OnDNzz5mL2cN+z8w9c4+5OyIiIvFIiboAERFp/BQmIiISN4WJiIjETWEiIiJxU5iIiEjcFCYiIhI3hYkIYGYvmtmUCNf/EzO79yDzrzSzt+uzJpHDoTCRSJjZJWY2y8z2mNnH4fPrzcyiqMfdz3b3B6JYd7j+37j71QBm1sfM3MzS6vp+ZrbGzDaZWVbMtKvNbMZhvk9nM5tqZhvMrMjM3jGzMVXaTDazteHP8hkz61DXuqXxUphIvTOz7wO3A7cARwBdgGuBE4GMCEtratKAb8f5Hq2BOcAooAPwAPCCmbUGMLNhwN3A5QQ/x73AnXGuUxohhYnUKzPLBm4Grnf3J919lwc+cPdL3b04bHeumX1gZjvNbL2Z/SLmPcabWUGV911jZqeHz0ebWV647CYzuy2c3sLMHjKzrWa2w8zmmFmXcN4MM6vcM+hvZv8N220xs4fNrF2Vdf3AzBaE39YfM7MWNfx715rZqPD5ZeEex9Dw9dVm9kz4/Bdm9lC42Jvh3zvMbLeZHR/zfn80s+1mttrMzj7E5r4F+EFs7YfL3Ve5+23uXuju5e5+D0HgDw6bXAo85+5vuvtu4KfAl82sTV3XKY2TwkTq2/FAJvDsIdrtAa4A2gHnAteZ2RdruY7bgdvdvS3QH3g8nD4FyAZ6Ah0J9ob2VbO8Ab8FugFHhu1/UaXNRcBZQF9gOHBlDbW8AYwPn48DVgGnxLx+o5plxoV/t3P31u7+Xvh6DLAM6AT8AfjHIQ4L5gEzgB9UNzMMwx01PKrduzCzEQRhkh9OGgbMr5zv7iuBEmDQQeqSJkhhIvWtE7DF3csqJ5jZu+EvsH1mNg7A3We4+0J3r3D3BcBUPvklfCilwAAz6+Tuu919Zsz0jsCA8Fv2XHffWXVhd89391fcvdjdNwO3VbPuv7j7BnffBjwHjKihljdilj2ZIKQqX59C9WFSk7Xu/nd3Lyc43NSV4NDSwfwM+KaZ5VSd4e7D3b1dDY/rq7Y3s7bAg8Av3b0onNwaKKrStAjQnkkzozCR+rYV6BTbuezuJ7h7u3BeCoCZjTGz181ss5kVEexFdKrlOq4i+Gb8YXgo67xw+oPAS8CjYYfyH8wsverCYafzo2b2kZntBB6qZt0bY57vJfilWp03gJPN7AggFXgMONHM+hDsJc2r5b/pU+t0973h05rWW9luEfA8cONhrOczzKwlQWjOdPffxszaDbSt0rwtsCue9UnjozCR+vYeUAxMPES7R4BpQE93zwbuIjj8BMEhsFaVDc0sFTjwzdvdV7j7JKAz8HvgSTPLcvdSd/+luw8FTgDOIziUVtVvAQeGh4fKLotZ92Fx93yCsPkW8Ka77yIIhWuAt929orrF6rKug/g58DWge+xEM1sc9slU97grpl0m8AzwEfD1Ku+9GDgmpm0/gsOYyxP8b5AGTmEi9crddwC/BO40swvMrLWZpYTH4rNimrYBtrn7fjMbDUyOmbccaBF20qcD/0PwCww40NGdE/6i3hFOLjezU83s6DB8dhIc9iqvpsw2BN+4d5hZd+D/xfnPfgO4gU8Oac2o8rqqzUAF0C/O9QIHAu0xgkCLnT4s7JOp7nEtQLh9nyToW7qimvB7GPiCmZ1swWnINwP/DkNTmhGFidQ7d/8D8D3gh8DHwCaC00t/BLwbNrseuNnMdhEc9388ZvmicP69BN+W9wCxZ3edBSw2s90EnfGXuPt+gtOQnyQIkqUEv8wf4rN+CYwkOPb/AvDvOP/JbxAE1Js1vP6U8BDWr4F3wr6ksXGuH4Jf8lmHbPVZlXtwZ/LJ2WW7zezksNbFBIcgHyb4WbYh+NlIM2O6OZaIiMRLeyYiIhK3pIaJmZ1lZsvMLN/MPnM2iZl9z8yWhOe7v2ZmvWPmlZvZvPAxLZl1iohIfJJ2mCvs5FwOnEFwPHsOMMndl8S0ORWY5e57zew6YLy7XxzO2+3uBz3tUUREGoZk7pmMBvLD4RhKgEepcjqou78ec778TKBHEusREZEkqfOopLXQHVgf87qAYDiImlwFvBjzuoWZ5QFlwO/c/ZmqC5jZNQTn65OVlTVqyJAhcRctItKczJ07d4u7f2aEhMOVzDCp7iKvao+pmdllQC6fHrKil7tvCC+C+q+ZLQzH/fnkzYJB5+4ByM3N9by8vMRULiLSTJjZ2kS8TzIPcxUQDJBXqQewoWojC0Z6vQmYUDliLIC7bwj/XkVwkdexSaxVRETikMwwmQMMNLO+ZpYBXEIwPMYBZnYswcVqE9z945jp7cMhHDCzTgT3uViCiIg0SEk7zOXuZWZ2A8HAeqnAfe6+2MxuBvLcfRrB/RZaA0+EI2mvc/cJBMN+321mFQSB97vYs8BERKRhaTJXwKvPRETk8JnZXHfPjfd9dAW8iIjETWEiIiJxU5iIiEjcFCYiIhI3hYmIiMRNYSIiInFTmIiISNwUJiIiEjeFiYiIxE1hIiIicVOYiIhI3BQmIiISN4WJiIjETWEiIiJxU5iIiEjcFCYiIhI3hYmIiMRNYSIiInFTmIiISNwUJiIiEjeFiYiIxE1hIiIicVOYiIhI3BQmIiISN4WJiIjETWEiIiJxU5iIiEjcFCYiIhI3hYmIiMRNYSIiInFTmIiISNwUJiIiEjeFiYiIxE1hIiIicVOYiIhI3JIaJmZ2lpktM7N8M7uxmvnfM7MlZrbAzF4zs94x86aY2YrwMSWZdYqISHySFiZmlgrcAZwNDAUmmdnQKs0+AHLdfTjwJPCHcNkOwM+BMcBo4Odm1j5ZtYqISHySuWcyGsh391XuXgI8CkyMbeDur7v73vDlTKBH+PzzwCvuvs3dtwOvAGclsVYREYlDMsOkO7A+5nVBOK0mVwEvHs6yZnaNmeWZWd7mzZvjLFdEROoqmWFi1UzzahuaXQbkArcczrLufo+757p7bk5OTp0LFRGR+CQzTAqAnjGvewAbqjYys9OBm4AJ7l58OMuKiEjDkMwwmQMMNLO+ZpYBXAJMi21gZscCdxMEyccxs14CzjSz9mHH+5nhNBERaYDSkvXG7l5mZjcQhEAqcJ+7Lzazm4E8d59GcFirNfCEmQGsc/cJ7r7NzH5FEEgAN7v7tmTVKiIi8TH3arsxGp3c3FzPy8uLugwRkUbFzOa6e26876Mr4EVEJG4KExERiZvCRERE4qYwERGRuClMREQkbgoTERGJm8JERETipjAREZG4KUxERCRuChMREYmbwkREROKmMBERkbglbdRgkeaqaF8pzy/YwJNzC1i8YeeB6ekpxvjBnTl/VHfGDcwhLVXf5aTpUJiIJMiOvSX86vmlPLdgAyVlFQzu0oYpx/cmNSUIjaJ9pby0eCMvLCwkp00m15zcj6tO6ktKSnU3FhVpXBQmIgkwZ802vj31AzbvLmby6F5cmNuTYd3aEt6n54BfThjG68s+5uFZ6/j19KW8uWIzt100gpw2mRFVLpIYup+JSBwqKpw7Z+Tzp1dX0KN9S/466ViG92h3yOXcnamz1/PL5xbTpkU6f754BCcN7FQPFYt8mu5nIhIxd+emZxbxx5eXc+7RXXn+myfVKkgAzIzJY3rx7A0n0q5VOlfcN4uXFm9McsUiyaMwEamjW15axtTZ67h+fH9uv2QEbVqkH/Z7DDmiLc9+40SG92jHN6d+wLsrtyShUpHkU5iI1MHf31zFnTNWMnlML/7f5wd/pm/kcGRlpnH/lcfRu0MrrvnXXBYWFCWwUpH6oTAROUxPzS3g19OXcu7RXfnVxKPiCpJK7bMyePCqMWS3TGfK/bNZvWVPAioVqT8KE5HDsHzTLn7y9EJO6N+R2y4+htQEntZ7RHYLHrp6DO7ODY+8T3FZecLeWyTZFCYitbS/tJxvTf2A1plp3H7JsWSmpSZ8HX07ZfH784ezeMNObn15ecLfXyRZFCYitfSH/yzjw427uOXC4Um9LuTMYUcweUwv7nlzFW+vUIe8NA4KE5FaeGP5Zu57ZzVTju/N54Z0Sfr6fnruUPrlZPH9J+axfU9J0tcnEi+Ficgh7Nhbwg+emM+gLq358TlH1ss6W2ak8pdLjmXbnhJ+8vTCelmnSDwUJiKHcNsry9m6u5jbLhpBi/TE95PU5Kju2Xz7tIG8uGgjbyzfXG/rFakLhYnIQSzZsJOHZq7l8rG9Oap7dr2v/2vj+tGnYyt+OW0xJWUV9b5+kdpSmIjUwN35+bRFtGuVwffOGBxJDZlpqfz8C8NYtWUP972zOpIaRGpDYSJSg2fnbWDOmu388PODyW51+EOlJMqpQzpz+pGd+etrK9i0c39kdYgcjMJEpBq7i8v4zfSlHNMjm4tye0ZdDj89byilFc5vpy+NuhSRailMRKrxfzPy+XhXMb+YMKxB3Lyqd8csvj6uH8/M28D767ZHXY7IZyhMRKrYsruY+99Zw3nDu3Jsr/ZRl3PAdeP706l1BrfpynhpgBQmIlXcNWMl+0vL+e4Zg6Iu5VNaZaRx3fgBvJ2/hfdWbo26HJFPUZiIxNhYtJ8HZ67lyyN70D+nddTlfMalY3pxRNsW3PbKMprKXVKlaUhqmJjZWWa2zMzyzezGauaPM7P3zazMzC6oMq/czOaFj2nJrFOk0h2v51Ne4Xz7tIFRl1KtFump3PC5AcxZs503NW6XNCBJCxMzSwXuAM4GhgKTzGxolWbrgCuBR6p5i33uPiJ8TEhWnSKV1m/by6Nz1nHxcT3p2aFV1OXU6KLcnvRo35JbX9beiTQcydwzGQ3ku/sqdy8BHgUmxjZw9zXuvgDQpb0Sub/+dwVmxg2fGxB1KQeVkZbCt08byIKCIl5esinqckSA5IZJd2B9zOuCcFpttTCzPDObaWZfrK6BmV0TtsnbvFljF0ndFWzfy1Pvf8Tk0b3omt0y6nIO6UvHdqdvpyz+9t987Z1Ig5DMMKnu5PzD+dT3cvdcYDLwZzPr/5k3c7/H3XPdPTcnJ6eudYpw71urMeCacf2iLqVW0lJT+Pq4fiz8qIh3dWaXNADJDJMCIPbS4R7Ahtou7O4bwr9XATOAYxNZnEilbXtKeHTOOiaO6E63dg1/r6TSl0Z2p3ObTP5vxsqoSxFJapjMAQaaWV8zywAuAWp1VpaZtTezzPB5J+BEYEnSKpVm7YF317C/tIJrT2kceyWVMtNSueqkvrydv4WFBUVRlyPNXNLCxN3LgBuAl4ClwOPuvtjMbjazCQBmdpyZFQAXAneb2eJw8SOBPDObD7wO/M7dFSaScHtLynjgvTWcfmQXBnZpE3U5h23ymF60aZHGXW9o70SilZbMN3f36cD0KtN+FvN8DsHhr6rLvQscnczaRAAenb2eHXtLuW5849orqdSmRTqXje3N3W+sZM2WPfTplBV1SdJM6Qp4abZKyyu4961VjO7TgVG9O0RdTp195cQ+pKWmcM9bq6IuRZoxhYk0Wy8sKGRD0X6+3sj6Sqrq3KYFF4zqwZNzC9i6uzjqcqSZUphIs+Tu3P/Oavp1yuLUwZ2jLiduXz2xLyVlFTwya13UpUgzpTCRZun9dTuYX1DElSf2aRD3K4nXgM6tGTcohwdnrtW94iUSChNplu5/ZzVtWqRx/sjPnP/RaH3lxD58vKuYFxcVRl2KNEO1ChMz62FmPzCzZ81sjpm9aWZ3mtm5ZqZAkkalsGgfLy7ayMW5PcnKTOoJjfXqlIE59OuUxX3vrIm6FGmGDhkEZnY/cB9QAvwemARcD7wKnAW8bWbjklmkSCI9+N5a3J0pJ/SJupSESkkxrjyxD/PX79CtfaXe1eZr2a3uvqia6YuAf4dXt/dKbFkiybG/tJyps9dx+pFdGvQw83V1/sge3PLSMu5/Zw0jG9Ath6XpO+SeSXVBEg53MjycX+Lu+ckoTiTRnvngI7bvLeUrJ/aNupSkyMpM4+Lcnry4sJDCon1RlyPNSK37O8xshpm1NbMOwHzgfjO7LXmliSSWu/PPd9cw5Ig2jO3XeC9SPJQpJ/Sh3J2pOk1Y6tHhdJ5nu/tO4MvA/e4+Cjg9OWWJJN7767bz4cZdXHF8H8wa/+nANenZoRWnDu7M1DnrdZqw1JvDCZM0M+sKXAQ8n6R6RJLmwffW0jozjYkjukVdStJdPrY3m3cV8/KSjVGXIs3E4YTJzQQjAOe7+xwz6wesSE5ZIom1dXcx0xdu5PyR3ZvU6cA1GTcoh54dWvLQzLVRlyLNRK3DxN2fcPfh7n59+HqVu5+fvNJEEufxvAJKyiu4bGzvqEupF6kpxqVjejNz1TZWbNoVdTnSDNT5gkMzu97MLjazpv81Txq18grn4VlrGduvQ6O8Z0ldXTiqBxmpKdo7kXoRz9XrBpwE/DtBtYgkxRvLP6Zg+z4uH9sn6lLqVcfWmZw7vCtPvf8Re4rLoi5Hmrg6h4m73+Hu33T3CYksSCTRHnxvLTltMjlzWJeoS6l3l43tze7iMp6Z91HUpUgTd8hDVGZ2xUFmu7s/mMB6RBJq/ba9zFi+mRtOHUB6avMbRm5kr3Yc2bUtD81cx+TRvZr0KdESrdr0dxxXzTQDvgB0BxQm0mA9Nmc9BlwyunmO+GNmTB7Ti58+s4j5BUWM6Nku6pKkiarNcCrfrHwA3wJmAacAM4GRSa5PpM5Kyyt4LG894wd3pnu7llGXE5kvjuhGq4xUHpmljnhJntoOQZ9mZlcDSwiuer/A3S929wVJrU4kDq8t3cTmXcVMbqZ7JZXatEhnwjHdeG5+ITv3l0ZdjjRRtRmC/hsEITIKOMvdr3T3ZUmvTCROD89aR9fsFowfnBN1KZGbPKYX+0rLeeYDdcRLctRmz+SvQFuC04CfM7MF4WOhmWnPRBqkdVv38taKLVx8XE/SmmHHe1XDe7TjqO5teWTWOtw96nKkCapNB3zTHKtbmrSpc9aRYnDxcT2jLqXBmDy6Nz95eiHvr9vBqN6614kkVm2+sq1z97U1PQBM5xtKA1JSVsETeev53JAudM1uvh3vVU0Y0Y2sjFSmztbQ9JJ4tQmT183sm2b2qV5MM8sws8+Z2QPAlOSUJ3L4Xl26iS27S5g8RnslsVpnpjHx2O48v2ADRXvVES+JVZswOQsoB6aaWaGZLTGzVQQjBk8C/uTu/0xijSKHZersdXTLbsEpgzpHXUqDM3l0L/aXVuiKeEm4Q/aZuPt+4E7gTjNLBzoB+9x9R7KLEzlclR3v3z19EKkpOvpa1VHdsxneI5ups9dxxfG9dUW8JExtTg1uYWbfMbO/AV8BNitIpKF6NOx4v+i4HlGX0mBNGt2LDzfu4oP1+m8siVObw1wPALnAQuAc4NakViRSR6XlFTyeV8DnhnRWx/tBfOGYsCNe94iXBKpNmAx198vc/W7gAuDkJNckUievLd3Elt3FTGrmV7wfSuvMNCaM6MZzCzboinhJmNqEyYFPm7vrpgjSYD0yez1ds1twyiBd8X4ok8KO+Gd1RbwkSG3C5Bgz2xk+dgHDK5+b2c5kFyhSG+u37eWtFZu5KFdXvNfG0d2zGdatLY/MXq8r4iUhajNqcKq7tw0fbdw9LeZ52/ooUuRQHs8Lhpq/SFe814qZMWl0L5YW7mRBQVHU5UgToK9w0uiVllfw2Jz1nDIop1kPNX+4Jo7oRsv0VB5RR7wkQFLDxMzOMrNlZpZvZjdWM3+cmb1vZmVmdkGVeVPMbEX40BX2UqP/fvgxH+8qZvKY3lGX0qhUDk0/bf4GdqkjXuKUtDAxs1TgDuBsYCgwycyGVmm2DrgSeKTKsh2AnwNjgNHAz81MI9NJtR6ZtY4j2rbgVA01f9gmVQ5NP29D1KVII5fMPZPRQL67r3L3EuBRYGJsA3dfE95gq6LKsp8HXnH3be6+HXiFYFgXkU9Zv20vb67YzEUaar5OjumRzdCuGppe4pfM/33dgfUxrwvCaQlb1syuMbM8M8vbvHlznQuVxqvyHu8aar5uKu8Rv7RwJ/PVES9xSGaYVDfoT22/+tRqWXe/x91z3T03J0eHOJqb4Ip33eM9XhPDe8TriniJRzLDpACI/brYA6jtgdl4lpVm4rWlYce7rniPS2xHvK6Il7pKZpjMAQaaWV8zywAuAabVctmXgDPNrH3Y8X5mOE3kgKmzg4533eM9fpNGBx3xz6ojXuooaWESDr1yA0EILAUed/fFZnazmU0AMLPjzKwAuBC428wWh8tuA35FEEhzgJvDaSLAJx3vusd7YgzvEVwR//DMteqIlzqpzT3g68zdpwPTq0z7WczzOQSHsKpb9j7gvmTWJ43XI7PXqeM9gSo74m96epHuES91oq900ugUl5Xz+Jz1nHZkF7qp4z1hvjiiO60z03h45tqoS5FGSGEijc5/Fm1k654SLhurK94TKSszjS8d253nFxayfU9J1OVII6MwkUbn4Znr6N2xFScP6BR1KU3OZWN7U1JWwRNz1x+6sUgMhYk0Kss27mL2mm1MHt2LFN3jPeEGH9GG4/q05+FZ66ioUEe81J7CRBqVh2etJSMthQtz1fGeLJeN7c3arXt5O39L1KVII6IwkUZjT3EZ/37/I849uisdsjKiLqfJOuuoI+iYlcHDs9QRL7WnMJFGY9r8DewuLuOysbriPZky01K5MLcnry79mMKifVGXI42EwkQaBXfngXfXMOSINozspWsgku3SMb2ocNeNs6TWFCbSKMxevY0PN+7iyhP6YKaO92Tr2aEVpw3pzNTZ6yguK4+6HGkEFCbSKPzrvbVkt0xn4oja3sVA4jXlhD5s2V3C9IWFUZcijYDCRBq8wqJ9/GfxRi4+rictM1KjLqfZOLF/J/rlZPHPd9URL4emMJEG7+GZ66hw53Jd8V6vUlKMKcf3Yf76HcxbvyPqcqSBU5hIg7a/tJyps9dx2pAu9OzQKupymp3zR/WgdWYaD7y7JupSpIFTmEiDNn1hIVv3lHDlCX2iLqVZap2ZxgWjevDCgkI27yqOuhxpwBQm0mBVng7cPyeLEwd0jLqcZuvy43tTUl7B1Nk6TVhqpjCRBmvu2u3MLyjS6cAR65/TmnGDcnhw5lqdJiw1UphIg3XvW6vJbpnO+aOqvX+a1KOrT+rL5l3FPDdfpwlL9RQm0iCt3bqHl5Zs5LKxvWiVkdQbgkotnDywE4O7tOHet1bptr5SLYWJNEj3v7OGtBTjiuP7RF2KENzW96qT+/Lhxl28u3Jr1OVIA6QwkQanaF8pj+et5wvHdKNL2xZRlyOhiSO60al1Jve+tSrqUqQBUphIgzN19jr2lpRz1Ul9oy5FYmSmpXLF8b15fdlm8j/eFXU50sAoTKRBKS2v4J/vrOGE/h0Z1i076nKkikvH9CIzLYV/vL066lKkgVGYSIPy/IINbNy5n6tP1l5JQ9SxdSZfHtmDp97/SBcxyqcoTKTBqKhw/m/GSgZ3acP4QZ2jLkdqcM24fpSVV3DfO9o7kU8oTKTBeHXpJpZv2s114/uTkqKLFBuqvp2yOPvorjz43lqK9pVGXY40EAoTaRDcnTtnrKRnh5acN7xr1OXIIVw/vj+7i8t4aKaGp5eAwkQahPdWbWXe+h18fVx/0lL1sWzohnXLZvzgHP7x9mr2lWiIFVGYSANx5+sryWmTyQUaOqXR+MapA9i2p4TH5mgASFGYSAMwf/0O3s7fwtUn9aVFuu6k2Fgc16cDx/Vpzz1vrqKkrCLqciRiChOJ3B2v59O2RRqX6k6Kjc71pw5gQ9F+nv6gIOpSJGIKE4nUwoIiXl6yiatO6kfrTA3o2NiMH5TD8B7Z/OW1fO2dNHMKE4nUn15dTrtW6Xz1pD5RlyJ1YGZ894xBfLRjH0/MXR91ORIhhYlE5v112/nvhx9zzbh+tGmRHnU5UkfjB+Uwslc7/vbffPaX6syu5kphIpH50yvL6ZiVwRQNM9+omRnfP3MwhUX7eVS39m22khomZnaWmS0zs3wzu7Ga+Zlm9lg4f5aZ9Qmn9zGzfWY2L3zclcw6pf7NXr2Nt1Zs4brx/clSX0mjd0L/jozp24E7ZqzUdSfNVNLCxMxSgTuAs4GhwCQzG1ql2VXAdncfAPwJ+H3MvJXuPiJ8XJusOqX+uTu3vryMzm0yuUxncDUJlXsnm3cV66r4ZiqZeyajgXx3X+XuJcCjwMQqbSYCD4TPnwROMzMNytTEvbF8M7NWb+Mbpw7QdSVNyOi+HTh5YCfunJGvMbuaoWSGSXcg9vSOgnBatW3cvQwoAjqG8/qa2Qdm9oaZnVzdCszsGjPLM7O8zZs3J7Z6SYqy8gp+M30pfTq2YtLoXlGXIwl249lD2LGvlDtfz4+6FKlnyQyT6vYwvJZtCoFe7n4s8D3gETNr+5mG7ve4e6675+bk5MRdsCTfE3MLWL5pNzeePYSMNJ3/0dQM65bNl4/twf3vrGH9tr1RlyP1KJn/mwuAnjGvewAbampjZmlANrDN3YvdfSuAu88FVgKDklir1IPdxWXc+vJyjuvTns8POyLqciRJfvD5QaSkwB9eWhZ1KVKPkhkmc4CBZtbXzDKAS4BpVdpMA6aEzy8A/uvubmY5YQc+ZtYPGAisSmKtUg/ueWMlW3YXc9O5Q1HXWNPVNbslXzu5H8/N38AH67ZHXY7Uk6SFSdgHcgPwErAUeNzdF5vZzWY2IWz2D6CjmeUTHM6qPH14HLDAzOYTdMxf6+7bklWrJF9h0T7ueWsVE47pxoie7aIuR5Ls66f0p1PrTH79wlLcqx7dlqYoqSf4u/t0YHqVaT+Leb4fuLCa5Z4CnkpmbVK/fjP9Qyoc/t/nB0dditSD1plpfO+MQfzk6YVMm7+BiSOqnnsjTY16QCXp3lqxmefeweuhAAAS40lEQVTmb+D68f3p2aFV1OVIPbn4uJ4M75HN/76wlJ37dapwU6cwkaTaX1rOz55dTJ+Orbj2lP5RlyP1KDXF+PUXj2br7mJuVWd8k6cwkaS6+41VrN6yh5snHqULFJuho3tkc/nY3jw4cy0LC4qiLkeSSGEiSbNmyx7umJHPucO7Mm6QrgNqrr535mA6ZGVy0zMLKa9QZ3xTpTCRpHB3fvrsIjJSU/jZeVWHZJPmJLtlOj8970gWFBRp3K4mTGEiSTF19nreWrGFH541mC5tW0RdjkRswjHdGDcoh9+9+CFrtuyJuhxJAoWJJNy6rXv53xeWcOKAjlw2RqMCSzCq8O/PP5q0VOMHT8zX4a4mSGEiCVVe4fzgifmkmnHLBceQkqIr3SXQNbslN08cRt7a7fz9LQ1o0dQoTCSh7nt7NbPXbOPnE4bRrV3LqMuRBuaLI7pz1rAjuO3l5Xy4cWfU5UgCKUwkYZYW7uSWl5dxxtAunD9SVzzLZ5kZv/7SUbRtmcZ3H5uve8Y3IQoTSYid+0u5/uH3adcynd9++WgN5Cg16tg6kz9cMJylhTv55XOLoy5HEkRhInFzd3705ALWbdvL3yaPpFPrzKhLkgbuc0O6cP34/kydvZ4n5xZEXY4kgMJE4vaPt1fz4qKN/OiswYzu2yHqcqSR+N4Zgzi+X0duenohSwvVf9LYKUwkLnlrtvG7Fz/kzKFd+NrJ/aIuRxqRtNQU/jLpWLJbpnPdQ3M1GGQjpzCROlu7dQ/XPDiXHu1bcsuFx6ifRA5bTptM7rh0JAXb93H9Q+9TUlYRdUlSRwoTqZNte0q48v45uDv3XXkc2S3Toy5JGqnj+nTgd+cP5+38Lfz43wt1M61GKqk3x5KmaX9pOVc/MIePduxj6tfG0C+nddQlSSN3wagefLR9H396dTk92rfku2cMirokOUwKEzksZeUVfOfReXywfgd3Th7JqN7qcJfE+NZpAyjYvpfbX1tBt3YtuPi4XlGXJIdBYSK1VlZewXcfn89/Fm/kp+cN5eyju0ZdkjQhZsZvvnw0m3YVc+O/F5JixoW5PaMuS2pJfSZSK2XlFXznsXk8N38DPz57CFed1DfqkqQJSk9N4Z7LR3HSgE788KkFPJ63PuqSpJYUJnJIZeUVfPuxeTy/oJCfnDOEr+v2u5JELdJT+fsVuZw0oBM/emoBj89RoDQGChM5qF37S/nqA3m8sKCQm845kmvGKUgk+SoD5eSBOfzwqQXc8Xq+zvJq4BQmUqOPduzjwrve4938Lfzuy0fztXG6KFHqT4v0VO65fBQTR3TjlpeW8aOnFug6lAZMHfBSrfnrd3D1v/LYX1LOP78ympMGdoq6JGmGWqSn8ueLR9C7YxZ/eW0FBdv3ceelI2nXKiPq0qQK7ZnIp7g79729mgvveo+M1BSeuv4EBYlEysz43hmDuPXCY5izZhvn/uVt8tZsi7osqUJhIgds3V3MVQ/kcfPzSxg3qBPPffMkBnVpE3VZIgCcP6oHT1x7AqkpxsX3zOSvr63Q7X8bEIWJ4O78Z1EhZ9/+Fm+v2MIvvjCUv1+RS4csHUqQhmVEz3a88K2TOPfortz6ynIm3TOT/I93RV2WoDBp9tZv28tVD+Rx7UPv07F1Jk9/4wSuPLGvBm2UBqtNi3Ruv2QEf7zwGD7cuJOzb3+LP760THdtjJg64JupnftLufet1dzz5kpSzPifc4/kyhP6kJaq7xfS8JkZF4zqwfjBOfzmhaX87fV8ps3fwPfPHMQXhncjJUVfhuqbNZVzt3Nzcz0vLy/qMhq8vSVl/PPdNdz9xiqK9pVyztFH8D/nDqVbu5ZRlyZSZ+/mb+GXzy1h2aZdDO7Shu+dOYgzh3bRHnYtmNlcd8+N+30UJs3DRzv28dDMtTw6ex3b95Zy6uAcvn/mYI7qnh11aSIJUVHhPL+wkD+/spxVW/YwuEsbppzQhy8e241WGToIUxOFSRUKk88qKavgrRWbeSKvgJeXbATgjKFduGZcP432K01WWXkFz8zbwH1vr2ZJ4U7atkjjwtyefOnY7gzr1lZ7K1UoTKpQmASKy8qZs3o70xcVMn1hITv2ltK+VToXHdeTy8f2pkf7VlGXKFIv3J28tdv55ztreGnxRsoqnAGdWzPhmG6cMbQLQ45oo2BBYfIZzTVMyiuc5Zt2kbdmG28s38y7K7eyt6SclumpnDmsCxNHdOOkATlkpKljXZqv7XtKmL6okGfnbWD26uCCxyPatuCUQTmcNLATuX3a0zW7efYbKkyqaA5hUl7hrNm6h6WFO1lauJOFH+3kg7Xb2VVcBkDPDi0ZP6gz4wfncHz/jjpOLFKNTTv388ayzcxY/jFvLd9y4P9P93YtObZXO4Z1y+bIrm0Y2rUtOW0ym/zeS6MIEzM7C7gdSAXudfffVZmfCfwLGAVsBS529zXhvB8DVwHlwLfc/aWDrasphElpeQXb9pSwaed+NuzYT2HRPjbs2MfqLXtZs3UP67bupaQ8GOguLcUY0Lk1I3u3J7d3e3J7d6Bnh5ZN/oMvkkhl5RUsKdzJ3LXbyVu7nXnrdvDRjn0H5rdpkUbfTln07ZRF7w6t6NquJV2zW9CtXUtyWmeS3TK90Z+G3ODDxMxSgeXAGUABMAeY5O5LYtpcDwx392vN7BLgS+5+sZkNBaYCo4FuwKvAIHev8aqkRISJu+MODlS4UxG+rnCnvMKpqIDy8Hl5hVNWUUF5hVNa7pSWV1BW7pSUl1NcVhE8SisoLitnX0k5e0vK2VtSxu7icvYUl7FrfylF+4LHjn2lbNtTwo69pZ+pKSMthT4dW9GnY/CBHtC5NUd2bcvALq3JTEuN698rIp9VtLeUpRuDvf9Vm/ewZuseVm3eQ2HRPqqO3pKWYnTIyqBDVgZtW6bTrmU6bVum0zozjdaZaWRlptEqI5WWGam0TE+lRXoqGWkpZKalkJGWQkZqCumpKaSlGukpKaSmGmkpRmqKkWpGSvg8xSDFjBQz7MBzEvLlMVFhkszjIKOBfHdfBWBmjwITgSUxbSYCvwifPwn8zYKtMxF41N2LgdVmlh++33s1rWzRhiKG/PRFAKrmox/4A5xPAqMySL2aZZIlPdXICj9o7Vqlk90ynSOyW9AxK5OOrTPo1DqTnDaZdMtuSdd2LeiYlaG9DZF6lN0qnbH9OjK2X8dPTS8tr+DjXcUU7tjHhqL9bNlVzNY9xWzZVcL2vSXs2FfK2q172bm/lN3FZewpLvtM+CSLGRhBuFj4GsAIZ3DgrwPtEy2ZYdIdiL1FWgEwpqY27l5mZkVAx3D6zCrLdq+6AjO7BrgmfFm87H/PWZSY0pOqE7Al6iJqQXUmlupMrMZQZ2OoEWBwIt4kmWFSXfZVzema2tRmWdz9HuAeADPLS8SuWrKpzsRSnYmlOhOnMdQIQZ2JeJ9kni9aAPSMed0D2FBTGzNLA7KBbbVcVkREGohkhskcYKCZ9TWzDOASYFqVNtOAKeHzC4D/etCRMQ24xMwyzawvMBCYncRaRUQkDkk7zBX2gdwAvERwavB97r7YzG4G8tx9GvAP4MGwg30bQeAQtnucoLO+DPjGwc7kCt2TrH9LgqnOxFKdiaU6E6cx1AgJqrPJXLQoIiLR0RgbIiISN4WJiIjErVGFiZldaGaLzazCzHKrzPuxmeWb2TIz+3wNy/c1s1lmtsLMHgtPDEh2zY+Z2bzwscbM5tXQbo2ZLQzb1fu4MGb2CzP7KKbWc2pod1a4jfPN7MYI6rzFzD40swVm9rSZtauhXSTb81DbJzyp5LFw/iwz61NftYXr72lmr5vZ0vD/0reraTPezIpiPgs/q88aY+o46M/QAn8Jt+UCMxsZQY2DY7bTPDPbaWbfqdImku1pZveZ2cdmtihmWgczeyX8HfiKmbWvYdkpYZsVZjalujafEQwh0jgewJEEF9jMAHJjpg8F5gOZQF9gJZBazfKPA5eEz+8Crqvn+m8FflbDvDVApwi37S+AHxyiTWq4bfsBGeE2H1rPdZ4JpIXPfw/8vqFsz9psH+B64K7w+SXAY/VcY1dgZPi8DcGQR1VrHA88X5911eVnCJwDvEhwXdpYYFbE9aYCG4HeDWF7AuOAkcCimGl/AG4Mn99Y3f8foAOwKvy7ffi8/aHW16j2TNx9qbsvq2bWgeFX3H01UDn8ygHhMC2fIxi2BeAB4IvJrLea9V9EMOZYY3VgiBx3LwEqh8ipN+7+sruXhS9nElyD1FDUZvtMJPjsQfBZPM3qcbwcdy909/fD57uApVQzukQjMRH4lwdmAu3MrGuE9ZwGrHT3tRHWcIC7v0lwlmys2M9fTb8DPw+84u7b3H078Apw1qHW16jC5CCqG7ql6n+QjsCOmF9E1Q7RkkQnA5vcfUUN8x142czmhsPEROGG8HDBfTXs/tZmO9enrxJ8M61OFNuzNtvnU0MIAZVDCNW78BDbscCsamYfb2bzzexFMxtWr4V94lA/w4b2ebyEmr8sNoTtCdDF3Qsh+GIBdK6mTZ22a4O74YWZvQocUc2sm9z92ZoWq2ZabYduiVsta57EwfdKTnT3DWbWGXjFzD4Mv1kkzMHqBP4P+BXBNvkVwSG5r1Z9i2qWTfi55bXZnmZ2E8E1SA/X8DZJ357ViPRzeDjMrDXwFPAdd99ZZfb7BIdqdod9Z88QXDhc3w71M2wQ2xIg7H+dAPy4mtkNZXvWVp22a4MLE3c/vQ6L1Wb4lS0Eu8Fp4TfChA3RcqiaLRgq5ssE922p6T02hH9/bGZPExwySegvv9puWzP7O/B8NbPqZZibWmzPKcB5wGkeHuSt5j2Svj2rcThDCBXYp4cQqjdmlk4QJA+7+7+rzo8NF3efbmZ3mlknd6/XQQtr8TNsSMMunQ287+6bqs5oKNsztMnMurp7YXhI8ONq2hQQ9PNU6kHQT31QTeUw1yGHXwl/6bxOMGwLBMO41LSnk2inAx+6e0F1M80sy8zaVD4n6GSu1xGQqxxr/lIN66/NEDlJZcEN134ETHD3vTW0iWp7xjOEUL0I+2f+ASx199tqaHNEZT+OmY0m+D2xtb5qDNdbm5/hNOCK8KyusUBR5SGcCNR45KEhbM8YsZ+/mn4HvgScaWbtw8PdZ4bTDq6+zzCI50HwS64AKAY2AS/FzLuJ4EyaZcDZMdOnA93C5/0IQiYfeALIrKe6/wlcW2VaN2B6TF3zw8digsM59b1tHwQWAgvCD1zXqnWGr88hOANoZUR15hMcz50XPu6qWmeU27O67QPcTBB+AC3Cz15++FnsV8/b7ySCQxYLYrbhOcC1lZ9R4IZwu80nOMnhhAh+ztX+DKvUacAd4bZeSMwZnvVcayuCcMiOmRb59iQIt0KgNPy9eRVB/9xrwIrw7w5h21yCu+FWLvvV8DOaD3ylNuvTcCoiIhK3pnKYS0REIqQwERGRuClMREQkbgoTERGJm8JERETipjARAczsz2Y2Lnx+r5kNPYxlrzSzvx3m+taYWafDrbMuzOzVmkaHFUkUhYk0e2bWARjr4VAd7n61uy+JuKxEepBgtGKRpFGYSJMWXkn9QjjI3iIzu7iaZhcA/4lZZoaF98sxs91m9vtwsMFXzWx0OH+VmU2IeY+eZvYfC+5l8vOY93omXHZxTQNO1tQmXPevw9pnmlmXcHoXC+7lMj98nBBOv8zMZltwz4y7zSw1fKtpBFdoiySNwkSaurOADe5+jLsfRUxoxDgRmFvD8lnADHcfBewC/hc4g2A0hptj2o0GLgVGABfaJzdv+2q4bC7wLTOrboTgmtpkATPd/RiCMam+Fk7/C/BGOH0ksNjMjgQuJhgccQRQHtaDB8OIZ9awbpGEUJhIU7cQOD3cuzjZ3YuqadMV2FzD8iV8EkALCX6Jl4bP+8S0e8Xdt7r7PuDfBMOWQBAOlcNo9KT60WJralPCJwNuzo1Z3+cIRnnG3cvDf9NpBAOJzrHgbp6nEQxJUuljgiFnRJKiwY0aLJJI7r7czEYRjD/1WzN72d1vrtJsH8GYWdUp9U/GHKogGBcOd68IR/09sKqqqzaz8QSDfB7v7nvNbEbV9RyiTey6yzn4/1cDHnD36oZAJ3zPfQdZXiQu2jORJs3MugF73f0h4I8Eh4WqWgoMiHNVZ1hwf+2WBHeve4dgePntYUgMIbi1bFW1aVPVa8B1AGaWamZtw2kXWHDvj8p7ffcOnxvB/WHWxPUvFDkIhYk0dUcDs8NDPzcR9HlU9QKfvn9DXbxNcNbUPOApd88jODyWZmYLCG44NrOa5WrTpqpvA6ea2UKCw1/DwrPP/ofgzoQLCG61WnlbgVEEfS9l1b6bSAJo1GARwMzeBs5z9x1R15JoZnY7MM3dX4u6Fmm6tGciEvg+0CvqIpJkkYJEkk17JiIiEjftmYiISNwUJiIiEjeFiYiIxE1hIiIicVOYiIhI3P4/ZHQae1BFCH8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(ins,data20)\n", "plt.ylabel('P(N,s)')\n", "plt.xlabel('s (imbalance)')\n", "plt.title('Gaussian with N=20')\n", "plt.axis([-10,10,0,0.25])\n", "plt.savefig('Gauss20.png')" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEWCAYAAACjYXoKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xl8VXed//HXhwQChCVAwhYChKWlUKBASldaShdpq22traVarVqtWqnOQ/3NOOP8Rgd/s6gzjjpjp61ara3aTavYxUpboKWFlrBvBUIIJARIAiQsgayf3x/n0N5ebyBwcrO+n4/HfXDvOd9zzyfnhvvO+Z5zvsfcHRERkSi6tXUBIiLS8SlMREQkMoWJiIhEpjAREZHIFCYiIhKZwkRERCJTmEiXZWYvmtndbbj+fzCzn51i/qfMbFlr1iRythQm0mrMbJ6ZvWVmx8ysLHx+n5lZW9Tj7te7+6Ntse5w/f/q7p8FMLPRZuZmlnq272dmRWa238zSY6Z91syWnMV7fcfMNphZvZl9O8H8j5nZrvCz/IOZDYyZN9DMng3n7TKzj53tzyQdh8JEWoWZfQ34EfB9YCgwBPgCcBnQow1L62xSga+0wPsUAH8LPB8/w8wmAQ8BnyD4HKuBB2Ka/ASoDed9HPjfcBnpxBQmknRm1h9YANzn7s+4+xEPrHH3j7t7TdjuRjNbY2aHzaw49i9iM5ttZiVx71tkZteEz2eaWX647H4z+0E4vaeZPW5mB8ys0sxWmtmQcN4SMzu5ZzDWzF4N21WY2a/NLCNuXV83s/VmVmVmT5pZzyZ+3l1mNiN8fle4xzExfP1ZM/tD+PzbZvZ4uNhr4b+VZnbUzC6Jeb//MLNDZrbTzK4/zeb+PvD12NrPhrs/6u4vAkcSzP448Cd3f83djwL/F7jVzPqGe0UfAf6vux9192XAQoLgkU5MYSKt4RIgDfjjadodAz4JZAA3Al80s1uauY4fAT9y937AWOCpcPrdQH8gBxhEsDd0PMHyBvwbMBw4L2z/7bg2HwXmArnAFOBTTdSyFJgdPr8CKASujHm9NMEyV4T/Zrh7H3dfHr6+CNgKZALfA35+mm7BfGAJ8PVEM8MwrGzi8UCiZRKYBKw7+cLddxDsiZwTPhrcfVtM+3XhMtKJKUykNWQCFe5ef3KCmb0ZfoEdN7MrANx9ibtvcPdGd18P/Jb3voRPpw4YZ2aZ4V/EK2KmDwLGuXuDu69y98PxC7t7gbsvcvcady8HfpBg3T9291J3Pwj8CbigiVqWxiw7iyCkTr6+ksRh0pRd7v5Td28AHgWGEXQfnco/AfebWVb8DHef4u4ZTTzua2ZNfYCquGlVQN/TzJNOTGEireEAkBl7cNndL3X3jHBeNwAzu8jMFptZuZlVEexFZDZzHfcQ/FX8TtiV9cFw+mPAS8ATZlZqZt8zs+7xC5vZYDN7wsz2mNlh4PEE694X87ya4IszkaXALDMbCqQATwKXmdlogr2ktc38md63TnevDp82td6T7TYCzwHfOIP1nImjQL+4af0IusRONU86MYWJtIblQA1w82na/Yagfz3H3fsDDxJ0P0HQBdb7ZEMzSwHe/cvb3be7+53AYOC7wDNmlu7ude7+z+4+EbgU+CBBV1q8fwMcmBJ2ld0Vs+4z4u4FBGHzZeA1dz9CEAr3AsvcvTHRYmezrlP4FvA5IDt2opltCo/JJHo82Mz33gRMjXnPMQTdmNvCR6qZjY9pPzVcRjoxhYkknbtXAv8MPGBmt5lZHzPrZmYXAOkxTfsCB939hJnNBGJPKd0G9AwP0ncH/pHgCwx490B3VvhFXRlObjCzq8xschg+hwm6vRoSlNmX4K/qSjPLBv5PxB97KTCf97q0lsS9jlcONAJjIq4XeDfQniQItNjpk8JjMokeXzjZzsy6hycYdCMIh57hNgT4NfAhM5sVHnBfAPw+PLHiGPB7YIGZpZvZZQR/RDzWEj+XtF8KE2kV7v494KsEp5uWAfsJTi/9O+DNsNl9BF9CRwj6/Z+KWb4qnP8zYA/Bnkrs2V1zgU1mdpTgYPw8dz9BcBryMwRBsoXgy/xx/to/A9MJ+vefJ/hCjGIpQUC91sTr9wm7sP4FeCM8lnRxxPVD8CWfftpWif2U4ESFO4Fvhs8/Eda6iaAL8tcEn2Vfgs/mpPuAXuG83wJfDJeRTsx0cywREYlKeyYiIhJZUsPEzOaa2VYzKzCzvzqzxMy+amabw3PfXzGzUTHzGsxsbfhYmMw6RUQkmqR1c4UH67YB1xL0ba8E7nT3zTFtrgLecvdqM/siMNvd7wjnHXX3U54CKSIi7UMy90xmAgXuXujutcATxJ0a6u6LY86dXwGMSGI9IiKSJGc9QmkzZAPFMa9LCIaGaMo9wIsxr3uaWT5QD/y7u/8hfgEzu5fg3H3S09NnTJgwIXLRIiJdyapVqyrc/a9GSzhTyQyTRBd8JexTM7O7gDzeP3zFSHcvDS+IetXMNoRjAL33Zu4PAw8D5OXleX5+fstULiLSRZjZrpZ4n2R2c5UQDJZ30gigNL6RBaO+fhO46eTosQDuXhr+W0hwwde0JNYqIiIRJDNMVgLjzSzXzHoA8wiGyniXmU0juHDtJncvi5k+wMzSwueZBPe82IyIiLRLSevmcvd6M5tPMMheCvCIu28yswVAvrsvJLj3Qh/g6XBU7d3ufhPBEOAPmVkjQeD9e+xZYCIi0r50mivgdcxEROTMmdkqd8+L+j66Al5ERCJTmIiISGQKExERiUxhIiIikSlMREQkMoWJiIhEpjAREZHIFCYiIhKZwkRERCJTmIiISGQKExERiUxhIiIikSlMREQkMoWJiIhEpjAREZHIFCYiIhKZwkRERCJTmIiISGQKExERiUxhIiIikSlMREQkMoWJiIhEpjAREZHIFCYiIhKZwkRERCJTmIiISGQKExERiUxhIiIikSlMREQkMoWJiIhEpjAREZHIFCYiIhKZwkRERCJTmIiISGQKExERiSypYWJmc81sq5kVmNk3Esz/qpltNrP1ZvaKmY2KmXe3mW0PH3cns04REYkmaWFiZinAT4DrgYnAnWY2Ma7ZGiDP3acAzwDfC5cdCHwLuAiYCXzLzAYkq1YREYkmmXsmM4ECdy9091rgCeDm2Abuvtjdq8OXK4AR4fMPAIvc/aC7HwIWAXOTWKuIiESQzDDJBopjXpeE05pyD/DimSxrZveaWb6Z5ZeXl0csV0REzlYyw8QSTPOEDc3uAvKA75/Jsu7+sLvnuXteVlbWWRcqIiLRJDNMSoCcmNcjgNL4RmZ2DfBN4CZ3rzmTZUVEpH1IZpisBMabWa6Z9QDmAQtjG5jZNOAhgiApi5n1EnCdmQ0ID7xfF04TEZF2KDVZb+zu9WY2nyAEUoBH3H2TmS0A8t19IUG3Vh/gaTMD2O3uN7n7QTP7DkEgASxw94PJqlVERKIx94SHMTqcvLw8z8/Pb+syREQ6FDNb5e55Ud9HV8CLiEhkChMREYlMYSIiIpEpTEREJDKFiYiIRKYwERGRyBQmIiISmcJEREQiU5iIiEhkChMREYlMYSIiIpEpTEREJDKFiYiIRKYwERGRyBQmIiISmcJEREQiU5iIiEhkChMREYlMYSIiIpEpTEREJDKFiYiIRKYwERGRyBQmIiISmcJEREQiU5iIiEhkChMREYlMYSIiIpEpTEREJDKFiYiIRKYwERGRyBQmIiISmcJEREQiU5iIiEhkChMREYksqWFiZnPNbKuZFZjZNxLMv8LMVptZvZndFjevwczWho+FyaxTRESiSU3WG5tZCvAT4FqgBFhpZgvdfXNMs93Ap4CvJ3iL4+5+QbLqExGRlpO0MAFmAgXuXghgZk8ANwPvhom7F4XzGpNYh4iIJFkyu7mygeKY1yXhtObqaWb5ZrbCzG5J1MDM7g3b5JeXl0epVUREIkhmmFiCaX4Gy4909zzgY8APzWzsX72Z+8PunufueVlZWWdbp4iIRJTMMCkBcmJejwBKm7uwu5eG/xYCS4BpLVmciIi0nGSGyUpgvJnlmlkPYB7QrLOyzGyAmaWFzzOBy4g51iIiIu1L0sLE3euB+cBLwBbgKXffZGYLzOwmADO70MxKgNuBh8xsU7j4eUC+ma0DFgP/HncWmIiItCPmfiaHMdqvvLw8z8/Pb+syREQ6FDNbFR6fjkRXwIuISGQKExERiUxhIiIikSlMREQkMoWJiIhE1qyxucxsBMF1IrOA4cBxYCPwPPCiu2tsLRGRLuy0YWJmvyAYU+s54LtAGdATOAeYC3zTzL7h7q8ls1AREWm/mrNn8p/uvjHB9I3A78Or20e2bFkiItKRnPaYSaIgCYc7mRLOr3X3gmQUJyIiHUOzD8Cb2RIz62dmA4F1wC/M7AfJK01ERDqKMzmbq7+7HwZuBX7h7jOAa5JTloiIdCRnEiapZjYM+CjBwXgRERHgzMJkAcEIwAXuvtLMxgDbk1OWiIh0JM2+B7y7Pw08HfO6EPhIMooSEZGO5ayvgDez+8zsDjNrdiCJiEjnFGU4FQMuB37fQrWIiEgHddZ7Fe7+k5YsREREOq7mDKfyyVPMdnd/rAXrERGRDqg5eyYXJphmwIcIxuxSmIiIdHGnDRN3v//kczMz4OPA3wErgH9JXmkiItJRNHcI+lTgU8DXgLeA29x9axLrEhGRDqQ5x0y+BHwFeAWY6+67kl6ViIh0KM3ZM/lvgnuYXA78KejpAoLjJu7uU5JUm4iIdBDNCZPcpFchIiIdWnPCZLe7+6kamJmdro2IiHRezbkCfrGZ3W9m77ubopn1MLM5ZvYocHdyyhMRkY6gOXsmc4HPAL8NRwo+RHAP+BTgL8B/ufva5JUoIiLtXXOuMzkBPAA8YGbdgUzguLtXJrs4ERHpGJpzanBP4AvAOGA98Ii71ye7MBER6Tiac8zkUSAP2ADcAPxnUisSEZEOpznHTCa6+2QAM/s58HZySxIRkY6mOXsmdSefqHtLREQSac6eyVQzOxw+N6BX+PrkFfD9kladiIh0CM05myulNQoREZGOK8pte0VERIAkh4mZzTWzrWZWYGbfSDD/CjNbbWb1ZnZb3Ly7zWx7+NAV9iIi7VjSwsTMUoCfANcDE4E7zWxiXLPdBPdJ+U3csgOBbwEXATOBb5nZgGTVKiIi0SRzz2QmUODuhe5eCzwB3BzbwN2L3H090Bi37AeARe5+0N0PAYsIhnUREZF2KJlhkg0Ux7wuCae12LJmdq+Z5ZtZfnl5+VkXKiIi0SQzTCzBtOYOU9+sZd39YXfPc/e8rKysMypORERaTjLDpATIiXk9AihthWVFRKSVJTNMVgLjzSzXzHoA84CFzVz2JeA6MxsQHni/LpwmIiLtUNLCJBx6ZT5BCGwBnnL3TWa2wMxuAjCzC82sBLgdeMjMNoXLHgS+QxBIK4EF4TQREWmHrLPcbTcvL8/z8/PbugwRkQ7FzFa5e17U99EV8CIiEpnCREREIlOYiIhIZAoTERGJTGEiIiKRKUxERCQyhYmIiESmMBERkcgUJiIiEpnCREREIktt6wJEOrqGRufw8Tqqjtdx+EQd1bUNnKhroKa+kfqG9w9X1CO1Gz27dyMtNYX0tBT69+pO/17d6ZOWilmiOy+IdAwKE5HTqK6tZ0fZMXYdPEbxwePsPljN3qrj7D9cQ/mRExw4VkvUIe66pxhZfdIY3K8nQ/qlkZ3Rm5EDe5EzsDejM9MZNbA3qSnqSJD2S2EiEnJ3ig8eZ2NpFRv3VLFl72G2lx2l5NDx97UbmN6D4Rk9Gd6/JxfkZJDVN42McA+jX6/upPdIIa17Cmmp3eie0o2TOxzuUFvfyIn6YM/lWE09VeEezcFjdZQdOUHZ4RoKy4/x2rYKjtc1vLvOHindGJOVzrjBfZg0vD/nZ/fj/OH9GZDeozU3kUiTFCbSZZ2oa2D17kOs3nWIVbsOsaa4ksrqOgBSuxnjBvdh+sgB3JGXw7jBfRidmU7OwN70SUv+fxt3p+JoLbsPVrOz4hjby45QsP8oa3ZX8tz6ve+2GzmwNzNGDWD6yAxmjBrIhKF96dZN3WXS+hQm0mXUNzSypriSZdsrWF54gLW7K6ltaATgnCF9mDtpKFNGZDA5uz/nDO1DWmpKm9VqZmT1TSOrbxozRg1437zK6lo2lx5mw54q1hZXsqyggmfX7AEgo3d3LsodyKVjM5k1PpPczHQdi5FWoTCRTq38SA2L3yljybYyXt9ewZET9XQzmDS8P5+6bDSXjBnE9JED6N+7e1uX2mwZvXtw6bhMLh2XCQR7MSWHjrOy6CDLdxxgeeEBXtq0H4Ccgb248pws5kwYzKVjM+nZve0CUjo33RxLOp2iimP8edM+Fm3ez+rdh3CHof16cuU5Wcw+N4tLx2Z2qPA4G7sPVLN0ezlLt5bz5o4Kqmsb6NU9hSvOyeS6iUO5ZuIQ+vfq3NtAmqelbo6lMJFOYfeBap7bUMrz6/eyqfQwAOdn9+Pa84Zy7cQhnDesb5ft7qmpb2BF4UEWbd7Hy5vL2Hf4BN1TjFnjs7hx8jCunTSEfj0VLF2VwiSOwqTrOXC0huc37OUPa/awenclABfkZPDBKcOYe/5QRgzo3cYVtj/uztriSl7YsJcXNuxjT+VxeqR249rzhnDzBcOZfe5geqTqFOSuRGESR2HSNdQ1NLJkazlP5Rez+J0y6hudc4f05ZZp2Xxo6jAFyBlwd9YUV/LHNXt4bv1eDhyrJaN3d265IJvbZozg/Oz+bV2itAKFSRyFSee2s+IYT6zcze9WlVBxtJbMPj348LRsbp0+gvOG9Wvr8jq8uoZGlm2v4JnVJSzatJ/ahkbOG9aPO2fmcMu0bHWDdWIKkzgKk86ntr6Rv2zex2/e2s2bOw6Q0s2YM2EwH83LYfa5WXTXFeFJUVldy8J1pTy5sphNpYfp2b0bH5oynI9fPIoLcjLaujxpYQqTOAqTzmNf1Ql+8/Zufvv2bsqP1JCd0Ys7Z+Zwe14OQ/r1bOvyupT1JZX89u3d/HFtKdW1DUwZ0Z9PXDyKD00drtOMOwmFSRyFSce3evchHlm2kxc37qPRndnnZPGJS0Zx5TmDSdFV3W3qyIk6nl2zh18t30VB2VEyenfnYzNH8olLRjGsf6+2Lk8iUJjEUZh0THUNjby4cR8/X7aTdcWV9O2ZyrwLc7jr4lGMGpTe1uVJHHdneeEBHn2ziEWb99PNjOsnD+Oey3PVBdZBtVSY6Ap4aRNHTtTxxNvF/OKNnZRWnWBMZjoLbp7ER6aPIL0Vxr6Ss2NmXDo2k0vHZlJ8sJpfLS/iibeL+dO6Ui4cPYDPzhrDNecN0Z5kF6Q9E2lV+w+f4JE3dvKbFbs5UlPPRbkD+dysMcyZMFgDFHZQR2vqeXJlMY8s28meyuPkZqbzuVljuHV6to6rdADq5oqjMGnfdpQf5aGlO3h2zR4aGp0bJg/j81eMZfIIXcvQWdQ3NPLnTft4aGkhG/ZUkdknjU9fNpq7Lh6loVvaMYVJHIVJ+7ShpIoHlhTw50376JHSjY/m5fC5WWMYOUgXF3ZW7s7yHQd48LVCXttWTt+0VO66ZBSfuSyXrL5pbV2exFGYxFGYtC9v7zzIf7+6nde3V9C3ZyqfvGQUn74sl8w++jLpSjbuqeJ/l+zghY176ZHSjTsuzOHzV44lO0NngLUXCpM4CpO25+4sK6jgv18t4O2dBxmU3oN7ZuVy18WjdAV1F1dYfpQHl+7g96v3YAa3ThvBF2ePZXSmzthrawqTOAqTtuPuLN5axo9eKWBdcSVD+/Xk81eOYd6FI+nVQwdg5T0lh6p5+LVCnlhZTH1DI7dckM2X5oxjbFafti6ty1KYxFGYtD53Z9Hm/fz41e1s3HOY7Ixe3HfVWG6bMaJN71Io7V/Z4RM8/Fohj7+1i5r6Rj44ZTj3zxnHOUP6tnVpXY7CJI7CpPU0Njp/2byfH7+ync17DzNqUG++NHscH56erfGy5IxUHK3hZ6/v5LHlRVTXNXDD5GF8ec54zh2qUGktHSJMzGwu8CMgBfiZu/973Pw04FfADOAAcIe7F5nZaGALsDVsusLdv3CqdSlMku9kiPzole1s2XuY3Mx05l81jpsvGE6qQkQiOHSslp8tK+TRN3dxtKaeGyYP5ctXj2fCUI0InWztPkzMLAXYBlwLlAArgTvdfXNMm/uAKe7+BTObB3zY3e8Iw+Q5dz+/uetTmCRPfIiMyUzn/qvH8aEpChFpWZXVtfx82U5+8UaRQqWVdIThVGYCBe5eCGBmTwA3A5tj2twMfDt8/gzwP9ZV763aDrk7L216f4j81x1TFSKSNBm9e/C1687lnstzeWTZTh55o4gXNuzjxsnD+PLV6v5qz5IZJtlAcczrEuCiptq4e72ZVQGDwnm5ZrYGOAz8o7u/Hr8CM7sXuBdg5MiRLVt9F3bywPoPXw6OieRmpvODj07lpqkKEWkdGb178NXrzuUzl+e+u6fywsa93DB5GF+5erwO1LdDyQyTRHsY8X1qTbXZC4x09wNmNgP4g5lNcvfD72vo/jDwMATdXC1Qc5fm7ry8pYwfvryNTaWHGT2ot0JE2tTJPZXPXHYyVHbywoa93BiGyniFSruRzDApAXJiXo8ASptoU2JmqUB/4KAHB3JqANx9lZntAM4BdFAkCdydV7aU8cNXtrFxT3B21n/cPpVbdGBd2okB6T34+geC7q+fvl7Io28W8fyGvXxoynC+fPU4xg1WqLS1ZIbJSmC8meUCe4B5wMfi2iwE7gaWA7cBr7q7m1kWQag0mNkYYDxQmMRau6T4EBk5sDffv20KH56WrRCRdmlAeg/+du4EPjtrzLuh8qf1pQqVdiDZpwbfAPyQ4NTgR9z9X8xsAZDv7gvNrCfwGDANOAjMc/dCM/sIsACoBxqAb7n7n061Lp3N1XzuzqvvlPHDl7ezYU8VOQN7cf+c8Xx4mq4TkY7l4LFaHn6tkF8tL+J4XQM3TR3O/XPGM26wrqhvrnZ/anBrU5ic3sk9kR+98l6IzL9qHLdOH6EQkQ7twNEaHn69kF+9uYsT9Q3aUzkDCpM4CpOmub93xfqm0qA7a/6ccdoTkU7nwNEafvr6znf3VE6eUqyzv5qmMImjMPlrjY3OS5v28eNXC9gSDnsy/6px3KIQkU7uwNEafr5sJ4++WcSx2gZumDyU+VeNZ+JwXfwYT2ESR2HynoZG57n1pfxkcQHb9h9lTGY68+eM0ym+0uUcOlb7bqgcqann2olD+PKc8brDZwyFSRyFCdQ1NPLsmj08uGQHhRXHGD+4D/dfPZ4bJw8jRfdXly6s6ngdv3yjiEfe2EnV8TquPCeL+XPGceHogW1dWptTmMTpymFyoq6Bp/OLeXBpIXsqjzNpeD/mXzWOD0waSjeFiMi7jpyo41fLd/HIsp0cOFbLzNyBzL9qHLPGZ9JVR3JSmMTpimFy+EQdj6/YxSPLiqg4WsP0kRncP2c8s8/N6rL/MUSa43htA799ezcPvbaD/YdrmJzdn/tmj+2Sf4ApTOJ0pTApP1LDL97YyWMrdnHkRD2zxmdy3+xxXDxmoEJE5AzU1Dfw7Oo9PLh0B0UHqhmblc7nrxjLzdOGd5kbvClM4nSFMCmqOMbDrxfyzKoS6hoamTtpKPfNHqeDiSIRNTQ6L2zYywNLdrBl72GG9EvjnstzuXPmSPr27N7W5SWVwiROZw6T1bsP8fDSQl7avI/u3brxkRkjuPeKMeRmprd1aSKdirvz2vYKHlyyg+WFB+iblsrHLhrJpy/LZWj/nm1dXlIoTOJ0tjBpaAyGgf/Z64Xk7zpE/17duevikdx9yWgG9+ucv9Qi7cm64kp++nohL2zYSzczbrpgOPdcnsuk4Z2rJ0BhEqezhMnRmnqeWlnML98sYvfBakYM6MU9l+fy0bwc0tOSOS6niCRSfLCany/byVP5xVTXNnDJmEHcc3kucyYM7hQH6xUmcTp6mBRVHOPR5UU8nV/C0Zp68kYN4J7Lc7lu0lBdIyLSDlRV1/HEyt08+mYRpVUnGD2oN5+8ZDS3543o0MdVFCZxOmKYNDY6r20v51fLd7F4axmp3YwbJw/j05flMjUno63LE5EE6hoa+fPGffzyzSJW7TpEeo8Ubpsxgk9cMqpDDiypMInTkcKksrqWZ1aV8PiKXRQdqCazTxofm5nDXReP0vEQkQ5kXXElv3yziOfX76W2oZFLxw7ik5eM4przhnSYoYsUJnHae5i4O6t3V/Lrt3bx/Pq91NQ3kjdqAJ+4ZBTXnz+MHqkd4xdPRP5axdEanlxZzK9X7KK06gRD+qVxx4UjmXdhDsMzerV1eaekMInTXsOksrqWZ9fs4cmVxbyz7wh90lK5ZdpwPjZzlEYwFelk6hsaWby1nN+8tYsl28oxYPa5g7njwhzmTBjcLkfrVpjEaU9h0tDovLmjgqfzS/jzpn3U1jcydUR/7rhwJDddMJw+OitLpNMrPljNEyt383R+CWVHasjqm8at07O5fUZOu7oTpMIkTnsIk50Vx/jdqhJ+v7qE0qoT9OuZyoenZXPHhSO1FyLSRdU3NLJkazlPrNzN4q3lNDQ600ZmcPuMHG6cPIz+vdv2TDCFSZy2CpNDx2p5bn0pv1u9h7XFlXQzmDU+i9vzRnDNeUPo2b1rjO8jIqdXduQEf1izh6fzS9hedpQeKd24+rzBfHhaNrPPHdwmx04VJnFaM0yO1tSzaPM+Fq4t5fXtFdQ3OhOG9uXW6dncfEE2Q3RGloicgruzYU8Vz67Zw5/WlVJxtJaM3t25/vyhfGjqcC7KHdRq15cpTOIkO0yO1dTz6jtlPL9+L4u3llFT30h2Ri8+OHUYN0/NVjeWiJyVuoZGlm2v4I9r9/CXzfuprm1gcN80rj9/KDdMHkbe6IFJDRaFSZxkhEnV8TpefWc/f964jyVby6mpbyQr/JBvmjqc6SMHdIrhFESkfThe28DLW/bz3PrSd79zBvdN4wOThvKBSUO5aMzAFj8jTGESp6XCpORQNa9sKePlLftZvuMA9Y3OkH7Bh3ljK/yVICICQXf6q++U8cL6vSzdVs7NpFZAAAALw0lEQVTxugb69+rO1RMGc83EIVxxTlaLnBmqMIlztmFS39DI2uJKFm8t49V3ytmy9zAAYzLTuXbSEOZOGsrUERnaAxGRNnO8toHXtpfz0sZ9vLq1jMrqOnqkdOOiMQO56tzBXDVh8FnfkkJhEudMwmRv1XFe31bB0u3lLNteQdXxOlK6GTNGDeCa8wZz9XlDGJvVfs4DFxE5qb6hkVW7DvHKO0EPSmH5MQBGD+rNledkMWt8FpeMHdTsUcYVJnFOFSZV1XWs2HmANwsqeGPHAQrKjgIwuG8as8ZnMWfCYC4fn0n/Xh135E8R6Zp2H6hmybYyFr9TxorCgxyva6B7ijFt5AAuHTuIy8ZlMnVERpOnHStM4sSGycFjtby98yBv7zzIWzsPsGXvYRodenVP4cLcgcwal8kV52RxzpA+ume6iHQaNfUNrCo6xNLt5bxZcICNpVV4+N03Y9QALsodyMzcgUzNyXj3GjiFSZyxE6f4bQseJ3/XQXaEu31pqd2YPnIAF40ZyKVjM7kgp+l0FhHpbE72yizfcYAVhQd4Z98RAHqkdGPyiP7kjR7AP9wwUWESK23YeJ/whZ8wY9QAZowawMzcgUwZ0Z+0VF2BLiICwcCzb+88SP6uQ+QXHWTDnioK/vVGhUmsyVOn+7o1q3TWlYhIM52oa6BXj9QWCZNO0+eT1r2bgkRE5Ay05NiBnSZMRESk7ShMREQkMoWJiIhEpjAREZHIkhomZjbXzLaaWYGZfSPB/DQzezKc/5aZjY6Z9/fh9K1m9oFk1ikiItEkLUzMLAX4CXA9MBG408wmxjW7Bzjk7uOA/wK+Gy47EZgHTALmAg+E7yciIu1QMvdMZgIF7l7o7rXAE8DNcW1uBh4Nnz8DXG3B+CY3A0+4e4277wQKwvcTEZF2KPpg+E3LBopjXpcAFzXVxt3rzawKGBROXxG3bHb8CszsXuDe8GWNmW1smdKTKhOoaOsimkF1tizV2bI6Qp0doUaAc1viTZIZJomuIIy/3L6pNs1ZFnd/GHgYwMzyW+IqzmRTnS1LdbYs1dlyOkKNENTZEu+TzG6uEiAn5vUIoLSpNmaWCvQHDjZzWRERaSeSGSYrgfFmlmtmPQgOqC+Ma7MQuDt8fhvwqgeDhS0E5oVne+UC44G3k1iriIhEkLRurvAYyHzgJSAFeMTdN5nZAiDf3RcCPwceM7MCgj2SeeGym8zsKWAzUA98yd0bTrPKh5P1s7Qw1dmyVGfLUp0tpyPUCC1UZ6cZNVhERNqOroAXEZHIFCYiIhJZhwoTM7vdzDaZWaOZ5cXNO+3wK+HJAG+Z2fZwGJcerVDzk2a2NnwUmdnaJtoVmdmGsF2LnKp3hnV+28z2xNR6QxPtTjlETivU+X0ze8fM1pvZs2aW0US7NtmeUYYQaqX6csxssZltCf8vfSVBm9lmVhXzu/BPrVljTB2n/Awt8ONwW643s+ltUOO5MdtprZkdNrO/iWvTJtvTzB4xs7LY6+/MbKCZLQq/AxeZ2YAmlr07bLPdzO5O1OavuHuHeQDnEVxgswTIi5k+EVgHpAG5wA4gJcHyTwHzwucPAl9s5fr/E/inJuYVAZltuG2/DXz9NG1Swm07BugRbvOJrVzndUBq+Py7wHfby/ZszvYB7gMeDJ/PA55s5RqHAdPD532BbQlqnA0815p1nc1nCNwAvEhwXdrFwFttXG8KsA8Y1R62J3AFMB3YGDPte8A3wuffSPT/BxgIFIb/DgifDzjd+jrUnom7b3H3rQlmnXb4lXCYljkEw7ZAMIzLLcmsN8H6Pwr8trXWmQTNGSInqdz9L+5eH75cQXANUnsRZQihVuHue919dfj8CLCFBKNLdBA3A7/ywAogw8yGtWE9VwM73H1XG9bwLnd/jeAs2Vixv39NfQd+AFjk7gfd/RCwiGCMxFPqUGFyComGbon/DzIIqIz5Iko4REsSzQL2u/v2JuY78BczWxUOE9MW5ofdBY80sfvbnO3cmj5D8JdpIm2xPZuzfd43hBBwcgihVhd2sU0D3kow+xIzW2dmL5rZpFYt7D2n+wzb2+/jPJr+Y7E9bE+AIe6+F4I/LIDBCdqc1XZN5nAqZ8XMXgaGJpj1TXf/Y1OLJZjW3KFbImtmzXdy6r2Sy9y91MwGA4vM7J3wL4sWc6o6gf8FvkOwTb5D0CX3mfi3SLBsi59b3pztaWbfJLgG6ddNvE3St2cCbfp7eCbMrA/wO+Bv3P1w3OzVBF01R8NjZ38guHC4tZ3uM2wX2xIgPP56E/D3CWa3l+3ZXGe1XdtdmLj7NWexWHOGX6kg2A1ODf8ibLEhWk5XswVDxdwKzDjFe5SG/5aZ2bMEXSYt+uXX3G1rZj8Fnkswq1WGuWnG9rwb+CBwtYedvAneI+nbM4EzGUKoxN4/hFCrMbPuBEHya3f/ffz82HBx9xfM7AEzy3T3Vh20sBmfYXsadul6YLW774+f0V62Z2i/mQ1z971hl2BZgjYlBMd5ThpBcJz6lDpLN9dph18Jv3QWEwzbAsEwLk3t6bS0a4B33L0k0UwzSzezviefExxkbtURkOP6mj/cxPqbM0ROUpnZXODvgJvcvbqJNm21PaMMIdQqwuMzPwe2uPsPmmgz9ORxHDObSfA9caC1agzX25zPcCHwyfCsrouBqpNdOG2gyZ6H9rA9Y8T+/jX1HfgScJ2ZDQi7u68Lp51aa59hEOVB8CVXAtQA+4GXYuZ9k+BMmq3A9THTXwCGh8/HEIRMAfA0kNZKdf8S+ELctOHACzF1rQsfmwi6c1p72z4GbADWh79ww+LrDF/fQHAG0I42qrOAoD93bfh4ML7OttyeibYPsIAg/AB6hr97BeHv4phW3n6XE3RZrI/ZhjcAXzj5OwrMD7fbOoKTHC5tg8854WcYV6cR3IBvR/i7m9fadYZ19CYIh/4x09p8exKE216gLvzevIfg+NwrwPbw34Fh2zzgZzHLfib8HS0APt2c9Wk4FRERiayzdHOJiEgbUpiIiEhkChMREYlMYSIiIpEpTEREJDKFiQhgZj80syvC5z8zs4lnsOynzOx/znB9RWaWeaZ1ng0ze7mp0WFFWorCRLo8MxsIXOzhUB3u/ll339zGZbWkxwhGKxZJGoWJdGrhldTPh4PsbTSzOxI0uw34c8wySyy8X46ZHTWz74aDDb5sZjPD+YVmdlPMe+SY2Z8tuJfJt2Le6w/hspuaGnCyqTbhuv8lrH2FmQ0Jpw+x4F4u68LHpeH0u8zsbQvumfGQmaWEb7WQ4AptkaRRmEhnNxcodfep7n4+MaER4zJgVRPLpwNL3H0GcAT4f8C1BKMxLIhpNxP4OHABcLu9d/O2z4TL5gFfNrNEIwQ31SYdWOHuUwnGpPpcOP3HwNJw+nRgk5mdB9xBMDjiBUBDWA8eDCOe1sS6RVqEwkQ6uw3ANeHexSx3r0rQZhhQ3sTytbwXQBsIvsTrwuejY9otcvcD7n4c+D3BsCUQhMPJYTRySDxabFNtanlvwM1VMeubQzDKM+7eEP5MVxMMJLrSgrt5Xk0wJMlJZQRDzogkRbsbNVikJbn7NjObQTD+1L+Z2V/cfUFcs+MEY2YlUufvjTnUSDAuHO7eGI76++6q4ldtZrMJBvm8xN2rzWxJ/HpO0yZ23Q2c+v+rAY+6e6Ih0Anf8/gplheJRHsm0qmZ2XCg2t0fB/6DoFso3hZgXMRVXWvB/bV7Edy97g2C4eUPhSExgeDWsvGa0ybeK8AXAcwsxcz6hdNus+DeHyfv9T0qfG4E94cpivQTipyCwkQ6u8nA22HXzzcJjnnEe57337/hbCwjOGtqLfA7d88n6B5LNbP1BDccW5Fguea0ifcV4Coz20DQ/TUpPPvsHwnuTLie4FarJ28rMIPg2Et9wncTaQEaNVgEMLNlwAfdvbKta2lpZvYjYKG7v9LWtUjnpT0TkcDXgJFtXUSSbFSQSLJpz0RERCLTnomIiESmMBERkcgUJiIiEpnCREREIlOYiIhIZP8fkNWJ8qYGT3IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(ins,data100)\n", "plt.ylabel('P(N,s)')\n", "plt.xlabel('s (imbalance)')\n", "plt.title('Gaussian with N=100')\n", "plt.axis([-10,10,0,0.25])\n", "plt.savefig('Gauss100.png')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.5" } }, "nbformat": 4, "nbformat_minor": 2 }