Skip to content

Instantly share code, notes, and snippets.

@kimitoboku
Created November 13, 2017 14:21
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 kimitoboku/d4a86ae041c7395b1b200b9c1ae1e2d2 to your computer and use it in GitHub Desktop.
Save kimitoboku/d4a86ae041c7395b1b200b9c1ae1e2d2 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG3ZJREFUeJzt3X2sZVdZx/Hfw7QVExASZ4xNZ8Zp4phY8YV6U0pITCMQ\np62Z+QPUNkEtApOgVQxEU9QUrH9YJEFFq9gA4UWhVDRkhMGGhBoSk9aZ8ibTWjKpxU4l6VCgaBDq\n6OMf55x2u2atvdfeZ7/v7yeZzD3nrLPP2vfc+5x1n/Wstc3dBQCYl2cM3QEAQPsI7gAwQwR3AJgh\ngjsAzBDBHQBmiOAOADNEcAeAGSK4A8AMEdwBYIYuGOqFd+/e7QcOHBjq5QFgku67776vuPueqnaD\nBfcDBw7o5MmTQ708AEySmX0ppx1pGQCYIYI7AMwQwR0AZojgDgAzRHAHgBmqDO5m9m4ze8zMvpB4\n3Mzs7WZ22sw+b2aXt99NAEAdOSP390g6VPL41ZIOrv8dlfTn23cLALCNyuDu7p+S9NWSJkckvc9X\n7pH0XDO7uK0OAsCcHLjpYzpw08c6f502FjFdIumRwu0z6/u+HDY0s6Naje61f//+Fl4aAKahGNAf\nvvXazl+v1wlVd7/d3XfcfWfPnsrVswAwebGR+lRG7o9K2le4vXd9HwAsVlkA72Pk3kZwPybpRjO7\nQ9ILJD3h7uelZABgCfoYleeoDO5m9kFJV0nabWZnJL1J0oWS5O7vkHRc0jWSTkv6pqRXdtVZABir\nOkH9wE0f63z0Xhnc3f36isdd0q+01iMAmJAmI/XZTagCwFxsU9I4lQlVAFiMNgLzVCZUAWD2xjJR\nmou0DACUCNMvD996be2Rdx8j9RDBHQAitsmphx8Am+P0GeRJywBAQRvpl9QxNvePohQSAJZgajn1\nKqRlACxa0/TLNiNvqmUAoANdjdIfvvXa0fwFQHAHsBhdB96xBHaJ4A5gAboIumMK5DHk3AHMVpdX\nPaqTNy+27St1w8gdwOyMaVQdBvM+yiAlgjuAGRlTUN8Yqk+kZQBMXl8Xnc4VjszZfgAAamg7qKeC\ncN38+hDXTA0R3AFMTjGoNx0Vx55XtW1AeF9sE7HY/UNMqNrqQkr929nZ8ZMnTw7y2gCmKdydsSpI\n1gmkuW1jG4LVsW2Kxszuc/edqnaM3AGMXjhS3wTIzdepgFkWfNvKgxf7kvMafaVoGLkDGK1wpF6n\nfV25fwnkvmbV8Zp+uDByBzBpsZF628pG3WV9St2OPRb766KP6hlG7gBGp8lkaV/pjtgE6rbHqIOR\nO4BJa2N0W3WMbUfTdSZrN/+z/QCAxem7HjzcFqCu3EBdvALT5nldY+QOYFCbSphtAl+d+vSmr1Gl\nbFVq3/l2iZE7gAGEQTdVhdJVIOxqC+BNfzcfVsXb0tMjfTYOAzAbOQG9+FjuaH6oidSUMJDHHusD\n1TIAOlMnoOcco41KlbrCIF21QrUsD9/GyD23WobgDqB1dRcf1T1mURfVJ1UBum7fYumZpgjuAHrV\nxii97muFQbNK3Q+CqqAcO+eqVanbBvhW69zN7JCZPWhmp83spsjj+83sbjP7jJl93syuadJpANMT\n2/el6wnDpvXiqTRKbvuq16sK3H1NpkoZwd3Mdkm6TdLVki6TdL2ZXRY0+x1Jd7r78yVdJ+nP2u4o\ngHFJbeY1RD9S6uxHU7X5WGrUnpOmCfP0Y6lzv0LSaXd/yN2flHSHpCNBG5f0XeuvnyPp39vrIoAx\nGUtQb/v1q+YJUkE5Z36hzxH7Rk5wv0TSI4XbZ9b3Fb1Z0ivM7Iyk45J+tZXeARiNMQX1UJ0tf6tK\nMMseSwXyqtx87CLZXWtrher1kt7j7nslXSPp/WZ23rHN7KiZnTSzk2fPnm3ppQF0baxBvSh3FWjO\nTo6xY9cJyLEJ33Af+q5VVsuY2Qslvdndf2p9+42S5O6/X2hzStIhd39kffshSVe6+2Op41ItA4xf\nFyWNXatbe14mp7olt3KnjUoZqd1qmROSDprZpWZ2kVYTpseCNv8m6cXrF/5BSc+UxNAcmLApjNZj\n2rr6USyVUlU9U5Zb73tTtMrg7u7nJN0o6S5JD2hVFXPKzG4xs8PrZm+Q9Boz+5ykD0q6wYcqoAew\nlTYuPj20sk28cp+X8z2IVcGEz+szFVPEIiYAT5l6UC+qu7CpaluBTdvwNXL/Umjre8oKVQC1zCmw\nb9RdjRrOMVQF6Zzjx44zlpw7gJmbY2CXtjufMNA3PVbfufYNtvwFFmyK1TBdaVIimWuI7y0jd2Dh\nplYNU1eTydSurpwUq7jpCjl3YKGGWBI/pLr58bIVqeExc/eX6bPOnbQMsEBD5YHHpmzhUaqmvcl+\n7mVtukJaBliYuU6eVkkF5dieL7l18rEtBXKu1NQHgjuwIEsN7Bt1FiQVVU22xkboqb1k+hrBE9yB\nhVlqYK/SJOgWA3s4CVv8q2CIqiSCO7AQ5NnryQ3CsRF7atFS0w3MmiC4AwvCqD1fk2utpm4PgeAO\nLACj9pW6l+SrCtpV2/zGXp+0DIBWjWE0OWapvHnYpux2mGsPkZYB0BpG7flyyiJDm6qYsq19czYk\naxvBHVgARu35K1Rjz8sN3LGRfN094ttCcAeAtVTJYu7+7pu2xVz8EKN2ieAOAFGxKyqlcu7h1avK\nLuLBhCqArZFvj8tdqSrlBeScQF63smZbbBwGzBz59vrbCdS5ryzVMuSHKyN3AChRtWlYbk69bOK1\nC4zcgZkiJbO9qgnQWGBPpV/63jyMkTswA5sJveI/rFSVH5alraq+j6mAXVY+2dekKldiAiYod5/w\npV1tKabuFZiq2pUdv2rStI33gisxATMy5EUfpixn0rPYbiN1ub3ww7Iq3546Th/vHcEdGJm2R32M\n3vNG100viZfaTyYW2Ittu35PCO7AwLoclfe5InIq6gbWqgttlAXwsA0jd2DmyiopMIycD8GcD4aq\nbYVZxATMCDnz/oWj5FQ1S6qcsSoIl+Xby3L2faFaBujAmIJ5n6mAMapbBVO3/DH1gVB2ub1tUC0D\n9GhMwTy0xLx7Tm177orRpt+/1JYFbBwGjFy4YKh40YaxBPaipQX4jbLNu4oLjcL3LXejr9w0zuaY\no9ry18wOmdmDZnbazG5KtPlZM7vfzE6Z2Qfa7SYwDmUBfczG3r+upKpWwtuxnHhVSiXWNlyVOuT3\nvTItY2a7JN0m6aWSzkg6YWbH3P3+QpuDkt4o6UXu/jUz+56uOgz0aczpliaWVvOemtgs3o6VLtbN\nn8cWOQ2dDsvJuV8h6bS7PyRJZnaHpCOS7i+0eY2k29z9a5Lk7o+13VGgL3MtUxw62PQllW8P0yc5\ntefhoqNYxU1YeVO2YKrP739OcL9E0iOF22ckvSBo8wOSZGb/KGmXpDe7+9+HBzKzo5KOStL+/fub\n9BfoxFwDeszSRu9lYsE+NfLOWdFaZ6K2a21Vy1wg6aCkqyTtlfQpM/thd/96sZG73y7pdmlVCtnS\nawONLCmgb8x99J6qZQ8f32YnyBxlk7HhaL8rOcH9UUn7Crf3ru8rOiPpXnf/b0n/amZf1CrYn2il\nl0BLlhjQY5Yyek/VoZft3xNb4FT1oZibfukrsEt51TInJB00s0vN7CJJ10k6FrT5iFajdpnZbq3S\nNA+12E+gsalWuHQlt8RvLsoCbupnIFVVE95XlWsvq7DpWmVwd/dzkm6UdJekByTd6e6nzOwWMzu8\nbnaXpMfN7H5Jd0v6DXd/vKtOA1UI6OXmGOBzg3fOlr052/PGKmpi91c9rytZOXd3Py7peHDfzYWv\nXdLr1/+AQZByqWfu+feY2CIl6fyRd7FN6ntUVSq5aVP8IO0zLcP2A5g0Avr25pB/b7I6tCwv3vS4\nZbX0m+Ox/QBQgpRLO+aWnsmtggm3HSjeF/5s1Xnt8Gcytrgp7EtX2BUSk8EovTtNgtlY1Mlvl6VZ\nUh90qdx62Wun9rBp4/ucuyskI3eMHqP07s1tBL+RSoPERu5Vq0urjh0ea2gEd4xSqtoF3RlTYMpV\n1ddYPjy2p0wsPVP8Oqdksqxap/g6fSG4Y1QYpQ9rigG+qtxxI1W7XgzyqfY5k6/hfbHcfZ8BnmoZ\nDI5c+riEI9qpSa1EDYN46oOsznnn5tCHqJhh5I7BMEofr1Qeeqxik5xlP085E6ubY5X9TMY+IMIP\nlKG+lwR39I5c+jSMbYIwFI6GpbzVp2Vi5Yu5I+3UIqawb339vFMKid6kyskwfmMrlawK5mWqthaI\nTcDGXrMY/OvqoxSS4I5OkU+fj7EE+NTPVG69eyxYl43yU1U2qdfJ6es2qHPHoMinz88Y8vCpdF6q\nDDGVE68qWyzel5vOifUjLLHs83tHcEeryKfP25B5+Jy/HHIDafizmVqFWqeOPdbXWB/IuWNSyKcv\nT5/veSqwt/UBk5uaKXtuLEefel4fOXdG7tgKI/Xl6nsUn7tYKef5uR8SYbuyCdjYaH3IKiOCOxoh\nqGOj61x8W9UoTVI0YVlk3dcd8veC4I5aCOqI6WoUXyfPHnteWZtU+7pVQWUTu6kPhT5G9AR3ZCGo\nI0cXo/imqZewP1UfEOEiozrljmXHTq1Y7RoTqijFRCma2vZnp84Ius4HSWrVac6kaJ2/TsomZZlQ\nxaAYqWMbYaqmz5WcqaX+YT49lUbZPFb1+sXXCb9OlUf2hV0hcR5G62hTLMD3VToZuz8cfZcF3aoc\neuy5VSWTfWHkjqeQV0eX6ky61g2CVVUqZaWPOatdU8dJPadqkRQTqugFQR19yk3VbJOSqZMSqfrZ\nj6V4tp0H6ON3jOC+cAR1DCGWj992NLvNXwWxVGTd0s6ylNMmndPn7xjBfaH6zH8CKZuA10YJZVnQ\nrhp1p9IyYf9SUimn3HRPF5hQXRgmSzFWXSyCCo8ZliOm0jexdjm17DmLpXKO1wbq3BeEkTqmomnw\nq1t/HlvklFrdmlMpk3qt2P9NUeeOp5CCwVKlft7DIB2rbY+1KTt2KvUT1tb3VRJJWmbGSMFgqYrB\nOhRLvaTq38sCcVlFTlV1DtUyaIwqGEzdNqPbOgOb3BF5Kl1T1Y++J1I3soK7mR0yswfN7LSZ3VTS\n7mVm5mZWmQ9CN0jBYA7a+tlNbUMgnT8yrzPSb9KH1O2uVAZ3M9sl6TZJV0u6TNL1ZnZZpN2zJb1O\n0r1tdxJ5GK0D50vtExML1lWrSqsmVMM9ZsJ2Y6tzv0LSaXd/yN2flHSHpCORdr8n6S2SvtVi/5CJ\n0TrmqM2URk7Ne+6GYbHnhs+PBfWx7S1ziaRHCrfPrO97ipldLmmfuw93TamFIg2Dueo6IJaN2nOD\nfOyDoljumPoA6MPWE6pm9gxJb5P0hoy2R83spJmdPHv27LYvvXikYTB3dX6uc7cVKN4OA3hsJF9H\naoOwIVIzOcH9UUn7Crf3ru/beLak50n6BzN7WNKVko7FJlXd/XZ333H3nT179jTv9cIxWsfStLXv\nzOZY4URq8f+c2vjY/WWrXcsmdruSE9xPSDpoZpea2UWSrpN0bPOguz/h7rvd/YC7H5B0j6TD7s7y\n0w4wWsfS5KRn2kp1lE2S5r5G+LsZ5uFHk5Zx93OSbpR0l6QHJN3p7qfM7BYzO9x1B/E0RutYqqb5\n93BDstheM2UDpXDlaur1w9dJpXv6/N1lb5kJYKUpsBLbE6ZKrHIlFaxTpY6x/WhyXp+9ZZBEGgZ4\nWhu142WlipvHw2qX1F/NVSPytgJ6EwT3ESMNA5yvyai9GIRzqmZixy/7PcwZwW/ajalaBgMgsANp\nub8XsYqYVP49Z+uBsoqYzf9lo/jU63SB4D4ylDkCecrq03NU1bgX7yv7nYyldMLXSX0odIktf0eE\noA7U10agDCdMUyWRTY670XdhBMF9JAjsQH1Nd2oMhWmaWBqmLF+e6sOQlW6UQo4AgR3YTln9edko\nPFbNEhvB5/xu5pZV5h4vhVLIiSCwA9trUg5ZnOAMA3nsQ6FsdWlZtU1YrdMXgvtAmDgF2rXNDoxV\n+8LkHCeWghkqsEvk3AdBUAe6kVqNWta+TNN8fixX3/fvPTn3nhHYge6lKl/KRvc5i6NiqZc62xGU\nHStXbs6dkXuPCOxAv8rqy8vSJbmpmKbVOn3EAHLuPSGwA/2pU91SXLXaZBK07PGclatdIS3TAwI7\nMIxUGWSdssXY88N2dUfwlELOAIEdGE5VBU3Z/jJVWw1sRuPFOvkx/Z4T3DtEYAeGV7ZNb2rxUu6+\nNWFVTO4eNH0guHeEwA6MU2ynyDDgpx4PpdIxsbRPX7tBbhDcO0BgB8albj17G7Xy4YRs34uZCO4t\nI7AD41SngiacYC3bMCznd75s64KuENxbRGAHxq3O72asgqYsQKcWNPUd1DcI7i0hsAPzEY60UymV\nWB17kx0ku0BwbwGBHZiOpqN3KV03v3ms7EMhtplYlwjuWyKwA9NTp1SxbNOvqtRLuGlYn1ihugUC\nOzBdORUwuZfbK9uorOxCIk2wcVjHCOzAtJUF3tw69Vggr/vB0BXSMg0Q2IF5yP0dzg3suc/rI01D\nWqaBviZEAPSj6aZfTYM0G4eN0FA1qwC6k7vKtOz5bbRpE8G9BtIxwDKU7RYZa7dpU9wpMvZ4n4ND\n0jKZCOzA/DXZQ6ZOuWMb11IlLdMiAjuA1ORoaj/32GrWPkfuWaWQZnZI0h9L2iXpne5+a/D46yW9\nWtI5SWcl/ZK7f6nlvg6CwA4sR1mdetUCptxY0VcsqQzuZrZL0m2SXirpjKQTZnbM3e8vNPuMpB13\n/6aZvVbSH0j6uS463CcCO7A8VfXvqf1lim1ix6v7IbCtnLTMFZJOu/tD7v6kpDskHSk2cPe73f2b\n65v3SNrbbjeHQ2AHliu80lKqTWyytGwvmbHsLXOJpEcKt8+s70t5laSPb9OpMaDkEViuJpuLhSP+\nnA+GLrU6oWpmr5C0I+mticePmtlJMzt59uzZNl+6VaRjANRZvRrL1YfHSH0IdCUnuD8qaV/h9t71\nff+Pmb1E0m9LOuzu344dyN1vd/cdd9/Zs2dPk/52jsAOoK6wWqaqLn4saZkTkg6a2aVmdpGk6yQd\nKzYws+dL+gutAvtj7XezXwR2AFLzWFBWBtlXiqYyuLv7OUk3SrpL0gOS7nT3U2Z2i5kdXjd7q6Rn\nSfprM/usmR1LHG7UyLMDaCrMsQ8R0Iuy6tzd/bik48F9Nxe+fknL/eod6RgAMXVz5GWpmXDRU5dY\noSoCO4BydSZXq/Zx72tXWS7WsUZgB1CmagQfpmFiWwP3GWcWH9zJswOoI/fqTbGA3meQX3RahnQM\ngDrqxopU/XsfFhvcCewA2la1z0yfV3FbbHCXCOwA6iuLG1Wj9HATsS4tMriTZwfQtdwLZndlcROq\npGMAbCtVOZNzpaW+UjOLHLkT2AFsKxZHynaC7PsaqosK7qRjAPRpyO0HFhPcSccAaFtOPBlqULmY\n4C4R2AG0r+qi2OGFs/uyiAlV0jEA+pBapTqE2Y/cSccA6Fqqpn1Isw/u0vDfZADzF1ud2ueVl0Kz\nDu6kYwAMKYxBm/w7K1RbwKgdwNBiV2nq2myDO6N2AH2ryr33Wfc+y+DOJCqApZtlcJcI7ACGUbXt\nb6xNF2YX3EnHABhabDFT32YX3CVG7QDGhb1ltsSoHcBYDL2waTbBnUlUAGNHtUxDBHYAY5KKSUyo\nZiIdA2BKGLnXwKgdwBgNlXuffHBn1A5gahi5Z2LUDmDMwgt2MHKvwKgdwFSUXTy7C1nB3cwOmdmD\nZnbazG6KPP4dZvah9eP3mtmBtjuawqgdwBSMbuRuZrsk3SbpakmXSbrezC4Lmr1K0tfc/fsl/aGk\nt7Td0RCjdgBTM7ZFTFdIOu3uD7n7k5LukHQkaHNE0nvXX39Y0ovNzNrrZhyjdgBTNJa0zCWSHinc\nPrO+L9rG3c9JekLSd7fRwRhG7QCmqq9B6QW9vMqamR2VdFSS9u/f3/g4jNgBTNkocu6SHpW0r3B7\n7/q+aBszu0DScyQ9Hh7I3W939x1339mzZ0+zHgMAKuUE9xOSDprZpWZ2kaTrJB0L2hyT9Ivrr18u\n6ZPu7u11EwBQR2Vaxt3PmdmNku6StEvSu939lJndIumkux+T9C5J7zez05K+qtUHAABgIFk5d3c/\nLul4cN/Nha+/Jeln2u0aAKCpSa9QBQDEEdwBYIYI7gAwQwR3AJghgjsAzJANVY5uZmclfanh03dL\n+kqL3RkS5zI+czkPiXMZq23O5fvcvXIV6GDBfRtmdtLdd4buRxs4l/GZy3lInMtY9XEupGUAYIYI\n7gAwQ1MN7rcP3YEWcS7jM5fzkDiXser8XCaZcwcAlJvqyB0AUGK0wd3M3m1mj5nZFxKPm5m9fX1R\n7s+b2eV99zFXxrlcZWZPmNln1/9ujrUbAzPbZ2Z3m9n9ZnbKzF4XaTP69ybzPCbxvpjZM83sn8zs\nc+tz+d1Im8EuYl9H5rncYGZnC+/Lq4foaw4z22VmnzGzj0Ye6/Y9cfdR/pP0E5Iul/SFxOPXSPq4\nJJN0paR7h+7zFudylaSPDt3PzHO5WNLl66+fLemLki6b2nuTeR6TeF/W3+dnrb++UNK9kq4M2vyy\npHesv75O0oeG7vcW53KDpD8duq+Z5/N6SR+I/Rx1/Z6MduTu7p/Sam/4lCOS3ucr90h6rpld3E/v\n6sk4l8lw9y+7+6fXX/+HpAd0/jV1R//eZJ7HJKy/z/+5vnnh+l84mTbIRezryjyXSTCzvZKulfTO\nRJNO35PRBvcMORfunpIXrv8U/biZ/dDQncmx/jPy+VqNroom9d6UnIc0kfdl/ef/ZyU9JukT7p58\nT7yHi9hvI+NcJOll65Tfh81sX+TxMfgjSb8p6X8Tj3f6nkw5uM/Jp7VaUvyjkv5E0kcG7k8lM3uW\npL+R9Ovu/o2h+9NUxXlM5n1x9/9x9x/T6hrHV5jZ84buU1MZ5/J3kg64+49I+oSeHv2Ohpn9tKTH\n3P2+ofow5eCec+HuSXD3b2z+FPXVVa8uNLPdA3crycwu1Cog/pW7/22kySTem6rzmNr7Iknu/nVJ\nd0s6FDyUdRH7MUmdi7s/7u7fXt98p6Qf77tvGV4k6bCZPSzpDkk/aWZ/GbTp9D2ZcnA/JukX1pUZ\nV0p6wt2/PHSnmjCz793k2szsCq3el1H+4q37+S5JD7j72xLNRv/e5JzHVN4XM9tjZs9df/2dkl4q\n6V+CZpO4iH3OuQTzN4e1mi8ZFXd/o7vvdfcDWk2WftLdXxE06/Q9ybqG6hDM7INaVSvsNrMzkt6k\n1eSK3P0dWl3T9RpJpyV9U9Irh+lptYxzebmk15rZOUn/Jem6Mf7irb1I0s9L+ud1XlSSfkvSfmlS\n703OeUzlfblY0nvNbJdWH0B3uvtHbZoXsc85l18zs8OSzml1LjcM1tua+nxPWKEKADM05bQMACCB\n4A4AM0RwB4AZIrgDwAwR3AFghgjuADBDBHcAmCGCOwDM0P8BDv/PHbYjQPoAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x110f3db38>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"import random\n",
"\n",
"a1 = 1.0\n",
"a2 = 4.0\n",
"\n",
"xs = []\n",
"a = []\n",
"for i in range(100000):\n",
"\tan = a1 + (a2 - a1) * float(i) / 100000\n",
"\tx = random.random()\n",
"\tfor j in range(1000):\n",
" \tx = an * x * (1 - x)\n",
"\ta.append(an)\n",
"\txs.append(x)\n",
"\n",
"plt.plot(a,xs,\",\")\n",
"plt.show()"
]
}
],
"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.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment