28:null 29:["$","title",null,{"children":"Taylor approximation | Mattias Villani | Observable"}] 2f:T864,f = { if (funcname == "square"){ return [function(x){return x**2}, function(x){return 2*x}, function(x){return 2}, function(x){return 0}, function(x){return 0}, function(x){return 0} ] }else if(funcname == "exp"){ return [function(x){return Math.exp(x)}, function(x){return Math.exp(x)}, function(x){return Math.exp(x)}, function(x){return Math.exp(x)}, function(x){return Math.exp(x)}, function(x){return Math.exp(x)} ] }else if(funcname == "cube"){ return [function(x){return x**3}, function(x){return 3*x**2}, function(x){return 6*x}, function(x){return 6}, function(x){return 0}, function(x){return 0} ] } else if(funcname == "log"){ return [function(x){return Math.log(x)}, function(x){return 1/x}, function(x){return -1/(x**2)}, function(x){return 2/(x**3)}, function(x){return -6/(x**4)}, function(x){return 24/(x**5)} ] } else if(funcname == "sqrt"){ return [function(x){return Math.sqrt(x)}, function(x){return (1/2)/x**(1/2)}, function(x){return -(1/4)/x**(3/2)}, function(x){return (3/8)/x**(5/2)}, function(x){return -(15/16)/x**(7/2)}, function(x){return (105/32)/x**(9/2)} ] } else if(funcname == "sin"){ return [function(x){return Math.sin(x)}, function(x){return Math.cos(x)}, function(x){return -Math.sin(x)}, function(x){return -Math.cos(x)}, function(x){return Math.sin(x)}, function(x){return Math.cos(x)} ] } else if(funcname == "cos"){ return [function(x){return Math.cos(x)}, function(x){return -Math.sin(x)}, function(x){return -Math.cos(x)}, function(x){return Math.sin(x)}, function(x){return Math.cos(x)}, function(x){return -Math.sin(x)} ] } }30:[] 2c:{"id":"0df4323495d58088","slug":null,"title":"Taylor approximation","thumbnail":"ce1868650de09803b47215707e918ebd8617272f9170eef42fa7444866dd3827","access_level":"public","create_time":"$D2025-03-13T00:12:52.418Z","update_time":"$D2025-03-13T00:13:08.803Z","view_time":null,"first_public_time":"$D2025-03-13T00:12:52.418Z","fork":{"id":"4a49ee2a3be496f7","slug":"taylor-approximation","title":"The Taylor approximation","owner":{"id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani"}},"stars":0,"starred":false,"owner":{"type":"team","id":"7592a7ab121f74a0","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":"pro_2024"},"version":426,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"56499be8cf43f0ef","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":null}],"members":[{"type":"individual","id":"56499be8cf43f0ef","login":"mattiasvillani","name":"Mattias Villani","avatar_url":"https://avatars-next.observableusercontent.com/a440aaf39a3abefb86cbcc758fd7a5ab42555ae68d6825b52ec006f385a94d1d","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"Taylor approximation","cells":[{"id":0,"value":"# Taylor approximation","pinned":false,"hidden":false,"mode":"md"},{"id":185,"value":"> *This notebook shows how a smooth function ${tex`f(x)`} can be approximated locally by polynomials tailored to ${tex`f(x)`}.*","pinned":false,"hidden":false,"mode":"md"},{"id":360,"value":"The *n*th order Taylor approximation of the function ${tex`f(x)`} around ${tex`x=x_0`} is\n${tex.block`f(x) \\approx f(x_0) + \\sum_{k=1}^n \\frac{f^{(k)}(x_0)}{k!}(x-x_0)^k,`}\nwhere \n${tex.block`f^{(k)}(x_0) = \\frac{\\mathrm{d}^k}{\\mathrm{d} x^k}f(x)\\vert_{x=x_0}`} \nis the *k*th derivative of ${tex`f(x)`} evaluated at ${tex`x=x_0`}.\n\nThe function ${tex`f(x)`} must have *n* continuous derivatives. The approximation is local around ${tex`x=x_0`} and may not be at all accurate far from that point. \n\nExperiment below with different functions ${tex`f(x)`}, approximation orders *n*, and expansions points ${tex`x_0`}.","pinned":false,"hidden":false,"mode":"md"},{"id":297,"value":"viewof funcname = Inputs.select([\"exp\", \"log\", \"square\", \"cube\", \"sqrt\", \"sin\", \"cos\"], {value: \"exp\", label: \"Function f(x)\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"viewof inputs = Inputs.form({\n order: Inputs.range([1, 5], {label: tex`\\text{Taylor order, }n`, value: 1, step: 1}),\n x0: Inputs.range(xdomain, {label: tex`\\text{Expansion point, }x_0`, value: x0init, step: 0.1})\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":143,"value":"Plot.plot({\n title: latexfunc,\n color: {\n legend: true,\n domain: plotlabels.labels,\n range: [colors[0],colors[1]]\n },\n x: {\n label: \"x\",\n axis: true,\n },\n y: {\n label: \"f(x)\",\n axis: true,\n domain: ydomain\n },\n grid: false,\n marks: [\n Plot.ruleY([ydomain[0]]),\n Plot.ruleX([0]),\n Plot.line([{x: inputs.x0, y: ydomain[0]}, {x: inputs.x0, y: f[0](inputs.x0)}], \n {x: \"x\", y: \"y\", stroke: \"gray\", strokeDasharray: \"4,4\", label: \"x\"}),\n Plot.text([{x: inputs.x0, y: ydomain[0]}], \n {x: \"x\", y: \"y\", text: [\"x₀\"], dx: 0, dy: 20, fill: \"gray\"}),\n Plot.lineY(funcdata, {x: \"x\", y: \"f\", stroke: colors[0], strokeWidth: 3}),\n Plot.lineY(funcdata, {x: \"x\", y: \"fapprox\", stroke: colors[1], strokeWidth: 3}),\n Plot.dot([{x: inputs.x0, y: f[0](inputs.x0)}], {x: \"x\", y: \"y\", stroke: colors[1], fill: colors[1]})\n \n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":31,"value":"
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n\n---\n## 🙈 Appendix 🙈\n\nHere is only code and details that you probably don't want to read.","pinned":false,"hidden":false,"mode":"md"},{"id":28,"value":"#### Imports","pinned":false,"hidden":false,"mode":"md"},{"id":96,"value":"mathjs = require('mathjs')","pinned":false,"hidden":false,"mode":"ojs"},{"id":166,"value":"import {colors} from \"@mattiasvillani/statistics-tools\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":16,"value":"#### Data","pinned":false,"hidden":false,"mode":"md"},{"id":19,"value":"latexfunc = {\n if (funcname == \"linear\"){\n return tex`f(x) = x`\n }else if(funcname == \"exp\"){\n return tex`f(x) = \\exp(x)`\n }else if(funcname == \"log\"){\n return tex`f(x) = \\log(x)`\n }else if(funcname == \"square\"){\n return tex`f(x) = x^2`\n }else if(funcname == \"cube\"){\n return tex`f(x) = x^3`\n }else if(funcname == \"sqrt\"){\n return tex`f(x) = \\sqrt{x}`\n }else if(funcname == \"cos\"){\n return tex`f(x) = \\cos(x)`\n }else if(funcname == \"tan\"){\n return tex`f(x) = \\tan(x)`\n }else if(funcname == \"gaussian\"){\n return tex`f(x) = \\frac{1}{\\sqrt{2\\pi}}\\exp(-\\frac{1}{2}x^2)`\n }else{\n return tex`f(x) = \\sin(x)`\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":21,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs"},{"id":36,"value":"xdomain = {\n if (funcname == \"log\"){\n return [0.1, 5]\n }\n if (funcname == \"exp\"){\n return [-3,3]\n }\n if (funcname == \"square\"){\n return [-2,2]\n }\n if (funcname == \"cube\"){\n return [-2,2]\n }\n if (funcname == \"sqrt\"){\n return [0.01,5]\n }\n if (funcname == \"cos\"){\n return [-2*Math.PI,2*Math.PI]\n }\n if (funcname == \"sin\"){\n return [-2*Math.PI,2*Math.PI]\n }\n else{\n return [-5,5]\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":385,"value":"x0init = {\n if (funcname == \"log\"){\n return 0.5\n }\n if (funcname == \"exp\"){\n return 0\n }\n if (funcname == \"square\"){\n return 0.5\n }\n if (funcname == \"cube\"){\n return 0.5\n }\n if (funcname == \"sqrt\"){\n return 0.5\n }\n if (funcname == \"cos\"){\n return 0.5\n }\n if (funcname == \"sin\"){\n return 0.5\n }\n else{\n return 0\n }\n}\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":145,"value":"funcdata = {\n const ngrid = 500;\n let fval = [];\n const xgrid = d3.range(xdomain[0], xdomain[1], (xdomain[1] - xdomain[0])/ngrid)\n for (let i = 0; i < ngrid; ++i) {\n let fapprox = 0;\n for (let j = 0; j < inputs.order + 1; j++){\n fapprox += (f[j](inputs.x0)/mathjs.factorial(j))*(xgrid[i] - inputs.x0)**j\n }\n fval.push({x: xgrid[i], f: f[0](xgrid[i]), fapprox: fapprox})\n }\n return fval\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":173,"value":"plotlabels = {\n let labels = [`f(x)`, `Taylor approximation of order ${inputs.order}`];\n let labelColors = [colors[0],colors[inputs.order + 1]]\n return {labels: labels, colors: labelColors}\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":202,"value":"ydomain = d3.extent(funcdata.map(d => d.f))","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W30","stdlib":"1","titling":"auto"}} 2d:"$2c:body"