{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "%pylab inline\n",
    "from network.network import (is_valid_network,\n",
    "                             get_number_of_activated_antenna,\n",
    "                             can_antenna_be_turned_on,\n",
    "                             pp_network,\n",
    "                             turn_off_antenna,\n",
    "                             turn_on_antenna)\n",
    "from network.gibbs import get_random_networks, estimate_number_of_activated_antenna\n",
    "from numpy import mean\n",
    "from random import Random\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "random = Random()\n",
    "antennas = []\n",
    "box_size = 8\n",
    "N = 500\n",
    "for _, network in zip(range(N), get_random_networks(box_size, random)):\n",
    "    antennas.append(get_number_of_activated_antenna(network))\n",
    "    \n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEZCAYAAACHCd7XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8XFV99/HPNyREFAhBTQ5ySQTLRVuBVhFFZRRRQCU8\nPpUKaokoPvYitvVCog/leAdtq1YfaVXEiKJGLAStkjSS8YYoahAkGG8QIpADhJBA8REkv/6x1iT7\nTGbOmTlnJjNnz/f9ep3X2bMva/9mz96/WbP23msrIjAzs6lvWq8DMDOzznBCNzMrCSd0M7OScEI3\nMysJJ3Qzs5JwQjczKwkn9C6T9GxJN+/E9f1M0nN3wnrOk3TJBJe9RdLzOxTHtrIkLZb0iU6Um8u7\nX9L8PHyxpHd1sOwLJb2jU+WNsZ5jJa3v9nomS9IqSWf2Oo6pbnqvA+hXkm4F5gB/AAQE8JmIOHuc\n5bYCT4qI3wBExHeBw7oU48XA+oj4x9q4iPjjbqyrib66iSEi3t/KfJJWAZdExKfHKW+PTsQl6Qzg\ndRHxnELZf9WJslvUV5+TdY8TenMBvDgiVk1gOZvCJO0SEY90ski8X3RMFz6f0nCTy9jUcKR0kKSq\npPsk3SXpC3n8t/IyN0jaIunl9T95cxPBWyTV5vmUpDmSvi5ps6QVkmYV5l8q6U5Jm/I6D8vjzwJe\nCbwtl7OsUH6tCWJXSR+WdLuk30r6kKQZedqxktZL+gdJI3mehU03hDQ/r3+zpOXA4+qmHy3peznO\n1ZKObWkDS0+XdE1e7nZJH5XUtKIh6dWSbpV0t6S3103b1gwkaaakSyTdk8v+gaTHS3oP8BzgY3m7\n/Wuef6ukv5b0C+AXhXEHFlbx+Pz5bMlNBAfk+ebleacVYlkl6UxJhwIXAs/MTTj35umjmnAknSXp\nlzneKyTtU5i2VdL/kfQLSfdK+tgY2+dRkj6T5/sZ8PS66aPe01hNSZLOkPQdSR/M5f1a0gmF6ftI\nWiZpY47tdXWfxdL8GWyR9FNJfyRpUd7fbpV0fN0qn5Q/p82SLpe0V932PVPSOuCbefyE9rlSiwj/\nNfgDbgGe32TapcDiPLwr8KzCtK3AEwuvjwVuqyv3GlJC3AcYAX4MPDWX9U3g3ML8C4FHAzOAfwFW\nF6ZdDLyrWdzAu/K6Hpv/vge8sxDXw8B5wC7AicB/A7OavOdrgA/mOJ4DbAE+m6ftC9wDvCi/Pi6/\nfux42xb4U+Ao0hfhAcBNwNlNlnsycD9wTI7jn4GHCmWdV4jp9cAyYGYu+0hg9zxtFXBmXdlbgeXA\nXsDMPO4R4MDCtt5cWPeHge/kafPyvNMK5W1bB3AG8O269W377IDnA3cDh+ey/xX4Vl1sVwJ7APsD\ndwEvbLKNzge+BczKn8uNjN7/tr2nZvtQYdoZwO+BM/M2fANwe2H6t4GP5pgPz3FVCp/Fg8ALSBXH\nJcCtwOK8v70O+E3d9lpPap7cDbiM1CxW275bgc/kaTOBJ9DGPjcof66hj+2KXDPZlP+/No9/GJgn\nad+IeCgirqlbrmHNvuCjEXFPRNwJfAe4NiJuiIiHgMtJyQeAiPhMRDwYEQ+TEvThklpt2z2dlMA3\nRsRG4J3AqwvTHwLeHRGPRMQ3gAeAQ+oLkbQ/8DTgHyPi4Yj4DvDVwiyvBP4zIpbnmL8J/Ag4abwA\nI+InEfHDSG4DPkH6smnkfwNfjYjv5e1xLs2bMh4mfYkdnMteHREPjBPO+yLivoj4fX5d/zn+Z2Hd\n7yDVuvcd7z224HTgooj4aS57cS77gMI874+I+yNiPSn5HdGkrJcD74mIzRFxO+nLoWi8fbPeuoj4\ndKSsuQTYR+kX5X7AM4Fz8j7xU+BTwF8Wlv1ORKyMiK3Al4G9gfMjNZd8EZgvac/C/JdExM0R8TvS\nZ3uqpFq8AZwXEb/Ln8+rmOA+V2ZO6GNbEBF7R8Ts/P+iPP6tpG33Q0k3SnpNm+WOFIZ/1+D17gCS\npkk6X9KvJN1HqtkGdc0dY3gCcFvh9bo8rmZjPthqHqytu0E5m/KBViyrZh7p4Lu39gVIqsnuwzjy\nz/CvKjUr3Qe8l+bv7wmkWhwAEfEgsLHJvJeQatxfVGpuukDSLuOE89txphfX/d/AvYzenhP1BArb\nM5e9kVTDrinuI80+p1pZxfexrsl8rdpQiKv2+e+e13Nv/gyK62oW8++Ae/IXQ+11raya4tU460g1\n/+K+UHxfE97nyswJfWwNazMRcVdEvD4i9iX9DP14XVtrp7wSeCmpSWEvYH6OqVhrGcsdpB2/Zl4e\n1647gdmSdiuMK9Ye15OaOvYufAHuEREfaKHsC4GbgYPye3wHzWuRd5KaHACQ9GhSLXwHEfGHiHh3\nRDwFeBbwErbXHpttt/G2Z3Hdu5NqnLeTmqogNY3VDLVR7qjPSdJjSO9rvC+YZmXtX3g9r276gzSP\ns9317J1jrTmAtD0mqj7uh0jNKDXF7TiZfa60nNAnQNKfF35q30dq36vVdDcAnUruu5PaMDflA+f9\njN6pR8ZZ1xeA/yvpcZIeR/oZ2/a147kp5EfAOyXNkPRs0hdNzeeAl0p6Yf5V8Silk66t1F73ALZE\nxIP5BOJYl/NdBrxE0rOUTu6+i+YnriuS/jifqHyA1ARTuzJivO3WzEl53bsC7wa+HxF3RMQ9pET2\nqvz+zwQOKiw3AuyXY27kC8BrJD1V0kzgfaRmuIlcP/5lYLGkvXKzyN/WTV8NnJ7jPIHmzVtjiojf\nks6rvF/pBPRTgdcygf2r4FWSDs1f1O8Evlyo0dd/zpPZ50rLCX1sX81n6Gt/X8njnw78QNIW4ArS\nSbxb87Rh4LP5Z+CfNyizvrY2Vu3ts6Qmk9uBn5EOoKKLgKfkdf1Hg/LeQ0rENwA/zcPvHWN9Y8Vy\nOnA0qSngXFJ7alooHdwLgLeTTu6tA95C8/2ruJ63AK/M2/LfSW2rjReKWAP8DSkB3pFjaVaLHSJ9\nAWwmnWhdRUoCAB8BXp6vzvhwg5gaxRmkk+HDeb1Hktpxa84C3kaqUR5GOgFdc3WOYYOkuxq8r2+S\ntul/kD7rJwKvaBJHs1hr3knaZ24BriLtQ0V/B5wMbAJOI52zaUdx3aflWO8AvkI6md/OZb712/cS\n0n51B+kCgTc1mXci+9xA0PYvQDMzm8oG+tvMzKxMnNDNzErCCd3MrCSc0M3MSsIJvYvkLkFbotF9\nsOzQL0qX1vlyScvzJYjdXM9OeT+9IumY3P/KXr2OxZzQdxqljo4u7lBZ9Z1GlUH9JWxdI+kIUv8k\nC3J3C91WykvJ8nXu7wFOioj7eh2Pufvcna1TB3YpE8TOEhHXkzojszap0HVtvhb8eT0OyQpcQ+8g\nScdLujl35vVRxugISTt2i/vkwrSLJX1M0tfyDU3fl/TEPG2HLnrz+JcodSG6SdJ3Jf1JobxbJL1Z\nqQvTTZK+UGtq0Djd6Cp1wftPktbleD+e72ZE0mOV+mHZlG/S+dYY7/fDkm5T6hr1uny3aSvbdNTT\njZo0zyzMZd8j6Q2Snpbf6735cyiWd6akNXnaVZKKt9wfqtQ97sb8Ob68MO0kSTflbb5e0j80iXda\n3l53S/oV8OK66Qvz+rco9dHz+sK08T6LVmMY9TSp+mYfpabAd+X9ZEveDns3KasW09sk3Ql8Oo8v\n7m/fk3R4YZlzlPrP2ZK34/PyeCl1n/urvH2+KGl2YbmxjolRzZfKXfs2inmgdbr7xkH9I/W9sQX4\nX6TuQf+OdLv5mU3mX8jY3eLeDfwZ6Uv3c8Clhen1XfQeSbq9/GmkZP9q0p2CM/L0W4Brgbmk7mHX\nAK/P08bsRhf4EOlu2FnAY0hd0r43T3sf8PEc4y7AMWNsn9PzuqcBf0/ql2XXPK3Y7e2ormip68a4\nwbxbcwy7AscD/5909+NjSR1IjQDPyfMvAH4JHJrjOBf4QZ72aNIdln+Zt+Hh+TM4NE+/g9xNct4W\nRzR5n2/I2/cJ+f1eXfd+TgTm5+Hn5G19RIufRasxbNtGTbbpqrwdDiJ1RbuK1NNko7JqMb2PtK/O\nJO1vd7G92+OFpDs1dwUOzttxbl7+APK+Srrz8xpSB1ozSP34FPfrhTQ/JkZ1eUyDLon9F07oHduQ\nKYleUzduPU0Set18e+XEtEd+fTHwicL0E4E1hddbGd2n9cfJ/ZwXxv2c7YnsFuC0wrQLgI/n4WNz\n0ij25T0CHJWHH2D0l8czyf1Yk24zv5zUsVa72+te4E/y8GQS+iPAUGH6PcCfF15fRu5fHfg68NrC\ntF1InVUdAJxKoQ/yPP3fyH3Tk/ryPqv2GY3xvr5J/rLMr4+nrq/0uvkvB97Y4mfRagytJPS3F6b/\nFfD1JmUdS/qSnFG3v727br61wHNJXxIbSP2TT6+bZw3wvMLrfUgdcO2wbdjxmHBCb+HPTS6dM6pr\n16xh50pqrVvcDYXhsbpLhXTAvlmjuxLdj9Fdu47V/WrDbnQlPZ5UY/pxrWzgG2zv4fCDwK+BFfm9\nnNMsQKWnNK3JP6c3AXvSejfA4yn2j/K7Bq9r73UecG6OYw2pf5z7SP2+zAOOrtuGp5N+1UDqi/3F\nwLr88//oJrHU7wejuq+VdKJSE9rGvI4TGb0dxurSuNUYWtHO/nV3pH7aa+YBtaarNUoPQd8dmBMR\nvyb9Oh0GRiRdKmmosNzlhX1pDan2P7fFY8LG4ZOinXMno7uUhdHdgRadzvZucW9TeuTcJtp/+EDN\nelIzSEsPSW7DPaSD/SmRHsYxSqQHRrwFeEtu71wl6YdR10GTUnv5W0m1szV53L209n7/m8509wpp\nO10SEZ+vn6B01VA1Il7UaMGI+DFwilKf6m8ElrLj5w11XfwyulvcXUm/GF4FLIuIrZIup8XPvY0Y\n6rfZZPsIrz8Jv570a7RhR28R8UVSP/S7kx5YcgGpRn0bqZb9/fplJL2KsY+JTu4HpeUaeuf8J/Bk\nSadI2kXSm9heu6u3B2N3izue+i56Pwm8QdJRkPrTzifQHtNw6RZF+m37SeDDubaOpH0lvTAPv1hS\nrZvY+4E/sL0b4aI9SDWxjUonWf8xj2ummOCuB14habqkpwH1PVi28yX4b8DbJT0lxz9L23vE/Bpw\nsKRX5XXNUDq5emgePl3SnpGu8Lif7V3x1lsKnJ2302yg+Ktl1/x3T07mJwIvbCXwNmO4HniupP1z\nYlzUyjraUNvfjs4nOrftb5IOlvS8/OX1EOkXUm2f+Hfgfdr+LNbHSzo5TxvvmLgeeJmk3SQ9idRV\nr9VxQu+QSI94ezmpNnIPqS3xe01mH69b3PEMU+iiN9fcziI9+Phe0kOOzyiG12b5xfkXAb8Crs0/\nhVeQTnwB/BGwUtL9pPf6/yKi0ZUuy/PfL0g/pR+kSXNUg/WfCzyJ1OZ+HlBfu65/b01fR8QVpJN7\nX8jv5QbghDztAVJyfQXp5OMdpOdz1m48ejVwS17u9aRfWY18kvRea90V17pcrq3jbODL+XN6Bekk\n81iK76elGCJiJfCl/P6uY/TjAuvLbFve315LerzdRkbvbzNJ2+1u0jZ8POmRepC6LV5GaqLbTNrv\nj8rTxjsmPkSqFGwgnWP6HLaDlrrPzbXN2hO9PxkR/5prH18i/aS8FTg1IjZ3K1AzMxvbuDX0/PP0\ntaRL4o4gPTHmIFLNbWVEHEK6NGtx81LMzKzbWmlyOYx0re7vc9vdt4GXkZ56UntqzRLglO6EaGZm\nrWglof8MeI6k2UrP+juJdBZ/bkSMAETEBmBO98I0M7PxjHvZYkT8XNIFwH+RbjJZTeOz65M60WJm\nZpPT0nXoEXEx6cwykt5LukJhRNLciBjJNw7s8PDbPL8TvZnZBEREW/emtHTZYuEa5ANIfZVcClxJ\n6nsB0iVLTS+/6vXtsK38nXfeeT2PwXHu3L+8d1L7cdnreKbytnScnf+biFbvFP2KUm9sDwN/HRFb\ncjPMUqUe0NaR+sIwM7MeabXJ5bkNxt0LvKDjEZmZ2YT4TtGsUqn0OoSWOM7BM1W2pePsvZbuFJ3U\nCqTo9jrMJkIS2y/O0oTbLc26QRLRjZOiZmbW/5zQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ\n3cysJJzQzcxKwgnd+trQ0HwkIYmhofm9Dsesr/lOUetr3byb03eKWj/znaJmZgPMCd3MrCSc0M3M\nSsIJ3cysJFp9BN1iSTdJukHS5yXtKmm2pBWS1kpaLmlWt4M1M7Pmxk3okuYBZwFHRsRTSU85Og1Y\nBKyMiEOAq4HF3QzUzMzG1koNfQvwEPAYSdOB3YDbgQXAkjzPEuCUrkRoZmYtGTehR8Qm4J+B20iJ\nfHNErATmRsRInmcDMKebgZqZ2djGfUi0pAOBvwfmAZuBL0t6JdvvyKhpelfG8PDwtuFKpVLqZ/qZ\nmU1EtVqlWq1Oqoxx7xSVdCpwfESclV+/GjgaeD5QiYgRSUPAqog4rMHyvlPUJsx3itqg6tadomuB\noyU9SukIOA5YA1wJLMzznAEsa2fFZmbWWS315SLpraTk/QiwGngdsAewFNgfWAecGhH3NVjWNXSb\nMNfQbVBNpIbuzrmsrzmh26By51xmZgPMCd3MrCSc0M3MSsIJ3cysJJzQzcxKwgndzKwknNDNzErC\nCd3MrCSc0M3MSsIJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJMZN6JIOlrRa\n0k/y/82SzpY0W9IKSWslLZc0a2cEbGZmjbX1xCJJ04DfAs8A/hbYGBEfkHQOMDsiFjVYxk8ssgnz\nE4tsUO2MJxa9APh1RKwHFgBL8vglwCltlmVmZh3UbkL/C+DSPDw3IkYAImIDMKeTgZmZWXumtzqj\npBnAycA5eVT979Omv1eHh4e3DVcqFSqVSssBmpkNgmq1SrVanVQZLbehSzoZ+OuIOCG/vhmoRMSI\npCFgVUQc1mA5t6HbhLkN3QZVt9vQTwO+UHh9JbAwD58BLGtnxWZm1lkt1dAlPRpYBxwYEffncXsD\nS4H987RTI+K+Bsu6hm4T5hq6DaqJ1NDbumxxIpzQbTKc0G1Q7YzLFs3MrE85oZuZlYQTuplZSTih\nm5mVhBO6mVlJOKFbxwwNzUcSkhgamt/rcMwGji9btI7pxmWAvmzRBpUvWzQzG2BO6GZmJeGEbmZW\nEk7oZmYl4YRuZlYSTuhmZiXhhG5mVhJO6GZmJeGEbmZWEi0ldEmzJH1Z0s2SbpL0DEmzJa2QtFbS\nckmzuh2smZk112oN/SPA1/NDoA8Hfg4sAlZGxCHA1cDi7oRoZmatGLcvF0l7Aqsj4qC68T8Hjo2I\nEUlDQDUiDm2wvPtyGRDuy8Wsc7rVl8sTgXskXSzpJ5I+kR8aPTciRgAiYgMwp/2QzcysU6a3OM+f\nAn8TET+S9CFSc0t9daZp9WZ4eHjbcKVSoVKptB2omVmZVatVqtXqpMpopcllLvD9iDgwv342KaEf\nBFQKTS6rcht7/fJuchkQbnIx65yuNLnkZpX1kg7Oo44DbgKuBBbmcWcAy9pZsVkZ+SEf1kstPeBC\n0uHAp4AZwG+A1wC7AEuB/YF1wKkRcV+DZV1DHxCuobvWb50zkRq6n1hkHTPVEuRUi9cGi59YZGY2\nwJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJJzQB5jvajQrF99YNMA6fRPMVLtRZ6rF\na4PFNxaZmQ0wJ3Qzs5JwQjczKwkndDOzknBCNzMrCSd0M7OSaOWZoki6FdgMbAUejoijJM0GvgTM\nA24lPeBic5fiNDOzcbRaQ99Ken7okRFxVB63CFgZEYcAVwOLuxGgmZm1ptWErgbzLgCW5OElwCmd\nCsrMzNrXakIP4L8kXSfpdXnc3PwAaSJiAzCnGwGamVlrWmpDB46JiDslPR5YIWkt2+9vrvE9zmZm\nPdRSQo+IO/P/uyVdARwFjEiaGxEjkoaAu5otPzw8vG24UqlQqVQmE7OZWelUq1Wq1eqkyhi3cy5J\njwamRcQDkh4DrADeCRwH3BsRF0g6B5gdEYsaLO/OufqUO+eaWvHaYJlI51yt1NDnApdLijz/5yNi\nhaQfAUslnQmsA05tO2IzM+sYd587wFxDn1rx2mBx97lmZgPMCd3MrCSc0M3MSsIJ3cysJJzQzcxK\nwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ\n3cysJFpO6JKmSfqJpCvz69mSVkhaK2m5pFndC9PMzMbTTg39TcCawutFwMqIOAS4GljcycDMzKw9\nLSV0SfsBJwGfKoxeACzJw0uAUzobmpkNDc1HEpIYGprf63Csz7XykGiADwFvBYrNKnMjYgQgIjZI\nmtPp4MwG3cjIOmrPKB0ZaevxkjaAxq2hS3oxMBIR1wNj7VF+Gq6ZWQ+1UkM/BjhZ0knAbsAeki4B\nNkiaGxEjkoaAu5oVMDw8vG24UqlQqVQmFbSZWdlUq1Wq1eqkylBE6xVrSccCb46IkyV9ANgYERdI\nOgeYHRGLGiwT7azDdh5JbP9hJSb7OXW6vG6V2c2yp8I2talBEhHRVjvbZK5DPx84XtJa4Lj82szM\neqStGvqEVuAaet+aCrVJ19BdQx9UO7uGbmZmfcQJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc\n0M3MSsIJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJJzQzcxKwgndzKwkWnlI\n9ExJP5C0WtJNkt6Xx8+WtELSWknLJc3qfrhmZtbMuAk9In4PPC8ijgSeCjxf0jHAImBlRBwCXA0s\n7mqkZmY2ppaaXCLiwTw4My+zCVgALMnjlwCndDw6M+u4oaH5SEISQ0Pzex2OdVBLCV3SNEmrgQ1A\nNSLWAHMjYgQgIjYAc7oXppl1ysjIOtJzSiMPW1lMb2WmiNgKHClpT2C5pArbn1y7bbZmyw8PD28b\nrlQqVCqVduM0Myu1arVKtVqdVBlq9yniks4Ffge8FqhExIikIWBVRBzWYP7wk8r701R4Qn03n3o/\nFeKdCjFad0giItTOMq1c5fK42hUsknYDjgdWA1cCC/NsZwDL2orWzMw6qpUml32AJUpf69OASyLi\nm7lNfamkM4F1wKldjNPMzMbRdpNL2ytwk0vfGvTmgakQ71SI0bqjK00u1h98qZmZjcc19CliKtTU\npkKM3S570LepdY5r6GZmA8wJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJJzQ\nzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJFp5BN1+kq6WdJOkGyWdncfPlrRC0lpJy2uP\nqTMzs95opYb+B+AfIuIpwDOBv5F0KLAIWBkRhwBXA4u7F6aZmY1n3IQeERsi4vo8/ABwM7AfsABY\nkmdbApzSrSDNzGx8bbWhS5oPHAFcC8yNiBFISR+Y0+ngzMysddNbnVHS7sBlwJsi4gFJ9c+tavoc\nq+Hh4W3DlUqFSqXSXpRmZiVXrVapVquTKqOlZ4pKmg58DfhGRHwkj7sZqETEiKQhYFVEHNZgWT9T\ntAOmwrMlp0KM3S570LepdU43nyn6aWBNLZlnVwIL8/AZwLJ2VmxmZp01bg1d0jHAt4EbSV/rAbwd\n+CGwFNgfWAecGhH3NVjeNfQOmAo1takQY7fLHvRtap0zkRp6S00uk+GE3hlT4cCeCjF2u+xB36bW\nOd1scjEzsz7nhG5mVhJO6GZmJeGEbmZWEk7oZmYl4YRuZlYSTuhmZiXhhG5mVhJO6GZmJeGEbmZW\nEk7oZmYl4YRuZlYSTuhmZiXhhN4lQ0PzkYQkhobm9zocMxsA7j63SwaxG9WpEGO3yx70bWqd4+5z\nzcwG2LgJXdJFkkYk3VAYN1vSCklrJS2XNKu7YZqZ2XhaqaFfDLyobtwiYGVEHAJcDSzudGBmZtae\ncRN6RHwX2FQ3egGwJA8vAU7pcFxmNkX4AoD+MX2Cy82JiBGAiNggaU4HYzKzKWRkZB21k6wjI22d\nw7MOm2hCrzfmafLh4eFtw5VKhUql0qHVmpmVQ7VapVqtTqqMli5blDQP+GpEPDW/vhmoRMSIpCFg\nVUQc1mRZX7Y4IJevTYUYu132IG5TXwbZHd28bFH5r+ZKYGEePgNY1s5Kzcys81q5bPFS4BrgYEm3\nSXoNcD5wvKS1wHH5tZmZ9ZDvFO2SqfCzdhBj7HbZg7hN3eTSHb5T1MxsgDmhm5mVhBM6vjHCzMrB\nbegMbjvlIMbY7bIHcZu6Db073IZuZjbAnNDNrO90uhl0UJpV3eTC4P6sHcQYu132IG7TQY2x29zk\nYmY2wJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSsIJ3cysJDr1CLqd5sorv86FF6bnU8+cOZ0L\nL/wn9tlnnx5HZWbWe5OqoUs6QdLPJf1C0jmdCmosn/3sUq66am+uuuplrFjxc6677rqdsVozs743\n4YQuaRrwMeBFwFOA0yQd2qnAxvYM4C+YMWPfDpZZ7WBZ3VTtdQC201V7HUCLqr0OYOBNpoZ+FPDL\niFgXEQ8DXwQWdCasXqj2OoAWVXsdgO101V4H0KJqrwPYqfqxf5jJtKHvC6wvvP4tKcmbmZXeyMg6\nav3DjIy01eVK10y5q1xmzpzBbrt9mD33fCkPPfQDZsyY0euQzMz6woR7W5R0NDAcESfk14uAiIgL\n6ubr/27NzMz6ULu9LU4moe8CrAWOA+4EfgicFhE3T6hAMzOblAm3oUfEI5L+FlhBarq5yMnczKx3\nuv6ACzMz2zm6dlJU0mJJN0m6QdLnJe3arXW1S9JFkkYk3VAYN1vSCklrJS2XNKsPY/yApJslXS/p\nK5L27GWMOaYd4ixMe7OkrZL27kVsdbE0jFPSG/M2vVHS+b2KrxBPo8/96ZJ+KGl1/v+0Hse4n6Sr\n8/F9o6Sz8/h+O4bq43xjHt9Xx1Gz7VmY3vpxFBEd/wPmAb8Bds2vvwT8ZTfWNcH4ng0cAdxQGHcB\n8LY8fA67FsdUAAADNklEQVRwfh/G+AJgWh4+H3h/P27LPH4/4CrgFmDvfowTqJCaDKfn14/r0zhX\nAS/MwycCq3oc4xBwRB7enXQu7dA+PIaaxdlXx1GzOPPrto6jbtXQtwAPAY+RNB14NHBHl9bVtoj4\nLrCpbvQCYEkeXgKcslODqtMoxohYGRFb88trSR92TzXZlgAfAt66k8Npqkmcf0VKOn/I89yz0wOr\n0yTOO4FabXcv4PadGlSdiNgQEdfn4QeAm0n7Yr8dQ43i3LffjqNmcebJbR1HXUnoEbEJ+GfgNtLO\nd19ErOzGujpoTkSMQNrAwJwexzOeM4Fv9DqIRiSdDKyPiBt7Hcs4DgaeK+laSat63ZQxhkXAv0i6\nDfgAsLjH8WwjaT7pF8W1wNx+PYYKcf6gblJfHUfFOCdyHHUloUs6EPh7UtPLE4DdJZ3ejXV1Ud+e\nLZb0DuDhiLi017HUk7Qb8HbgvOLoHoUznunA7Ig4GngbsLTH8TRzEfDGiDiAdFx9usfxACBpd+Ay\n4E25Zll/zPTFMdQgztr4vjqOinECjzCB46hbTS5PA74XEfdGxCPAfwDP6tK6OmVE0lwASUPAXT2O\npyFJC4GTgH79gjwImA/8VNItpJ+zP5bUN7W1gvWkfZOIuA7YKumxvQ2poWdExBUAEXEZfdDFRm5K\nvQy4JCKW5dF9dww1ibPvjqMGcU7oOOpWQl8LHC3pUZJEuvmo365RF6O/8a4EFubhM4Bl9Qv0wKgY\nJZ1Aak87OSJ+37OodrQtzoj4WUQMRcSBEfFEUh8/R0ZEzw9udvzMrwCeDyDpYGBGRGzsRWB16uP8\npaRjASQdB/yiJ1GN9mlgTUR8pDCuH4+hHeLs0+NoVJwTPo66eOb2rcBNwA2kEyQzenkmuS62S0kn\naX9Paud/DTAbWEn6MloB7NWHMf4SWAf8JP99vB+3Zd3039AfV7k02p7TgUuAG4EfAcf2aZx/Rmr7\nXQ18Px/YvYzxGFKTwPU5pp8AJwB799kx1CjOE/vtOGq2Pevmaek48o1FZmYlMeV6WzQzs8ac0M3M\nSsIJ3cysJJzQzcxKwgndzKwknNDNzErCCd3MrCSc0M3MSuJ/AL1lzC1zT2xRAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f88746b0208>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "hist(antennas, bins=100)\n",
    "title(\"Estimation de la distribution du nombre \\n d'antennes allumées dans un réseau\")\n",
    "savefig(\"/home/didier/Downloads/gibbs.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.4.3+"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}