{ "cells": [ { "cell_type": "markdown", "metadata": { "toc": "true" }, "source": [ "# Table of Contents\n", "

1  課題解答例
1.1  Gaussian(正規分布)へのフィット
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# 課題解答例\n", "\n", "## Gaussian(正規分布)へのフィット\n", "\n", "正規分布で知られる,ガウス関数\n", "$$\n", "f(x)= \\frac{1}{\\sqrt{2\\pi\\sigma}}\n", "\\exp \\left(\\frac{- (x-\\mu)^2}{2\\sigma^2} \\right)\n", "$$\n", "フィットをやってみましょう.\n", "\n", "例えば,平均値($\\mu$)が60点,偏差値($\\sigma$)が15点として,ピークの人数が20人としましょう." ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPXZ//H3PVkhQFgSMKxhiSEJu2GVIqICgksrWsHdqlTFrZtibftof/XRat21RSu2LtQVrRubUhGtoAl7ICyBsIQtYUuAELLdvz9m8IkxkZCZyZnlfl3XXMw5c+bM56vhw+Q7Z84RVcUYY0z4cDkdwBhjTNOy4jfGmDBjxW+MMWHGit8YY8KMFb8xxoQZK35jjAkzVvzGGBNmrPiNMSbMWPEbY0yYiXQ6QF0SEhI0OTnZ6RjGGBM0li1btk9VExuybUAWf3JyMtnZ2U7HMMaYoCEi2xq6rU31GGNMmLHiN8aYMGPFb4wxYcaK3xhjwowVvzHGhJmTFr+IdBGRz0RknYisFZE7PevbisgnIrLJ82ebep4/XkQ2iEieiEz39QCMMcacmoa8468EfqWq6cAwYJqIpAPTgYWqmgIs9Cx/h4hEAM8B5wPpwBTPc40xxjjkpMWvqrtVdbnn/mEgF+gEXAy87NnsZeDHdTx9CJCnqltUtRx4w/M8Y8xJHNx8gK/ueIPPL32GvA/WodV2mVTjG6f0BS4RSQYGAl8DHVR1t+ehPUCHOp7SCdhRY7kAGFrPvqcCUwG6du16KrGMCSlLfvk2Lf7+OOlHvmEE1e6Vs6Egoht5Ay5l6IL/R7O2zZwNaYJagz/cFZEWwGzgLlUtqfmYuq/Y7tXbEVV9QVUzVTUzMbFB3zo2JqRotbJo1B8Y/sRPia04zBejfkfO35dQ8EU+i698np2JAxi97DHyuo2hKGev03FNEGvQO34RicJd+rNU9V3P6r0ikqSqu0UkCSis46k7gS41ljt71hljaig7VMby/tczevsbfJFyPUOXzyClRfS3j3ce6f6FeOnd79Lv0as4MGAom97+iJSf9HEutAlaDTmqR4CZQK6qPl7joQ+Aaz33rwXer+PpWUCKiHQXkWhgsud5xhiP6spq1px+CSO2v8Gi8Q8zcv1MomuUfk3DHrmErS8vJlLLaX/JmWxbmNfEaU0oaMhUz5nA1cAYEVnpuU0AHgbOE5FNwLmeZUSko4jMAVDVSuA2YD7uD4XfUtW1fhiHMUFr8U+eYHDRXBZP/iuj596DuOQHt0+/JpPKz7+iWiI4etEUyo+UN1FSEyrEPT0fWDIzM9XOzmnCQe5ry+h59XCWd7yQoTveOWnp17T07ncZ9ugkFg29h9FLH/ZjShMMRGSZqmY2ZFv75q4xDjmy5wixP5vCfld7Uhf//ZRKH9zTPot7T2XU14+w/NGFfkppQpEVvzEOWX72r+hWkcfex16jTc+2jdrHGZ8/Tn50Kh2nX03xtkM+TmhClRW/MQ7Y9F4OI9f/ncWD7mLAXaMbvZ+49nGUz3yN06p3s+KKR30X0IQ0K35jHHDg1t9xmJb0e+M+r/eVdtUZfNV1MoO/epLC1Xt8kM6EOit+Y5pYzotLGbrnfVac8xvaprTzyT47vvhHYjhO7lUP+mR/JrRZ8RvThLRaqfjNvRRJezJfu8tn+00+L4Wv0m5k+Jrn2bE432f7NaHJit+YJrT8z58w8NAi1l3yO1qc1sKn+0555fdUEcG26//Hp/s1oceK35gmotVK7IO/oyCiG8Nemurz/SdlduLrwbczYstrbP4o1+f7N6HDit+YJpLz9yVkHM1iy2X3ENMqxi+v0efl31BONDunP+OX/ZvQYMVvTBM5/L/PUEw8g5642m+vkZCWSFavKQxa+wrF24v99jomuFnxG9ME9izfxeDt77Bi0M98PrdfW8L9t9OCo6y88x9+fR0TvKz4jWkC63/xPBFU0f3RaX5/rbQrB7G65Qi6ffwc1ZXVfn89E3ys+I3xs/Ij5aR/+TzZ7SfQbUzPJnnNI9fdTnJFHssenNckr2eCixW/MX6WdffbtK/ei+uO25vsNQc/PIk9riR41j7kNd9nxW+Mn7V69Vnyo05n0D3nNdlrRjWPYv1ZNzN43zzy529sstc1waEhV+B6SUQKRSSnxro3a1yUZauIrKznuVtFZI1nOzvBvgk7eR+so++RpWybcDOuyKZ9n5X+5FQqiWDbH19u0tc1ga8hP4n/BMbXXKGql6vqAFUdgPtavO/W9USPsz3bNugCAcaEkoKHXqWSCDL+dEWTv3b7fqexInEsvb5+zT7kNd9x0uJX1cXAgboe81yP96fA6z7OZUzQq66sJiVrFisSx5HYp4MjGSouv5rOVdtZ/exiR17fBCZvf/f8EbBXVTfV87gCn4rIMhHx/XfUjQlgq57+nE5VO6iY7L8vbJ3MgP+5mBJacvivrzqWwQQeb4t/Cj/8bn+kZzrofGCaiIyqb0MRmSoi2SKSXVRU5GUsY5x3ZMarlNCSgfdf7FiG5gnNWdVrEv03vc2xA8ccy2ECS6OLX0QigUuAN+vbRlV3ev4sBN4DhvzAti+oaqaqZiYmJjY2ljEBoXRfKQM2vc2qlEtp1raZo1la3Hw1rTjMivvfdzSHCRzevOM/F1ivqgV1PSgicSLS8sR9YCyQU9e2xoSalQ+8T0uO0PLWa5yOQv87R7MzogtRb9h0j3FryOGcrwNLgFQRKRCRGzwPTabWNI+IdBSROZ7FDsCXIrIK+Ab4WFXta4QmLES9+SoFEV3pd1u9s5tNxhXpYtPgKxlYNJ+inL1OxzEBoCFH9UxR1SRVjVLVzqo607P+OlWdUWvbXao6wXN/i6r299wyVNWuCWfCwr7cIgYWLSBv8BVNfux+fTrfezWRVLHu/recjmICQGD8VBoTQnIffp9Iqki663Kno3yr10XpbIrJIP7T2U5HMQHAit8YH4udM5ttkT04/bL+Tkf5jp1DJ9G3+AuK1hY6HcU4zIrfGB86lH+Q/vsWkj9wEuISp+N8R8fbJxFBNbkP/dvpKMZhVvzG+FDOQx8STQUJP5/kdJTvSbmkL1ujehE39x2noxiHWfEb40NRH85mZ0QXMq6v9ysrjhGXsPWMS+l34DMObq7zLCwmTFjxG+Mjh3cdpv+e+eT1vSTgpnlOaH/LJKKoZO1DHzgdxTjIit8YH1nz5znEcpzWNwTeNM8JaVedQUFEN6I/tKN7wpkVvzE+IrPfYa/rNPpMHeF0lHqJS8gbMIn+hQsoKShxOo5xiBW/MT5w7MAx+u2cw4a0nxARHeF0nB/U7qZJxFDOmoc+cjqKcYgVvzE+sObJhcRRStxVP3E6ykll3DCMva7TkA9tnj9cWfEb4wNlb3/IYVrQ97aznI5yUq5IFxt7TSRjxzwqSiucjmMcYMVvjJe0WknZ+BE5ncYR3SLa6TgNEn3JBcRTTM6ML52OYhxgxW+Ml9a/voKk6l1UnX+h01EarM8vzqOMGIpnfeh0FOMAK35jvLT3xQ+pRuj9ywlOR2mwuPZxrEkcQ/KaD9FqdTqOaWJW/MZ4qf03H5LTcjgJacF15bhj51xIckUe+fM2OB3FNDErfmO8sGf5LtJLl3Fg+AVORzllPe+YCMD2v9phneGmIVfgeklECkUkp8a6+0Vkp4is9Nzq/B1XRMaLyAYRyROR6b4Mbkwg2PjkxwB0ujl45vdP6DS8Kxti+9P6S5vnDzcNecf/T2B8HeufUNUBntuc2g+KSATwHHA+kA5MEZF0b8IaE2hiFnzIjshkel2c4XSURtmdeSF9iv9rJ20LMw259OJioDE/FUOAPM8lGMuBN4CLG7EfYwLSsQPH6Lv3U7akXxiwJ2U7mcSfXei+JONjc52OYpqQN3P8t4vIas9UUJs6Hu8E7KixXOBZZ0xIyHl2Ec05RovJwTe/f0La1ZkUujrAnO/90m5CWGOL/29AD2AAsBt4zNsgIjJVRLJFJLuoqMjb3Rnjd6Wz51JKMzJuGeV0lEZzRbrYmDyO3tvnU1Ve5XQc00QaVfyquldVq1S1Gvg77mmd2nYCXWosd/asq2+fL6hqpqpmJiYG12FxJjx1y53L2sSziW0d63QUr7gmjKed7if31Wyno5gm0qjiF5GkGos/AXLq2CwLSBGR7iISDUwG7KxQJiRsW5hHckUepWed73QUr/W+YyxVuNj32jyno5gm0pDDOV8HlgCpIlIgIjcAj4jIGhFZDZwN/MKzbUcRmQOgqpXAbcB8IBd4S1XX+mkcxjSprTPcJZl8S/AXf9uUdqxrMYSELPuAN1xEnmwDVZ1Sx+qZ9Wy7C5hQY3kOYJ8amZDT/PO5bI3qRfKYnk5H8Yn9meMZtegB9m/YR7vUBKfjGD+zb+4ac4rKDpWRUfQZ29KC/93+CYnXnI8LZf3TC5yOYpqAFb8xp2jt3xbTnGM0nxQ6xZ92dSb7JAGdZ/P84cCK35hTdOTtuZQRQ8atgX/RlYZyRbrY0HUsqfnzqK6sdjqO8TMrfmNOUZe1c8lpN5rmCc2djuJb559Pohax/l/LnU5i/MyK35hTsGNxPj3KN3B0VOhM85yQesc4qhEKX7aje0KdFb8xpyB/xnwAutw4zuEkvpeQlsj65mfQNmu+01GMn1nxG3MKohctoCCiK93HpzodxS8KB4wl/fBSircXOx3F+JEVvzENVFFaQdruhWxJGRe0Z+M8mTaTxxFJFbnP/cfpKMaPrPiNaaDcl78hnhKiJox1OorfpF0/jMO0oPwjO54/lFnxG9NAB16fTxUu0m47x+kofhPdIpp1p40heeN8uwh7CLPiN6aBEpYvYF2LIbTuXtflJ0JH2Y/G0rUyn+2fbXY6ivETK35jGuDg5gOkHc1i/xmhdzRPbck/d49x6/N2dE+osuI3pgHWP7eQCKppNyV05/dP6Hp2T7ZHdif2C5vnD1VW/MY0QMWcBRQTT9q1dV1zKLSIS9h6+jjS9/yHitIKp+MYP7DiN+YktFrpmTef3I7nEBl70jOZh4ToC8bSkiOsm7nE6SjGD6z4jTmJ/Hkb6FS1g/KzQ39+/4S0aWOoJIKDb9g8fyhqyBW4XhKRQhHJqbHuURFZLyKrReQ9EWldz3O3eq7UtVJE7IKeJijtmOme6+4+9TyHkzSd+K7xrGs5jMRVnzgdxfhBQ97x/xMYX2vdJ0AfVe0HbATu/YHnn62qA1Q1s3ERjXFWsy8XkB+VQpdR3Z2O0qQOZI4l7Wg2BzbtdzqK8bGTFr+qLgYO1Fq3wHNNXYClQGc/ZDPGceVHykkvXMT23qF/NE9tCVeMdV+V668LnY5ifMwXc/w/A+o7j6sCn4rIMhGZ+kM7EZGpIpItItlFRUU+iGWM99bNXEILjhJ7QfhM85zQ+6pMiomncq5N94Qar4pfRO4DKoFZ9WwyUlUHAOcD00RkVH37UtUXVDVTVTMTExO9iWWMzxx6awGVRJB269lOR2lykbGR5HY8hx55C+z0DSGm0cUvItcBFwBXqmqdPxWqutPzZyHwHhD6B0GbkNJ+5QLWthpOq86tnI7iiPLRY+lctZ38+RudjmJ8qFHFLyLjgbuBi1S1tJ5t4kSk5Yn7wFggp65tjQlEBzbtp3fpMg6eEX7TPCck3+Qe+46XbLonlDTkcM7XgSVAqogUiMgNwLNAS+ATz6GaMzzbdhSROZ6ndgC+FJFVwDfAx6o6zy+jMMYP1v91IS6UhCvC74PdE7qO7sG2yJ52+oYQc9KvIarqlDpWz6xn213ABM/9LUB/r9IZ46CqOQs4JK3pfVV4H4m8LXUsA9e+SkVpBVHNo5yOY3zAvrlrTB20Wum++ZOwOk1DfWIu9Jy+4aWlTkcxPmLFb0wd8udvpHPVdirOCt/5/RN633K2+/QNb9p0T6iw4jemDidO03Diw81w5j59w1ASV9h5e0KFFb8xdWj2xQK2RvWi6+geTkcJCAfOcJ++4eDmAyff2AQ8K35janGfpuEztqWG79E8tbWb4jl9w3N2+oZQYMVvTC3fnqbhQiv+E9KuGUwx8VTMsXn+UGDFb0wth95aQAWRYXmahvrY6RtCixW/MbW0XzmfdWF8mob6lJ89zn36hnkbnI5ivGTFb0wN+3KL6F26nIODbZqnthMXojlxxJMJXlb8xtSw8W/u0zQkXmnFX1uXUd3Jj0qh2ZdW/MHOit+YGqrmLeCgtKH3lWc4HSUgbe89lozCzzhectzpKMYLVvzGeGi10mvzAtZ3OpeI6Ain4wSk2AvHEkcpuS8tcTqK8YIVvzEemz9cR1L1TirH2DRPfdJuGU0FkRx6077FG8ys+I3xKPiHe+76xIeY5vtadW7FulbD6bDKij+YWfEb49Hiy3lsju5N5zO7OR0loB0cMo60YysoytnrdBTTSA25EMtLIlIoIjk11rUVkU9EZJPnzzb1PHe8iGwQkTwRme7L4Mb40rEDx8jYv5iCjPFORwl4Ha51/zfa+JxdlStYNeQd/z+B2n8bpgMLVTUFWOhZ/g4RiQCew32h9XRgioike5XWGD9Z+9fPaUYZcZOs+E8mdfJAiiQR5tkF9YLVSYtfVRcDtU/JdzHwsuf+y8CP63jqECBPVbeoajnwhud5xgSco7PncYxYMm4Z5XSUgOeKdLEheRyp2+ZTXVntdBzTCI2d4++gqrs99/fgvr5ubZ2AHTWWCzzrjAk4XdbNIydhNM3aNnM6SlCQ8eNJ0H2s/9dyp6OYRvD6w11VVcDrszaJyFQRyRaR7KKiIm93Z0yDFXy5lR7lGygdOc7pKEHj9GnuI58KX7HpnmDU2OLfKyJJAJ4/C+vYZifQpcZyZ8+6OqnqC6qaqaqZiYmJjYxlzKnb8jf3oYldbrL5/YZKzGjPuuZn0CbLDusMRo0t/g+Aaz33rwXer2ObLCBFRLqLSDQw2fM8YwJK9GfzKIjoRvfxqU5HCSqFA8eTUbKE4m2HnI5iTlFDDud8HVgCpIpIgYjcADwMnCcim4BzPcuISEcRmQOgqpXAbcB8IBd4S1XX+mcYxjRO+ZFy0ncvZMvp4xGXOB0nqLS9YjyRVLHuWbsqV7CJPNkGqjqlnofOqWPbXcCEGstzgDmNTmeMn62buYQBHCb6QpvfP1Vp1w2leForKj6cB49OcjqOOQX2zV0T1g69Mc99ta1pY5yOEnSimkeR2/FceuXNs6tyBRkrfhPWOq78mJzWI4nvGu90lKBUcd4EOlYVsOndNU5HMafAit+ErV1f7+D0sjUUj5hw8o1NnVLuOB+AXS/ajG4wseI3YSvvaXdZdf75RIeTBK/TBnUkt9lA2iz52Oko5hRY8ZuwFbPwY3ZEJtPzgjSnowS1vWdMpE/JVxzcXPvMLiZQWfGbsFR2qIw+exeypfcEO4zTS+2unkAE1eQ+ZdfiDRZW/CYs5Tz3OXGU0vxSm+bxVvp1Q9gnCVR/ZNM9wcKK34Slo2/P4Rix9Ln9bKejBL2I6AjWJ4+n99Z5VJVXOR3HNIAVvwk7Wq0kr/2YnMQxdjZOH3FNnECC7iP3lSyno5gGsOI3YSd//ka6VW6mdIxN8/hK2l3jqMLFvpdtuicYWPGbsLP9efdhnD1vs+P3faVNz7bktBpBh2VW/MHAit+EnfjPPyQvJp3OI5OdjhJSDg6fSNqxFezOrvfs6yZAWPGbsHJw8wH6HlpMwSC7Cqivdb71IgA2PWZnXw90VvwmrKx7bC6RVJFwgxW/r/W8II2tUb1o9qkVf6Cz4jdhxfXh++x1nUb6tYOdjhJyxCVs7XcR/fb9h8O7Djsdx/wAK34TNo6XHKdPwVw2nH4hrkj70feHNtdeTAzl5PzFrsUbyBr90y8iqSKyssatRETuqrXNaBEprrHNH7yPbEzjrHlmES05QrPLbZrHXzJuGsF+aUfVv226J5Cd9Apc9VHVDcAAABGJwH0h9ffq2PQLVb2gsa9jjK+Uvv4+R2lOnzvsoiv+EhkbSW73iWTkf0hFaQVRzaOcjmTq4Kvfd88BNqvqNh/tzxif0molJfcD1iSNs2/r+lnkpItpowfJmfGl01FMPXxV/JOB1+t5bISIrBaRuSKS4aPXM+aUrP/XcpKqd1I50aZ5/K3PL8dSRgwlr9l0T6DyuvhFJBq4CHi7joeXA11VtR/wDPDvH9jPVBHJFpHsoqIib2MZ8x17X3ifKlyk/dpO0+BvLU5rwZrEc+i+5n27Fm+A8sU7/vOB5aq6t/YDqlqiqkc89+cAUSKSUNdOVPUFVc1U1czExEQfxDLm/3T65j3WxI+kXWqdP37Gx46N/zFdK/PZ+M5qp6OYOvii+KdQzzSPiJwmIuK5P8Tzevt98JrGNNiWOetJOZ5D8bmXOh0lbKTf+2OqcLH7mXecjmLq4FXxi0gccB7wbo11N4vIzZ7FS4EcEVkFPA1MVlX73c80qe1PzAbg9OmXOJwkfCSkJbKqzWi6fvO2TfcEIK+KX1WPqmo7VS2usW6Gqs7w3H9WVTNUtb+qDlPVr7wNbMypSvrybVa3HEFSZieno4SVI+MupUf5BvLeX+t0FFOLfX3RhLStn2witWwVB8bYNE9T633vT6hG2PmUTfcEGit+E9K2Puae5km5Z5LDScJP+36nsTr+R3RaasUfaKz4TUhr/8U75MQNodPwrk5HCUvF515KyvG1bP4o1+kopgYrfhOydizOJ710GfvOsmkep5z4QH3Hk7MdTmJqsuI3IWvzI+6y6fEbm+ZxSlJmJ/cH6/+16Z5AYsVvQlbiorfIbTaIrqN7OB0lrB0YcympZavYMneD01GMhxW/CUlbP9lExtEs9o6Z4nSUsJf6u59SjbD94X85HcV4WPGbkLT1wVlUI6Teb8XvtKTMTqxsczbdv5plX+YKEFb8JuRotZL81SxWtR5tX9oKEEcvvpJulZtZ+49vnI5isOI3IWjdy1kkV+Rx5OIrnY5iPPo9MIkyYtj/9Cynoxis+E0I2vfULI4TTb8H7GieQBHfNZ4VnS4gfc2bVJZVOh0n7Fnxm5BSWVZJ+uo3WNHxAuK7tXY6jqlBrrySRC1k5V8+dTpK2LPiNyFl5WMLSdRCuNKmeQLNwPsmcEhaU/aSTfc4zYrfhJSyl2ZxSFoz4LcTnI5iaolpFcPq1MsYkP8eRwuPOh0nrFnxm5BRUlDCwC2zWZ16GbGtY52OY+oQf+tVtOAoK35np3BwkhW/CRkrp79BHKW0/fUNTkcx9eg37UdsjepFyzdfdDpKWPP2ClxbRWSNiKwUkew6HhcReVpE8kRktYgM8ub1jPkh7f79Ipti+pBx/RCno5h6iEvYes6N9C/5wk7h4CBfvOM/W1UHqGpmHY+dD6R4blOBv/ng9Yz5no3vrCbjaBa7JtyIuMTpOOYHpP/5WiqJYPv9LzkdJWz5e6rnYuAVdVsKtBaRJD+/pglDu/80k+NE0/fPVzkdxZxE+36nkZ10IRlZ/6SitMLpOGHJ2+JX4FMRWSYiU+t4vBOwo8ZygWfd94jIVBHJFpHsoqIiL2OZcFJ2qIx+q19lWddLaJvSzuk4pgEipt5Iohay7IGPnI4Slrwt/pGqOgD3lM40ERnV2B2p6guqmqmqmYmJiV7GMuFk+e/fo40eJPa2G52OYhpo4PRx7HZ1IuIf9iGvE7wqflXd6fmzEHgPqP2p2k6gS43lzp51xvhM7KwX2R7ZnQG/ONvpKKaBImMj2TDiegYVzWPX1ztO/gTjU40ufhGJE5GWJ+4DY4GcWpt9AFzjObpnGFCsqrsbndaYWvI+WMegg/9hy9k34oq0o5ODSa+H3IfdbvzlDIeThB9v/qZ0AL4UkVXAN8DHqjpPRG4WkZs928wBtgB5wN+BW71Ka0wtu6Y/TRkx9H2mro+YTCDrPDKZrKSL6LvkeY4dOOZ0nLDS6OJX1S2q2t9zy1DVBz3rZ6jqDM99VdVpqtpTVfuq6veO9TemsQ7lHyQz9xWyUq6kXWqC03FMI8T8+g7a6X6yf/W601HCiv1ubILWyttepDnH6PC/dzodxTTSgLtGsyG2H+3feMquztWErPhNUKosq6TX/GdZ0Xo0p1/az+k4ppHEJRRefgepZatZ9fTnTscJG1b8Jihl/f59Oldtp/xme7cf7DIfv4L90o6yR55yOkrYsOI3Qan5C0+xPbI7mQ9c6HQU46VmbZuxZvjPGbL7fXYsznc6Tliw4jdBZ83zX7lP8jXxdiKiI5yOY3wg9albqSSSLdP+4nSUsGDFb4JO2e//xD5JYPALdghnqEjK7MTS3tcxNGcme5bvcjpOyLPiN0El97VlDC6aS855vySufZzTcYwPJf9tOpFUsv4me9fvb1b8JqgU3/0nDklrBs2c5nQU42NdR/dgaY8rGbJ8BkVrC52OE9Ks+E3Q2Dh7DcN2/5uVo+6kVedWTscxftDx2d8SSxlrb3zC6SghzYrfBI19v3iQw7Sg/8w7nI5i/KTH+aks7fJTzlj6LAc3H3A6Tsiy4jdBYePsNQzb8RbLhk6jTc+2TscxfpTwxH205Airrra5fn+x4jdBofjmuymRePrPutvpKMbPTp/Ul/8mX8nQJU/YKZv9xIrfBLxlD3/C4H3zWDnxd/ZuP0x0e+1BBGXLFb9zOkpIsuI3Aa2qvIoWD/yGHZHJDJ91m9NxTBPpfGY3lg6+gxFbXmXDmyudjhNyrPhNQFsy7TVSy1ax45aHiGkV43Qc04QGvv1bDkkbjt7yaztzp49Z8ZuAdbTwKD1fuo+1cYMZ/uTlTscxTSy+W2vW/PgPDDq4kOz/N9fpOCHFm0svdhGRz0RknYisFZHvnSZRREaLSLGIrPTc/uBdXBNOsibeT1L1TqoeeRxxidNxjAOGv3IL+VGn0+FPt1G6r9TpOCHDm3f8lcCvVDUdGAZME5H0Orb7QlUHeG5/9OL1TBjJnbWcH2U/zuLeN9Hv1pFOxzEOiW4RTfGjL9C1Mp9vJtzvdJyQ4c2lF3er6nLP/cNALtDJV8FM+Kosq4SbbmK/K5F+c/7sdBzjsAF3nsUXp9/AyKzHWf/6CqfjhASfzPGLSDIwEPi6jodHiMhqEZkrIhk/sI+pIpItItlFRUW+iGWC1JeXPUXaseVsvvMZWndv43QcEwD6znuUA65U05oRAAAMSUlEQVQEqm+8yf3GwHjF6+IXkRbAbOAuVS2p9fByoKuq9gOeAf5d335U9QVVzVTVzMTERG9jmSC1bWEegz/6A193uJBhf7nU6TgmQLTu3oa8O54mvXQZX06y8/h4y6viF5Eo3KU/S1Xfrf24qpao6hHP/TlAlIgkePOaJnQdLznO0YsmUy4xdHn/OftA13zH8McuY2nSTxgx5z7WvZLtdJyg5s1RPQLMBHJV9fF6tjnNsx0iMsTzevsb+5omtC0dPZ300mVsvOclOg7t4nQcE2DEJaR+8SKFEUnE3XA5xduLnY4UtLx5x38mcDUwpsbhmhNE5GYRudmzzaVAjoisAp4GJquqfRPDfM83v/+Qs1Y8yef9bmfoQz92Oo4JUG16tuXAM/+iU+U21o78uX2xq5EkEHs4MzNTs7PtV7lwUfDfbTT/0SAKY7uSvGsJsa1jnY5kAtyisf/L6E/uY/EVMxg16+dOxwkIIrJMVTMbsq19c9c4qnh7McfOmUiEVhLz3ptW+qZBRs2ZTna7cQz/120sf+RTp+MEHSt+45iK0gryBl5K8vENbH7kXbqPO93pSCZIuCJdpCx/k60xvel5zyQ2vZfjdKSgYsVvHKHVytIBN3PGgU/5+oa/M+g35zgdyQSZ+K7xNP/sY4654mh22UT2rtztdKSgYcVvmpxWK58PvZsfbXqJRaP+wMgXr3M6kglSnYZ35eArH9G6aj/Fw8fZRdobyIrfNCmtVhYPuovR2X/h877TOOuz+52OZIJc2pWD2PjI+3Quy6Nk0Gh7598AVvymyVRXVvNFn1s4a9XTLDrjl4xa+Yx9Scv4xKDfnMPGJ+fSoXw7pUPOYndWgdORApoVv2kSR/Yc4evknzIq93kWDb+Xs775i5W+8akBd55F/owFtK3YS/XwEeTOWu50pIBlxW/8bvuiLezqPoIhO99j0cWPc9aXD1rpG7/o+/MR7HnjcxQh+aoz+er2152OFJCs+I1fZT0whxZjBpN4vICVD81j9L9/YaVv/Cr18gHErMpiQ/wQRjx7BYsyf8XxkuNOxwooVvzGL4q3HeKLlOsZfP9E9sV0ouTTLM6Yfp7TsUyYSMxoT8auT/m87zRGL3ucHe0Hse7lLKdjBQwrfuNTWq0svftdSnv0YXjeqywa8Vu67c2i25ieTkczYSaqeRRnrX6WrAfmEFdRzOnXDWfR0HsoKah99vjwY8VvfGbN81+xpvVIhj06iSNRbdj48lJG//dBYlrFOB3NhLHBfzif5ltyWJJyLaO/eYTjXXvx+WXPUn6k3OlojrHiN17RamX5I5+S1X4ifW8+k/ZH81l89Qt0P7SC9GsadL4oY/wuvltrfrRxJutezqIgPoOz3rmd3W3S+PzSZziy54jT8ZqcFb9plAOb9rP4ihlsiuvPoHvOI3lfNovO/RMtd29i1Cs3ERkb6XREY74n/ZpMBuz/D1n/8xElse05a/YdVCV1ZlHmr9j4zmqn4zUZOy2zabCitYVseGY+0e+9xcDCeURRyYbYfhROuYvBj0+xM2uaoJPz4lJK/vgkQ3a8QyRVbIrpw85RU+g49QJSLukbVEegncppmb0qfhEZDzwFRAAvqurDtR4Xz+MTgFLgOlU96bcqrPgDQ+HqPeS/vpRj/1lChzWfknbM/b9ut6sTGwZNIenXV3L6Zf2D6i+HMXXZl1vE2gfeps2cWfQ7/BUAe1xJ5CWfR/WZI0n6yXB6TEwjIjrC4aT1a5LiF5EIYCNwHlAAZAFTVHVdjW0mALfjLv6hwFOqOvRk+7bibzpV5VXsXbGLom/yKVmxmaqcXOK25tDpwBo6Vrm/9l5OFLmthnJwyHg6XDue1MkDcUXaLKEJTbuzd5L31wVEfDqf1IJPaafuq8UepgVbW/TlYOc+aFoGzfv1om1mDzqOSKZZ22YOp2664h8O3K+q4zzL9wKo6kM1tnkeWKSqr3uWNwCjVfUHz6JkxX9y1ZXVVJRWUH6knIqj5ZQfPk55SRnlxccoP1RK+aFSKg4dpWJfMZUHSqg+WAz79xNxaD/RJfuIO7ybtmW7SKzaQxSV3+63jBi2xfZmX4cMKgZk0m7CMFJ+OtCmcUxY0mpl28I8CmYvpeqrr4nfnkOXkpxv/zE4Yb+0Y39MR0rikihr1Z6q+HZo23ZI2zZEtGlFVEI8UW1aEBXfnJi2cUTHNyO6Vaz71jKGqOZRRLeI9uo3ilMpfm8+gesE7KixXID7Xf3JtukE+OX0eeviMomuOnZKzxGvPuOo+7lSY/3/7V+/Xe9e514WtI7lagTFRTUurSKCqm//jKSSKCqIoJoY4FQOlDxONAdcCRyObkdJXBKbk9LZmJiEq0cycX26kzC4O53OTCY1NpLUxvznMCbEiEtIPi+F5PNScF9i3P2PQVFuEXuXbKF4xRYqNmzBtWcXMft3EXd4N+2LNxGfv594Tv37Antdp9Ghyv9nFw2YQy9EZCowFaBr166N2sfB9r1xVTTiq9nS+Dlqrfe5NdZ7tnFvK9+uUxHPY/LdZXGhLheIC1wu1BUBERFoRARERkFkpPsWHQ3R0Uh0FNK8Ga7msbiaxRAZH0dU6zii28TRrEMrmp/Wipad42me0Jwkl5DU6NEaY8QlJGa0JzGjPTCs3u3Kj5RzeGcJR3YWU7q7mPKDR92/hRcfpfrIMaqPHae6tAw9VgYVFe5bbCwdmmAM3hT/TqBLjeXOnnWnug0AqvoC8AK4p3oaE+jM/Nca8zRjjPG56BbRtEtNoF1qgtNRvsebT+iygBQR6S4i0cBk4INa23wAXCNuw4Dik83vG2OM8a9Gv+NX1UoRuQ2Yj/twzpdUda2I3Ox5fAYwB/cRPXm4D+e83vvIxhhjvOHVHL+qzsFd7jXXzahxX4Fp3ryGMcYY37KDsY0xJsxY8RtjTJix4jfGmDBjxW+MMWHGit8YY8JMQJ6WWUSKgG2n8JQEYJ+f4gSqcBwzhOe4w3HMEJ7j9mbM3VQ1sSEbBmTxnyoRyW7oyYlCRTiOGcJz3OE4ZgjPcTfVmG2qxxhjwowVvzHGhJlQKf4XnA7ggHAcM4TnuMNxzBCe426SMYfEHL8xxpiGC5V3/MYYYxooqItfRMaLyAYRyROR6U7n8RcR6SIin4nIOhFZKyJ3eta3FZFPRGST5882Tmf1NRGJEJEVIvKRZzkcxtxaRN4RkfUikisiw0N93CLyC8/Pdo6IvC4isaE4ZhF5SUQKRSSnxrp6xyki93r6bYOIjPNVjqAtfs/F3p8DzgfSgSkiku5sKr+pBH6lqum4L/kzzTPW6cBCVU0BFnqWQ82dQG6N5XAY81PAPFXtDfTHPf6QHbeIdALuADJVtQ/u07xPJjTH/E9gfK11dY7T83d8MpDhec5fPb3ntaAtfmAIkKeqW1S1HHgDuNjhTH6hqrtVdbnn/mHcRdAJ93hf9mz2MvBjZxL6h4h0BiYCL9ZYHepjjgdGATMBVLVcVQ8R4uPGfYr4ZiISCTQHdhGCY1bVxcCBWqvrG+fFwBuqelxV83Ff12SIL3IEc/HXdyH3kCYiycBA4GugQ40rmu2BJrlcZ1N6ErgbqK6xLtTH3B0oAv7hmeJ6UUTiCOFxq+pO4C/AdmA37iv1LSCEx1xLfeP0W8cFc/GHHRFpAcwG7lLVkpqPeS56EzKHaInIBUChqi6rb5tQG7NHJDAI+JuqDgSOUmuKI9TG7ZnTvhj3P3odgTgRuarmNqE25vo01TiDufgbfCH3UCAiUbhLf5aqvutZvVdEkjyPJwGFTuXzgzOBi0RkK+5pvDEi8hqhPWZwv6srUNWvPcvv4P6HIJTHfS6Qr6pFqloBvAuMILTHXFN94/RbxwVz8TfkYu8hQUQE95xvrqo+XuOhD4BrPfevBd5v6mz+oqr3qmpnVU3G/f/2P6p6FSE8ZgBV3QPsEJFUz6pzgHWE9ri3A8NEpLnnZ/0c3J9jhfKYa6pvnB8Ak0UkRkS6AynANz55RVUN2hvuC7lvBDYD9zmdx4/jHIn717/VwErPbQLQDvdRAJuAT4G2Tmf10/hHAx957of8mIEBQLbn//e/gTahPm7gAWA9kAO8CsSE4piB13F/jlGB+7e7G35onMB9nn7bAJzvqxz2zV1jjAkzwTzVY4wxphGs+I0xJsxY8RtjTJix4jfGmDBjxW+MMWHGit8YY8KMFb8xxoQZK35jjAkz/x+8gnCHEzihyAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "def func(x, a1, a2, a3):\n", " return a1*np.exp(-(x-a2)**2/a3**2)\n", "\n", "ndata = 100\n", "xdata = np.linspace(1, ndata, ndata)\n", "y = func(xdata, 20, 60, 15)\n", "ydata = y\n", "plt.plot(xdata, ydata, 'b-', label='data')\n", "\n", "popt, pcov = curve_fit(func, xdata, ydata)\n", "plt.plot(xdata, func(xdata, *popt), 'r-', label='fit')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "from pprint import pprint\n", "import scipy.linalg as linalg\n", "\n", "def dfda1(x,a1,a2,a3):\n", " return np.exp(-(x - a2) ** 2 / a3 ** 2 / 2)\n", "def dfda2(x,a1,a2,a3):\n", " return a1 * (x - a2) / a3 ** 2 * np.exp(-(x - a2) ** 2 / a3 ** 2 / 2)\n", "def dfda3(x,a1,a2,a3):\n", " return a1 * (x - a2) ** 2 / a3 ** 3 * np.exp(-(x - a2) ** 2 / a3 ** 2 / 2)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "nparam = 3\n", "guess1 = [10,50,10]" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "array([ 19.866, 59.841, 14.878])\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VHX6/vH3k0oSOgkdDSgCQRAlgihFkYRQpOki4CoqLuqqq65+d3WLq7u6u+66doVFRVEUC6L0Joq9EJASSugdQiBAIL08vz9m2F+MCYTMJGfK87quXJlz5syc+6N4O3zmFFFVjDHGBI8QpwMYY4ypXVb8xhgTZKz4jTEmyFjxG2NMkLHiN8aYIGPFb4wxQcaK3xhjgowVvzHGBBkrfmOMCTJhTgeoSGxsrMbHxzsdwxhj/MbKlSsPq2pcVbb1yeKPj48nNTXV6RjGGOM3RGRXVbe1qR5jjAkyVvzGGBNkrPiNMSbIWPEbY0yQseI3xpggc8biF5E2IvKZiGwQkfUicq97fWMRWSoiW9y/G1Xy+hQRSReRrSLykLcHYIwx5uxU5RN/MfCAqiYAlwF3iUgC8BCwTFXbA8vcyz8hIqHAS8AgIAEY636tMcYYh5yx+FX1gKqucj8+AWwEWgHDgWnuzaYBIyp4eQ9gq6puV9VC4F3364wxZ5CVBe++Cy+8ABs2gN0l1XjLWc3xi0g8cDHwPdBMVQ+4nzoINKvgJa2APWWW97rXVfTeE0UkVURSMzMzzyaWMQHlgw+gVy+Ii4OJY7N57DeH6dwZ2raFBx+EvDynExp/V+XiF5G6wIfAfaqaXfY5dd2x3aPPI6o6RVUTVTUxLq5KZx0bE1BU4ZFH4JXRS/jzhrFkNOpANg3IDGnG7oQU7op7jxf/k0///pCR4XRa48+qVPwiEo6r9N9W1Vnu1Rki0sL9fAvgUAUv3Qe0KbPc2r3OGFNGfj78ckwxdf72B5YwkEHRy4ntkwB//Svy8MO0ObGB/0sdw9HG51Hy41p69oS0NKdTG39VlaN6BHgN2KiqT5d5ag4w3v14PDC7gpevANqLSFsRiQDGuF9njHErLYVbh2Qw4f1k/sA/0Nt+hezYAR99BH/+Mzz+OOzYAYsWERUlfBN5JZ1Pfs8VV8DWrU6nN/6oKp/4rwBuBPqLyGr3z2Dgn0CSiGwBBriXEZGWIrIAQFWLgbuBxbi+FH5fVdfXwDiM8Vsv/DOHP3x6NX3Cv4M33kBemQJ16vx0o9BQGDgQvvySsNhGzM0fQN/S5YwdC4WFzuQ2/kvUBw8VSExMVLs6pwkGK1OVjT3GM06nI4sXI8lJZ37R/v2QlETJ1u10KVzJsN8n8M9/1nxW49tEZKWqJlZlWztz1xiHnDwJHw15lV/qW+T//i9VK32Ali1h2TJC68UwP+5m/vNkMcuW1WxWE1is+I1xyDPjV/OnQ/eQlZhE9BN/OrsXN28OL79M28wV/Cv2X9x4Ixw7VjM5TeCx4jfGAWmrixk565cUxDSh8fzprjn8szV6NIwezb3HHqXJgXX8+9/ez2kCkxW/MQ5YfssbXMh6Ql58Hpo2rf4bvfQSIY0bMafReF58poiDB72X0QQuK35jatkPn+UwavUj7GnTi3rjR3n2ZrGxrimfoz8yOv9NnnjCOxlNYLPiN6YWqcLaCc/SkgM0ee3fIOL5m44aBZdeyj/qPs5rk4vYscPztzSBzYrfmFr0+QeHGL3jSbZ1HUl00hXeeVMRePRRYk/sZDzT+MtfvPO2JnBZ8RtTS1Th0D1/I5pc2kz/h3fffNAg6NmTJ6If5723Ctm40btvbwKLFb8xtSR1wSGGH5rC5ituJaJLB+++uftTf+PsXdwW9gYvvODdtzeBxYrfmFqy86HJRFLIuc8/UDM7GDgQLruMxyKfYMa0Qo4fr5ndGP9nxW9MLdi/s5DeaZPYGJ9CzCVe/rR/yqm5/pzdDMl9n9dfr5ndGP9nxW9MLfj2/vdpwUEaPHJfze4oORkuuIDf1ZvESy+5rvxpTHlW/MbUsMIC5fx5z7InpiMtb06u2Z2JwB130PXEN0RvXcOiRTW7O+OfrPiNqWHL//4NFxWv5MQtv/HOcftnMn48WqcOD0RPsi95TYWs+I2pYaEvPsfxkIZ0/PtNtbPDxo2RMWMYUzydrxdls3lz7ezW+I+q3IFrqogcEpG0MuveK3NTlp0isrqS1+4UkXXu7ewC+ybopH+RQb+sWWzpO4GQejG1t+Nf/5qIwhxukulMm1Z7uzX+oSqf+N8AUsquUNXrVbWbqnbDdS/eWRW90O0q97ZVukGAMYFky19nEEYJ8Y/eUrs7vvRS6N6dB2NeZvpbal/ymp84Y/Gr6hdAVkXPue/HOxqY4eVcxvi90lKI/2IaW+p3J7Zf59oPcOedxJ9cT5s9X/PFF7W/e+O7PJ3j7wNkqOqWSp5X4BMRWSkiEz3clzF+JXXqWi4sWk32yPHOBLj+ejQ6mlvCp/PWW85EML7J0+Ify+k/7fd2TwcNAu4Skb6VbSgiE0UkVURSMzMzPYxljPOynn2TIsJI+NtYZwLUrYuMGMH1IR/w8fuF5OU5E8P4nmoXv4iEAaOA9yrbRlX3uX8fAj4Cepxm2ymqmqiqiXFxcdWNZYxPyM0u5uIN01l3zhCi2sQ6F2TcOOoWZHHFyUXMnu1cDONbPPnEPwDYpKp7K3pSRGJEpN6px0AykFbRtsYEmh+eWEozzSDiNoemeU5JTkZjY7kt6h2b7jH/U5XDOWcA3wIdRGSviExwPzWGctM8ItJSRBa4F5sBX4nIGuAHYL6q2nmEJijIW9M4GtKYhAcHOxskPBwZPZqUojl8vegEGRnOxjG+oSpH9YxV1RaqGq6qrVX1Nff6m1V1crlt96vqYPfj7ap6kfuns6raTeFMUDi84wQ9DsxmU7cxhERFOh0HbriBiOI8hpV+xPvvOx3G+AI7c9cYL1v793lEkU/s3Q59qVter14QH8/Euu/w4YdOhzG+wIrfGC+LmDuTjNAWnH/T5U5HcRGBceO4PGcp6V9kcOiQ04GM06z4jfGio3tO0j1jAVu7jkJCfeg/rxtuIERLuVY/4OOPnQ5jnOZDfzKN8X9rn1xIFPk0uu06p6P8VEIC2qkT46I+YuZMp8MYp1nxG+NFoR/N5HBIHJ0m9nE6ys/IyJH0zP+cNZ8eIavCi7CYYGHFb4yXnMjIpdv++aQnjELCQp2O83MjRxKqJQwqmcucOU6HMU6y4jfGS9b8azF1yaH+hF84HaVi3bujbdpwQ9RHdnRPkLPiN8ZL9IOZZEkTEu7s53SUiokgI0dyZeFivl58kuxspwMZp1jxG+MFeccKuGjPXDZ0HEloZJjTcSo3ahThJQX0L1rEvHlOhzFOseI3xgvWPvcZ9TlBzA0jnY5yer17o7GxjKszy+b5g5gVvzFekPf+XHKIpvM9/Z2OcnqhociwYaSUzGfZwkKKipwOZJxgxW+Mh7RUaZ8+l7SWyUTUr+N0nDMbNYroomwSs5fx1VdOhzFOsOI3xkOb3l9Lq5I9FKdc43SUqrn6ajQmhhEhc5k71+kwxglW/MZ46MCUuZQidPjtEKejVE2dOkhSEiMj5jF3jqLqdCBT26z4jfFQs+/nsLFuD2I7N3M6StUNHUrT/D1EbVtHerrTYUxts+I3xgMHfzxA59wVHO7lJ9M8pwx23SBmCPPtsM4gVJU7cE0VkUMiklZm3aMisk9EVrt/KrzNkIikiEi6iGwVkYe8GdwYX7D5mfkAtLzdz4q/RQvo3p3ro+fZPH8Qqson/jeAlArWP6Oq3dw/C8o/KSKhwEvAICABGCsiCZ6ENcbXRC6Zy97Qczl/ZBeno5y9IUPomvcdm746bBdtCzJVufXiF0B1/lj0ALa6b8FYCLwLDK/G+xjjk/Ky8uiSsZRtna9BQsTpOGdv6FBCtJSk0kUsXOh0GFObPJnjv0dE1rqnghpV8HwrYE+Z5b3udcYEhLQXlxNNHnWvH+p0lOrp3h1t1oxrI+ez4Gd/ZzeBrLrFPwloB3QDDgD/8TSIiEwUkVQRSc3MzPT07YypcbkfLiSXKDr/2kcvynYmISHI4MEk6yKWLSqipMTpQKa2VKv4VTVDVUtUtRR4Bde0Tnn7gDZlllu711X2nlNUNVFVE+Pi4qoTy5hade7GhayPu4o6Df3gbN3KDB1KTOExOmR9Q2qq02FMbalW8YtIizKLI4G0CjZbAbQXkbYiEgGMAeyyUCYg7Fq2lfiireT2G+R0FM8kJaHh4QxlPosWOR3G1JaqHM45A/gW6CAie0VkAvAvEVknImuBq4D73du2FJEFAKpaDNwNLAY2Au+r6voaGocxtWrXZNe3ofF3+nnx16uH9O7NiOjF9gVvEDnjhcNVdWwFq1+rZNv9wOAyywsA+9rIBJyozxeyI7w9bfuf53QUz6Wk0P6z37Pn+/0cPtyS2FinA5maZmfuGnOW8o/mcWHmZ+xK8PNP+6ekuE7TSWYxS5Y4nMXUCit+Y87S+pc/J4p8YkZVdF6jH+rSBW3RguERi2yeP0hY8RtzlnJmLiSPOnS+60qno3iHCJKSwgCWsnRhMaWlTgcyNc2K35izdM76haTFXkl0kyino3hPSgp1C48Sf3gFq1Y5HcbUNCt+Y87Cns+3E1+0hZw+ATK/f8qAAWhICIOwyzcEAyt+Y87CjsmLAWjzqwCZ3z+lcWOkZ09GxSxi8WKnw5iaZsVvzFmI/Hwxe0LjaTewvdNRvC8lhYScFWz59jDHjzsdxtQkK35jqqgot4hOBz5lR/tk/7wa55mkpBCC0r90KZ9+6nQYU5Os+I2pok1vfEd9ThA2dKDTUWpG9+5okyZcE7bQjucPcFb8xlTRkRlLKCaUhLv6Ox2lZoSGIklJpIQuZfEiuwl7ILPiN6aK4n5czIa6PWkY39DpKDUnOZnGBQeJ2ZnGtm1OhzE1xYrfmCo4uvUInXJSOZIYoNM8pyQlAZDMEju6J4BZ8RtTBZte+oQQlCZjk52OUrNat0YTEhgetcTm+QOYFb8xVVCyYAnHaEjC+EudjlLjJDmZywq/4JtleRQVOZ3G1AQrfmPOQEuV87YtZkOrAYRFhjodp+YlJxNRkk+3nK/49lunw5iaYMVvzBnsWLCRFiX7KOof4PP7p/Tti0ZEkCI2zx+oqnIHrqkickhE0sqs+7eIbBKRtSLykYhUeJiDiOx036lrtYjYHT2NX9o71TXZ3W5iksNJaklMDNK7N8OjlrB0qdNhTE2oyif+N4DyFyZZClyoql2BzcDDp3n9VaraTVUTqxfRGGdFf7WE7eEdaNP7XKej1J7kZM7PXcveFQc4csTpMMbbzlj8qvoFkFVu3RL3PXUBvgNa10A2YxxXeKKATpmfs6djkHzaP2Wga1rraj5h2TKHsxiv88Yc/61AZRdyVeATEVkpIhNP9yYiMlFEUkUkNTMz0wuxjPHcxte+IYZcIq8J8MM4y+vaFW3alKHhi226JwB5VPwi8kegGHi7kk16q2o3YBBwl4j0rey9VHWKqiaqamJcXJwnsYzxmqMfLKWIMBJ+faXTUWpXSAgyYADJIZ+wdHGpXb4hwFS7+EXkZmAocINqxX8sVHWf+/ch4COgR3X3Z4wTmq5ewvr6vajfqp7TUWpfcjKNCjJosGcdmzc7HcZ4U7WKX0RSgN8Bw1Q1t5JtYkSk3qnHQDKQVtG2xviirM2H6Zi7iqOJQTbNc8qAAQAksdSmewJMVQ7nnAF8C3QQkb0iMgF4EagHLHUfqjnZvW1LEVngfmkz4CsRWQP8AMxX1UU1MgpjakD6y8sIQYkbF2Rf7J7SqhV07swIu3xDwAk70waqOraC1a9Vsu1+YLD78XbgIo/SGeOg4oVLOCYN6fjLID4SOTmZHs9P4ttP8ygqiiI83OlAxhvszF1jKqClSrttS9nYMkgu01CZpKT/Xb7hu++cDmO8xYrfmArsXJxOq5I9FF0ZpNM8p5S5fINN9wQOK35jKrD7NVfLxf8qyIs/Jga54gqGRS216/YEECt+YyoQ9eVSdoS355x+bZ2O4rzkZNrnrmH3igyyss68ufF9VvzGlFN4spCEQ5+xu0OQHsZZXrLrn4NdviFwWPEbU86G176lLjnUGWbFD0C3bmhsLNeEL7Z5/gBhxW9MOcfeX0IRYXS680qno/gG9+UbBoYsZclitcs3BAArfmPKiVu9hA31e1G/dX2no/iO5GQaFRykwZ51pKc7HcZ4yorfmDIObzpMp9yVHL3Upnl+wj3Pn4wd1hkIrPiNKSN9kvsyDTdY8f/Eqcs3RFvxBwIrfmPKKF24hKPSiI43dHc6iu9JTqZnwRd892kuBQVOhzGesOI3xk1LlfO2L2FTqwGERgTxZRoqM3Ag4SUFdM/7km+/dTqM8YQVvzFuW+dtomXJXor72zRPhfr0QSMjGSSL7SxeP2fFb4zbvqmuNms7Mcgv01CZ6GikTx+GRS2x4vdzVvzGuMV8tZjtER1ofcW5TkfxXQMH0i53PRk/7iMjw+kwprqqciOWqSJySETSyqxrLCJLRWSL+3ejSl6bIiLpIrJVRB7yZnBjvCkvK48LjyxnT0KK01F8W5nDOu2uXP6rKp/43wDK/9fwELBMVdsDy9zLPyEiocBLuG60ngCMFZEEj9IaU0PSJn1JFPnUvXag01F8W5cuaPPmXBO5hEV2Pz2/dcbiV9UvgPLX5BsOTHM/ngaMqOClPYCtqrpdVQuBd92vM8bn5Hy4iHwiSbizn9NRfJsIkpxMki5h6aISSkudDmSqo7pz/M1U9YD78UFc99ctrxWwp8zyXvc6Y3xO6/WL2NCkL1FNop2O4vtSUqhXmEX8kVRWrXI6jKkOj7/cVVUFPL5sk4hMFJFUEUnNzMz09O2MqbI9X+/m/MKNnOht8/tVkpyMipDCIpvu8VPVLf4MEWkB4P59qIJt9gFtyiy3dq+rkKpOUdVEVU2Mi4urZixjzt72Sa5jE9tMsPn9KmnSBOnRg+tiFtphnX6qusU/BxjvfjwemF3BNiuA9iLSVkQigDHu1xnjU8I/W8z+0Na0HWLHHlRZSgqdc34g/ZsjHDvmdBhztqpyOOcM4Fugg4jsFZEJwD+BJBHZAgxwLyMiLUVkAYCqFgN3A4uBjcD7qrq+ZoZhTPUU5hTRef9StrcfiISI03H8x6BBhKD0L11qd+XyQ2Fn2kBVx1by1NUVbLsfGFxmeQGwoNrpjKlh66d+z8VkEz7U5vfPSmIi2rgxw04sZNGiMVx7rdOBzNmwM3dNUDv67iJKCKHj3T/7HGNOJzQUSU5mUMhiFi8stbty+RkrfhPUWqxawPoGl9Pg3ApPPjenM2gQjQoyaLJvDevWOR3GnA0rfhO09q3YT6f8HznWa/CZNzY/N9B1FNQgFrLAJnT9ihW/CVpbXnAdhN5q4hCHk/ipZs3gkkv4RcxC5s93Oow5G1b8JmhFLJ3PgdBWtBvexeko/islha6537Lh66Nklb+wi/FZVvwmKOVnF3LhwaVs6zDYDuP0xNChhGoJSbrY7sXrR6z4TVBaN/lr6nOCOqNsmscjPXqgsbGMipxv0z1+xIrfBKUT786nkHA6/8YO4/RIaCgyaBCDWMiShSWUlDgdyFSFFb8JOqpwTtoC1jfpR1RcXafj+L8hQ6hXcITzjnzPihVOhzFVYcVvgs6OT3dwftFGcq60aR6vGDgQDQ3lGubZdI+fsOI3QWfny66DzuN/bcfve0XDhkjv3vwi2ub5/YUVvwk6McvnsSvifFr3v8DpKIFjyBDOz13L4R93s6/Si68bX2HFb4JK1q4TdMv6lN3dhjkdJbAMHQrAYBYwxy6+7vOs+E1QSXt6CZEU0uRmK36v6tgRbduW0dHzrfj9gBW/CSo6ezZZ0piOE65wOkpgEUGGDqVPwSd8uyyXEyecDmROx4rfBI2CnGK67JpP+nlDCIk4460ozNkaNozwknz6FS21e/H6uGoXv4h0EJHVZX6yReS+cttcKSLHy2zziOeRjameNZO+oTFZRFxn0zw1ol8/tEEDro/82KZ7fFy1P/aoajrQDUBEQnHdSP2jCjb9UlWHVnc/xnhL9vTZFBBBwv12U/UaER6ODB3KNTPncu+8YoqKwggPdzqUqYi3pnquBrap6i4vvZ8xXqWlyvlps1nfrD9RTes5HSdwjRhBvYIjJBz7mq++cjqMqYy3in8MMKOS5y4XkbUislBEOntpf8aclQ2zNhFfso2CgcOdjhLYBg5EIyO5NtSme3yZx8UvIhHAMOCDCp5eBZyjql2BF4CPT/M+E0UkVURSMzMzPY1lzE/snzQbgA4P2KxjjapXDxkwgNERHzP7Y7V78foob3ziHwSsUtWM8k+oaraqnnQ/XgCEi0hsRW+iqlNUNVFVE+Pi4rwQy5j/r+m3H7Opbncad23tdJTAN2IEzfN2Um/nWtaudTqMqYg3in8slUzziEhzERH34x7u/R3xwj6NqbKtn+7morzvybrqOqejBIdrrkFFGCkfM3Om02FMRTwqfhGJAZKAWWXW3SEid7gXrwPSRGQN8DwwRtX+8mdq1/b/uP54tvu/ax1OEiSaNUMuv5wbYj7mgw+w6R4f5FHxq2qOqjZR1eNl1k1W1cnuxy+qamdVvUhVL1PVbzwNbMzZils+ky3RXWnep73TUYLHiBG0P7mawvTtrF/vdBhTnp25awLa9q/2c3Hu12T2tWmeWnWd65/3aD6w6R4fZMVvAtrWf7vOKYx/0Iq/VsXHQ48e3FzXit8XWfGbgNb405lsj0qg5dWdnI4SfEaPpuPJleSt38bGjU6HMWVZ8ZuAteuHDC4++QUHL7cvdR3hnu75BR/w4YcOZzE/YcVvAlb6kx8TSinn/NameRxx7rnQsyc3x9h0j6+x4jcBq+EnH7A78nxaD+ridJTgNXo0HXNWcXLNVtLTnQ5jTrHiNwFp+9cH6J79Gft6jwHXOYTGCWWO7nnnHYezmP+x4jcBafPf3iOUUtr9+QanowS3c86BXr24pe77vP22nczlK6z4TcBRhZbL32Zz3Uto1q+j03HM6NG0P7ma0G3p/PCD02EMWPGbALTuw810LUjl6KBxTkcxANdfj4aEcHPodN5+2+kwBqz4TQDa9+93KEXo+OgYp6MYgBYtkKQkJkS+xfvvllJc7HQgY8VvAkpxkdJh5dukxV1Fg4RWTscxp9x4I01zd3FB5ld88onTYYwVvwkoqZNW0K5kK0XX2TSPTxkxAq1bl9si3rTpHh9gxW8CyvFJ71BABJ0fsbN1fUpMDHLttfxCP2DhrDxycpwOFNys+E3AyD5cyMWbZrA+fih1mjd0Oo4p78YbiSrK5urcOXYJB4dZ8ZuA8d2f5tGUQ8Tcc6vTUUxFrrwSbd2aO6Lf4tVXnQ4T3Dy9A9dOEVknIqtFJLWC50VEnheRrSKyVkQu8WR/xpxO9LuvkRHWigvuGeh0FFOR0FDkhhvol7+I9C8z7BIODvLGJ/6rVLWbqiZW8NwgoL37ZyIwyQv7M+ZnNi7ZQ6/ji9jZ/xYkPMzpOKYy48cTUlrCLTKNqVOdDhO8anqqZzjwprp8BzQUkRY1vE8ThLb/5Q1CKeWCv9/idBRzOp06QZ8+3Bs9hWmvl1JU5HSg4ORp8SvwiYisFJGJFTzfCthTZnmve93PiMhEEUkVkdTMzEwPY5lgkp9byoU/TGVds6tp1L2d03HMmdx+Oy1ytnFh5qfMm+d0mODkafH3VtVuuKZ07hKRvtV9I1WdoqqJqpoYFxfnYSwTTL55/FPOLd2JTJjgdBRTFddeizZuzL1RU+xLXod4VPyqus/9+xDwEdCj3Cb7gDZlllu71xnjPa+9xrGQRiT8caTTSUxV1KmDjB/P4IKPWLUwgz17zvwS413VLn4RiRGReqceA8lAWrnN5gA3uY/uuQw4rqoHqp3WmHLSv8jgikOz2NLzRkKi6zgdx1TVxImElhZzM68zebLTYYKPJ5/4mwFficga4AdgvqouEpE7ROQO9zYLgO3AVuAV4NcepTWmnM2/nUwkhZz3zF1ORzFno2NH6NePe6NeYcrkUvLynA4UXER98M4IiYmJmpr6s9MCjPmJowcLKGxxLhmtu9N1z3yn45iz9c47cMMNpLCQ0a+lcKudd+cREVlZyWH1P2Nn7hq/9c39H9CMDOr+4TdORzHVcd11aIsW/KXu0zz3nN2dqzZZ8Ru/VFyktJ71HLuiO9LujmSn45jqiIhA7r6bXieXUrp2HZ9/7nSg4GHFb/zSl//+josKUzl642/sZur+7I470OhoHo50feo3tcOK3/il0meeIzukAV2evNHpKMYTjRsjt9zC6KK3+f7jA+zY4XSg4GDFb/zOypk76Hd4Jpv73kZog7pOxzGeuu8+QrWYe0Je4qmnnA4THKz4jd/J+O2TlBBKwpT7nY5ivOH885Hhw7knfBIzXs1h/36nAwU+K37jV9Yu2MvVe15nfY9biW5v99QNGA88QN2CLMYXvWqf+muBFb/xK7vueYoQSmn/6u+djmK8qXdv6NePv9T5J29MyuPQIacDBTYrfuM3Ni7PYMD2/5LW7UbqdYl3Oo7xtsceo2HeQW7On8wzzzgdJrBZ8Ru/seXOp4mgkLb/fcjpKKYm9OsHV13FI3We5LUXcsnKcjpQ4LLiN35hw/JDXLXpZdZ1up6GPS5wOo6pKY89RsP8DG7MmWRz/TXIit/4ha3j/0oUecRPfcTpKKYm9ekDAwbwSOST/PfpHLtkcw2x4jc+7+up6QzePZm0XhNpeFlHp+OYmvbYYzQoyOTu4mf505+cDhOYrPiNTyspgYL7HyJPoun03qNOxzG14fLLYcQI/hDyDz55cz+rVzsdKPBY8RuftuTPX9I/+2N2/OL3RLZp6nQcU1ueeooIKeI/EQ/z4IN25U5vs+I3PivnpNL8qQfJCG9Fl6l2lm5QOe885Le/ZUzhm5xY9j0LFzodKLB4cuvFNiLymYhsEJH1InJvBdtcKSLHRWS1+8e+mTNVNv+6qVxc9APHH3wciYl2Oo6pbX/4A9q8Of+tcy/33FVKbq7TgQKHJ5/4i4EHVDUBuAy4S0QSKtjuS1Xt5v6vDQdGAAANIElEQVT5qwf7M0Fk7ZKDJC1+kM3N+3LB4zc5Hcc4oV495B//oFv+9/Te+RaPPup0oMBR7eJX1QOqusr9+ASwEbCLpxiPFRfDgTH3EU0uzWZPgRCbkQxaN90EvXrxcsT9vPOfA/z4o9OBAoNX/osSkXjgYuD7Cp6+XETWishCEel8mveYKCKpIpKamZnpjVjGT825fT4Dj77HltF/okGPDk7HMU4KCYHXXydacnk1/A5+dZtSXOx0KP/ncfGLSF3gQ+A+Vc0u9/Qq4BxV7Qq8AHxc2fuo6hRVTVTVxLi4OE9jGT+1bfUJLn39TnbV60znN+1CbAbo0AF5/HFSCubQYdU7dh0fL/Co+EUkHFfpv62qs8o/r6rZqnrS/XgBEC4isZ7s0wSugnxl89V30FL3UefNV5DICKcjGV9x//3oZZcxOfwenv/DQVJTnQ7k3zw5qkeA14CNqvp0Jds0d2+HiPRw7+9IdfdpAtusYa8zKOsd0sc9RrMRvZyOY3xJaCjy+uvUDclleth4xo4u4fhxp0P5L08+8V8B3Aj0L3O45mARuUNE7nBvcx2QJiJrgOeBMap2Kob5uc9e2sDwpXeT3ro/CW8+7HQc44s6dkSef55++UsYv/Mxbr/dTuyqLvHFHk5MTNRU+7tc0Ni9KZeTF/akuWQQvXkNddq2cDqS8VWqcOut8MYbDGEewyYP4fbbnQ7lG0RkpaomVmVbO07OOOp4VgnpPW6kY8l6cie/ZaVvTk8EXn4Z7daNd8N+yX/u2s4nnzgdyv9Y8RvHFBXB0ot/R9KJWWy982laTxjodCTjD6KikA8/pG4MLAwbym0jj5CW5nQo/2LFbxyhCu/3fZHrdj/NhgG/4YKX73M6kvEn7dohsz+mHdv5sGAIvxh0kgMHnA7lP6z4Ta1ThbeHvceY7+5lY4fhJCyq8KAwY06vXz/k3Xe5pGQFzx+4jiFJhXaT9iqy4je1ShXeHvgmY+eNY0eLK+iY+jaEhjody/irESOQV14hqWQxf944jqS+BfbJvwqs+E2tKS2Ft/tNYdzSm9l2Tn/O27wQqRvjdCzj7269FZ59lpGlH/LsliEM7nOCvXudDuXbrPhNrTh5Qnm36xP88svb2XzeINpvmmulb7zn3nth2jSulOVM3XEVQ3tmsmqV06F8lxW/qXE71p3km9ajGbf+T2y8ZBwd0mYhUXWcjmUCzU03IbNnc1HEBuZm9ODuXiuZMcPpUL7Jit/UqM9f2Uzexb24OnsWmyc+RafU6UidSKdjmUA1ZAghny+nZfMSlhddzpfjXuaB3yoFBU4H8y1W/KZGHDtSwowez9Bj4kW0Yh8Zry/kgv8+4DoBx5ia1KMHoWt+JCxlAC9zF72e+QXJXQ+yYoXTwXyHFb/xKlVY+ux6trTsy9gVv2Vn+yTqbEmj5c3JTkczwaRJE0LmzYUnn2RU2FzmbOnItJ4v89D/lZBd/uLxQciK33hN6kd7mNd8Av3v78oFJRvZ/tfpdEqfTWTblk5HM8EoJAR+9ztC1q8juk93XtS7uO6pnvzqnEW8+IJSWOh0QOdY8RuPqMLX07Yyu+29dB51AcmHprNhwL3E7N1Muz/fYFM7xnkXXED48k9g+nS6NM/kveOD6P6by7n93EW88FwpJ086HbD22dU5TbUcOVjEd08sI3raJPqdmEsxYWy8eBzt336M6E7nOh3PmIoVFqJTXyf/z48TdXgvWzifaXVuJ+SWm7nujli6dnU6YPWdzdU5rfhNlR3aU8Cal7+maMZMLt31AXEcJis0lt2D76Tjc3falTWN/ygogJkzyf73ZOqv+YoiwljG1XzX8loajh9O/zFN6dLFv/7CWmvFLyIpwHNAKPCqqv6z3PPifn4wkAvcrKpnPK3Cit83HNyey+b3fiR7yXc0XPUpF2cvJ4Zc8iSK9AuGUf/2sbS9M8UOzzT+LS2N3EnTKHp/Fg0Ob6cUYS1d+SHmKnIvvZLYwT3oNqgFnTr59tVFaqX4RSQU2AwkAXuBFcBYVd1QZpvBwD24ir8n8Jyq9jzTe1vx156SEjiw8RgZ3+/k+I/bKVyzkchtG2h6eD0ditIIowSAPXXac7DbQBqNTqbdhKsIqV/X4eTGeJkqrFtH9vQ5nJz3GbHpXxNR6joBYD8tWBt6CYebJVByXgeiL+5A40viadm9BfHnhRIV5XB2aq/4ewGPqupA9/LDAKr6jzLb/BdYrqoz3MvpwJWqetrLKFnxn1lpiVKUX0JhThFFOYUUniigIDufwuP5FBzLo+hYDkXHcig4fILirGxKjx5Hj2QRcvQIYdlHiMk+SP28gzQt3k8Dfnp82/7wc8iMTSC/8yU0SOpJ/Oge1Ilv7tBIjXFIfj76wwoOL1lF9vKVRG38kdijm4nQ/384UBFh7KMVh8OacyKqGfn1m1LcsAk0bIQ0aURY4waEN6pLZGw9IhrFENEgisiGUUQ2qENk/Ugi6kYQWS+C8OhwIiLFo79RnE3xh1V/N7QC9pRZ3ovrU/2ZtmkF1Mj18zbEJBJRkndWrxEPprrKvlLKLP3ksfv9tcx6Qd3r9eePUUK0FFBCKCVESwilhBAtJYQSwij+308kytlOshyVRpwIb0x2dHOOtbqQw7EDkPhziUpoS+NL4mnR7wJaNqqHHYBpgl6dOkjfPsT17UPcqXUlJejOXRz/IZ2s1bvJTd+N7txNzJEMmmTvpP6h76m3L4sIis56d8WEsjekFa1Ldnl1GBXxpPi9SkQmAhMBzjnnnGq9x9GmHQkpqsa52V76Bkd/8j5lHsvPH7u2FdeyyE+WNSQEJOR/vwkNRUNDXcclh4VDWJhrsjEiAomMgPBwQqLrEBIVSUhUJGENYghvEE14wxiimtYjunl96rasT3TLhjQKC6WRV0ZrTBAKDUXOa0fD89rRcGwl26hCbi6FGUfJOXiCnIMnyDt0gsJjuRQdz6U4O5eS3AI0v4DS/EJK8wtdt6MrKkKjomhdC8PwpPj3AW3KLLd2rzvbbQBQ1SnAFHBN9VQn0BU7plfnZcYY4z0iEBNDRLsYItrhkx+0PDmBawXQXkTaikgEMAaYU26bOcBN4nIZcPxM8/vGGGNqVrU/8atqsYjcDSzGdTjnVFVdLyJ3uJ+fDCzAdUTPVlyHc97ieWRjjDGe8GiOX1UX4Cr3susml3mswF2e7MMYY4x32bV6jDEmyFjxG2NMkLHiN8aYIGPFb4wxQcaK3xhjgoxPXpZZRDKBszlvORY4XENxfFUwjhmCc9zBOGYIznF7MuZzVTXuzJv5aPGfLRFJrerFiQJFMI4ZgnPcwThmCM5x19aYbarHGGOCjBW/McYEmUAp/ilOB3BAMI4ZgnPcwThmCM5x18qYA2KO3xhjTNUFyid+Y4wxVeTXxS8iKSKSLiJbReQhp/PUFBFpIyKficgGEVkvIve61zcWkaUissX92xcv/e0REQkVkR9FZJ57ORjG3FBEZorIJhHZKCK9An3cInK/+892mojMEJE6gThmEZkqIodEJK3MukrHKSIPu/stXUQGeiuH3xa/+2bvLwGDgARgrIgkOJuqxhQDD6hqAnAZcJd7rA8By1S1PbDMvRxo7gU2llkOhjE/ByxS1Y7ARbjGH7DjFpFWwG+ARFW9ENdl3scQmGN+A0gpt67Ccbr/Gx8DdHa/5mV373nMb4sf6AFsVdXtqloIvAsMdzhTjVDVA6q6yv34BK4iaIVrvNPcm00DRjiTsGaISGtgCPBqmdWBPuYGQF/gNQBVLVTVYwT4uHFdIj5KRMKAaGA/AThmVf0CyCq3urJxDgfeVdUCVd2B674mPbyRw5+Lv7IbuQc0EYkHLga+B5qVuaPZQaCZQ7FqyrPA74DSMusCfcxtgUzgdfcU16siEkMAj1tV9wFPAbuBA7ju1LeEAB5zOZWNs8Y6zp+LP+iISF3gQ+A+Vc0u+5z7pjcBc4iWiAwFDqnqysq2CbQxu4UBlwCTVPViIIdyUxyBNm73nPZwXP/TawnEiMgvy24TaGOuTG2N05+Lv8o3cg8EIhKOq/TfVtVZ7tUZItLC/XwL4JBT+WrAFcAwEdmJaxqvv4hMJ7DHDK5PdXtV9Xv38kxc/yMI5HEPAHaoaqaqFgGzgMsJ7DGXVdk4a6zj/Ln4q3Kz94AgIoJrznejqj5d5qk5wHj34/HA7NrOVlNU9WFVba2q8bj+3X6qqr8kgMcMoKoHgT0i0sG96mpgA4E97t3AZSIS7f6zfjWu77ECecxlVTbOOcAYEYkUkbZAe+AHr+xRVf32B9eN3DcD24A/Op2nBsfZG9df/9YCq90/g4EmuI4C2AJ8AjR2OmsNjf9KYJ77ccCPGegGpLr/fX8MNAr0cQOPAZuANOAtIDIQxwzMwPU9RhGuv91NON04gT+6+y0dGOStHHbmrjHGBBl/nuoxxhhTDVb8xhgTZKz4jTEmyFjxG2NMkLHiN8aYIGPFb4wxQcaK3xhjgowVvzHGBJn/B5Xzx84CkY4PAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "df=np.zeros([ndata])\n", "for i in range(0,ndata):\n", " dy = ydata[i]-func(xdata[i], *guess1)\n", " df[i]=dy\n", "#pprint(df)\n", "Jac=np.zeros([ndata,nparam])\n", "for i in range(0,ndata):\n", " Jac[i,0] = dfda1(xdata[i], *guess1)\n", " Jac[i,1] = dfda2(xdata[i], *guess1)\n", " Jac[i,2] = dfda3(xdata[i], *guess1)\n", "# pprint(Jac)\n", "iJac = linalg.inv(np.dot(np.transpose(Jac),Jac))\n", "# print(iJac)\n", "Jdf = np.dot(np.transpose(Jac),df)\n", "# pprint(Jdf)\n", "guess1 = guess1 + np.dot(iJac, Jdf)\n", "pprint(guess1)\n", "plt.plot(xdata, ydata, 'b-', label='data')\n", "\n", "popt, pcov = curve_fit(func, xdata, ydata)\n", "plt.plot(xdata, func(xdata, *guess1), 'r-', label='fit')\n", "plt.show()" ] }, { "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.6.1" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "toc": { "colors": { "hover_highlight": "#DAA520", "navigate_num": "#000000", "navigate_text": "#333333", "running_highlight": "#FF0000", "selected_highlight": "#FFD700", "sidebar_border": "#EEEEEE", "wrapper_background": "#FFFFFF" }, "moveMenuLeft": true, "nav_menu": { "height": "12px", "width": "252px" }, "navigate_menu": true, "number_sections": true, "sideBar": true, "threshold": 4, "toc_cell": true, "toc_section_display": "block", "toc_window_display": true, "widenNotebook": false } }, "nbformat": 4, "nbformat_minor": 2 }