{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Secant Method" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as pt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here's a function: (Look here! No derivative.)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHv5JREFUeJzt3Xl8FeW9x/HPDxQFt5SXigpo6oJaq6ZYqXWpaRWLtC6t\n3uv1esXorcW1BquCS9HbDaxwRWptXbChV60LdtFbF9RLlOLKEgERERRcMEBRKCIaIM/94zkhISbh\n5MycM/Oc832/XvNK5mSY8309wO9MfjPzjDnnEBGR4tAl6QAiIhIfFXURkSKioi4iUkRU1EVEioiK\nuohIEVFRFxEpIlkVdTO728yWmdmcFq/1NLOnzGyBmU02s7L8xRQRkWxke6T+e2BQq9dGAE855/oB\nz2TWRUQkQZbtzUdmVg486pw7OLM+HzjWObfMzHYDap1zB+QrqIiIbFmUnnov59yyzPfLgF4x5BER\nkQhiOVHq/OG+5hsQEUnYVhH+7DIz2805V29muwPL29rIzFTsRURy4Jyzzv6ZKEfqjwDnZL4/B/hL\nexs651K/XH/99YlnKJacIWRUTuVM+5KrbC9p/CPwPLC/mb1rZucCo4GBZrYA+FZmPViLFy9OOkJW\nQsgZQkZQzrgpZzpk1X5xzp3Zzo+OjzGLiIhEpDtKM6qqqpKOkJUQcoaQEZQzbsqZDllfp57zG5i5\nfL+HiEixMTNcgU+UFpXa2tqkI2QlhJwhZATljJtypoOKuohIEVH7RUQkhdR+ERERFfUmofTZQsgZ\nQkZQzrgpZzqoqIuIFBH11EVEUkg9dRERUVFvEkqfLYScIWQE5YybcqaDirqISBFRT11EJIXUUxcR\nERX1JqH02ULIGUJGUM64KWc6qKiLiBQR9dRFRFJIPXUREVFRbxJKny2EnCFkBOWMm3Kmg4q6iEgR\nUU9dRCSF1FMXEREV9Sah9NlCyBlCRlDOuClnOqioi4gUEfXURURSSD11ERGJXtTN7Goze83M5pjZ\nfWa2TRzBCi2UPlsIOUPICMoZN+VMh0hF3czKgfOB/s65g4GuwL9FjyUiIrmI1FM3s57AC8ARwBrg\nz8AtzrmnW2yjnrqISCcl0lN3zn0IjAXeAZYCq1oWdBERKaytovxhM9sHqAbKgdXAQ2Z2lnPu3pbb\nVVVVUV5eDkBZWRkVFRVUVlYCzf2tpNebXktLnvbWx40bl8rxa7leV1dHdXV1avK0t9767z7pPO2t\nazxLYzxra2upqakB2FQvc+Kcy3kBzgDuarF+NvCbVtu4EEyZMiXpCFkJIWcIGZ1TzrgpZ7wytbPT\ndTlqT/1Q4F7gcOBToAZ42Tn3mxbbuCjvISJSipLqqb8K/AGYDszOvHxHlH2KiEjuIl+n7pz7lXPu\nIOfcwc65c5xz6+MIVmgt+4FpFkLOEDKCcsZNOdNBd5SKiBQRzf0iIpJCmvtFRERU1JuE0mcLIWcI\nGUE546ac6aCiLiJSRNRTFxFJIfXURURERb1JKH22EHKGkBGUM27KmQ4q6iIiRUQ9dRGRFFJPXURE\nVNSbhNJnCyFnCBlBOeOmnPGZNCn3P6uiLiKSIosXw9Chuf959dRFRFJk2DDYaisYMya3nrqKuohI\nSnz0EeyzD8yeDX376kRpJCH02SCMnCFkBOWMm3JGd8cd8N3vQp8+ue8j0oOnRUQkHg0NMH48PPZY\ntP2o/SIikgITJ8K998LkyX491+vUdaQuIpIw52DMGBg7Nvq+1FPPSHOfraUQcoaQEZQzbsqZuyef\nBDMYODD6vlTURUQSNmoUDB/uC3tU6qmLiCRo2jQ4+2xYsMBfn95Ec7+IiARo1Ci46qrNC3oUKuoZ\naeyztSWEnCFkBOWMm3J23quvwsyZUFUV3z5V1EVEEjJ6NFRXw7bbxrfPyD11MysD7gIOAhxwnnPu\nxRY/V09dRKSVhQvhiCPgrbdgxx0///Mkr1O/BXjMOXe6mW0FbBfDPkVEitqvfgUXXNB2QY8iUvvF\nzHYCjnHO3Q3gnNvgnFsdS7ICS1OfrSMh5AwhIyhn3JQze0uW+DnTq6vj33fUnvoXgRVm9nszm2lm\nd5pZjziCiYgUq1Gj4Ic/hJ13jn/fUdsvWwH9gUucc6+Y2ThgBDCy5UZVVVWUl5cDUFZWRkVFBZWV\nlUDzp6bWs1tvei0tedpbb5k1DXnaWq+srExVno7Wm6Qlj8Yz9/Vly+DBByt5443Nf15bW0tNTQ3A\npnqZi0gnSs1sN+AF59wXM+tHAyOcc99tsY1OlIqIZFx0EeywA9x4Y8fbJXLzkXOuHnjXzPplXjoe\neC3KPpPS+hM8rULIGUJGUM64KeeWvfce3H8/XHFF/t4jjqtfLgXuNbNuwCLg3Bj2KSJSdEaPhv/8\nT9hll/y9h+Z+EREpgPfeg0MOgddfh169trx9ru0XFXURkQIYOhS+8AV/tJ4NTegVkfqB8QkhIyhn\n3JSzfQsXwsMP+4m78k1FXUQkz264AS67DHr2zP97qf0iIpJHc+fCccf5o/Uddsj+z6n9IiKSQj/5\niW+7dKagR6GinqF+YHxCyAjKGTfl/LxXXvHLRRcV7C1V1EVE8sE5f4Q+ciR0716491VPXUQkDx57\nDH78Y5gzJ7dH1amnLiKSEhs3wvDhfn6XuJ49mi0V9Qz1A+MTQkZQzrgpZ7M//MHfaHTSSXl/q88p\n8GeIiEhx++QT30efNAms082T6NRTFxGJ0ahRMGsWPPhgtP1o7hcRkYTV18OXvwwvvgj77httXzpR\nGpH6gfEJISMoZ9yUE669Fs49N3pBj0I9dRGRGMyc6S9jnD8/2Rxqv4iIROQcfOMbMGQInH9+PPtU\n+0VEJCEPPQRr1sB55yWdREV9E/UD4xNCRlDOuJVqznXr4Mor4ZZboGvXWHedExV1EZEIRo2CAQPg\n2GOTTuKppy4ikqM334Svfx3q6qBPn3j3rZ66iEgBOQeXXurneIm7oEehop5Rqv3AfAghIyhn3Eot\n55/+BO++C9XVsewuNrpOXUSkkz7+GIYN8xN3bb110mk2p566iEgnXXUVLF0K99yTv/fItaeuI3UR\nkU6YNQtqavzDL9Iolp66mXU1s1lm9mgc+0tCqfUD8ymEjKCccSuFnBs2wA9+4B9+0atXfJniFNeJ\n0suAeYD6LCJStG65BcrKoKoq6STti9xTN7M+QA3wC+By59xJrX6unrqIBO+tt/xNRnFMq5uNJK9T\nvxm4EmiMYV8iIqnjHFxwgZ8OIMlpdbMRqaib2XeB5c65WUACD26KTyn0AwslhIygnHEr5px33QUr\nV8Lll8efJ25Rr345EjjZzAYD2wI7mtkfnHNDWm5UVVVFeXk5AGVlZVRUVFBZWQk0D3DS603Skqe9\n9bq6ulTlaWu9rq4uVXlCX9d4Jjue9fVwzTWVTJkC06blL19tbS01NTUAm+plLmK7Tt3MjgWuUE9d\nRIpFYyMMHOiXESMK+95pmftF1VtEisbvfgdr18IVVySdJHuxFXXn3LPOuZPj2l+hNf0alHYh5Awh\nIyhn3Iot56JFcP31MHEibBXQbZqa0EtEpJX16+Gss/yDpPffP+k0naO5X0REWhk5El5+2T9IuktC\nh76a+0VEJAZTp8Kdd8LMmckV9CgCjJwfxdYPTFIIGUE541YMOVetgrPPhjvugN13L1ymOKmoi4jQ\nfNfo4MFw0klb3j6t1FMXEcFfvnjbbX5ulx49kk6Te09dRV1ESt7MmfDtb8O0adCvX9JpvLTcfBSs\nYugHpkUIGUE54xZqztWr4V//FW69NT0FPQoVdREpWc75h16ccAKccUbSaeKh9ouIlKwxY+D+++Hv\nf4dtt006zeZ0nbqISCdMngxjx/qbjNJW0KNQ+yUj1H5gGoWQEZQzbiHlXLTIX4/+wAPQt2/SieKl\noi4iJWXdOjj1VD8VwDe+kXSa+KmnLiIlY+NGOP10+MIXYMIEsBQ/r009dRGRLRg+HD76yLdd0lzQ\no1D7JSOkfmDahZARlDNuac95++3wyCMwbFgt3bolnSZ/dKQuIkVv8mT/wIupU+H995NOk1/qqYtI\nUZs1y08B8PDDcMwxSafJnqYJEBFpZeFC+M534Le/DaugR6GinpH2fmCTEHKGkBGUM25py1lf74/Q\nr78eTjut+fW05YybirqIFJ1Vq2DQIKiqgqFDk05TWOqpi0hRWbPGT9A1YACMGxfupYuaT11ESt7a\ntf4I/ctf9g+8CLWgg06URhZKny2EnCFkBOWMW9I5P/nEP4auXz/4zW/aL+hJ58w3FXURCd7atXDK\nKdC7t39odJcSrmxqv4hI0Fav9pct9usHd94JXbsmnSgeibRfzKyvmU0xs9fMbK6Z/SjK/kREOmPl\nSjj+eKiogLvuKp6CHkXUX1LWA8OccwcBRwAXm9mB0WMVXih9thByhpARlDNuhc75/vtQWQnf/Cb8\n+tfZt1xCGc9cRSrqzrl651xd5vuPgdeBPeIIJiLSnnnz4Mgj4ayz4MYbw77KJW6x9dTNrBx4Fjgo\nU+CbXldPXURiM3WqnxN9zBj/9KJileh86ma2PTAJuKxlQW9SVVVFeXk5AGVlZVRUVFBZWQk0/yqk\nda1rXetbWh85spbx4+GhhyoZODD5PHGu19bWUlNTA7CpXubEORdpAbYGngSq2/m5C8GUKVOSjpCV\nEHKGkNE55YxbPnNu3Ojcddc5t9deztXVRdtXKOOZqZ2drsmRjtTNzIAJwDzn3Lgo+xIRacvHH8OQ\nIbBiBbz8Muy6a9KJ0i1ST93MjgaeA2YDTTu62jn3RIttXJT3EJHS9cYbvn8+YIC/7X+bbZJOVDia\n+0VEisoDD8All8AvfgHnn196V7ho7peImk5YpF0IOUPICMoZt7hyfvopXHwxXHutfwzdD38Yb0EP\nZTxzpaIuIqnx6qvw1a/C8uUwYwZ85StJJwqP2i8ikriNG2HsWLjpJv/17LNLr93SWqLXqYuI5Or1\n133PvEsXeOUViHKJtqj9skkofbYQcoaQEZQzbp3N2dAAP/uZfyD0mWdCbW1hCnoo45krHamXuE8/\nhXfegSVL/LJ0qb8eePly//Wf//SPB1uzxm+7fr3/z7hxoz+y6tLFz4y3zTbQowd07w7OQZ8+UFbm\nl113hV128V93393Ped27N+y0k37FLlXPPAM/+hHsvTfMmgV9+yadqHiop14iNmyA+fNh5kx/Mmr+\nfP9r79KlvsDutZdf9tgDevXyBXjnnX3h3WEHv3TvDltv7ZeuXX3xbmz0+/70U1i3zi9r1vg5rlev\nhg8/9B8OTR8UH3zgZ9dbutR/MJSX+/ctL4d9921e9t67tK5JLhVvvw1XXOH/HY4dC9/7nj7Y26Pr\n1GUzq1fDtGl+8qOpU/3RUO/e0L+/n3v6S1+CAw6AL37RF+mkMi5ZAosX+//sixbBwoXw5pvw7rv+\n6O2AA/xy0EH+uZMHHgjbbZdMXsndypV+NsUJE2DYMPjxj/1BgrRPRT2i2traTZPspFl7ORsb/SVg\nTzzhl9mz/aVhxxzjlwED/FF3khk7o6HBF/k33vDTrL72ml8WLPAfTocc4pdDD/UfUnvt1fkjvtD/\nztOmrZwffwzjxvnl9NPhJz/xf39JCmU8dfVLCdqwwR+FT5oEf/6z71+feCLccIMv5Ntum3TC3HXr\n5o/KDzwQTj21+fUNG/yR/OzZvo10551QV+cfOlxR4X8TOeww/3W//fQknKR89BHceqtfjjsOXnzR\nt9Uk/3SkHqDZs2HiRLj3Xn9C8rTT/NKvX9LJkrN8uW8xzZrlf2OZMQP+8Q9f6L/61eZl331L+6HE\n+bZkiS/kd98NJ58Mw4f79pl0ntovRW7tWl/Ef/c7X6yGDPFLKRfyLfnwQ39CbsYMf/3zjBn+tcMO\nay7yhx/uT9LqZF3unIMpU/wj5Z57zv+7HDYM9twz6WRhU1GPKK19trff9kc+EyfCUUfB0UfXcvnl\nlaluK6R1LMF/IE6f7pfHH6/l7bcraWjY/Gj+sMP8b0BpKfRpHc933oH/+R//b7NbNzjhhFp++tNK\ntt8+6WQdS+t4tqaeepGZOxdGj/YnPc87zxeh8nJ/g0aaC3ra7bwzDBrkl6OP9g8uXrrUH8VPnw53\n3OG/mvni3tSf79/fH3mmpdAnpb7en7+ZNMmfyzjjDF/YBwyAZ58l9QW9FOhIPWVmz4aRI/2Jpepq\nuPDCwl21Ip5z8N57vrjPmtXcwmlo8D36puXQQ32/uFu3pBPnj3MwZ44/uPjf//XfDx7sz+EMHhz2\nyfi0U/slcIsW+WL+zDMwYgQMHarreNOmvt4fndbV+WI/Z45vj+23Hxx8sL+WvmkpL4etAvw92Dl/\nGWnT/Q1PP+3vFB40yF9ZdfzxuimsUFTUI0qqz7Zqlb8E8Z57/G3Tw4b5uzfbE0I/MISMEE/Odev8\ndfRz5/rr6OfO9evLlvm7Yg84wBf9lnfK9u7duRZavsZzwwZ46y3/4dT028iMGf7mrqb7G775TZ8/\nyZxxCyWneuqB2bjR3103ciSccoq/ZX+XXZJOJZ3VvXtz772lTz7x19PPn+/vkp02DWpq/N2zK1b4\nwr7nns3z4PTu7adn2GUXvzRN0dCjR265nPN37C5f7pf6+s3n+FmwwBf03r393cWHHQYXXeTPHfTp\nE3VUJEk6Uk9AXR384Ae+Hzl+vP+PJKXjs8/8NAjvvOPnwWlamgrwihX+Kp3Vq/3R9I47Nk+W1r27\n7+E3TaYGzZOsNTT4OzjXrPFft9uueR6fXXf1HyJNc/zst59f1OJLL7VfArBuHfzXf/kbM268Eaqq\ndDWFdKyhwRf3psnSPvnEF/rGRr8454t8t25+Dp/tt/ftu+23D7OnL830jNKI8j3H8gsv+LlK3n7b\nX+Fy7rm5FfQQ5oIOISOEkbNbN3jttVr23BP2398/3u3ww+FrX4Ovfx2OPNJfW3/IIX5Khb59/XQR\nSRT0EMYTwsmZK32W59n69fDTn/o5Sm67Db7//aQTiUgxU/sljxYsgLPO8ie+JkzwD4gQEcmG2i8p\n8+CD/rb+qir4299U0EWkMFTUM+Lqs332GVx6KVx9NTz5JFx8cbwnQ0PoB4aQEZQzbsqZDpGLupkN\nMrP5ZvammQ2PI1Soli6FY4/1l6fNmKFLFUWk8CL11M2sK/AGcDzwPvAKcKZz7vUW25RET/2VV/xJ\n0Asv9EfpulRRRKJIqqc+AFjonFvsnFsP3A+cEnGfwbnvPvjOd/wUuddco4IuIsmJWtR7A++2WH8v\n81pwcumzOecvV7z2Wj8R1ykF+DgLoR8YQkZQzrgpZzpEvU49q75KVVUV5eXlAJSVlVFRUbFpQp2m\nAU56vUm22x99dCUXXgjPPVfL2LFw8MGFyVtXV5fX/cexXldXl6o8oa9rPEtjPGtra6mpqQHYVC9z\nEbWnfgRwg3NuUGb9aqDROXdji22Krqe+dq1/OMCGDf5hAXowgIjELame+nRgPzMrN7NuwBnAIxH3\nmWqrVvk5pXfeGR59VAVdRNIlUlF3zm0ALgGeBOYBD7S88iUkTb8GdeQf/4BvfcvPu/H73/sJlAot\nm5xJCyEjKGfclDMdIs/94px7HHg8hiyp9sEH/gj91FPh5z/XFS4ikk6a+yUL77/vnwBTVeUvWRQR\nyTc9+ShP6uvhuOP8Qy2uuirpNCIiHdPcLxlt9dlWrPAF/T/+Iz0FPYR+YAgZQTnjppzpoKLejg8/\n9D30738frrsu6TQiItlRT70Na9f6gn7UUXDTTTopKiKFp2eUxmT9en+7f69e/lmiKugikgQ9JCOi\n2tpaGhvhvPOga1f/+Lk0FvQQ+oEhZATljJtypoOufmnhqqv8g6EnT9aT2EUkTGq/ZPz2t3DLLfD8\n89CzZ9JpRKTUqacewRNPwLnnwt//Dvvsk3QaERH11HM2Zw4MGQLXXFMbREEPoR8YQkZQzrgpZzqU\ndOd4+XI46SQYNw722CPpNCIi0ZVs+2X9en+36LHHws9+lnQaEZHNqafeSZdcAkuWwF//Cl1Kvgkl\nImmjnnon3H03PP003HNPc0EPpc8WQs4QMoJyxk0506HkeuovvQQjRsBzz8FOOyWdRkQkXiXVflm5\nEvr3h/Hj/VQAIiJppZ76FjQ2+itdDjwQxoxJOo2ISMfUU9+Cm26Cjz6CUaPa/nkofbYQcoaQEZQz\nbsqZDiXRU586Ff77v2H69GQeFi0iUihF335ZuRIqKuD222Hw4MRiiIh0inrqbXAOTjsNysv9kbqI\nSCjUU2/DhAmwaBH88pdb3jaUPlsIOUPICMoZN+VMh6LtqS9YAFdfDc8+C9tum3QaEZHCKMr2S0MD\nHHmkf4rRRRcV9K1FRGJR8PaLmd1kZq+b2atm9iczS839mT//uX/G6IUXJp1ERKSwovTUJwMHOecO\nBRYAV8cTKZrp0/2VLnfd1blnjIbSZwshZwgZQTnjppzpkHNRd8495ZxrzKy+BPSJJ1LuPvsMzjkH\nbr4Zdt896TQiIoUXS0/dzB4F/uicu6+NnxWspz5ihD9B+vDDnTtKFxFJm1x76h1e/WJmTwG7tfGj\na5xzj2a2uRZoaKugF9KLL0JNDcyerYIuIqWrw6LunBvY0c/NrAoYDBzX0XZVVVWUl5cDUFZWRkVF\nBZWVlUBzfyvKekMDVFdXMn48zJtXy7x5nd9f02tx5Mnn+rhx42Ifv7jX6+rqqK6uTk2e9tZb/90n\nnae9dY1naYxnbW0tNTU1AJvqZU6cczktwCDgNWDnLWzn8u366507+WTnGhtz38eUKVPiipNXIeQM\nIaNzyhk35YxXpnZ2ujbn3FM3szeBbsCHmZdecM597qrwfPfU583zzxmdNQv6JH6qVkQkHiU590tj\nIxxzDJx1lm4yEpHiUpJzv9x+u/96wQXR99WyH5hmIeQMISMoZ9yUMx2Cnfvlgw9g5Eg/t0uXoD+a\nRETiE2z75d//3U+pm80MjCIiocnLdepp9X//B88/D3femXQSEZF0Ca5x0dAAF18M48bBdtvFt99Q\n+mwh5AwhIyhn3JQzHYIr6jffDHvvDaecknQSEZH0Caqn/s470L8/vPQS7LNPLLsUEUmlkrik8cor\nfetFBV1EpG3BFPWpU+GFF2D48PzsP5Q+Wwg5Q8gIyhk35UyHIIp6YyNUV8Po0dCjR37eo66uLj87\njlkIOUPICMoZN+VMhyCKek0NbLMNnHlm/t5j1apV+dt5jELIGUJGUM64KWc6pP469X/+E667Dv76\nV82TLiKyJak/Uv/lL2HgQDj88Py+z+LFi/P7BjEJIWcIGUE546ac6VCQSxrz+gYiIkUqlVPviohI\n4aS+/SIiItlTURcRKSKxF3Uzu8nMXjezV83sT2a2UzvbDTKz+Wb2ppnl6ZaidjP+i5m9ZmYbzax/\nB9stNrPZZjbLzF4uZMbM+2ebM7GxzLx/TzN7yswWmNlkMytrZ7tExjOb8TGz8Zmfv2pmXylUtlYZ\nOsxpZpVmtjozfrPM7LoEMt5tZsvMbE4H26RhLDvMmZKx7GtmUzL/x+ea2Y/a2a5z45nLg007WoCB\nQJfM96OB0W1s0xVYCJQDWwN1wIFxZ+kg4wFAP2AK0L+D7d4GehYqVy45kx7LTIZfAVdlvh/e1t95\nUuOZzfgAg4HHMt9/DXgxgb/rbHJWAo8UOlurDMcAXwHmtPPzxMcyy5xpGMvdgIrM99sDb8TxbzP2\nI3Xn3FPOucbM6ktAW4+DHgAsdM4tds6tB+4HCjbvonNuvnNuQZabJ3Z1fJY5Ex3LjJOBiZnvJwKn\ndrBtocczm/HZlN859xJQZma9Chsz67/HRO/WcM5NBT7qYJM0jGU2OSH5sax3ztVlvv8YeB3Yo9Vm\nnR7PfPfUzwMea+P13sC7Ldbfy7yWNg542symm9n5SYdpRxrGspdzblnm+2VAe//okhjPbManrW3a\nOhjJp2xyOuDIzK/hj5nZlwqWLntpGMtspGoszawc/5vFS61+1OnxzOmOUjN7Cv+rQ2vXOOcezWxz\nLdDgnLuvje3yfh1lNhmzcJRz7gMz2wV4yszmZ44AYhNDzoJck9pBzms3C+Oc6+DehLyPZxuyHZ/W\nR22FvtY3m/ebCfR1zn1iZicCf8G359Im6bHMRmrG0sy2ByYBl2WO2D+3Sav1Dsczp6LunBvY0c/N\nrArfCzqunU3eB/q2WO+L/wSKzZYyZrmPDzJfV5jZn/G/IsdahGLImfexhI5zZk5I7eacqzez3YHl\n7ewj7+PZhmzGp/U2fTKvFdIWczrn1rT4/nEzu83MejrnPixQxmykYSy3KC1jaWZbAw8D9zjn/tLG\nJp0ez3xc/TIIuBI4xTn3aTubTQf2M7NyM+sGnAE8EneWLLXZVzOzHma2Q+b77YATgHbP+BdAe/2/\nNIzlI8A5me/PwR/1bCbB8cxmfB4BhmSyHQGsatFOKpQt5jSzXmZ+BiQzG4C/eTBNBR3SMZZblIax\nzLz/BGCec25cO5t1fjzzcEb3TWAJMCuz3JZ5fQ/gby22OxF/tnchcHWBzzp/D9+nWgfUA4+3zgjs\njb8CoQ6YW+iM2eZMeiwz798TeBpYAEwGytI0nm2NDzAUGNpim1szP3+VDq6ISjIncHFm7OqA54Ej\nEsj4R2Ap0JD5t3leSseyw5wpGcujgcZMhqZ6eWLU8dQ0ASIiRUR3lIqIFBEVdRGRIqKiLiJSRFTU\nRUSKiIq6iEgRUVEXESkiKuoiIkVERV1EpIj8P8rJF+RaQ07MAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def f(x):\n", " return x**3 - x +1\n", "\n", "xmesh = np.linspace(-2, 2, 100)\n", "pt.ylim([-3, 10])\n", "pt.grid()\n", "pt.plot(xmesh, f(xmesh))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [], "source": [ "guesses = [2, 1.5]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Evaluate this cell many times in-place (using Ctrl-Enter)" ] }, { "cell_type": "code", "execution_count": 30, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-1.3247179572447458\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPoSlNkZWiwIooihUFGyASf4pSdRepIkVY\nu6tgAxuwlgVXWQGVXZCqFIEoVVipARQsQEioUlSKELBAABEC5Pz+OBkMkIRJptwp3/frdV9MZm7u\nPBzCkzPnPve5xlqLiIjEhkJeByAiIsGjpC4iEkOU1EVEYoiSuohIDFFSFxGJIUrqIiIxxK+kbowZ\nYYzZZYxZle25ssaYOcaYDcaY2caYMqELU0RE/OHvTH0k0Oik53oCc6y1lwDzsr4WEREPGX8vPjLG\nVAWmW2uvyvp6PdDAWrvLGFMRSLLW1ghVoCIicnqBrKlXsNbuynq8C6gQhHhERCQAQTlRat10X/0G\nREQ8ViSA791ljKlorU0zxpwH7M5pJ2OMkr2ISAFYa01+vyeQmfo0oFPW407AlNx2tNZG/Na7d2/P\nY4iVOKMhRsUZvG1s6ljK/ascHZ7o4HkssTCevq2g/C1pHA8sAS41xmwzxtwP9AMaGmM2AP+X9bWI\nxAlrLX0X9+X5ec8zv9N8qp1TzeuQBD+XX6y17XJ56fYgxiIiUeJo5lEe+/Qxvt7xNUu7LuX80ueT\nSKLXYQmBranHlISEBK9D8Es0xBkNMYLiLKgDGQdok9iGY5nHWNR5EaXPKA1EXpy5iZY4C8rvOvUC\nv4ExNtTvISLhsXP/TpqNb0atirUY3HQwRQsX9TqkmGWMwYb5RKmIxJE1u9dQZ3gdWtRowdDmQ5XQ\nI5SWX0TktBZ8v4A2iW3of0d/OtTs4HU4kgcldRHJ09jUsXT/rDsftfyI/7vw/7wOR05DSV1EcmSt\npe/nfRmyfAjzO83nyvJXeh2S+EFJXUROkVPJokQHJXUROcGBjAO0ntSaTJt5QsmiRAdVv4jIcTv3\n76TBqAZUKl2J6e2mK6FHISV1EQFUshgrtPwiIipZjCFK6iJxTiWLsUVJXSROqWQxNimpi8Sho5lH\nefTTR/lmxzcqWYwxSuoicWb/4f20SWyDxapkMQap+kUkjmQvWZzWdpoSegxSUheJE76SxZaXt1TJ\nYgzT8otIHFjw/QLaftyW/nf0576r7/M6HAkhJXWRGDc2dSxPzX6Kj+75iFsvvNXrcCTElNRFYpSv\nZHHo8qHM7zifK8pf4XVIEgZK6iIxKHvJ4pKuS1SyGEeU1EVijK9kUV0W45OqX0RiyI79O9RlMc4p\nqYvECJUsCmj5RSQm+Los/vvOf6tkMc4pqYtEuTGpY3h69tNMaDlBJYsSeFI3xjwP3AdkAquA+621\nhwM9rojkzVrLPxf/k/dXvK+SRTkuoDV1Y0xV4AGglrX2KqAw0DbwsEQkL0czj/LQjIf4eN3HLOm6\nRAldjgt0pr4POAKUMMYcA0oAPwYclYjkyleyCLDo/kWUKlbK44gkkgQ0U7fW/gr0B7YCO4C91tq5\nwQhMRE7l67JY+azKTGs3TQldThHQTN0YcxHQDagKpAOTjDHtrbVjs+/Xp0+f448TEhJISEgI5G1F\n4tKa3WtoMq4JD9V+iOdvfh5jjNchSRAlJSWRlJQU8HGMtbbg32xMG6ChtfZvWV93AG6y1j6WbR8b\nyHuICMz/fj5tE9vy9p1v0/7q9l6HI2FgjMFam+/f3IFefLQeuMkYU9y4acPtwNoAjyki2YxJHUPb\nxLZMaDlBCV1OK6DlF2ttijHmA2AZrqRxBTA0GIGJxDtfyeLQFUNZ0GmBKlzELwEtv/j1Blp+Ecm3\nI8eO8Oinj7J853Jm3DtDXRbjUEGXX3RFqUiE2X94P60TWwOwsPNCNeWSfFFDL5EIsmP/Dm4ZdQuV\nS1dWl0UpECV1kQjh67LY6vJWDG0+lCKF9EFa8k8/NSIRwNdlUSWLEigldRGPqcuiBJOSuohH1GVR\nQkFJXcQDvhtDL9+5nKVdl3Je6fO8DklihJK6SJj5ShYNhoWdF6oplwSVql9EwshXsljlrCrqsigh\noaQuEibZSxaHNBuikkUJCf1UiYSBr8uibgwtoaakLhJiH6Z8qJJFyZfExIJ/rxp6iYSItZbXF7/O\nsBXD+PTeT1WyKH754QeoXRt+/bVgDb2U1EVCIHuXxU/v/VQli+K37t2hSBF46y0ldZGIsP/wflpN\nakXhQoWZ0HKCKlzEb3v2wEUXQWoqVKnizZ2PRCQbX8niBWdfwNS2U5XQJV+GDoVmzaBy5YIfQzN1\nkSBZvXs1Tcc15eHaD9Pz5p66MbTkS0YGXHghzJwJNWvqJhkinvKVLA5oNIB7r7rX63AkCo0fD1dc\n4RJ6IJTURQL0YcqHPDPnGSa2mkhC1QSvw5EoZC289Rb07x/4sZTURQooe8nigk4LuLzc5V6HJFHq\ns8/AGGjYMPBjKamLFMCRY0d45NNHSE5LVpdFCVjfvtCjh0vsgVJSF8mn7CWL6rIogfriC9i2Ddq0\nCc7xVNIokg879u+g/sj6KlmUoOnbF557zl1wFAxK6iJ+Wr17NXWG16HtlW35b7P/qsuiBCwlBVas\ngM6dg3dM/VSK+EElixIK/fpBt25w5pnBO2bAFx8ZY8oAw4ArAAt0sdZ+me11XXwkUc1Xsjih5QSV\nLErQbNoEN90E330HZ5116uteXnw0EJhprW1pjCkClAzCMUU85ytZHJ48XCWLEnT/+hc8/HDOCT0Q\nAc3UjTFnA8nW2mp57KOZukSd7CWLM9rNUMmiBNWWLXDttbBhA5x7bs77eDVTvxD4yRgzEqgJLAee\ntNYeDPC4Ip7Zd3gfrSe1VsmihEzfvvDgg7kn9EAEmtSLALWAx6213xhjBgA9gV7Zd+rTp8/xxwkJ\nCSQkJAT4tiKh8eO+H2k6rik3Vb6Jd5u8qwoXCbqtW2HiRPj22xOfT0pKIikpKeDjB7r8UhFYaq29\nMOvrm4Ge1tpm2fbR8otEBV+XxUeve5Tn6j2nLosSEo8+CqVLwxtv5L2fJ8sv1to0Y8w2Y8wl1toN\nwO3AmkCOKeKFed/No93H7VSyKCG1fTt89NGps/RgCkZJY01cSWMxYDNwv7U2PdvrmqlLRPsg5QOe\nnfMsE1tOpEHVBl6HIzHs8ceheHF4883T71vQmbpukiFxy1rLa4teY8TKEcy8dyaXlbvM65Akhm3f\nDldfDevWQYUKp99fN8kQyYcjx47w8IyHWblrJUu6LFHJooTcq6+6ihd/EnoglNQl7uw7vI9Wk1pR\npFARlSxKWGzaBB9/7OrSQ00NvSSu/LjvR24ZeQsXlrlQXRYlbPr0gSefhLJlQ/9eWlOXuLFq1yqa\njW/GI9c9Qo96PVSyKGGxejXcdpubrZcu7f/3aU1dJA++ksWBjQbS7qp2XocjceTll12/9Pwk9EBo\npi4xTyWL4pVvvoG//hU2bnSljPmhmbrISbKXLCZ1SlLJooSVtW6G3qtX/hN6IJTUJSb5ShZTdqWw\ntOtSKpaq6HVIEmdmzYK0NOjSJbzvq6QuMcdXsli0UFGSOiepwkXC7tgx6NHD9XcJ1r1H/aWSRokp\nP+77kfoj61OtTDWmtJ2ihC6e+OADOOccaN48/O+tE6USM1btWkXTcU157PrH1GVRPHPwIFx6KSQm\nwo03Fvw4OlEqcU0lixIpBg6EOnUCS+iB0Exdot7olaN5bu5zKlkUz6WlwZVXwpdfwsUXB3YsdWmU\nuGOt5dVFrzJy5Uh1WZSI0LWrawXgT2vd09Hyi8SVI8eO8NCMh0jdlaqSRYkIK1bAzJmwfr23cSip\nS9TZd3gfLSe25IwiZ6hkUSKCta5h1yuvwNlnexuLSholqvhKFi8uezGT20xWQpeIMGkS7N8f/guN\ncqKkLlFj1a5V1B1Rl/ZXtee9Ju9RpJA+aIr3fv8dnn3WVb0ULux1NFp+kSjhK1kc1HgQba9s63U4\nIsf17Qs33AANIqTwStUvEvF8JYuTWk3ilgtu8TockeM2bnQ16StXQuXKwT22ql8k5mQvWVSXRYk0\n1sLf/+56vAQ7oQdCSV0ikkoWJdJ98gls2wbdunkdyYmU1CXiqGRRIt2BA9C9u2vcVbSo19GcSNUv\nElF8JYsXnXORShYlYr3yCtxyCyQkeB3JqTRTl4jh67L4+A2P82zdZ9VlUSJScjKMGgWrVnkdSc6C\nktSNMYWBZcB2a60HHYQl2s39bi73fnyvShYloh09Cn/7m7v5RYUKXkeTs2DN1J8E1gJhul+2xBJf\nyWJi60SVLEpEGzgQypSBzp29jiR3ASd1Y0xloAnwOvBUwBFJ3LDW8srCVxiVMkolixLxvvvOXWj0\n5ZcQySuDwZipvw08C5wVhGNJnDhy7AgPzniQVbtWqWRRIp618PDDrh1AoH3SQy2g6hdjTDNgt7U2\nGYjg310SSfYd3kfTcU35+eDPLOy8UAldIt6wYfDLL/BUFKxFBDpTrwvcZYxpApwJnGWM+cBa2zH7\nTn369Dn+OCEhgYRIrAOSsNi+bztNxzWlXpV6DGo8SE25JOJt2QIvvAALFoS2Jj0pKYmkpKSAjxO0\n3i/GmAbAMydXv6j3i/ioZFGiTWYmNGzotp49w/vekdL7RdlbcqSSRYlG//0v/PYbPPOM15H4T10a\nJeRGrRxFj7k91GVRosrmzXDTTfD553DppeF//0iZqYscp5JFiVZHjkD79vDii94k9EAoqUtIqGRR\notmrr7qLjJ54wutI8k9JXYIue5fFhZ0XUrJYSa9DEvHb4sXw/vuwYgUUisKWh1EYskSy7fu2n3Bj\naCV0iSZ790KHDjB0KJx3ntfRFIxOlErQpO5Kpdm4ZipZlKhkLbRrB2XLwuDBXkejE6XisTmb59D+\nk/a80/gd2lzZxutwRPJtyBBYu9b1dolmmqlLwHwli4mtEql/QX2vwxHJtxUr4M474Ysv4JJLvI7G\n0Uxdws5Xsjg6ZTQLOy+kxrk1vA5JJN/S06F1a3j33chJ6IHQTF0KxFeyuHr3ama0m0GFUhF6xwCR\nPFjrEnq5cpGxjp6dZuoSNifcGLpTkipcJGr17w/ffw8ffuh1JMGjkkbJF5UsSqyYPdsl9cmT4cwz\nvY4meJTUxW8paSnUGV6H9le1570m76ltrkStzZtdPfqECVClitfRBJf+V4pfVLIoseLAAfjLX6BX\nL7glBvvL6USpnNbI5JH0nNdTXRYl6h07Bi1bwjnnwPDhkX2vUZ0olaCz1vKPhf/gg5QPVLIoMaFH\nD9izxy27RHJCD4SSuuTIV7K4ZvcalnZdqpJFiXpDhsC0abB0KRQr5nU0oaOkLqfwlSyeWeRMFnRa\noAoXiXqzZ0Pv3q4D45/+5HU0oaXqFzmBr2SxetnqKlmUmJCcDPfdB5MmQfXqXkcTekrqclzqrlTq\nDq9Lh6s78G6TdylcqLDXIYkEZNMmaNoU/vMfqB8nbYm0/CIAzN48m/s+uU8lixIz0tJck67eveGe\ne7yOJnyU1OV4yeLHrT9Wl0WJCXv3QqNG0LkzPPSQ19GEl+rU45i1lj5Jffgw9UNmtp+pkkWJCfv3\nwx13wA03wIAB0Vu6qDp1yZeMYxk8OP1B1vykkkWJHb/9Bk2awDXXRHdCD4SSehxKP5ROy0ktKV6k\nuLosSsw4eBCaN3c90d97Lz4TOqj6Je74ShYvKXuJShYlZvz2G9x9N1Sq5G4aXSiOM1sc/9Xjj6/L\nokoWJZakp7sqlypVYNQoKBznP9YBJXVjTBVjzAJjzBpjzGpjzBPBCkyCa/bm2TT8sCH97+jPs/We\nxcTrZ1OJKb/8Arff7tbQhw1TQocAq1+MMRWBitbalcaYUsBy4C/W2nXZ9lH1i8dGJo/k+XnPM6nV\nJJUsSsz48UdXtti4MbzxRuytoXtS/WKtTQPSsh4fMMasA84H1uX5jRIW2bssJnVOUsmixIy1a10y\nf+QR13kx1hJ6IIJW/WKMqQpcC3wVrGNKwflKFtf+tFYlixJTFi92PdHfesvdvUhOFJSknrX0kgg8\naa09cPLrffr0Of44ISGBhISEYLyt5CJ7yaK6LEosmTAB/v53GDsWGjb0OprgSkpKIikpKeDjBHxF\nqTGmKDADmGWtHZDD61pTD6Nt6dtoMq4JDS5owMBGA1XhIjEhM9P1cPnwQ5g6FWrW9Dqi0PNkTd24\nEorhwNqcErqEV0paCs3GN6Pbjd14qs5TqnCRmHDgAHTsCD/9BF9/DeXLex1RZAu0Tr0ecB9wqzEm\nOWtrFIS4JJ+ylyw+XfdpJXSJCd9+C3XquHuKzp2rhO4PNfSKASOSR/D8vOdJbJWokkWJGRMmwOOP\nw+uvwwMPxF+Fixp6xSFfl8Uxq8awqPMiLj33Uq9DEgnYoUPw9NPw2WfuNnTXXut1RNFFST1KZRzL\n4IHpD7Dup3Us6bJEJYsSE1JSoH17uOwyWL4czj7b64iij3q/RKH0Q+k0GduEPb/vYUGnBUroEvWO\nHYN//ctd8v/cczBxohJ6QWmmHmVUsiixZt06t2ZeqBB88w1Urep1RNFNM/UokpKWQt0RdelcszPv\nNH5HCV2iWkYGvPqquyF0u3aQlKSEHgyaqUcJ342h323yLq2vaB204x46BFu3wpYtbtuxw9UD797t\n/ty3z90ebP9+t++RI+4/47FjbmZVqJDrjHfGGVCiBBQvDqVKuY/OZcq4rXx5KFfO/Xneea7ndaVK\nbp94q2gQZ948eOIJqFYNkpNd21wJDpU0RoERySN4Yd4LJLZO5OY/31ygYxw9CuvXw4oV7mTU+vXu\nY++OHS7BXnCB284/HypUcAn43HNd4i1d2m3Fi0PRom4rXBisdVf6HT3qEv7vv7tt/37X4zo9HX79\n1f1y8P2i2LnTddfbscP9Yqha1b1v1apw8cV/bNWquV8UElu+/x6eecb9HPbvD3/9q36x56agJY1K\n6hEse8nizHtn5qtkMT0dvvjCNT9avNjNhipVglq1XO/pyy+HGjXgwgtdkvZCerr7dPDDD+4/++bN\nsGkTbNwI27a52VuNGm674gq48kpXFVFSrWyizi+/uPa4w4dD9+6uZLF4ca+jimxK6jEme8ni9HbT\nT1vhkpnpSsD+9z+3pabCdde59cr69d2d1aOpmiAjwyX5b791bVbXrHHbhg3ul9PVV7utZk33S+qC\nCzTji0QHDrgbQA8Y4Dorvvyy+/eT01NSjyHph9JpMbEFpYqVYlyLcbl2WTx61M3CExNh8mS3ft24\nsbtxQP36cOaZYQ48DI4edTP51FS3jJSSAitXupsOX3ON+yRSu7b7s3p13QnHK3v2wLvvuu222+CV\nV9yymvhPST1G+FOymJoKo0e79qOVK8M997jtkks8CDhC7N7tlpiSk90nluXL4eefXaK/7ro/tosv\nju+bEofali0ukY8YAXfd5W5gUUP3ZikQJfUYsDJtJc3HN+fJG5/k6TonNuX67TeXxP/7X5esOnZ0\nWzwn8tP59Vd3Qm75clf/vHy5e6527T+S/PXXu5O0WropOGthwQJ45x1YtMj9XHbvDn/+s9eRRTcl\n9Sj32abP6DC5A+81eY9WV7Q6/vz337uZz+jRUK8ePPaY+zirZYWC+flnWLbMbd9847aMjBNn87Vr\nu09ASvR527rV9TcfPRqKFXPNt+67z5W0SuCU1KNYTiWLq1dDv37upGeXLvDoo7owI1R27HCzeF+i\nX7bMJfTatf9Yn69Vy8084z3Rp6W58zeJie5cRps20KmTOxEf72MTbErqUchaS++k3oxdNfZ4yWJq\nKvTqBV9+Cd26uRvrRlPVSiywFrZvd8k9OfmPJZyMDLdG79tq1nTrxcWKeR1x6FgLq1a5ycWMGe5x\nkybuHE6TJrF5Mj5SKKlHmewlizPuncH+tPL06uWutOvZEx56SHW8kSYtzc1OV650yX7VKrc8Vr06\nXHWVq6X3bVWrQpEovF7bWldG6ru+Ye5cd6Vwo0ausur223VRWLgoqUeR9EPp3DPxHkoWK8ng28bx\n5uslGTPGXTbdvbu7elOiw++/uzr61atdHf3q1e7rXbvcVbE1arikn/1K2UqVIuOcyNGj8N137peT\n79PI8uXu4i7f9Q233uril/BTUo8Sx0sW/5zAldsH0Kd3Ye6+G157zfVHkdhw8KCrp1+/3l0l67tS\n9ocfXMuESpXcGr2vD06lSq49Q7lybvO1aChRomBr1da6K3Z373ZbWtqJPX42bHAJvVIld3Vx9nMH\nlSsHfTikAJTUo4CvZLF1le4k9etO8TMNgwa5/0gSPw4fdm0Qtm51fXB8my8B//STq9JJT3ez6bPO\n+qNZWvHibg3f10wN/miylpHhruDcv9/9WbLkH318ypd3v0R8PX6qV3eblvgil5J6hPOVLNbd8x5L\nhrXijTegc2dVDEjeMjJccvc1Szt40CX6zEy3WeuSfLFirodPqVJu+a5Uqehc05c/6B6lEWz4iuE8\n99mLFJ8+mTOq1CM1FSpW9DoqiQbFimlZTvJHST2ErLW8NK83gz8fS+HxixjU9xJatPA6KhGJZUrq\nIZJxLIM2Yx5gTvJ6bvxuKWOSynPeeV5HJSKxTkk9BPYe2kv9d+5hw+rS9Lt+Ad3eLqG1cxEJCyX1\nINv00zZuHNSYo5tu5fNnBnB97QgoSBaRuBFwE1JjTCNjzHpjzEZjTI9gBBWtZq9ayRX963L+7vv5\nYfAgJXQRCbuAShqNMYWBb4HbgR+Bb4B21tp12faJyZLGN954gyHTppFZtCiFjhzhmssuY+q5U2ld\najDjXmyl5RYRCYhXJY03AJustT9kBfERcDewLq9vinZvvPEG/ebNY+/rrx9/bver/6DR760YP7ZV\nHt8pIhJagS6/VAK2Zft6e9ZzMW3ItGnsfeGFE5777eXerP8hxaOIREScQGfqfq2rmFhci7j11lOe\n+o4Y/buKSNQINKn/CFTJ9nUV3Gz9BLG2pl61bl22/POfpzxf7cUX2fzFFx5EJCKxpqATxECXX5YB\n1Y0xVY0xxYA2wLQAjxnRVm/bxi/Fd1Dq1VdPeL7M66/z4F13eRSViIgTcEMvY0xjYABQGBhure17\n0usxU/2yYN1K7hzdnJvsUzQ5J4P3p0/jWNGiFD5yhAfvuosePeK6olNEgkhdGkNs3Nf/o+PkjtxV\neDAfv9pSJYsiElIFTeoBX3wUD96cN4yOn3Sma8nJfPKaErqIRC61CciDtZbu019m8KLxdD93EW/2\nvMTrkERE8qSknouMYxm0n9iVT5du5KlyS+nXs7zXIYmInJaSeg72HtpL87EtSP36LJ4oP59+L5fw\nOiQREb9oTf0kW9O3UnfYzWz6/Cq6lv6Yvq8ooYtI9FBSz2Zl2krqDq+LXdGVRgyk/5uFdVJURKKK\nll+y/G+TK1mssXkwZ//Ukvcn66bQIhJ9NFMHhq0Yxv1T7+e2n6aQubolEyboTuwiEp3iOnVZa+m1\noBfjV4/n4WKLGD+tOkuWQAkto4tIlIrbpJ5xLIOu07qy8ZeNvHbhEro/WJ7PP4eyZb2OTESk4OIy\nqe89tJcWE1pw9plnM6j2fJrdWYLJk+Gii7yOTEQkMHG3pr41fSs3j7iZq8pfxeAGibT+awkGDIB6\n9byOTEQkcHGV1JN3JlN3eF3+VutvvHX7QNq0LkyHDnDvvV5HJiISHHHTpdFXsvifpv/hnsvv4fHH\nYcsWmDoVCsXVrzYRiQZe3Xg6Kry//H1eXvAyU9pOoW6VuowYAXPnwldfKaGLSGyJ6Zm6tZaXF7zM\nR6s/Ylb7WVT/U3W++gqaN4dFi6BGDU/CEhE5Lc3UT5K9ZHFp16WUK1mOX36B1q3h/feV0EUkNsVk\nUveVLJY5swzzO82nRNESZGZCx47QqhXcfbfXEYqIhEbMrShnL1mc1GoSJYq6y0PffBP27IG+fU9z\nABGRKBZTM/Xknck0H9+cZ+o+Q7ebuh1/fvFi+Pe/YdkyKFrUwwBFREIsZpL6rI2z6DSl0/GSRZ9f\nfnF16CNHQpUqHgYoIhIGMVH9MmzFMF6a/xKftPmEulXqHn/eWrjnHqha1c3URUSiRVxWv1hreWn+\nS0xYM4HF9y+m+p+qn/D68OGweTOMG+dRgCIiYRa1M/WMYxl0mdqFTb9uYnq76ZQrWe6E1zdscP1c\nFi6Eyy8P+tuLiIRUQWfqUVn9svfQXhqNacTBIweZ32n+KQk9I8Oto//jH0roIhJfCpzUjTFvGmPW\nGWNSjDGfGGPODmZgudmydwv1RtQ7pWQxu9degwoV4JFHwhGRiEjkCGSmPhu4wlpbE9gAPB+ckHKX\nvDOZeiPq8UCtBxjYeCCFCxU+ZZ9ly2DIEBg2TPcYFZH4U+Ckbq2dY63NzPryK6BycELK2ayNs7hz\nzJ0MbDTwhBr07A4fhk6d4O234bzzQhmNiEhkClb1SxdgfJCOdYr3l79Pr6ReTG07lTpV6uS6X+/e\ncOml0K5dqCIREYlseSZ1Y8wcoGIOL71grZ2etc+LQIa1NuiFg74uixPWTGBR50WnlCxm9+WXMGoU\npKZq2UVE4leeSd1a2zCv140xnYEmwG157denT5/jjxMSEkhISDhtYL4ui5t+3cSSLktOqXDJ7vBh\n6NIFBg2C8uVPe2gRkYiTlJREUlJSwMcpcJ26MaYR0B9oYK39OY/98l2nnr3L4tgWYyletHie+/fp\nA8nJMGWKZukiEhsKWqceSFLfCBQDfs16aqm19tEc9stXUt+ydwtNxjWhYbWG9L+jf44VLtmtXQsN\nGrikXjmkp2pFRMIn7End7zfIR1JfsXMFd42/65Qui7nJzIT69aF9e3j0lF8nIiLRK+p7v8zcOJNO\nUzoxpNkQWlzWwq/vGTLE/fnwwyEMTEQkikTETH3o8qH0WtCLyW0m51mymN3OnXD11ertIiKxKSpn\n6r4uixPXTsyxy2Jenn4aHnhACV1EJDvPkvrho4fpMq0L3+357rQliyebPx+WLHE3kBYRkT940qVx\nz+97aDS2Eb8f+Z35HU/tspiXjAx47DEYMABKlgxhkCIiUSjsSX3L3i3cPPJmalaoyaRWk05bg36y\nt9+GatVnBdnOAAAFFUlEQVTg7rtDFKCISBQL64nS/JYsnmzrVqhVC776Ci66KNiRiohEjog/UTpr\n4yw6TumYr5LFkz37rFt6UUIXEclZWGbqQ5YNyXfJ4skWL3YXGa1fDyVOvS+GiEhMiejb2b255E0W\n37+4wAk9MxO6dYN+/UKX0IPRSCccoiHOaIgRFGewKc7IEJakvqTLknzVoJ9s1Cg444zQ9kmPln/o\naIgzGmIExRlsijMyhGVNPT8liyfbtw9eegmmTlUHRhGR0/GkTj0//vlPaNgQrr/e60hERCJfWE6U\nhvQNRERiVES23hURkfCJ+OUXERHxn5K6iEgMCXpSN8a8aYxZZ4xJMcZ8Yow5O5f9Ghlj1htjNhpj\negQ7jtPE2MoYs8YYc8wYUyuP/X4wxqQaY5KNMV+HM8as9/c3Ts/GMuv9yxpj5hhjNhhjZhtjyuSy\nnyfj6c/4GGMGZb2eYoy5NlyxnRRDnnEaYxKMMelZ45dsjHnJgxhHGGN2GWNW5bFPJIxlnnFGyFhW\nMcYsyPo/vtoY80Qu++VvPK21Qd2AhkChrMf9gH457FMY2ARUBYoCK4HLgh1LHjHWAC4BFgC18tjv\ne6BsuOIqSJxej2VWDP8Cnst63COnf3OvxtOf8QGaADOzHt8IfOnBv7U/cSYA08Id20kx1AeuBVbl\n8rrnY+lnnJEwlhWBa7IelwK+DcbPZtBn6tbaOdbazKwvvwJyuh30DcAma+0P1tojwEdA2PouWmvX\nW2s3+Lm7Z9Xxfsbp6VhmuQsYnfV4NPCXPPYN93j6Mz7H47fWfgWUMcZUCG+Yfv87enq1hrV2MbAn\nj10iYSz9iRO8H8s0a+3KrMcHgHXA+Sftlu/xDPWaehdgZg7PVwK2Zft6e9ZzkcYCc40xy4wxD3gd\nTC4iYSwrWGt3ZT3eBeT2Q+fFePozPjntk9NkJJT8idMCdbM+hs80xkTifb8iYSz9EVFjaYypivtk\n8dVJL+V7PAt0RakxZg7uo8PJXrDWTs/a50Ugw1o7Lof9Ql5H6U+Mfqhnrd1pjCkHzDHGrM+aAQRN\nEOIMS01qHnG+eEIw1to8rk0I+XjmwN/xOXnWFu5aX3/ebwVQxVp70BjTGJiCW56LNF6PpT8iZiyN\nMaWARODJrBn7Kbuc9HWe41mgpG6tbZjX68aYzri1oNty2eVHoEq2r6vgfgMFzeli9PMYO7P+/MkY\nMxn3ETmoSSgIcYZ8LCHvOLNOSFW01qYZY84DdudyjJCPZw78GZ+T96mc9Vw4nTZOa+3+bI9nGWMG\nG2PKWmt/DVOM/oiEsTytSBlLY0xR4GNgjLV2Sg675Hs8Q1H90gh4FrjbWnsol92WAdWNMVWNMcWA\nNsC0YMfipxzX1YwxJYwxpbMelwTuAHI94x8Gua3/RcJYTgM6ZT3uhJv1nMDD8fRnfKYBHbNiuwnY\nm205KVxOG6cxpoIxrgOSMeYG3MWDkZTQITLG8rQiYSyz3n84sNZaOyCX3fI/niE4o7sR2AIkZ22D\ns54/H/g0236NcWd7NwHPh/ms819x61S/A2nArJNjBKrhKhBWAqvDHaO/cXo9llnvXxaYC2wAZgNl\nImk8cxof4CHgoWz7vJv1egp5VER5GSfwWNbYrQSWADd5EON4YAeQkfWz2SVCxzLPOCNkLG8GMrNi\n8OXLxoGOp9oEiIjEEF1RKiISQ5TURURiiJK6iEgMUVIXEYkhSuoiIjFESV1EJIYoqYuIxBAldRGR\nGPL/piHkgymROi8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "# grab last two guesses\n", "x = guesses[-1]\n", "xbefore = guesses[-2]\n", "\n", "slope = (f(x)-f(xbefore))/(x-xbefore)\n", "\n", "# plot approximate function\n", "pt.plot(xmesh, f(xmesh))\n", "pt.plot(xmesh, f(x) + slope*(xmesh-x))\n", "pt.plot(x, f(x), \"o\")\n", "pt.plot(xbefore, f(xbefore), \"o\")\n", "pt.ylim([-3, 10])\n", "pt.axhline(0, color=\"black\")\n", "\n", "# Compute approximate root\n", "xnew = x - f(x) / slope\n", "guesses.append(xnew)\n", "print(xnew)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] } ], "metadata": {}, "nbformat": 4, "nbformat_minor": 0 }