Skip to content

Instantly share code, notes, and snippets.

@dragly
Created September 18, 2013 12:39
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save dragly/6608593 to your computer and use it in GitHub Desktop.
Save dragly/6608593 to your computer and use it in GitHub Desktop.
QMC Force notebook
{
"metadata": {
"name": ""
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"forces = loadtxt(\"/home/svenni/scratch/vmc-test-fast/forces.dat\", unpack=True)\n",
"plot(forces[0], forces[1])"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 79,
"text": [
"[<matplotlib.lines.Line2D at 0xbf65390>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD5CAYAAAAk7Y4VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHbNJREFUeJzt3X1UVOedB/DvKBjFFwQU1BkQFZA3RVyNXRMb1BAbakhi\naEJeKqvGtWY9mmxOd2u2qdomMcb0mLc2MTlqYk+ipklaaSS0NXbU1SJxxWjARHwBB1AiIhqDCo53\n//gpb4IMzAzP3Hu/n3M4OMOdmd949DsPv/vc57FomqaBiIgMr5vqAoiIqGsw8ImITIKBT0RkEgx8\nIiKTYOATEZkEA5+IyCT8VBdwncViUV0CEZHudGRmvU+N8DVNM+TXkiVLlNfA98f3x/dnvK+O8qnA\nJyIi72HgExGZBAO/C6SkpKguwav4/vSN7888LFpnGkFeYLFYOtWTIiIyq47mJkf4REQmwcAnIjIJ\nBj4RkUnoOvC3bQOcTtVVEBHpg65P2vbpA5w8CfTt66WiiIh8mKlO2t5yC3D5suoqiIj0QdeB37Mn\nA5+IyFW6DvxbbgEuXVJdBRGRPug+8DnCJyJyja4Dny0dIiLX6Trw2dIhInKdrgOfI3wiItfpOvA5\nwicicp3uA58jfCIi1+g68NnSISJynVuBP3v2bISFhWHUqFFtHrNw4UJER0cjKSkJBQUF7rzcDdjS\nISJynVuBP2vWLOTm5rb585ycHBw5cgTFxcV4++23MX/+fHde7gZs6RARuc6twJ80aRKCgoLa/Hl2\ndjaysrIAABMmTEBNTQ0qKyvdeclm2NIhInKdnzefvLy8HOHh4Q23bTYbysrKEBYW1urxS5cubfhz\nSkpKu3tRsqVDRGZit9tht9s7/XivBj6AG5butFgsbR7bNPBdwRE+EZlJy4HwsmXLOvR4r87SsVqt\ncDgcDbfLyspgtVo99vwc4RMRuc6rgZ+eno7169cDAPLy8tC/f/822zmdwZO2RESuc6ul8/DDD2P7\n9u2oqqpCeHg4li1bhvr6egDAvHnzkJaWhpycHERFRaF3795Yt26dR4q+ji0dIiLX6XqLw7feAgoK\ngNWrvVQUEZEPM9UWhxzhExG5TteBzx4+EZHrdB/4nKVDROQaXQc+WzpERK7TdeBzhE9E5DpdBz5H\n+ERErtN14POkLRGR63Qf+GzpEBG5RteBz5YOEZHrdB34HOETEblO14HPET4Rket0Hfg8aUtE5Drd\nBz5bOkRErtF94NfXA76x3icRkW/TdeBbLIC/P1BXp7oSIiLfp+vAB+TELds6RETt033g88QtEZFr\nDBH4HOETEbVP94HPufhERK7RfeCzpUNE5BrdBz5P2hIRuUb3gc8RPhGRaxj4REQmofvAZ0uHiMg1\nug98jvCJiFxjiMDnCJ+IqH26D3zOwycico3uA58tHSIi1+g+8Hv1AmprVVdBROT7dB/4ISHAmTOq\nqyAi8n26D/zQUODbb1VXQUTk+xj4REQmofvAHzgQOH1adRVERL5P94HPET4RkWsY+EREJqH7wA8M\nBC5elLn4+/dziiYRUVt0H/gWS2MfPzMTePZZ1RUREfkm3Qc+IG2d8nKgpAT4wx+AAwdUV0RE5HsM\nEfgDBwL//CcQHg488wywcqXcv38/4HCorY2IyFcYIvBDQ4GdO4HYWOCHP2wc4S9dCmzYoLQ0IiKf\nYZjA37FDAj8uDiguBurrZYRfXq66OiIi32CIwB84EKiqksDv1Quw2YC9e4HSUgY+EdF1hgj80FD5\nHhcn3xMSgPffl7XyGfhERMJQgT9ypHxPTAQ2bQKmTAHKytTVRUTkSwwR+AMHyldIiNxOTJQWz913\nA5WVgNOptj4iIl9giMBPSgJ+//vG24mJ8v3WW4GgIC69QEQEGCTwe/UCMjIab8fEACNGSPBbrezj\nExEBBgn8lvz9gSNHgIAABj4R0XWGDPymrFaeuCUiAkwQ+DYbR/hERIAJAp8tHSIiwcAnIjIJwwd+\nXBzw5ZfAlSuqKyEiUsvtwM/NzUVsbCyio6OxYsWKG35ut9sRGBiI5ORkJCcn47nnnnP3JTskPBwY\nNgzYvr1LX5aIyOf4ufNgp9OJBQsWYOvWrbBarRg/fjzS09MRd31Rm2vuuOMOZGdnu1WoOzIygI8+\nAqZOVVYCEZFybo3w8/PzERUVhcjISPj7+yMzMxObN2++4ThN09x5Gbc98ADwySfA9OnA448rLYWI\nSBm3Rvjl5eUIDw9vuG2z2bBnz55mx1gsFuzevRtJSUmwWq14+eWXER8f3+rzLV26tOHPKSkpSElJ\ncae8BiNGAOPHA6NHA6tXA0uWSKuHiEhP7HY77HZ7px/vVuBbLJZ2jxk7diwcDgcCAgLw2Wef4b77\n7sPhw4dbPbZp4Hvap5/K90uXgNdea9wGkYhIL1oOhJctW9ahx7vV0rFarXA02TTW4XDAZrM1O6Zv\n374ICAgAANx9992or69HdXW1Oy/rlkWLgLVrZYMUIiIzcSvwx40bh+LiYpSUlKCurg6bNm1Cenp6\ns2MqKysbevj5+fnQNA3BwcHuvKxbhg6VwE9LA/73f5WVQUTU5dxq6fj5+eGNN97AtGnT4HQ6MWfO\nHMTFxWH16tUAgHnz5uGjjz7Cm2++CT8/PwQEBGDjxo0eKdwd994ru2FlZAC7dkmPn4jI6Cya6ik0\n11gsli6fzfPqq0B2NvD55136skREHtHR3DR14FdWyv63VVVd+rJERB7R0dw0/NIKNxMaCtTXAwrP\nIRMRdRlTB77FAkRHA8XFqishIvI+Uwc+INshtnFZABGRoTDwGfhEZBKmD3y2dIjILEwf+BzhE5FZ\nmHpaJgDU1MhCaufPy0lcIiK94LTMDurfH+jVCzh1SnUlRETeZfrAB2QbxKIi1VUQEXkXAx9AUpLs\ne0tEZGQMfEjg79+vugoiIu9i4IMjfCIyB9PP0gGAixeB4GDg3DmgRw8lJRARdRhn6XRCr17AsGHA\noUOqKyEi8h4G/jVs6xCR0THwrxk9GjhwQHUVRETew8C/ZuhQoKxMdRVERN7DwL9myBCgokJ1FURE\n3sPAv2bwYODkSdVVEBF5DwP/musjfN+YpEpE5HkM/Gv69pXVMr/7TnUlRETewcBvgn18IjIyBn4T\n7OMTkZEx8JvgCJ+IjIyB3wRH+ERkZAz8JjjCJyIjY+A3wRE+ERkZA78JjvCJyMgY+E0MGcIRPhEZ\nFwO/iSFDgPJy4OpV4MIFYOdO1RUREXkOA7+Jvn2lj19UBHz8MTBzpuqKiIg8h4Hfwm23Abt2Adu3\nAyUlgMOhuiIiIs9g4LfQNPBjY9nWISLjYOC3cNttwJYtsqH5v/87sGOH6oqIiDyDgd9CXJyctP3h\nD+WLI3wiMgoGfgvdugG33w5Mniwbm5eVAWfOqK6KiMh9fqoL8EXvvgv06QP4+QHDhwOlpUBIiOqq\niIjcw8BvRdNwDwkBqqrU1UJE5Cls6bQjJIQtHSIyBgZ+OwYMYOATkTEw8NvBET4RGQUDvx0MfCIy\nCgZ+Oxj4RGQUDPx2MPCJyCgY+O1g4BORUTDw28HAJyKjYOC3gxdeEZFRMPDbERgIXLwI1NWproSI\nyD0M/HZYLEBQEFBdrboSIiL3MPBdwD4+ERkBA98FDHwiMgIGvgu4ng4RGYHbgZ+bm4vY2FhER0dj\nxYoVrR6zcOFCREdHIykpCQUFBe6+ZJfjCJ+IjMCtwHc6nViwYAFyc3NRVFSEDRs24NChQ82OycnJ\nwZEjR1BcXIy3334b8+fPd6tgFRj4RGQEbgV+fn4+oqKiEBkZCX9/f2RmZmLz5s3NjsnOzkZWVhYA\nYMKECaipqUFlZaU7L9vlBg4EdFYyEdEN3Nrxqry8HOHh4Q23bTYb9uzZ0+4xZWVlCAsLu+H5li5d\n2vDnlJQUpKSkuFOex8TEAHa76iqIyOzsdjvsboSRW4FvsVhcOk7TNJce1zTwfUlCAvDVV6qrICKz\nazkQXrZsWYce71ZLx2q1wuFwNNx2OByw2Ww3PaasrAxWq9Wdl+1yw4bJ8grnz6uuhIio89wK/HHj\nxqG4uBglJSWoq6vDpk2bkJ6e3uyY9PR0rF+/HgCQl5eH/v37t9rO8WXduwNxcUBRkepKiIg6z62W\njp+fH9544w1MmzYNTqcTc+bMQVxcHFavXg0AmDdvHtLS0pCTk4OoqCj07t0b69at80jhXS0hASgs\nBH7wA9WVEBF1jkVr2WBXxGKx3NDr9yUrVwIVFcCqVaorISISHc1NXmnrosREoKAAeOkloMWlBkRE\nusARvoscDiAiQk7gpqQAa9eqroiIzK6jucnA74CiIrnqNjYWKC0F+vVTXRERmRlbOl4UHw+EhQFT\npwIbNqiuhoioYxj4nTB/PrB8OZdbICJ9YeB3wtSpwL/9GzB9Orc+JCL9YOB30pIlEvb/93+qKyEi\ncg0Dv5MsFmD8eODLL1VXQkTkGga+G5KSgP37VVdBROQaBr4bxozhCJ+I9IPz8N1w7hxgtcr37t1V\nV0NEZsN5+F0oMBAIDQWOHFFdCRFR+xj4bkpKYluHiPSBge+mW28FPvwQ0Fk3iohMiIHvpkWLgJIS\nWUWTiMiX8aStBzgcsiNWdTXQo4fqaojILHjSVoHwcCA6Gti3T3UlRERtY+B7yMSJwK5dqqsgImob\nA99DbrsN2L1bdRVERG1j4HvIbbfJCF+npyGIyAQY+B4SEQH4+QHHj6uuhIiodQx8D7FYgAkTgC++\nUF0JEVHrGPgeNHIkcPiw6iqIiFrHwPegmBgGPhH5Lga+B8XEAMXFqqsgImodA9+DoqOBb77hTB0i\n8k0MfA8aMEC+nzmjtg4iotYw8D3IYmEfn4h8l5/qAozmeuBPnKi6EiL9q6sDqqqAb78FTp8GamuB\nbt0av2prZbXaCxfk2Lo64PJl2YXuu++Aq1cbv5xOoL5ejtE04JZb5JizZ2Xhw5oaec2f/hRYt07p\n2/YaBr6HRUc3jvBfeAEIDgZ+9jO1NRF1lVOnGq9Fqa8Hjh2Trx49pOV59qwEeFWV/NnplPB1OiWo\nL12S4K2ulgsZNU0eN3CgfPXuLfddD/GePYHISKBfPyAgAAgKktcKDAT69JGtR5t+QPTo0bii7eXL\nQN++8pigIKB/fznG6VT21+d1XB7ZwzZtAlavBj74ABg+HEhNBTZvVl0VUeuuXgUqK2VE7OcHHDgg\nodu/P/D11/KzoiJg5065sLBfP+Crr4CKChnMDBggI+a+fYGyMhkxjx8P+PtLeI4YAQwb1jhSDwmR\nxwwYICHbvbu0Qrt3l/Du2VNeIzhYgvf681DrOpqbDHwPq60FpkyRX0Hj42X7Q4dDdVVkFk4nsGOH\ntDlOn5ZWx9ChEqIXLgClpfKz0lLgxAkJ6cBACfi6OiAhQUbGZ88CsbHAkCES2LffDuTlyah49Gi5\nv7q6cQ+I8+eBsDBpaVosqv8WzIOB7wPOnJE2ziuvyH+g4mL5dZSos2prgb/8RUa8kZFAYaHcl5cn\nv0FeugQMGiSBPGgQkJgIhIZKm6O0FPj+e6BXL3lsZKR8CEREADab3E/6xMD3MZMnA4sXA3fdpboS\n8lVNT0yuXi0j85AQ4OJF6YmXlclvibffLqPnEyck0Pv1kxH1ww/LKP3kSRnhx8aqfkfUVTqamzxp\n62Vjx8pOWAx8qqoCfv972SHtxAlg+3YJ9OPHJbB79wZmzpRBQnW19LMHDZLjIyKkT34z0dFd8z5I\nvxj4XjZ2LE/amsHBgxLYAQHAtm0yIs/PB/bvl9H6oEHA3r3A3XfL1dihocB//Zf0wqOj2VahrsGW\njpcdOgT86EcyiuNsA+MoKAA++0xG4g6HzGIBpK8+darM0Bo9Grj1VjkBeuqUjNLHjlVbNxkLWzo+\nJjZWpp/9/e/AtGmqq6GWLl+WC3BupqQEeOst4MoV+dD+4gv5AL//fmDwYBmlr1kjI/yrV2WKIZEv\n4gi/C6xbB/zxj0BOjupKqKn8fODOO+VE6b33An/4g8wjHzpUphrm5Mh0w7//HZg9W9owTqeM1DMy\nZMYMkUqcpeODLl2SEPnNb4Cf/ERG/KTO+fPAJ58Av/gF8PTTwEsvSY996FCZ9XL0qMxZv+cemU47\nYQIQFaW6aqIbMfB91LZtwMqVcll4bq7qasxD02RmTHKyXBi0Zo2M6CdNAhYskBF+drZMiZwzhxcN\nkb4w8H3YuXOA1SrT83r2VF2NsV25Ir3011+XwK+vl3ZNZibwn/8pFx8R6R1P2vqwwEC58jYvT2Zv\nWCycjucphw5JuJeWSntm82YJ+Z495e87MlJG+5wpRWbGf/5dbOpU4PPP5STh8uWqq9Gfd98FnnpK\n+vCATHmcOxdISZGTqnPmyHTIL74AjhyRxcCGDZMPV4Y9mR1H+F1syhTgwQdlvrYf//Zd8u23EvQn\nTgBbtki4x8VJ0H/wgfThv/lGFgAjorZxzNPFJk6UNVJee02WXDD4aQu37d8vs2QKC2VpgZ07ZZrr\np5/K0r3PPis9eoY9Uft40laBykppP4SGyvLJQ4aorkitqiqZ0x4YKLe3bgVefFHaMAcOAL/9LfDY\nY2prJPJFHc1NjvAVCAuTMLu+sBogHwL5+WrrUkHT5IrV5GT58HM6gYULgR//GJg3T3ZLYtgTeQa7\nyAqNHStrskyfLnPDP/9cVlA0g++/lw+5Y8dkOeDf/EZ68ffeKyP9J5/knHgiT+MIX6Hk5MYRfm6u\nhP/Vq7KByoULamvztieflHWGHnsM+NWvgEcflVZOXp7sBcywJ/I89vAVKimR+fh798r8/MBAGeUv\nWSJzyD/++MbHHD0qJy9DQ7u8XI/Zu1eWLdi1S9apefxxLjhG1Bm80lZnHn1UphSGhclFWOnpMvrt\n2RPYsAG4447mx2dmyrouzz2npl5X7NghF0EFBgJpaTIV9aWXZDnhigr57eXNN2WzDyLqPAa+zpSU\nACNHAi+/LEH45z/Llne/+AWwapW0OJqKi5NZPZ9/7r2aysrk5OnQoa4df/w48MADwBNPyIfXH/8o\nm3s4ncBHH8lI/qGH5IPMZpOVKENCvFc/kVl0WeBXV1fjoYceQmlpKSIjI/Hhhx+ifyuToSMjI9Gv\nXz90794d/v7+yG9jKopZAx+QZXj/9V/l6tBp04BnngF+/WtZd2fXLmDECDnu4kVZabNHD7nC1Btt\nkOJiubApPl7aLS2tWiWzaZ55BvjnP6WON96QfVX/8hcZ1X/ySfNAP3euccolEXlOl03LfPHFF5Ga\nmorDhw9j6tSpePHFF9ssyG63o6CgoM2wN7u0NAnyf/mXxtvdu8t0xaZ9/MJCOdE5eDDw1Vedf729\ne2XJ5uucTvl+9apMh3zmGXn+wsLmjzt2DHj+eaBPH2DcOFmvpqQE+J//kVk2+/YB//jHjaN3hj2R\nb+j0CD82Nhbbt29HWFgYTp06hZSUFHz99dc3HDds2DDs3bsXIe38Dm/mEX5Tq1fLejB+ftK2Wby4\ncX7+mjUybbNbN+AHPwB+9rOOP//ly8CAAcDttwN/+pOcK0hPlw+ZxER5zoMH5TeM8nLg7bflcfX1\ncqJ18mTgv//bc++XiDqvy1o6QUFBOHv2LABA0zQEBwc33G5q+PDhCAwMRPfu3TFv3jzMnTu3zcKX\nLFnScDslJQUpKSmdKc0wrlyR0fyePbJH6qJFsttSnz7A7t3Ae+/JUgMHDgD/8R+uPedf/yqzgCIi\nZLbPL38pQR8aCqSmykJjixfLHPnYWOnJBwRIj97fX9o1PXp4930TUevsdjvsdnvD7WXLlnku8FNT\nU3Hq1Kkb7n/++eeRlZXVLOCDg4NRXV19w7EnT57E4MGDcfr0aaSmpuL111/HpEmTbiyEI/xWPfus\nTMX84APprf/yl7L8b3q63D93rpwYLSuTD4Lrxo+X3xZabpq9YAEQHg7Mny+BnpwsJ42/+EI+RI4e\nbVwrfv58IDhYlj44f162AOSCb0S+o8O5qXXSyJEjtZMnT2qapmkVFRXayJEj233M0qVLtZdffrnV\nn7lRiqFduKBpERGatmSJpvXtq2mVlZp29aqmWa2a9s03mjZihKbFx2vaW281PubMGU0DNO1Xv2r+\nXFevalp4uKYVFsrtdevkuKIiTcvN1bR7721+/NGj8prh4ZpWU+PNd0lEndHR3Oz0Sdv09HS89957\nAID33nsP99133w3H1NbW4rvvvgMAfP/99/jb3/6GUaNGdfYlTal3b5mzfvAg8P770nqxWIAf/Qh4\n5x2ZAfPKKzJ7ZscOOQG7Z4+cBN6ypfF5srJkZO/vL1M7AZkHv3Wr3J42TaaENjV8uOz5um4dT7wS\nGYFb0zIffPBBnDhxotm0zIqKCsydOxdbtmzBsWPHMGPGDADAlStX8Oijj2Lx4sWtF8KWTod8/LFc\ntDV9usx7//Wv5QNh1iyZvnnpknwgHDokM27mzZPpn7fcwu39iIyCF16ZxLlzMttm1SrpywMyP376\ndLkS9+mngfXr5Zht22R3rfvvV1szEXkWl0c2icBAWUZ4+vTG+5KS5OpYu12mbc6dK7NtFi4EWum4\nEZHJcIRvMBs2SHvn0CHVlRCRt7GlY3KaBlRXc60aIjNgS8cHNb1Qwtsslq4P+658fyrw/emb0d9f\nRzDwu4DR/8Hx/ekb3595MPCJiEyCgU9EZBI+ddKWiIg6piMR7jNLYfnI5w4RkWGxpUNEZBIMfCIi\nk1Ae+Lm5uYiNjUV0dDRWrFihuhyPcjgcmDx5MhISEpCYmIjXXntNdUle4XQ6kZycjHvuuUd1KR5V\nU1ODjIwMxMXFIT4+Hnktd5TXueXLlyMhIQGjRo3CI488gsuXL6suyS2zZ89GWFhYsxV5q6urkZqa\nipiYGNx1112oqalRWKF7Wnt/P//5zxEXF4ekpCTMmDED586du+lzKA18p9OJBQsWIDc3F0VFRdiw\nYQMOGWhNAH9/f6xatQqFhYXIy8vD7373O0O9v+teffVVxMfHG+7E+6JFi5CWloZDhw7hwIEDiLu+\nrrQBlJSU4J133sG+fftw8OBBOJ1ObNy4UXVZbpk1axZyc3Ob3efq3tt60Nr7u+uuu1BYWIgvv/wS\nMTExWL58+U2fQ2ng5+fnIyoqCpGRkfD390dmZiY2b96ssiSPGjRoEMaMGQMA6NOnD+Li4lBRUaG4\nKs8qKytDTk4OHn/8cUOdeD937hx27tyJ2bNnAwD8/PwQaKBNAfr16wd/f3/U1tbiypUrqK2thdVq\nVV2WWyZNmoSgoKBm92VnZyMrKwsAkJWVhT+33PRBR1p7f6mpqejWTWJ8woQJKCsru+lzKA388vJy\nhIeHN9y22WwoLy9XWJH3lJSUoKCgABMmTFBdikc99dRTWLlyZcM/OqM4fvw4Bg4ciFmzZmHs2LGY\nO3cuamtrVZflMcHBwXj66acRERGBIUOGoH///rjzzjtVl+VxlZWVCAsLAwCEhYWhsrJScUXes3bt\nWqSlpd30GKX/S43WAmjLhQsXkJGRgVdffRV9mm48q3OffvopQkNDkZycbKjRPSAb9uzbtw9PPPEE\n9u3bh969e+u6HdDS0aNH8corr6CkpAQVFRW4cOEC3n//fdVleZXFYjFs5jz//PPo0aMHHnnkkZse\npzTwrVYrHA5Hw22HwwGbzaawIs+rr6/HAw88gMcee6zVbSD1bPfu3cjOzsawYcPw8MMPY9u2bZg5\nc6bqsjzCZrPBZrNh/PjxAICMjAzs27dPcVWes3fvXkycOBEhISHw8/PDjBkzsHv3btVleVxYWBhO\nnToFADh58iRCQ0MVV+R57777LnJyclz6wFYa+OPGjUNxcTFKSkpQV1eHTZs2IT09XWVJHqVpGubM\nmYP4+Hg8+eSTqsvxuBdeeAEOhwPHjx/Hxo0bMWXKFKxfv151WR4xaNAghIeH4/DhwwCArVu3IiEh\nQXFVnhMbG4u8vDxcvHgRmqZh69atiI+PV12Wx7my97ae5ebmYuXKldi8eTN69uzZ/gPc3DTdbTk5\nOVpMTIw2YsQI7YUXXlBdjkft3LlTs1gsWlJSkjZmzBhtzJgx2meffaa6LK+w2+3aPffco7oMj9q/\nf782btw4bfTo0dr999+v1dTUqC7Jo1asWKHFx8driYmJ2syZM7W6ujrVJbklMzNTGzx4sObv76/Z\nbDZt7dq12pkzZ7SpU6dq0dHRWmpqqnb27FnVZXZay/e3Zs0aLSoqSouIiGjIl/nz59/0OXxmLR0i\nIvIuY02tICKiNjHwiYhMgoFPRGQSDHwiIpNg4BMRmQQDn4jIJP4feESw5msmHM8AAAAASUVORK5C\nYII=\n",
"text": [
"<matplotlib.figure.Figure at 0xbaa3b50>"
]
}
],
"prompt_number": 79
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"plot(forces[0], -cumsum(forces[1]))"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 80,
"text": [
"[<matplotlib.lines.Line2D at 0xc1a1190>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD9CAYAAABZVQdHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl0VGWax/FvCXHQUVlaU7EJTlQMIQFNEERtl5KQICAx\nLAZohMjm3i1uSNvTNrQKoXHaxn1DLRRZ9EjIICBELBYBgwYVBSc2TTRAEpckLLLEJHf+eJs0m0hS\nVblVt36fczyGIlX3qSP+ePLUe9/XZVmWhYiIONJJdhcgIiLBo5AXEXEwhbyIiIMp5EVEHEwhLyLi\nYAp5EREH8zvkq6qqGDRoEB07diQxMZEPP/yQiooK0tLSiI+PJz09naqqqkDUKiIiDeR3yN911130\n6dOHzZs389lnn5GQkEBOTg5paWkUFRWRmppKTk5OIGoVEZEGcvlzM9TOnTtJSUnhn//852GPJyQk\nsGLFCtxuN2VlZXg8Hr788ku/ixURkYbxq5PfunUrZ511FiNHjqRLly6MHTuWH3/8kfLyctxuNwBu\nt5vy8vKAFCsiIg3T3J8n19TUUFhYyFNPPUW3bt0YN27cUaMZl8uFy+U66rnHekxERH5ZQwYwfnXy\nsbGxxMbG0q1bNwAGDRpEYWEhMTExlJWVAVBaWkp0dPTPFurUf/785z/bXoPen95fJL4/J783y2r4\ndN2vkI+JiaFdu3YUFRUBkJ+fT1JSEv369cPr9QLg9XrJzMz05zIiItJIfo1rAJ588kmGDRtGdXU1\n559/Pq+88gq1tbVkZWUxY8YM4uLimDdvXiBqFRGRBvI75C+66CLWr19/1OP5+fn+vnRY83g8dpcQ\nVHp/4c3J78/J760x/FpC6deFXa5GzZdERCJZQ7NT2xqIiDiYQl5ExMEU8iIiDqaQFxFxMIW8iIiD\nKeRFRBxMIS8i4mAKeRGJWNXV8Mgj8NxzdlcSPAp5EYlIH3wAKSmwbh307m13NcHj97YGIiLhZOdO\nmDAB8vJg+nQYOBCcvPO5OnkRiRhvvw1JSWBZ8MUXMGiQswMe1MmLSATYvh3uvBO+/BJmz4Yrr7S7\noqajTl5EHKuuDp59FpKT4cIL4ZNPIivgQZ28iDjUpk1w881mNLNiBSQm2l2RPdTJi4ijHDgAEyfC\n1VfDsGGwalXkBjyokxcRB1m9GsaOhYQEM5pp29buiuynkBeRsLdzJzzwACxcCE88AQMG2F1R6NC4\nRkTC2sFlkQCff66AP5I6eREJS5G8LLIh1MmLSFipq4NnnjHLIi+6KDKXRTaEOnkRCRubNpkPViGy\nl0U2hDp5EQl5Bw7An/9slkXeeKOWRTaE3518XFwcZ5xxBs2aNSMqKoqCggIqKioYPHgwX3/9NXFx\nccybN49WrVoFol4RiTCrVpmbmjp21LLIxvC7k3e5XPh8PjZs2EBBQQEAOTk5pKWlUVRURGpqKjk5\nOX4XKiKRpaoKbrkFhg6FRx81q2gU8A0XkHGNZVmH/TovL4/s7GwAsrOzyc3NDcRlRCRC5OZCp05w\n0klmt0gti2w8v8c1LpeLnj170qxZM2655RbGjh1LeXk5brcbALfbTXl5+TGfO3HixPqvPR4PHo/H\n33JEJIx9+y387newYYOWRR7k8/nw+XyNfr7LOrINb6DS0lLOPvtsvvvuO9LS0njyySfJyMigsrKy\n/nvatGlDRUXF4Rd2uY76CUBEIpNlwaxZcO+9cNNNZu+ZU06xu6rQ1NDs9LuTP/vsswE466yz6N+/\nPwUFBbjdbsrKyoiJiaG0tJTo6Gh/LyMiDlVSArfeCtu2wTvvQNeudlfkLH7N5Pfu3cvu3bsB+PHH\nH1m6dCmdO3cmIyMDr9cLgNfrJTMz0/9KRcRR6urMAdpdusCll8L69Qr4YPCrky8vL6d///4A1NTU\nMGzYMNLT0+natStZWVnMmDGjfgmliMhB//gHjBkD+/eDz/fvvWck8PyeyTf6wprJi0Sc2lp4/HHI\nyYE//hF+/3to1szuqsJLk8/kRUROxOefw6hRcPrpUFAA551nd0WRQdsaiEhQVVeb1TLXXGP2ncnP\nV8A3JXXyIhI069eb7j0uzqx9j421u6LIo5AXkYDbuxceeghef93M4IcMAZfL7qoik8Y1IhJQK1aY\nfd63b4eNG83eMwp4+6iTF5GA2LXLnLP6v/9rDvXIyLC7IgF18iISAEuXQufOUFNjVtEo4EOHOnkR\nabRdu+C+++Ddd+HFFyE93e6K5Ejq5EWkUZYtM927ZcFnnyngQ5U6eRFpkF274P77YfFi07336mV3\nRXI86uRF5ITl58OFF5rtCTZuVMCHA3XyIvKLdu82s/fFi+GFF+Daa+2uSE6UOnkROa78/H+vnNm4\nUQEfbtTJi8gx7d4N48fDwoWme+/d2+6KpDHUyYvIUd57z3TvBw6Y7l0BH77UyYtIvUO79+efhz59\n7K5I/KVOXkQAWL7crJw52L0r4J1BnbxIhNu71+w5M3++mb0r3J1FnbxIBFu3DpKTobJS3btTqZMX\niUDV1TBpEsyYAU89BYMG2V2RBItCXiTCbNwIw4fDOefAJ59ATIzdFUkwaVwjEiFqa2HqVOjRA37/\ne1iwQAEfCdTJi0SALVsgOxuaNzfnrsbF2V2RNJWAdPK1tbWkpKTQr18/ACoqKkhLSyM+Pp709HSq\nqqoCcRkRaSDLgueeg+7dzdx9+XIFfKQJSMhPnz6dxMREXP86yDEnJ4e0tDSKiopITU0lJycnEJcR\nkQbYvt3cqfrSS7ByJYwbBydpQBtx/P5Pvm3bNhYtWsSYMWOwLAuAvLw8srOzAcjOziY3N9ffy4hI\nA8yeDSkpcNllsHYtJCbaXZHYxe+Z/N133820adPYtWtX/WPl5eW43W4A3G435eXlx3zuxIkT67/2\neDx4PB5/yxGJaFVVcPvtUFgIixZB1652VyT+8vl8+Hy+Rj/fr5BfuHAh0dHRpKSk/GwRLperfoxz\npENDXkT8s2IFjBgB111nQv7UU+2uSALhyAZ40qRJDXq+XyG/Zs0a8vLyWLRoEfv372fXrl0MHz4c\nt9tNWVkZMTExlJaWEh0d7c9lROQ4qqvhoYdg5kwzf9ddq3Iov2bykydPpqSkhK1btzJnzhx69OjB\na6+9RkZGBl6vFwCv10tmZmZAihWRw23eDJdeCps2mRubFPBypIB+1n5wLDNhwgSWLVtGfHw8y5cv\nZ8KECYG8jEjEsyx4+mm46iq49VZzY5N+YJZjcVkHl8Q09YVdLmy6tEhYKyuDUaPgu+9g1iyIj7e7\nImlKDc1OrZoVCSN5eWbXyC5dYM0aBbz8Mm1rIBIGfvwR7rkHli6Ft96CK66wuyIJF+rkRULc+vXm\nxqb9++HTTxXw0jDq5EVCVE0N5OTAE0+YPd+zsuyuSMKRQl4kBBUXw403QosW5sam2Fi7K5JwpXGN\nSIh580245BLIzDQzeAW8+EOdvEiI2LvX7BS5fLn2nZHAUScvEgI2boRu3cwqmsJCBbwEjkJexEaW\nBc8+a47kGz8eXn8dzjjD7qrESTSuEbFJRQWMHQv//CesXg0dOthdkTiROnkRG6xebda+t2sH69Yp\n4CV41MmLNKHaWpg82Wwu9tJLZu93kWBSyIs0kR07YNgwM4f/+GNo29buiiQSaFwj0gSWLYOLL4Zr\nroH33lPAS9NRJy8SRLW18PDD8MILZlvgHj3srkgijUJeJEi+/daMZ376yYxnzj7b7ookEmlcIxIE\nK1eaPd+7d4f8fAW82EedvEgA1dXBtGnw+OPwyivQu7fdFUmkU8iLBMgPP0B2trnJaf16swZexG4a\n14gEwIcfmtUzCQmwYoUCXkKHOnkRP1iWOdTj0UfNCprMTLsrEjmcQl6kkX78EUaNgn/8w2xNcN55\ndlckcjSNa0QaYetWuPxyOOUU+OADBbyELr9Cfv/+/XTv3p3k5GQSExP5wx/+AEBFRQVpaWnEx8eT\nnp5OVVVVQIoVCQXvvQeXXQZjxpgVNC1a2F2RyM9zWZZl+fMCe/fu5dRTT6WmpoYrrriCxx57jLy8\nPM4880zGjx/P1KlTqaysJCcn5/ALu1z4eWmRJmVZMH06TJ0Ks2eDx2N3RRKJGpqdfo9rTj31VACq\nq6upra2ldevW5OXlkZ2dDUB2dja5ubn+XkbEVvv2meWRXq+ZvyvgJVz4/cFrXV0dXbp0YcuWLdx2\n220kJSVRXl6O2+0GwO12U15efsznTpw4sf5rj8eDR//nSAgqKYH+/SE+3szf/9XXiDQJn8+Hz+dr\n9PP9HtcctHPnTnr16sWUKVMYMGAAlZWV9b/Xpk0bKioqDr+wxjUSBlatgsGD4Z574N57weWyuyKJ\ndE0+rjmoZcuW9O3bl48//hi3201ZWRkApaWlREdHB+oyIk3i4NmrgwaZD1fvu08BL+HJr5D//vvv\n61fO7Nu3j2XLlpGSkkJGRgZerxcAr9dLpu4QkTBy4ADcfLM5vemDD6BXL7srEmk8v2bypaWlZGdn\nU1dXR11dHcOHDyc1NZWUlBSysrKYMWMGcXFxzJs3L1D1igRVaSkMHGh2jVy7Fk4/3e6KRPwTsJl8\ngy+smbyEmHXrzHjm1lvhwQfhJN0qKCGoodmpbQ1EgJdfhgkTYMYM6NfP7mpEAkchLxHtp5/g7rvN\nGawrV5pdJEWcRCEvEevbb+GGG8zcvaAAWra0uyKRwNPUUSLSxx9Dt25w5ZWQl6eAF+dSJy8RZ9Ys\nGDfu3+vgRZxMIS8Ro6YGHngAcnNh+XLo3NnuikSCTyEvEWHXLrM9QW2tOX+1TRu7KxJpGprJi+N9\n/TX85jdw7rmwaJECXiKLQl4cbf16c4LT6NFmm4Lm+tlVIoz+yItjvf22uXv1pZcgI8PuakTsoZAX\nx7EseOwxeOIJWLIEunSxuyIR+yjkxVFqa83yyJUrzQZjsbF2VyRiL4W8OMaBAzB8OHz3nQl53eAk\nog9exSF27oRrrzWjmsWLFfAiBynkJeyVlsLVV0NSEsyZAy1a2F2RSOhQyEtYKyoya+CzsuDJJ6FZ\nM7srEgktmslL2CoogOuvh0cfhVGj7K5GJDQp5CUsLV4M2dnmsI/rrrO7GpHQpXGNhJ2ZM2HkSFiw\nQAEv8kvUyUvYsCyYNg2eeQbefx86drS7IpHQp5CXsFBXB/feC/n58MEH0Lat3RWJhAeFvIS8Awfg\npptgxw5YtQpatbK7IpHwoZm8hLRdu6BvX6iuhnffVcCLNJRfIV9SUsI111xDUlISnTp14oknngCg\noqKCtLQ04uPjSU9Pp6qqKiDFSmT57jvweOCCC2DePN3kJNIYLsuyrMY+uaysjLKyMpKTk9mzZw8X\nX3wxubm5vPLKK5x55pmMHz+eqVOnUllZSU5OzuEXdrnw49LicGVlkJoKmZnwyCPgctldkUhoaGh2\n+tXJx8TEkJycDMBpp51Gx44d2b59O3l5eWRnZwOQnZ1Nbm6uP5eRCLNtm9mmYOhQc6OTAl6k8QL2\nwWtxcTEbNmyge/fulJeX43a7AXC73ZSXlx/zORMnTqz/2uPx4PF4AlWOhKniYtPB3367WU0jEul8\nPh8+n6/Rz/drXHPQnj17uPrqq/nTn/5EZmYmrVu3prKysv7327RpQ0VFxeEX1rhGjvDVV9CzJ4wf\nD3fcYXc1IqGpScc1AD/99BMDBw5k+PDhZGZmAqZ7LysrA6C0tJTo6Gh/LyMOt2kTXHMN/OlPCniR\nQPIr5C3LYvTo0SQmJjJu3Lj6xzMyMvB6vQB4vd768Bc5lk8/NR18Tg6MGWN3NSLO4te4ZvXq1Vx1\n1VVceOGFuP716diUKVO45JJLyMrK4ptvviEuLo558+bR6ogFzhrXCMBHH5l18E89BTfcYHc1IqGv\nodkZkJl8YyjkZe1as1Xwiy+af4vIL2todmpbA7HFihWmc5850xzbJyLBoW0NpMktXWoCfs4cBbxI\nsKmTlya1bBnceCO8/TZccYXd1Yg4n2by0mRWrYIBA2D+fAW8SGM1+Tp5kRPx0UcwcCC88YYCXqQp\nKeQl6DZuNMf0zZgBaWl2VyMSWRTyElRFRdCrF0yfDv362V2NSORRyEvQFBebO1kfeQQGD7a7GpHI\npJCXoNixw+wmef/9MGqU3dWIRC6FvATcd9+ZDn7sWPjd7+yuRiSyKeQloCorIT3dLJWcMMHuakRE\n6+QlYHbvNgF/6aXwt7/pRCeRYNAGZWKL/fuhTx84/3x44QUFvEiwKOSlydXUmBudTjkFZs2CZs3s\nrkjEubQLpTSpujoYPRqqq+HNNxXwIqFGIS+NZllwzz2wZYvZWfLkk+2uSESOpJCXRnvkEXj/fbM3\n/Kmn2l2NiByLQl4a5emnweuF1avhiJMdRSSEKOSlwd54wxy6vXIlxMTYXY2IHI9CXhrknXfMHP69\n9+Dcc+2uRkR+iUJeTtiaNXDTTbBwISQl2V2NiJwIbWsgJ2TTJujfH157Dbp3t7saETlRCnn5Rdu2\nQe/e8NhjOnhbJNz4HfKjRo3C7XbTuXPn+scqKipIS0sjPj6e9PR0qqqq/L2M2KSy0gT7nXfC8OF2\nVyMiDeV3yI8cOZIlS5Yc9lhOTg5paWkUFRWRmppKTk6Ov5cRG+zbBxkZZtOx++6zuxoRaYyA7F1T\nXFxMv3792LhxIwAJCQmsWLECt9tNWVkZHo+HL7/88vALa++akFZTAzfcYPajef11OEmDPZGQEBJ7\n15SXl+N2uwFwu92Ul5cf8/smTpxY/7XH48Hj8QSjHGkgyzKHfezZA3PnKuBF7OTz+fD5fI1+flA6\n+datW1NZWVn/+23atKGiouLwC6uTD1lTp5obnlatgjPOsLsaETlUQ7MzKD3awTENQGlpKdHR0cG4\njATB7Nlmy4JFixTwIk4QlJDPyMjA6/UC4PV6yczMDMZlJMBWrIC77jJ3tbZta3c1IhIIfo9rhg4d\nyooVK/j+++9xu9385S9/4frrrycrK4tvvvmGuLg45s2bR6sjdrHSuCa0bN4MHo8Z06Sm2l2NiPwc\nnQwlDVZWBpddBpMmwYgRdlcjIscTEjN5CR979kDfvjBqlAJexInUyUewmhq4/nqzXfBLL+nwbZFw\noE5eTohlwR13QG0tPPecAl7EqbTVcISaOhUKCszBH1FRdlcjIsGikI9Ab7wBzz4La9fC6afbXY2I\nBJNm8hFm1SoYOBCWL4dOneyuRkQaSjN5+Vlbt0JWltlwTAEvEhkU8hFi1y7o1w8efNBsHSwikUHj\nmghQW2uWSsbGmlm8VtKIhC+Na+QoEybA3r3w5JMKeJFIo9U1Dvfqq5CbC+vWaamkSCTSuMbBVq+G\nAQPM7pIdO9pdjYgEgsY1AkBxsTm+b+ZMBbxIJFPIO9Du3WYlzYQJcO21dlcjInbSuMZhamuhf3+z\n6djzz+uDVhGnCYmDvMU+Dz5o1sS/9ZYCXkQU8o7i9ZpwLyiAk0+2uxoRCQUa1zjEmjWQmQk+HyQm\n2l2NiASLVtdEoK+/hkGDTCevgBeRQynkw9yePZCRAfffD717212NiIQajWvCWF2dudnpzDPhxRf1\nQatIJNDqmgjyxz9CZSXMm6eAF5FjU8iHqddeM+H+4YdaSSMiPy9oM/klS5aQkJDABRdcwNSpU4N1\nmYi0di3cey/k5ZlRjYjIzwnKTL62tpYOHTqQn59P27Zt6datG7Nnz6bjIZuoaCbfOF9/DZddZmbw\nffvaXY2INLWQWEJZUFBA+/btiYuLIyoqiiFDhrBgwYJgXCqi7NoF110H48cr4EXkxARlJr99+3ba\ntWtX/+vY2Fg+/PDDo75v4sSJ9V97PB48Hk8wynGE2loYOhR+8xu46y67qxGRpuLz+fD5fI1+flBC\n3nWCSz0ODXk5vvvugwMHdLqTSKQ5sgGeNGlSg54flJBv27YtJSUl9b8uKSkhNjY2IK+9axeccUZA\nXipsPPccLFqk051EpOGCMpPv2rUrX331FcXFxVRXVzN37lwyMjL8ft1vv4UOHcy/I8WcOfDww/DO\nO9C6td3ViEi4CUrIN2/enKeeeopevXqRmJjI4MGDD1tZ01jR0TBsGDzwQACKDAPLlpn5+5Il0L69\n3dWISDgKu20Ndu82x9m9/jo4+XPaL7+Eq6+GN9+Eq66yuxoRCRUhsYQymE4/3awRHz4cvv/e7mqC\n44cfzPF9U6Yo4EXEP2HXyR90//1QVAS5uc5abXLggDmXtWtXmDbN7mpEJNQ0NDvDNuSrq6FbN7O0\ncPjwABZmo7o6sxa+thbmzoVmzeyuSERCTcSEPMCGDdCrF3zyCfz61wEqzCaWBePGmffy7rvQooXd\nFYlIKHL8TP5QKSlw221w880mJMPZX/8Ky5fDggUKeBEJnLAOeTB7qm/bBq+8YncljffMM/D882ap\nZKtWdlcjIk4S1uOagz7/HK65BlasCL8zTr1e+O//hpUr4dxz7a5GREJdRI1rDurUCaZOhRtugH37\n7K7mxL31FkyYYG56UsCLSDA4opMHM5MfPBjatYP/+Z+AvWzQvPMOjBplPmRNTra7GhEJFxG1uuZI\n338PF10Er74KaWkBfemAWr4chgwxJztdeqnd1YhIOInIcc1BZ54Jb7xh9rfZvNnuao5t7VrzE8e8\neQp4EQk+R4U8mP1e/vpXc4JSqG17sGEDZGaaD1udvO+OiIQOx4U8wE03QVYW9O8P+/fbXY2xapXZ\nruCZZ6BPH7urEZFI4aiZ/KHq6szc2+WC2bPhJBv/Ops/H265BWbNCu3PCkQk9EX0TP5QJ50EM2dC\naak5+Nouzz0Hd9xhbnRSwItIUwvK8X+hokULs0vlFVeYpZVNeQB2XR089JDZaGzVKjj//Ka7tojI\nQY4OeYA2bWDxYhP0bdo0zY6Vu3fDiBHmmMIPPjAnWomI2MGx45pD/dd/mZuOxo83nX0wbdkCl10G\nZ51l1sMr4EXEThER8mD2tFm40OxYuWxZ4F/fsswmY5deambwzz8P//Efgb+OiEhDOHZ1zc9ZvRoG\nDIC//x1++9vAvOY338Do0bBzp7nbNtw2SROR8KHVNb/giivMGOXBB80hHQcONP61ampMx37xxdCj\nB6xZo4AXkdAScSEPZtfKwkLYvh3at4e//c18SHqiduww56926ABz5pi/NP7wB2ju+I+xRSTcNDrk\n33zzTZKSkmjWrBmFhYWH/d6UKVO44IILSEhIYOnSpX4XGQxt2sCbb5oblTZsgPh4s+fNRx8dfsrU\ngQNmH5y8PHj4YejZE5KS4P/+D157Dd5/Hzp3Pvr1fT5fk70XO+j9hTcnvz8nv7fGaHTId+7cmfnz\n53PVVVcd9vimTZuYO3cumzZtYsmSJdx+++3U1dX5XWiwdO1qwvrrr80OlllZprvv3h3OOw9atoTr\nrzdjmd27zYeq27fDSy/B5Zf//Os6/Q+a3l94c/L7c/J7a4xGDxgSEhKO+fiCBQsYOnQoUVFRxMXF\n0b59ewoKCrg0xLdcbNnSLLG8/3744gsT6L/6lTnMIyrK7upERBon4FPkHTt2HBbosbGxbN++PdCX\nCRqXy8zsRUSc4Lghn5aWRllZ2VGPT548mX79+p3wRVwuV4Med4pJkybZXUJQ6f2FNye/Pye/t4Y6\nbsgva8RdQ23btqWkpKT+19u2baNt27ZHfZ9Ny/NFRCJKQJZQHhrYGRkZzJkzh+rqarZu3cpXX33F\nJZdcEojLiIhIAzU65OfPn0+7du1Yt24dffv2pXfv3gAkJiaSlZVFYmIivXv35plnnnH8WEZEJGRZ\nNli8eLHVoUMHq3379lZOTo4dJQTNN998Y3k8HisxMdFKSkqypk+fbndJAVdTU2MlJydb1113nd2l\nBFxlZaU1cOBAKyEhwerYsaO1du1au0sKqMmTJ1uJiYlWp06drKFDh1r79++3uyS/jBw50oqOjrY6\ndepU/9gPP/xg9ezZ07rgggustLQ0q7Ky0sYK/XOs93ffffdZCQkJ1oUXXmj179/fqqqqOu5rNPkd\nr7W1tdx5550sWbKETZs2MXv2bDaH6qnbjRAVFcXjjz/OF198wbp163j66acd9f4Apk+fTmJioiN/\nQrvrrrvo06cPmzdv5rPPPqNjx452lxQwxcXFvPjiixQWFrJx40Zqa2uZM2eO3WX5ZeTIkSxZsuSw\nx3JyckhLS6OoqIjU1FRycnJsqs5/x3p/6enpfPHFF3z66afEx8czZcqU475Gk4d8QUEB7du3Jy4u\njqioKIYMGcKCBQuauoygiYmJITk5GYDTTjuNjh07smPHDpurCpxt27axaNEixowZ47gPz3fu3Mmq\nVasYNWoUAM2bN6dly5Y2VxU4Z5xxBlFRUezdu5eamhr27t17zEUR4eTKK6+kdevWhz2Wl5dHdnY2\nANnZ2eQGe3/xIDrW+0tLS+Okf51n2r17d7Zt23bc12jykN++fTvt2rWr/3W4raNviOLiYjZs2ED3\n7t3tLiVg7r77bqZNm1b/h8xJtm7dyllnncXIkSPp0qULY8eOZe/evXaXFTBt2rTh3nvv5ZxzzuHX\nv/41rVq1omfPnnaXFXDl5eW43W4A3G435eXlNlcUPC+//DJ9+vQ57vc0+f+pTvwR/1j27NnDoEGD\nmD59Oqeddprd5QTEwoULiY6OJiUlxXFdPEBNTQ2FhYXcfvvtFBYW8p//+Z9h/aP+kbZs2cLf//53\niouL2bFjB3v27GHWrFl2lxVULpfLsZnz6KOPcvLJJ/PbX9gzvclD/sh19CUlJcTGxjZ1GUH1008/\nMXDgQG688UYyMzPtLidg1qxZQ15eHueeey5Dhw5l+fLljBgxwu6yAiY2NpbY2Fi6desGwKBBg47a\nfC+cffTRR1x++eX86le/onnz5gwYMIA1a9bYXVbAud3u+ps4S0tLiXbg8WyvvvoqixYtOqG/pJs8\n5Lt27cpXX31FcXEx1dXVzJ07l4yMjKYuI2gsy2L06NEkJiYybtw4u8sJqMmTJ1NSUsLWrVuZM2cO\nPXr0YObMmXaXFTAxMTG0a9eOoqIiAPLz80lKSrK5qsBJSEhg3bp17Nu3D8uyyM/PJ9GBByBkZGTg\n9XoB8HpcY0VqAAAA2klEQVS9jmq0AJYsWcK0adNYsGABLVq0+OUnBHP5z89ZtGiRFR8fb51//vnW\n5MmT7SghaFatWmW5XC7roosuspKTk63k5GRr8eLFdpcVcD6fz+rXr5/dZQTcJ598YnXt2vWEl6eF\nm6lTp9YvoRwxYoRVXV1td0l+GTJkiHX22WdbUVFRVmxsrPXyyy9bP/zwg5WamuqIJZRHvr8ZM2ZY\n7du3t84555z6fLntttuO+xq2Hf8nIiLB57wlEiIiUk8hLyLiYAp5EREHU8iLiDiYQl5ExMEU8iIi\nDvb/hmJjQbNo01gAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0xbf69710>"
]
}
],
"prompt_number": 80
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"energies = loadtxt(\"/home/svenni/scratch/vmc-test-fast/energies.dat\", unpack=True)\n",
"plot(energies[0], energies[1])"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 81,
"text": [
"[<matplotlib.lines.Line2D at 0xc1ce8d0>]"
]
},
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD5CAYAAADLL+UrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl4VNX9/99DEhICISFAEiABlBAh7IKl2lbCEqxEKCpW\nXB5SqLu2arWWLv4EW2gstUq/2PZbFI3UivJtgRSRNhSjVESFILssIWFJSEBCyEZIMjm/Pz49zJ3J\nrHfunTt37uf1PHmY5c65Z8jkfd/zPp9zjk0IIcAwDMNYgi5Gd4BhGIYJHSz6DMMwFoJFn2EYxkKw\n6DMMw1gIFn2GYRgLwaLPMAxjIaKN7oDEZrMZ3QWGYRhTEkjlfVg5fSFExP4899xzhveB3xu/P35/\n3n8qKgQmTRKw243vt78/gRJWos8wDGMkTz8NTJkCdIlgZQybeIdhGMZISkqAnTuBN980uif6EsHX\ns/AiJyfH6C7oRiS/N4Dfn9nx9/0tXw784hdAt2769sdobEJNKKQDNptNVT7FMAyjlvp6YM8eYNQo\nYNAg4MQJICnJ6F4FRqDayU6fYRjLUlQE5OYCS5dSlm82wVcDZ/oMw1iWY8eAgQOBZcuAd94xujeh\ngUWfYRjLUlYG/PSnQHMzMHOm0b0JDRzvMAxjWY4dA4YOBR59NPIHcCUs+gzDhB3PPgv85z/6n6es\nDBgyRP/zhBNcvcMwTFhRVgZMmADExAB/+Qswfbo+56mvB/r1AxobATOvAhOodnKmzzBMWPHqq8CC\nBcDgwcC77+on+mVlwNVXm1vw1cDxDsNYDLsdqK42uhfuaW0FXn8duP9+ICEBaGvT71xlZUBmpn7t\nhyss+gxjMbZvB+691+heuGf/fqBPH2DYMIp39BZ9q+X5AIs+w1iOy5eBS5e0a2/lSmDXLm3aamoC\nevWi2zExQHu7Nu2649gxFn2GYSyA3a6tmL77LnDbbcD588G31dwMxMfT7ehofZ3+2bM0kGs1WPQZ\nxmJoLfrt7UDfvrRYWbAoRV/veKe9nS4sVoNFn2Eshh6i/81vajM4zKKvPyz6DGMx7HZtxdRuB2Jj\n6d9gcRV9PTN9u51Fn2EYC6CH0+/aVXvR1zvTZ6cfILW1tcjNzUVWVhamT5+Ouro6j8fa7XaMGzcO\nM62yohHDhDFmEX2Od/RBtegXFBQgNzcXR44cwdSpU1FQUODx2OXLlyM7Oxs2q019Y5gwRGvRl/FO\nR0fwbbHo649q0S8qKkJ+fj4AID8/H+vXr3d73OnTp7Fp0ybcd999vLYOw4QBHR3aiml7uzkz/fZ2\nICpKv/bDFdXXuZqaGqSmpgIAUlNTUVNT4/a4J598EsuWLUN9fb3PNhctWnTldk5OTsTv3ckwRmCW\neIczffeUlJSgpKRE9eu9vuXc3FxUu6nDWrJkidN9m83mNrrZuHEjUlJSMG7cOL86qRR9hmH0QY94\nx4yZvlmrd1wN8eLFiwN6vde3XFxc7PG51NRUVFdXIy0tDWfOnEFKSkqnY7Zv346ioiJs2rQJLS0t\nqK+vx7x58/Dmm28G1EmGYbRDD6fPmb55UJ3pz5o1C4WFhQCAwsJCzJ49u9MxS5cuxalTp1BeXo41\na9ZgypQpLPgMYzBmiXdCkemz6AfAwoULUVxcjKysLGzduhULFy4EAFRVVSEvL8/ta7h6h2GMxyyT\nszjT1wfVbzk5ORlbtmzp9Hj//v3x3nvvdXp80qRJmDRpktrTMQyjEXY7/QihzQYiZq7Tt2L1Ds/I\nZRiLIcVZC5EGHKLPmb45YNFnGIshxVmrvJzX3jEXLPoMYzGkOGslqFynby5Y9BnGYkhx1kpQ9a7e\n0WsiP4s+wzCWQGunr9XaO0LQdolS9G02GmjVK+Jh0WcYxhJoKfodHfSjhdNvayOhj4lxPKZnrs/V\nOwzDWAItRd9uJ+GMigpe9JXRjkTPXJ+dPsMwlkDLTF9WwHTpoo/o61W22dFBcVIXCyqgBd8yw1gb\nLZ2+jEiiooLP9EMp+vJiZcVFAlj0GcZiaFmnLyMSveIdvTJ9q0Y7AIs+w1gOrTN9PUVfr0zfqoO4\nAIs+w1gOPeIds2X67PQZhrEMWg7kKuMdM2b6VoRFn2EshpniHc70tYdFn2Eshl7VO2bL9Fn0GYax\nBFqLvt5On0VfW1j0GcZi6BHvmHEgl6t3GIaxBHLAVauBXL0nZ3Gmry0s+gxjMfSId+RyBsEIfygz\nfa7eYRjGMugR7wDB5/qc6YcGFn2GsRh6VO8Awef6ly/TuvxKWPS1h0WfYSyG3U7r32s5OQsIPtd3\nF7lwpq89LPoMYzHsdiAuLvziHXdCzGvvaI+pRf/PfwYOHTK6FwxjLuT2hlrHO3qIPsc72mNq0f/X\nv4D9+43uBcOYC61FX4pnsJm+O/fNa+9oj6lFv3dv4Px5o3vBMOaio0M70XeNd4LJ9D05fc70tcXU\not+nD/DVV0b3gmHMhXT6Wk7OAoKPd9y5b157R3tMLfrs9BkmcPSKdzjTNwemFn12+gwTOFqKvtKd\na5Hph1L0uXrHhLDTZ5jA0bN6J9hM31WIo6P1yfR5INeksOgzTOCEa7zjaXIWxzvaYmrR53iHMQOv\nvgq89ZbRvXCg5UCu3pOzWPS1x9Rvm50+Ywb27gUSE43uhQMtZ+RqufaOJ9Gvr1ffZiDnsgqmdvqJ\niUBTkz5OgGFc+eILYMeOwF/X0EArSLrj2DHg9Ong+hUoWtbpa7n2TigzfR7INSldugC9egG1tUb3\nhLECGzcCf/tb4K+rrwcuXXL/3CuvhD760at6h+Mdc2Bq0Qc44mFCR2sr/QRKQ4Nn0W9uVtdmMGi9\nyqaek7N4GQbtMb3o82AuEyra2tSLvqd4xyjR1zLTN2udPou+SWGnz4QKPZz+pUuhH5PSM97Reu0d\nPTN9Fn2Twk6fCRV6ib6yTSGC32DcF+Fapx/KVTZZ9E0MO30mVKiNd+rr/Y93NmwAHnlEXf/8JVzX\n0w/15Cyu3jEpffqw6DOhQY3TFyIwp19bC1y4oL6P/qDX5CzO9M2BatGvra1Fbm4usrKyMH36dNTV\n1bk9rq6uDnPmzMHw4cORnZ2NHWoKnb3QuzfHO0xoUCP6LS0khP5W77S16Z/xh3Odfqgyfa7eUUFB\nQQFyc3Nx5MgRTJ06FQUFBW6Pe/zxxzFjxgwcOnQIe/fuxfDhw1V31h0c7zChQo3oNzTQv57iHdeB\n3FCIfrjGO5zphwbVol9UVIT8/HwAQH5+PtavX9/pmIsXL2Lbtm1YsGABACA6OhqJGs9HT0wELl7U\ntEmGcYuaTL+hAeje3f94x2yiz5OzzIfqt11TU4PU1FQAQGpqKmpqajodU15ejr59+2L+/PnYs2cP\nxo8fj+XLlyM+Pt5tm4sWLbpyOycnBzk5OT77waLPhAo1Tr++HkhJAc6edf+8EfFOuO6Ry6ts+kdJ\nSQlKSkpUv97r287NzUV1dXWnx5csWeJ032azwWazdTquvb0dpaWlWLFiBa677jo88cQTKCgowPPP\nP+/2fErR9xcWfSZUqI13UlOBEydoUNf1z8RIp6/HjFyzZPpmrt5xNcSLFy8O6PVeRb+4uNjjc6mp\nqaiurkZaWhrOnDmDlJSUTsekp6cjPT0d1113HQBgzpw5HrN/tSQlsegzoUFtvNOrF7ngtjZa/kDS\n0UEDvUaIvlYzcjneMR+qM/1Zs2ahsLAQAFBYWIjZs2d3OiYtLQ0ZGRk4cuQIAGDLli0YMWKE2lO6\npWdP+gothKbNMkwn1Dr9hAQgPr7zYG5LC/1r9oFcM07O4uodFSxcuBDFxcXIysrC1q1bsXDhQgBA\nVVUV8vLyrhz3P//zP7jnnnswZswY7N27Fz/72c+C77WC6GhyLY2NmjbLMJ1Qm+knJADdunUezJX3\nzZ7pa7GefkcHRV9dXBSJnb72qH7bycnJ2LJlS6fH+/fvj/fee+/K/TFjxuDzzz9Xexq/kLl+QoKu\np2Esjtp4p2dP96Ivnb9RdfrhtHOWJxHmtXe0x/QzcgEezGW0Y9s2WgrBHVrHO5Hi9LWYnOVJhNnp\na09EiD4P5jJa8dFHwL//7f65YETfk9OPjQ2t6EtR7to1vCZnGSH6Zq3eCZaIEP3ERMDDKhAMExCN\njcDly+6fUxvveMv0ExNDO5Brt5PYaRWbaBnvuBNhHsjVnogRfXb6jBY0NXkW9mAGcj3FO4mJoXX6\nWou+VpOzPIkwZ/raw6LPMAoaG72LfkdHYMLmayA3Kck40Q+nyVmc6YcOFn2GUdDU5D3eAQJz+/7E\nO2Z2+npX77Doa09EiD4P5DJa4cvpuw68+kIp+q7xTnOzcaIfExNek7NCnenzQK7J4YFcRis8ZfpC\n0OPdu6sT/fh47QZyz53z//yudHRon+lrMTnLW6Zvt2u/hSQ7fZPD8Q6jFZ6qd6RIxMVpG+/07Emi\nJsWyrY0uMJ7EUwjgqqvUC7bdTuKsV7yjdaZvs1F5qdZun6t3TA6LPqMVnpx+WxtFDV27Bib6LS10\nofAU73Tr5ixqrv+660dTk2PdnkDRYyBXz0wfoP8fT+MsamGnb3I402e0wlOm39pK4hOo6MvXeYp3\n4uOd2/Ql+vI4tSIoRV+68mBjE70nZwE0jsKirx0RIfrs9Bmt8BTvqBF9Iait2FjP8Y50+v6Kvuxb\nsKJvszny8mDQamN0bwOrgQ6e+3s+Fn0TwwO5jFZoGe/I/DwqynO8Ex/vXKESKtEHtMn1tVp7x1vG\nrle8w9U7JoadPuNKa6tjU/JAXuNpqQU1Tv/yZcemKZ7inUCdvlbxDqBNrs/xjvmICNFPSKCv5VqX\ndTHmZfVq4MknA3tNUxP9q6Xox8bSbU8zco2KdwBtnL7ek7MAfeIdrt4xOVFRQI8egTs7JnI5fx6o\nqOj8eGWl513WGhspanEnqGriHXmhANzHO54GcqOi9HP6HR2OjUq0mKCl1do73uIWrt7RlogQfYDc\nfn290b1gwoX6euD06c6P5+UB+/a5f01TE5CcrI/TDyTeiY83j9PXe+0dgOMdrYko0ectExlJQwOJ\nvqurr631PP7T2Kit6MtlGwDP8Y67gVxvoq91ph8u8Y63uIWrd7QlYkSf4x1GSX09OXdXgW9ocGT3\nrkinf/ly54uFmnhHOZDrKd4x2ulrMZBr1slZXL1jctjpM0qkAVBGPEJ4F/3GRloWwV02HazTd7ee\nvqeBXH9EP9gZuYA2YqrV2juhjnd4IDcCSEhgp884qK8nMVKKfksL/bF7c/rdu7sX9mBLNtPSgOpq\n58zb00BuqOIdd5GTmva0WnuHJ2eFhogRfY53GCUNDcDQoc6iLz8fnr4RNjbS58ids2xt9R7vvPwy\n8NvfOj+mHMjt0YO+RZw543heGe/4m+lrGe9066b+G4NEq3jHiMlZLPomh+MdRklDAzB8uHvRV+P0\n29q8O/2qKuDkSefHlCWbADBkCHDsGN3u6KBvIz17OtqUa+HExZnH6fPkLPMRMaLPTp9RUl8fuOhL\np68m3mlqAi5ccH5M6fQBEv2yMrp97hwJfmwsfYOQs4FjYrxvHOKv029pAUpLOz+urNPXOt7RM9PX\nI97hgVyTw5m+NTh9GrjpJt/HNTQA2dnOoi/ncfjj9AONd5qaOq//pBzIBZxFv7ISGDCAbss2laLv\nqZTSm9Nvbweuv57afvFF4OGHOx+jh9PXcz19QPt4p6ODBvW7RIz6BUbEvO0ePTjesQJVVcCBA96P\nkVU6/jr9/fuB737XOdMPNN5pbnbv9F3jHX9FX43Tb2wEduwAHnoIWLasc7UQoK3ou4qnjHd++Uvg\nlVcCa8vXQK6Woi+/ndhs2rVpJiJG9NnpW4OGBt8X90uXSEgHD/ZP9A8fBoqKSLTVVu+4c/re4p2q\nKqB/f7otB3Kl6Hurn/fm9Bsbgd696SI2bpz+ou/qzqXonz1LlUqBEMrJWVbO84EIE312+pFPYyOJ\nt7v1cyorgZ/8hGKchASgVy8SYymQDQ204Y7r56S6mo7Zts2R6buLd5Si/9xzwGefOZ53l+m7G8jV\nwul7Wh9ITi774APgf/9Xf9F3FWqZ6be0uD+3N0IZ77DoRwg8kGsNGhooVnBXavj++7S6ZkMDDZLa\nbPSv/Fw0NFC9vKvTr66mY8vLvcc7ykx/61bgyy8dzzc3+3b6ffuS4Fy44Cz6gQ7kJiR4dvo9etA3\nnLS00Dh9ZSQjM/1LlzyPm3hrK1TVOyz6EQLHO9ZAunR3v+uSEhLw2lr6PAAk+nIA15Po19QA3/gG\n3fY33qmpcV7gr6mJRFb5Olenb7M5yjYrK53jHX9Fv7WV3pO7i54UfcD9sg9AaOKdlhZ1oh+qyVlt\nbSz6EQEP5FoDpWtXIgTw4YckmEeOkDACzqJfXw/06+fe6d96K932N96pqXHug2xT6fZdnT4AjBkD\n7NxJmb7aeKdnT8/xTvfudDsmhi4yru3I3bwA7eMdKfpqnH4oJ2fJmdBWJWJEn52+NZAXdtcL/PHj\nJBxjxgCHDgXm9KurgRtuoOcSE33HO/X1jh9JczOdSyn6riWbADB5MmXurpm+ciDXH6fvLd6RuFvv\np6PD4ajj4rSNd/TK9LWOd+TqplaFRZ/RhQ8/BLZs0b5dT07/ww+BSZOA9HQSfXdOv6HBvdOvqSHB\n//xz4JprfMc7p0517kNTE51bOZjrWrIJkOhv2eKotAG0c/ruRN/1vYYi3tEj09cy3mlqYtGPCDje\nCS/++U8aWNUapdO32x3i9/nn5Nal6Htz+srPiRDk9FNT6bU2m+94R4q+bFcOXvbv7zveGTiQKmz6\n9XPewSqQTN/TQK4y3gHcO30tRb+lhb4tSORArtpMP1TxDjv9CKF7d/qg8T654cHFi/pchJVO/913\nHbNOKyuBjAyKTI4dczj9xETv8U5dHYlft26Ox3zFOzU19Jhst6WFXpOc7Oz0XQdyJZMnO6IdQF+n\nr6foX7xI/7+SYJ1+qCZnNTc7XxytRsSMYUdFOSoWlB98xhguXvS8F20wyN2tGhpoH9zjx+nxM2cc\n0Y3d7tnpJyfTbSnI1dV0IVDiK94Rgi4wysle3bvTvABfTh8AZs1yvhi4E31PF0x/SjYleot+XZ2z\n6Csz/UA3Zwnl5CyrxzsRI/qAI9dn0Teeixf1WduksZHEvbER+OorcviAQ/SliHnK9Hv2pM9HU5PD\ntaemOp/DV7wD0LLNyrV8uneniV/+OP2ZM+lHeT4tBnKbmoA+fRz3jXT6gYo0xzuhI2LiHYAHc8MJ\nf+OdDz6gunp/kYOx0ulXVlKkd/YsOXYZm7hz+nKmrowCgcCcfkyMw7lnZjralXGBv07flUDjHS2d\nfjDr6bsTfT0yfY53tCWiRJ8Hc8OHixf9uwA/9xxQXOx/u9LpNzSQ0790yZHhd+3qXfQbGvwTfV8L\nrgEk+sp4Jz7evdP3R/QDHcj15vT9EX2t6vQvXqT3LHF1+oFsuh7KyVns9CMIdvrhg79O/+xZ50XR\nfCGdfmMjOX2AJjv160e34+PJcbvGO+3tJJTx8ST6sm9q4x2l01fGO65O312844qagVxPM3J9Ve8o\n6/T1iHdkTCUjNH8J5eQsq2f6qkW/trYWubm5yMrKwvTp01HnuvDIf/n1r3+NESNGYNSoUbj77rtx\nWestcBTw+jvhg7+if+5cYKLf2EilkdLp9+vnLPoAlUVKBypFX0YfNpuz0z9zxnu8U1kJ/PnPzuvp\nA7ScwqVLJFahjHe8lWwaXb3TpQudLy6O/j8CmaDFk7NCh2rRLygoQG5uLo4cOYKpU6eioKCg0zEV\nFRVYuXIlSktLsW/fPtjtdqxZsyaoDnuDV9oMD4RwCK032tooz/dH9P/0JxKsxkZHrf3588Do0Z1F\nv6gIuO46ui1FX0Y7gLPoV1QAgwY5n0sZJ+zaBTz7bGenn5bm+MbgLd7x1+n7O5DraxkGI6t3oqKo\nD1L0A3H6oZycxZm+SoqKipCfnw8AyM/Px/r16zsd07NnT8TExKC5uRnt7e1obm7GAGWBssaw0w8P\nGhtJAHyJvoxn5GQnT3R0AI8+Suved+1K4lpfTwI7ahSwe7ezWx840JFbK0VfRj7K6KG8HLjqKufz\nKeOECxcogjp40CH6XbrQbFrZdqAlm65oWacfyslZ7uKdpiZqV0vR5+odbVEt+jU1NUj9bxiampqK\nGjljRUFycjKeeuopDBw4EP3790dSUhKmTZumvrc+YKcfHtTXO0oHvTm0c+dINHw5/fp6Ev79++l3\nnJBAr4mPp2WE5eCuO3r2JHE6f94R+UhBamujgdyMDOfXKOMd6dy/+orEOCkJ+P73SeDkGJJ0jmqd\nvlYDuYHGO3FxNDagdj6FO9FvbFTv9EM1Ocvqmb7XOv3c3FxUu9kCZ8mSJU73bTYbbG72HisrK8PL\nL7+MiooKJCYm4o477sBbb72Fe+65x+35Fi1adOV2Tk4OcnJy/HgLDhIT6YPIGIsUg5YWx2Qqd5w9\nS07900+9Oz0ppPv3k6j16EGxTP/+juWJvYl+fT1w9CgNvgIOQTp5kl4XE+P8GmWcUFdHx5w5QwIe\nF0cZv7JtZbwjJ6XZbPo5fSlYrv9n7uKdqirn1ytFv0sXOndLi/OMZH9xrd6Rmb4apx/KyVlmd/ol\nJSUoKSlR/Xqvol/spZYuNTUV1dXVSEtLw5kzZ5CSktLpmJ07d+KGG25A7/+uLHXbbbdh+/btfom+\nGpKTHTsTMcYhRb+pybvonztHot23Lznu9HTn51euBGbMcIj+vn0Op9/SQhGLTAs9iX737nTsgQO0\nmJp8rLGRLhyDB3d+jWu8M3MmCb2ra3eNd2Qdf2Mj9dHfkk15Pn+dfteuDverFMpA4x3AEfGoFX1X\npy8EXRjdndsb/sQ78mIaLGbP9F0N8eLFiwN6vep4Z9asWSgsLAQAFBYWYvbs2Z2OGTZsGHbs2IFL\nly5BCIEtW7YgOztb7Sl9kpzsyImZ0FNTA/zf/znEwNcYy9mzJPjp6bScQlaWcya+ciVtSejq9OWA\nrD+ib7PR8Tt3Oot+U5P7PB/oHO9cfz2NGbjO9HaNdwDnsk1/SzZ796bPrb9OPza2c+TR1kaCrrzI\n+KrTB4LL9d2JvmxTy0w/Kor6HEjdvzfM7vSDRbXoL1y4EMXFxcjKysLWrVuxcOFCAEBVVRXy8vIA\nAGPGjMG8efMwYcIEjB49GgDwwAMPaNBt9/TuHdjsTsZ/zpzx/Ny+ffRvcTHwy186i763MZZz54CU\nFMrUV62iCEYZR5w/T8dcuEDtlZU54h2Axg1SU+m+J9EHyJHv2gUMG0b35UCuJ9FXxgkXLpCZ+PLL\nztm/a7wD0GCuvEj5G+/07u1ox5voC9HZ6UtktKN0wr7q9IHgRN9d9Q6gfaYPaBvxcKavkuTkZGxx\ns2B6//798d577125/8wzz+CZZ55Re5oA+8ROXw8OHwbGjnVU5Sg5cwYYP56eO3qUZsdKMfAl+mfP\nUrvp6cAf/uB4TH4ZrK0l0bfZgHHjaDvEhAQSgOhoEsvoaODECe/rLfXsSQO/MtMfPJiWfU5NBW6+\nufPxyninro7cu1LclO0q4x3A2en7O5DbpQt946mspPr/6Gj3oi+3+evSpbPou0Y7QGDxTqDIAWDl\n0srKmb5aZvqA43eiRSzDTj+CkF+TGW0pLKQ/ctdBQYCEqq2NxP7oUfqDOnyYRNJTNdXu3cAXX5Cg\ny3inrY1ilHPn6Bi7ncTz7FkS/9GjSfylm+3Rw1Eh5GnMQNKzJwm9FKhZs6gEs6TEv3inVy/37bqL\nd5Rlm/46fYDio9OnvTt9ZXuxsc6zcl0rdwB9RV9+m1N+s1A6fS0zfUDbCh6zZ/rBElGin5zM8Y7W\n2O3A6tUkzidPOh5/7TX645GrXB48SKLfrRvl53I1S3ei/9JLwK9+RYKekkLC+61vkes/e5aOkRGJ\njHdSU2nQV+b5CQmOnad80bOnI88HSNQfeoi+pXgayPVH9N3FO7JsUxnF+ENaGs1X8Cb6yvY8xTtK\n9BZ9ZeUOoF+mD2gb77DTjyCSksh5aTXgw9AqmH37AlOmUIwCkNg8+ii5dSn6Bw6Q6E+ZQk7eW7zz\nySe0s1ZlJbU9ezawcSNdAKToy4u3FP3kZBJoKWzBiD4APPggDRzLkk8lSkGV8Y47pNNXxjvS6be3\nU9zhLadWkprqW/Rdnb6R8Y7rIC4QfKbvT7yjBVbP9CNK9KOi6IPoYRkgJkCEAJYuBR55hGa5Sqe/\nezf9AR4/TpFPVhawbRt91Z84kQTIk+h/9RUJ+/DhVDKZkkK/t549O4t+XBzdl2578GCH009MpAuG\nP2RkANde6/xYWhrFUO7W/JdOX9bOe4oCpNN3rd65cMH/ck1lfy5d8t/px8V1Fn0j4h0lwWT6/gzk\naiH6HR3UjnIswmpElOgDPJirJcXF5Ma/9z1an0Y6/e3bSeDLy+n53FzamHzoUPoBPIv+p58CX/sa\n5epRUc7RiavoDx3qcPq9egH33gtMnUrPr15NkZA/vPgiMG+e/+9biv6FCyTinmrDExPp201dnXO8\nU1fnf7mmRK70qdbpBxLvuJZsqllT37VyB3COd9xtyu4NXwO5WsU7ck6CHhv8mIWIe+tctqkd/+//\nUfYeHe3s9D/5BMjJcTj9KVPojzYz01Eh41qn//bbwD330DeCr38d+M53KFpR/vEpRf/8eYpkvvqK\nfp+9egHf/jZ9kwAcVS56IAVVir4npk6lUtGjR53jnQsXAhvEBRxrB6nN9MMx3gnHgVyr5/lABIo+\nO31tOHKEnP1tt9H9gQPpvhDk9O+5x+H0hwwBrr6anLmr6Dc2Un7/xBOUWb/4IlXpjBoF7NnjfM6+\nfR3VO7W1JKhxcXQeT4OpeiCdfl2d9/PGxwPr1tGPXKlTOv1ABnEBdU5/+3bgzjvpfnm582brgGNt\nnY4Ox2Na1el7E31v8c4XX9BF0hWtMn1f78XqeT4QgaLPTj8wvv3tzuILAGvWAHfc4fhDlvHOyZPk\nFqdMIadfWUliM3YsMGIEiV5qKv0epOgvWkSza4uKgNtvB264gdp0FVTXeKd3b7oQeKug0QNlvOPr\nvF260ECLOPkQAAAWIklEQVS0/Maip9NXRkaxscDy5cDatfRtqKQEmDSpc9+k8EtcnX5cnDrRP3Cg\n80C6/D/wNpD7wgvAz3/e+XEtqnfsdhpf2rvX8zHs9CNQ9Nnp+48Qjkoa18fffhuYO9fxmKzJ/v3v\naT2cjAxaduHSJfrjX70amDOHji0tdQy6fvUV/RFOm0YXhDVrPEcmvXpRHNTaSr/D5GS6EADeYxat\nkZuUnDwZ+HmVTl+PeEfp9GNjgcmTgX/8gwbX5cVUiYx42ttJbMvLg3f6K1ZQTPfgg86P+5Pp798P\nbNjguLhLghnIXbkS+Ogj6tPp0zS+5Amr1+gDESj67PT9p7aWqk9cF+w7cID+OL7+dcdjNhu5/RUr\ngJ/8hFxZejrl8jab8+CYLIPs0QP4+GOaYeuPu+rShSZcyRw/OZmcfmKi/6WPWtC1Kwno3/8e+DcM\nWbIZ6EBuUhId72+8k5QE/OAH9C1jyRJgzBj3YiZF//hx+t0dOODYVwCg39ubb9Lg+qpV5JbdUVxM\ny2B8/jkttfGvfzkuyBJlpt+nD83dGDgQ+M9/6PHWVprEd/vtNOFPib8zct2xejXw+OPAO+/QZ237\nds/tsNMPYhmGcCU5mT5sjG/KyuiP8uOPnb9eFxfT8gSuFQ4DB1Kppax5v+oq75UfMt5x50A9ISMe\nZbwTymhH8u1v08Vt/PjAXqe2ZNNmo1hMir7dTm1ERTnWqVeOEyxbRo8fPQr88IfO38qUxMfTa8vK\n6PewcaPz7/X22+n/OCMD+OlP6RwPPeTcRkcHVXAJQd+Cli93P6lN6fQzMuji/e679P/4n//QONGg\nQTTH4+GHgR//2PHaYOKdI0fonK++St96vC3vxZk+O31Lc/w4Obz0dBpgk2zdSpm9Kz/+MfCb3zju\nX31158FDJbKEMBDR79vXIfrS6Rsl+nZ74PFOQgLFJU1NgTl9gCKe2Fi6AERF0QD5yy/Tc9/4Bjl1\neSHp2pWOueYauhhPnuy+zSFDSBQPH6YF56KinEtQhw+nmGbGDJqTsXp15za2b6ffxUsvkRmQg8eu\nKDN9gL6hLVhA33z++U9amG/UKPrMlZU5r8DqS/S7d3cs7KdE7sW8YgVVlN10E/3/e9qNjZ1+BIo+\nZ/r+c/w4CfekSbSswu7d9Mf30Uf0B+RKTo6zwxsxwlGt4w41oi+dvjLTN0L0s7PpYhjouW02Eruz\nZwNz+gDwxhtU2QQA3/0uufB//5sc89GjdNv1QmKzkSi7u0gDNCmttJRE33VWsivTp9N5jh93fnzt\nWurPnXfSRcjTvAXl5CxJVBSweDHtM7xvHzByJH2TGT2aPm8SX5n+woXA+vXA00877/R19CgN3ubl\nUeRks9HnTRnxnDxJRQQAZ/pABIo+O33PdHTQH9/hw3Rfiv5DD9EmJjfdBDz/PDlHN3vidOIHP6B8\n1xMpKSQUrksSeyMlhQbjZLwzYIDz/rehwmYDnnkGmDAh8NcmJZHTDNTpZ2c7XvPXvwI/+xmJ17Zt\nNBHtwQfdC9aAAZ6FWIr+l1/6Fv2YGPp9vfWW47Fz52iPhDvu8O89REV1nu162200RvHqqyT6AG1c\n//nnjmN8ZfpXX01FBx99BPzoRw7hP3KERB9w/B9Mmwa88go5/ro6+nYyfz71gZ0+ABEmaNWV8nIh\n0tM1aSri+Pe/hejXT4jevYUoKhIiJ0eI4mLH89u2CREVJcQPf2hcHz/6SIiMDOqH3S5EW5sQDQ3G\n9UcN48cLkZwsxKpVwbd1zTVCTJsmxC9/KcTly0LU1gb2+ooKIdLShOjbV4iqKt/H79lDx+7bJ8RT\nTwmRmCjEY4/5f76uXYXYtavz4xs2CAEIcfgw3S8sFGLuXMfzsbFCXLrku/3aWiG+9jUhbr1ViPPn\nhXjuOSGefdb5mPZ2anvsWCEGDRLiBz8Q4vrrhXj/fSFefFGIJ57w//2YgUC1M+JEv7VViJgY+gNh\nnJk3T4jf/U6If/xDiGuvJXEtK3M+5q236A/eSG69lUTTrEybRu+hoyP4tu67j8Tygw/Uvb6jg/4v\nExL8788bb5AI5+SQsAZCXJwQBw6478eqVSTIQtAxmZn0+G9/S2ZEPueLlhYhHnhAiBkzhLjrLiFW\nr+58TGurEGvWCLF7N53jpZfo8z9tGt2OJCwv+kLQ1f3YMc2aMz2PP05in5QkRE0NOeihQ4WIjiYn\nHW4cPSrEI48Y3Qv1HDggRF2dNm0VFpKJaWpS38a0aUJMmBDYaz74gMQ1ULp3F+L4cd/HtbfThSg3\nl/pWURHYeVpayLT06iXEp5/6Pv7UKSFsNiEmTw7Pz3wwBKqdEZfpA1RKWF5udC+M46GHHKWU7e00\nSPvXv1KFR0oKDbj98IeUteu1fk0wZGZSJmtWsrPd77SlhptuAp56Krgc+tprHVtF+ktOTuAD0YD7\nTN/TcffcQwPQ27c7lrHwl9hY4Be/oPJYucifN9LTqcLnnXfC8zMfSmz/vVIYjs1mg1ZdWbCAqiDu\nv1+T5kyFXHv+ww+BG2+kUsy77qJJOe3tjoHClhYaSPN3pUrGvJw4QeWRchBVTx59FPjd79RdMAKl\nrY0Mjeu8AqsRqHay048wZLmdnGX76ac0s7ZLF+dqkrg4FnyrMGhQaAQfoG9ooRB8gKqNrC74amDR\njzCOH6eSQSn6O3Y4L6fAMIy1YdGPMMrKqNb6s88owmHRZxhGScSKfkWF0b0whuPHaZnj7GyaPXn6\nNM2cZRiGASJU9NPSaE2OQHbuiRSOH6f1VvLzaeG03/yGqxUYhnEQkdU7AJWo/f3v5HitxFVXAVu2\nkPAzDBP5cPXOf1Hu6WoV2tpoz9qBA43uCcMw4UrEin7fvrQ6oZU4eZI2MImJMbonDMOEKyz6Jke5\n81VZGcc6DMN4J2KH+OS2e5HMxx9TeWZ9Pe07+tJLntdVZxiGASLY6VtB9FesoN2s1qwBbrmFlmBQ\nbkHHMAzjCou+idm9m5z9jBlUk79mDef5DMN4h+Mdk9LYSLszydUTrbi4HMMwgRPRTv/cOaN7oR97\n9tAcBHb2DMMEQkSLfiQ7/d27aZ10hmGYQIhY0ZcbpHd0GN0Tfdi9Gxg3zuheMAxjNiJW9GNigB49\naA2eSKS0lEWfYZjAiVjRByI34mlvB778Ehg1yuieMAxjNlj0Tcjx47TcQjD7pjIMY01Y9E3ApUvA\nH/8I2O10/+BBYPhwY/vEMIw5YdE3Ae+9RxtO5+dTtHPoEIs+wzDqiHjRN3Ot/s6d5O7//nfgxReB\no0eBtWvJ6VttnwCGYbRBteivXbsWI0aMQFRUFEpLSz0et3nzZgwbNgxDhw7FCy+8oPZ0qjDzSpt2\nO/DNb5LYv/8+cNddNOt2/Xp2+gzDqEe16I8aNQrr1q3DjTfe6PEYu92Oxx57DJs3b8bBgwfx9ttv\n49ChQ2pPGTBmdvrl5TRQu2gRVemkpQEzZwL//CdV7rDoMwyjBtVr7wyTi7544bPPPkNmZiYGDx4M\nAJg7dy42bNiA4SFSrDFjgCVLaIJWF5MFWYcOARMn0tLJSUn0WGoqxTonTgCJicb2j2EYc6LrgmuV\nlZXIyMi4cj89PR2ffvqpx+MXLVp05XZOTg5ycnKCOv+115Jb/ugjIMimQs6hQyTw3/ue8+Pf+Q5t\neM4wjDUpKSlBidw5SQVeRT83NxfV1dWdHl+6dClmzpzps3GbzRZQZ5SirwU2GzB/PvD66+YT/YMH\nKdN35ZFHgNmzQ98fhmHCA1dDvHjx4oBe71X0i4O0lAMGDMCpU6eu3D916hTS09ODajNQ7r0XWLwY\naGkB4uJCeuqgOHTI/XLJCQnANdeEvj8Mw0QGmiTdQgi3j0+YMAFHjx5FRUUFWltb8c4772DWrFla\nnNJvUlIoC6+oCOlpg0IIrtBhGEYfVIv+unXrkJGRgR07diAvLw8333wzAKCqqgp5eXkAgOjoaKxY\nsQI33XQTsrOzceedd4ZsEFfJoEE0+GkGVq8G3nyTxiKSk43uDcMwkYZNeLLpIcZms3n8xhAs998P\nTJgAPPigLs0HxYsvAg88QLFNQwMwcCDdzswEtm41uncMw4Q7gWpnxG6XqGTw4PCMdzo6gJ//HOje\nHXjoIeCvfwUmTwYKC4G6OqN7xzBMJGKy6nV1hKvoV1YCra3An/5EOf6f/gQ8/DA5fUWlK8MwjGZY\nRvTDMdMvKwOuv55inalTaeOXqVON7hXDMJGMJUR/0KDwdPrHj1N2/+yzwOjRwIcfmm/mMMMw5sIS\nmX6/fsD58+FXq19WBgwZ0nnWLcMwjF5YwldGRVFGfvKk0T1xpqwMuPpqo3vBMIyVsIToA+E1mHvi\nBO2GJZ0+wzBMqLBEvAOEl+gvWABcdx1l+iz6DMOEEsuIfmYm7TwVDuzbR7tidXTQRi8MwzChwjLx\nTnY2rVxpNGfP0j63eXmU5we4ECnDMExQWMbpDx8eHqJ/4AAwciTw/PPA9u1G94ZhGKthGdG/6iqg\nuhpoaqJlD0KJEMBPfkIrfsbGkuhnZtIPwzBMKLGM6EdHA1lZtL/s+PGhPfdPfwps3Ejr6dx8c+jP\nzzAMI7FMpg9Qrh/CfdkB0GDtK68AH3xAcwXWrCGnzzAMYwSWEn0jcv2KCqBXL9rIZf58oLkZGDEi\ntH1gGIaRWEr0jajg2bcPGDWKbs+dC+TnA717h7YPDMMwEkuJ/ujRjvr4UKEU/aQk4I03QnduhmEY\nVywl+llZQJ8+tJqlnrS0OG7v3esQfYZhGKOxlOgDtKKlnm67ooL2tn31VbqvdPoMwzBGY4k9cpWc\nPUuO/9Qp2qFKawoKaNLVoUPAHXcAL70EXLwIdO2q/bkYhmEC1U7LiT5AYvyNbwBPPKF926NGAX/8\nIzB0KJCbSxOz9u3T/jwMwzAAi75f7NoFfOc7tLRxbKx27e7bB9xyC1BeTjtgXbhAK2nyZCyGYfQi\nUO20XKYPkAiPHg0UFmrb7ssvA/PmObY87NWLBZ9hmPDCkqIPAA8/DKxdq117u3YBmzYBTz+tXZsM\nwzBaY8l4BwDOnaPcvbY2+M3IL18GbrwReOAB4Pvf16Z/DMMw/sDxjp/07UullYcPB9dORwctr5CR\nQf8yDMOEM5YVfQCYOBH49NPg2li/npZ2WL06+G8MDMMwemNpmZo4Efjss+DaeO014Ec/Arp106ZP\nDMMwemJ50Xd1+s8+S2WW/lBZCXzyCTBnjvZ9YxiG0QNLi/64ccCxY8Af/gC0tVFM86tfAcuXe37N\n888Dd95JWf6f/0wTveLjQ9dnhmGYYLBs9Y5k926amdunD/20tgJFRcDJk45tFTs6gP37gRMngAcf\nBNLTaX383btp8bYhQ0LebYZhGABcvRMw48YB//oXlW6+/jqwdCkt0SAXTANo56tp06i2//XXgb/9\njb4ZbN3Kgs8wjLmwvNOX1NUBH38M5OXRmvtz55Kbf+454N57abtD3vGKYZhwg9fe0Qi7ncown34a\nuPtu4Pe/N7pHDMMwnWHR15iGBiAuDoiJMbonDMMwnWHRZxiGsRA8kMswDMN4hEWfYRjGQrDoMwzD\nWAgWfYZhGAsRlOivXbsWI0aMQFRUFEpLS90ec+rUKUyePBkjRozAyJEj8XuL1j6WlJQY3QXdiOT3\nBvD7MzuR/v4CJSjRHzVqFNatW4cbb7zR4zExMTF46aWXcODAAezYsQOvvPIKDh06FMxpTUkkf/Ai\n+b0B/P7MTqS/v0CJDubFw4YN83lMWloa0tLSAAA9evTA8OHDUVVVheHDhwdzaoZhGEYFIc30Kyoq\nsHv3bkycODGUp2UYhmEkwgfTpk0TI0eO7PRTVFR05ZicnByxa9cur+00NDSI8ePHi3Xr1rl9HgD/\n8A//8A//qPgJBJ/xTnFxsa9DfNLW1obbb78d9957L2bPnu32GJ6NyzAMoz+axTueRFsIge9///vI\nzs7GE088odXpGIZhGBUEJfrr1q1DRkYGduzYgby8PNx8880AgKqqKuTl5QEAPv74Y/zlL3/BBx98\ngHHjxmHcuHHYvHlz8D1nGIZhAiegMEgH3n//fXHNNdeIzMxMUVBQYHR3NOXkyZMiJydHZGdnixEj\nRojly5cb3SVdaG9vF2PHjhW33HKL0V3RnAsXLojbb79dDBs2TAwfPlx88sknRndJU5YuXSqys7PF\nyJEjxV133SVaWlqM7pJq5s+fL1JSUsTIkSOvPHb+/Hkxbdo0MXToUJGbmysuXLhgYA+Dw937e/rp\np8WwYcPE6NGjxa233irq6up8tmPojFy73Y7HHnsMmzdvxsGDB/H2229HVA2/VeYoLF++HNnZ2bDZ\nbEZ3RXMef/xxzJgxA4cOHcLevXsjqtS4oqICK1euRGlpKfbt2we73Y41a9YY3S3VzJ8/v1OKUFBQ\ngNzcXBw5cgRTp05FQUGBQb0LHnfvb/r06Thw4AD27NmDrKws/PrXv/bZjqGi/9lnnyEzMxODBw9G\nTEwM5s6diw0bNhjZJU1JS0vD2LFjATjPUYgkTp8+jU2bNuG+++6LuMH4ixcvYtu2bViwYAEAIDo6\nGomJiQb3Sjt69uyJmJgYNDc3o729Hc3NzRgwYIDR3VLNt771LfTq1cvpsaKiIuTn5wMA8vPzsX79\neiO6pgnu3l9ubi66dCEZnzhxIk6fPu2zHUNFv7KyEhkZGVfup6eno7Ky0sAe6UekzlF48sknsWzZ\nsisfvEiivLwcffv2xfz583Httdfi/vvvR3Nzs9Hd0ozk5GQ89dRTGDhwIPr374+kpCRMmzbN6G5p\nSk1NDVJTUwEAqampqKmpMbhH+rFq1SrMmDHD53GG/qVGYhzgjsbGRsyZMwfLly9Hjx49jO6OZmzc\nuBEpKSkYN25cxLl8AGhvb0dpaSkeeeQRlJaWonv37qaOB1wpKyvDyy+/jIqKClRVVaGxsRFvvfWW\n0d3SDZvNFrGas2TJEnTt2hV33323z2MNFf0BAwbg1KlTV+6fOnUK6enpBvZIe/yZo2BWtm/fjqKi\nIlx11VW46667sHXrVsybN8/obmlGeno60tPTcd111wEA5syZ43FhQTOyc+dO3HDDDejduzeio6Nx\n2223Yfv27UZ3S1NSU1NRXV0NADhz5gxSUlIM7pH2vPHGG9i0aZPfF2xDRX/ChAk4evQoKioq0Nra\ninfeeQezZs0yskuaIiJ8jsLSpUtx6tQplJeXY82aNZgyZQrefPNNo7ulGWlpacjIyMCRI0cAAFu2\nbMGIESMM7pV2DBs2DDt27MClS5cghMCWLVuQnZ1tdLc0ZdasWSgsLAQAFBYWRpzx2rx5M5YtW4YN\nGzYgLi7OvxfpVV7kL5s2bRJZWVliyJAhYunSpUZ3R1O2bdsmbDabGDNmjBg7dqwYO3aseP/9943u\nli6UlJSImTNnGt0Nzfniiy/EhAkTAiqJMxMvvPDClZLNefPmidbWVqO7pJq5c+eKfv36iZiYGJGe\nni5WrVolzp8/L6ZOnRoRJZuu7++1114TmZmZYuDAgVf05eGHH/bZTthsjM4wDMPoT+SVXDAMwzAe\nYdFnGIaxECz6DMMwFoJFn2EYxkKw6DMMw1gIFn2GYRgL8f8BrQnZb60tDZUAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0xbf7e210>"
]
}
],
"prompt_number": 81
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment