opt_1d/Golden Section Proportions.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "# Proportions of the Golden Section"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 2,
      "metadata": {
        "collapsed": false
      },
      "outputs": [],
      "source": [
        "import matplotlib.pyplot as pt\n",
        "from math import sqrt"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 3,
      "metadata": {
        "collapsed": false
      },
      "outputs": [],
      "source": [
        "a = 0\n",
        "b = 1\n",
        "\n",
        "m1 = a + (1-(sqrt(5)-1)/2) * (b-a)\n",
        "m2 = a + (sqrt(5)-1)/2 * (b-a)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 7,
      "metadata": {
        "collapsed": false
      },
      "outputs": [
        {
          "data": {
            "text/plain": [
              "[<matplotlib.lines.Line2D at 0x7fa0717e64e0>]"
            ]
          },
          "execution_count": 7,
          "metadata": {},
          "output_type": "execute_result"
        },
        {
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD/CAYAAAD4xAEfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAElhJREFUeJzt3H2QXeVdwPFvQl54y2bRggVlurLQKeMU14ITaDBsWxIK\ngcrEUQuOZnGoFJ2Ewoy1JbyUOnRaHShNhoIiWCiOVqalE9lpEhzZRoLlpSOl6lCalZ0W5bU22Qgm\nkLD+8ZzNntzsZs+959l77j7n+5m5k/Pc59x7Hn6cPb+9v9+5C5IkSZIkSZIkSZIkSZIkSaqZucCd\nwGPAI0Bvw/xFwBPZ/OW55z+dPfcksHrmlylJim0VcE+2vQT4Zm5uPvBDYHG2/QRwHNAPbMz2OQq4\nqR0LlSQdbF6J1y4FNmXbjwNn5OZOBbYDO7Pxo8Ay4H3A9wnJogv44xLHlySVUCYBdAGjufE+Qlno\n7WxuZ25uF+HTwDuAdwErgZMInwbeU2INkqQWlUkAo8Ci3Hj84g/h4p+fWwTsAH4CPAvsBZ4DdhOS\nwmsl1iFJasHcEq/dBlyQbZ8JPJObexY4BTgGWEAo/zxGKAV9ONvnBEIf4CeNb3zCCSeMAT58+PDh\no/hjO200B7iDkAi2Ae8GLgE+ls1fSGj+PgVcmXvdF3LPL5/ivccUx4033lj1EpJiPOMynvEQkkBT\nypSAxjjwwg6hrDPuoezR6E9KHFNNGhkZqXoJSTGecRnPapUpAUmSZjETQOIGBgaqXkJSjGdcxrNa\nc6pewBSykpYkqYg5c+ZAk9d0PwEkbmhoqOolJMV4xmU8q2UCkKSasgQkSQmwBCRJKswEkDhrrHEZ\nz7iMZ7VMAJJUU/YAJCkB9gAkSYWZABJnjTUu4xmX8ayWCUCSasoegCQlwB6AJKkwE0DirLHGZTzj\nMp7VMgFIUk3ZA5CkBNgDkCQVZgJInDXWuIxnXMazWiYASaopewCSlAB7AJKkwkwAibPGGpfxjMt4\nVssEIEk1ZQ9AkhJgD0CSVJgJIHHWWOMynnEZz2qZACSppsr0AOYCXwZOA/YAlwPDufmLgOuBvcA9\nwF/l5o4Dvgt8CHhukve2ByBJTWilBzCvxPEuBhYA7weWALdkzwHMB24FzgDeALYBG4FXsrm/AF4v\ncWxJUkllSkBLgU3Z9uOEi/24U4HtwE7gLeBRYFk29+fAHcCLJY6tgqyxxmU84zKe1SqTALqA0dx4\nX+79uggX/3G7gMXAAPAqsCV7vlNvQ5Wk5JW5AN8CfAd4IBv/GDgx234v8HlgZTa+lVAGWguMZY8+\n4AfArwMvN7y3PQBJakK7ewDbCI3eB4AzgWdyc88CpwDHEGr9ywiln6/n9nkEuIKDL/4ADAwM0NPT\nA0B3dzd9fX309/cDEx8bHTt27Liu4/HtkZERWlXmE8AcJu4CArgMOB04GrgLuBC4gVAWuptQ988b\nTwDeBTSDhoaG9p84Ks94xmU842n3J4Ax4MqG5/IX84eyx1Q+UOLYkqSSOrUJ6ycASWqCfwtIklSY\nCSBx+YaRyjOecRnPapkAJKmm7AFIUgLsAUiSCjMBJM4aa1zGMy7jWS0TgCTVlD0ASUqAPQBJUmEm\ngMRZY43LeMZlPKtlApCkmrIHIEkJsAcgSSrMBJA4a6xxGc+4jGe1TACSVFP2ACQpAfYAJEmFmQAS\nZ401LuMZl/GslglAkmrKHoAkJcAegCSpMBNA4qyxxmU84zKe1TIBSFJN2QOQpATYA5AkFWYCSJw1\n1riMZ1zGs1omAEmqKXsAkpQAewCSpMLKJIC5wJ3AY8AjQG/D/EXAE9n85dlz84GvAluBx7N9NIOs\nscZlPOMyntWaV+K1FwMLgPcDS4BbsucgXOhvBc4A3gC2ARuBC4BXgd8FjgGeBv6hxBokSS0q0wO4\nhfBb/N9n4xeAX8i2TwO+AJyfjW8lfBL4VnbM/wV+lvAJofGTA9gDkKSmtNIDKPMJoAsYzY33EcpC\nb2dzO3Nzu4DFwOvZeBHwALCuxPElSSWUSQCjhAv5uPGLP4SLf35uEfDTbPtE4BvA7cDfTfXmAwMD\n9PT0ANDd3U1fXx/9/f3ARN3Q8fTjfI21E9Yz28fG03h2ynh8e2RkhFaVKQGtIjRxLwPOBK4HVmZz\n84F/J/QGXieUfy4iJIgh4A8JjeOpWAKKZGhoaP+Jo/KMZ1zGM55WSkBlEsAc4MuEej+ERHA6cDRw\nF3AhcAPhk8HdwB3Al4DfBH6Qe5/zgd0N720CkKQmtDsBzCQTgCQ1wS+C6SD5eqHKM55xGc9qmQAk\nqaYsAUlSAiwBSZIKMwEkzhprXMYzLuNZLROAJNWUPQBJSoA9AElSYSaAxFljjct4xmU8q2UCkKSa\nsgcgSQmwByBJKswEkDhrrHEZz7iMZ7VMAJJUU/YAJCkB9gAkSYWZABJnjTUu4xmX8ayWCUCSasoe\ngCQlwB6AJKkwE0DirLHGZTzjMp7VMgFIUk3ZA5CkBNgDkCQVZgJInDXWuIxnXMazWiYASaopewCS\nlAB7AJKkwkwAibPGGpfxjMt4VqtMApgL3Ak8BjwC9DbMXwQ8kc1fXvA1kqQ2KdMDWAVcCPw+sAT4\nNHBxNjcf+A/gDOANYFu279mExHDZJK/JswdQ0uDgVtav38KePfNYuHAva9euYOXKZVUvq222Dg6y\nZf165u3Zw96FC1mxdi3LVq6seln7dfr6Zlrdz8+Z0EoPYF6J4y0FNmXbjxMu9uNOBbYDO7Pxo8Ay\n4CzgW1O8RpEMDm7lqqs2Mzx88/7nhofXAdTih2zr4CCbr7qKm4eH9z+3LtvuhItsp69vptX9/Owk\nZUpAXcBobrwv935dTFz8AXYBi6d5jSJZv35L7odrCIDh4ZvZsOHhytbUTlvWrz/g4gpw8/AwD2/Y\nUPq9Y9SsZ3J9s0Hdz89OUuYTwCiwKDeeC7ydbe9smFsE7JjmNQcYGBigp6cHgO7ubvr6+ujv7wcm\nfggdTz5++eUXCD9YYTz+Q7Z792Edsb6ZHr/w8suT/NfDYbt3u74OGNf9/Iw1Ht8eGRmhCquAv862\nzwQGc3PzgeeAY4AFwFPA8dO8Jm9MrVuxYt0YjB30OO+866peWlusW7Hi4P94GLvuvPOqXtrY2Fjn\nr2+m1f38nClA043TMuWXB4HdhAbvLcDVwCXAx4C3gGuAzYQ7fu4GXpziNYps7doV9PauO+C53t5r\nWbNmeUUraq8Va9eyrvfAG8yu7e1l+Zo1Fa3oQJ2+vplW9/Ozk/hN4EQNDm5lw4aHeemlH/POd57I\nmjXLa9Vg2zo4yMMbNnDY7t3sO/xwlq9ZE6XBOjQ0tP+jeCeub7ao+/k5E1q5C8gEkLhYFywFxjMu\n4xmPCUCSasq/BSRJKswEkLj8LWMqz3jGZTyrZQKQpJqyByBJCbAHIEkqzASQOGuscRnPuIxntUwA\nklRT9gAkKQH2ACRJhZkAEmeNNS7jGZfxrJYJQJJqyh6AJCXAHoAkqTATQOKsscZlPOMyntUyAUhS\nTdkDkKQE2AOQJBVmAkicNda4jGdcxrNaJgBJqil7AJKUAHsAkqTCTACJs8Yal/GMy3hWywQgSTVl\nD0CSEmAPQJJUmAkgcdZY4zKecRnPapkAJKmmpqsXrQZ6Gp77CnA8cBuwF9gCfLZhnyOA+4FjgV3Z\n+7wGfAj4U+At4BXg94D/m+S49gAkqQmt9ACm27kXGG547mTgAWAV8DwwCKwDns7tcw1wNCEx/DZw\nFvAJ4Fng14BXgc8BLwIbJjmuCUCSmtCuJvCRwELCxR9gM3Buwz5LgU3Z9qbcfD/h4g8wn8l/+1dE\n1ljjMp5xGc9qtZIAuoDR3HgXsHiSfXZOMv9S9u8q4BzgvhaOL0mKYF4LrxkFFuXGXcCOSfbpyrYX\nNcxfTUgAHwbenOogAwMD9PT0ANDd3U1fXx/9/f3AxG8Njqcf9/f3d9R6ZvvYeBrPThmPb4+MjNCq\nMj2A3yCUgR4CPgM8mdvnGsKF/ybgo4S6/x8RegXvA34H2H2I49oDkKQmzEQT+GomfpMftxFYQLgL\n6DBCD+D6bG4zsJJQ37+XcLfQHuDS7Fg/Ar6bPQfwNeDOSY5rAohkaGho/28OKs94xmU842klAUxX\nAvriIebOmuS587J/9wK/Ncn8wiKLkiTNPP8WkCQlwL8FJEkqzASQuPwdAyrPeMZlPKtlApCkmrIH\nIEkJsAcgSSrMBJA4a6xxGc+4jGe1TACSVFP2ACQpAfYAJEmFmQASZ401LuMZl/GslglAkmrKHoAk\nJcAegCSpMBNA4qyxxmU84zKe1TIBSFJN2QOQpATYA5AkFWYCSJw11riMZ1zGs1omAEmqKXsAkpQA\newCSpMJMAImzxhqX8YzLeFbLBCBJNWUPQJISYA9AklSYCSBx1ljjMp5xGc9qmQAkqabsAUhSAlrp\nAcybZn410NPw3FeA44HbgL3AFuCzDfscAdwPHAvsyt7ntdz8tcB7gUuaWawkKZ7pSkCPAjc1POYD\ndxAu3mcDS4C+htddCXwPWAbcB1yXmzsfuADwV/w2sMYal/GMy3hWq5UewJHAQuD5bLwZOLdhn6XA\npmx7U27+ZOAPgBvp3PKTJNVCKwmgCxjNjXcBiyfZZ2fD/FHA7cAVwL4WjqsW9Pf3V72EpBjPuIxn\ntabrAUxmFFiUG3cBOybZpyvbXpTNrwB+Dvga0A2cAHwS+LMW1iBJKqmVBPAG8CZwEqEMtAL4TMM+\n2wh1/icJNf+twIPZA+Ac4OMc4uI/MDBAT08PAN3d3fT19e3/bWG8buh4+nG+xtoJ65ntY+NpPDtl\nPL49MjJCq6arw1/NxG/y4zYCCwh3AR1G6AFcn81tBlYSGsX3Eu4W2gNcCrySe49zCKWgS6c4rreB\nRjI0NLT/xFF5xjMu4xlPK7eBdmoj1gQgSU3wbwFJkgozASQuXy9UecYzLuNZLROAJNWUPQBJSoA9\nAElSYSaAxFljjct4xmU8q2UCkKSasgcgSQmwByBJKswEkDhrrHEZz7iMZ7VMAJJUU/YAJCkB9gAk\nSYWZABJnjTUu4xmX8ayWCUCSasoegCQlwB6AJKkwE0DirLHGZTzjMp7VMgFIUk3ZA5CkBNgDkCQV\nZgJInDXWuIxnXMazWiYASaopewCSlAB7AJKkwkwAibPGGpfxjMt4VssEIEk1ZQ9AkhJgD0CSVNh0\nCWA1cGPD413AmcB3gEeBGyZ53RHA14GtwCDwjuz5k4GHgW8Dm4Bjyi1f07HGGpfxjMt4Vmu6BPAo\ncFPDYz5wB3AJcDawBOhreN2VwPeAZcB9wHXZ838JrAPOybbfU/q/QIf09NNPV72EpBjPuIxntVop\nAR0JLASez8abgXMb9llK+A2f7N9zgcOBY4GPAI8AZwGPt3B8NWHHjh1VLyEpxjMu41mtVhJAFzCa\nG+8CFk+yz86G+Z8BfolQAvoAofyzuoXjS5IiaCUBjAKLcuMuoDGNj2bPk+27A/gfQjL4dvb8Q8AZ\nLRxfTRgZGal6CUkxnnEZz87WO8lzJwP/CpxEuOVoEPjVhn2uITSMAT4K3J5tP0noGwB8kdArmMx2\nYMyHDx8+fBR+bKdJ090zejUTv8mP2wgsAG4DDiP0AK7P5jYDKwmN4nuB44E9wKXAK8BphGQwD/hP\nQglob7OLliRJkiRJUidbTWtfLJsD/Bfh9tFHgM+1Y7GzzFzgTuAxQowa+zcXAU9k85e3d2mzznSx\nvBr4NybOx3e3dXWz1xJCvBp5bjZvqlh29Ll5qKbyL2bjQQ7+YtnJhN6DprYKuCfbXgJ8Mzc3H/gh\n4Xbc+YQftuPaurrZ5VCxBPgq8CttXdHs90ngGcJFPs9zs3lTxRKaPDc74W8BFfli2enAzwP/REgQ\nHZXVOkT+y3ePc+AttqcS7hDYCbxF+KS1rK2rm10OFUsI5+O1wD8Dn2rjumaz7YTE2njjiedm86aK\nJTR5bnZCAijyxbL/JpR9Ppj9e397ljarNMZxHxP/f/NfzIPJY6wJh4olwN8CVxDOx7MJd77p0L7B\n5Hf8eW42b6pYQpPnZickgCJfLHuKiRLQNuCENqxrtmmM41zg7Wx7Z8PcIuCnbVrXbHSoWAJ8ifDF\nxrcIn0gtB7XOczOups7NTkgAbwBvMvHFshWEvyKadwPwiWz7l4EftW11s8c24IJs+0xCjXDcs8Ap\nhD+/sYDwEftf2rq62eVQsVwMfB84inC+fpDwC4pa47kZT9Pn5rw2LCrvI0z+xbKPA3/DxBfLnszm\nxr9Y9nlC2ecCwkefgTasdbZ5EFhOuHgBXEb4i61HA3cRvp29mZD07wZerGCNs8V0sfwU4Q6LPcA/\nMtEv0PTGsn89N8ubLJaem5IkSZIkSZIkSZIkSZIkSZIkSZKUkv8Hv4PiiLPgGGYAAAAASUVORK5C\nYII=\n",
            "text/plain": [
              "<matplotlib.figure.Figure at 0x7fa0717e6a58>"
            ]
          },
          "metadata": {},
          "output_type": "display_data"
        }
      ],
      "source": [
        "pt.xlim([a-0.5, b+0.5])\n",
        "pt.grid()\n",
        "pt.plot([a,b], [0,0], \"ob\")\n",
        "pt.plot([m1, m2], [0,0], \"or\")"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "collapsed": false
      },
      "outputs": [],
      "source": []
    }
  ],
  "metadata": {},
  "nbformat": 4,
  "nbformat_minor": 0
}