The tidiness of elementary mathematics: Unterschied zwischen den Versionen
(→Applets) |
(→The awkward integral when k=-1) |
||
(10 dazwischenliegende Versionen von einem Benutzer werden nicht angezeigt) | |||
Zeile 10: | Zeile 10: | ||
== Main article == | == Main article == | ||
− | + | Exceptional cases such as this are particularly problematic when we attempt to automate knowledge, that is when we seek to devise a ''computer algebra system'' which will symbolically manipulate expressions. If we want to integrate <math>x^{k}</math> with respect to <math>x</math> we certainly have two choices | |
− | < | + | # Ignore the special case <math>k=-1</math>. |
+ | # Ask the user if <math>k=-1</math>. | ||
+ | Both of these approaches are problematic. <math>k=-1</math> is only one point out of an infinity of other values for <math>k</math>. However, in mathematics such exceptional cases occur rather more often than you might expect if you just stuck a pin in the number line, or even if you chose an integer <math>k</math> at random. Asking the user can also be problematic as they many not yet know the value of <math>k</math>, or the integral might occur deep within some other function and the user may have no idea what the question is asking! Automatically using the formula might result in an error creeping into a much longer calculation. Mistakes like this are very hard to track down and can be serious. These issues are addressed more fully by [2]. | ||
+ | |||
+ | What we would really like is a formula that works for all <math>k</math>. Here we are going to consider the comment in [2] that really, at a symbolic level, the answer to <math>\int x^k\mathrm{d}x </math> might better be written as <math> | ||
+ | \int x^k\mathrm{d}x = \frac{x^{k+1}-1}{k+1}+c. | ||
+ | </math> | ||
+ | |||
+ | First consider the (easy) case <math>k\neq -1</math>. If we fix <math>k</math> and differentiate then we obtain <math>\frac{\mathrm{d}}{\mathrm{d}x}\frac{x^{k+1}-1}{k+1} = x^k.</math> So <math>\frac{x^{k+1}-1}{k+1}=\frac{x^{k+1}}{k+1}-\frac{1}{k+1}</math> is an anti-derivative of <math>x^k</math>. Really we have chosen a slightly different constant of integration. | ||
+ | |||
+ | Whereas, for <math>k=-1</math> we originally had the form <math>\frac{x^{0}}{0}=\frac{1}{0}</math>, the new formula gives <math>\frac{x^{0}-1}{0}=\frac{0}{0}</math>. To understand what is going on we need to look at this expression for <math>k</math> close to <math>-1</math> and take a limit, i.e. <math>\lim_{k\rightarrow -1} \frac{x^{k+1}-1}{k+1}.</math> If we can make sense of this limit, then we will define the value of the expression at <math>k=-1</math> by the value of the limit, thus removing the exceptional case in a meaningful way. This sort of thing is a standard mathematician’s move. | ||
+ | |||
+ | = Euler’s formula for <math>e^x = \lim_{n\rightarrow\infty} \left( 1+\frac{x}{n}\right)^n</math> = | ||
+ | |||
+ | To consider the special case <math>k=-1</math> we shall need to make use of a special formula discovered by Euler. | ||
+ | |||
+ | <math> e^x = \lim_{n\rightarrow\infty} \left( 1+\frac{x}{n}\right)^n.</math> | ||
+ | |||
+ | If we consider the Taylor series for <math>e^x</math>, on the left hand side we have | ||
+ | |||
+ | <math>e^x = 1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\frac{x^4}{4!}+\cdots </math> | ||
+ | |||
+ | On the right, we can expand out the product using the Binomial Theorem | ||
+ | |||
+ | <math>\left( 1+\frac{x}{n}\right)^n | ||
+ | = \sum_{i=0}^n \frac{n(n-1)\cdots(n-i+1)}{i!}\left(\frac{x}{n}\right)^i | ||
+ | = \sum_{i=0}^n \frac{x^i}{i!} \frac{n(n-1)\cdots(n-i+1)}{n^i}.</math> | ||
+ | |||
+ | Consider the coefficients of <math>x^i</math>. This equals <math>\frac{1}{i!} \frac{n(n-1)\cdots(n-i+1)}{n^i} = \frac{1}{i!} \frac{n}{n}\frac{n-1}{n} \cdots \frac{n-i+1}{n}.</math> As <math>n\rightarrow \infty</math> each term here containing <math>n</math> converges to <math>1</math> giving <math>\lim_{n\rightarrow\infty} \left( 1+\frac{x}{n}\right)^n = \sum_{i=0}^\infty \frac{x^i}{i!},</math> or exactly the Taylor series for <math>e^x</math>. | ||
+ | |||
+ | You might think this argument has made somewhat cavalier use of limits. The following argument, adapted from [1], is close to Euler’s original. He was even more bold in his use of limits! | ||
+ | |||
+ | <blockquote>Euler unhesitatingly accepts the existence of both infinitely small and infinitely large numbers, and uses them to such effect that the modern reader’s own hesitation must be tinged with envy. [1] | ||
+ | </blockquote> | ||
+ | Thinking about <math>e^x</math> near <math>x=0</math>, and its gradient, we see that, for infinitely small <math>\epsilon</math> <math>e^\epsilon = 1 + \epsilon | ||
+ | </math> Let <math>x</math> be any given number then <math>N=x/\epsilon</math> is infinitely large. So <math>e^x = e^{N\epsilon} = (e^\epsilon)^N | ||
+ | = (1+\epsilon)^N | ||
+ | = \left( 1 +\frac{x}{N} \right)^N. | ||
+ | </math> Written in terms of limits, the infinitely large <math>N</math> would give us our formula. | ||
+ | |||
+ | <p>In this applet move the slider so that <math>n \rightarrow \infty</math>. | ||
<ggb_applet width="634" height="463" version="4.0" ggbBase64="UEsDBBQACAAIAJx0aD8AAAAAAAAAAAAAAAAWAAAAZ2VvZ2VicmFfamF2YXNjcmlwdC5qc0srzUsuyczPU0hPT/LP88zLLNHQVKiuBQBQSwcI1je9uRkAAAAXAAAAUEsDBBQACAAIAJx0aD8AAAAAAAAAAAAAAAAMAAAAZ2VvZ2VicmEueG1sxVjrbts2FP7dPsWBUAwJmtiibrY7u8VaoFiBtCuQbii2rgAl0TYbmVRFypde/u199k57kh2SkiynaZauw2ok4e3wXL5z4YmnD7arAtasUlyKmUcGvgdMZDLnYjHzaj0/HXsP7t+eLphcsLSiMJfViuqZFxlKnuMkz/NxnmanGQ2T0ygYZ6eThJHTbD4nScYiGqfEA9gqfk/IZ3TFVEkzdp4t2YqeyYxqK3ipdXlvONxsNoNW1EBWi+FikQ62KvcA1RRq5jWTe8ju4NImtOSB75Phy6dnjv0pF0pTkTEPjAk1v3/71nTDRS43sOG5Xs68JEYzlowvlmhT4uNiaIhKBKRkmeZrpvBqb2lt1qvSs2RUmPNbbgZFZ44HOV/znFUzzx+E/iiZBN3fkQey4kzohpY0Moctt+mas41ja2YOZQ+0lEVKDUf48AECP/DhxAzEDQEOSeKOfLfnh24I3BC5IXY0kbseOdLI0USOJgo9WHPF04LNvDktFCLIxbxC73VrpXcFs/o0G3vryQnapPg7JA4NpA5y3Pf9E/OLOJ9ELdY9I0lPqq7qLxTaikzC6OYig68yNGxlBleZGcSfMTO5Rqiz+yZ2krgnE0XZH/v7icTwOjMvS3TrrxOYRP+LidNhmyrTJjtALQ1tEz2arZTJl3AC8cSEPYEYcyMZYZTHQCY4jALAbAASQxTjkowhMeMIwhEeRBDCGAwdCcEmRzzGP9HIMksgRmZmd4Q5CQQFRRCHQGxORYCZBDYvMUeDECniGGK8ZMSTwLAIE4gSXIVjiFBHk5IjgoQhXsQ1ig8gJBCay2QEQQKJ4Ucik+rJ2KiOLANIfEiIYYhZjRntshnpxxAaa5IGLi7KWh9AlK3ydqpl2fkCqbEe7aueq08HRfHWtKApK/CdODeeBFjTwmSEFTSXQkPrxMDtLSpaLnmmzpnWeEvBG7qmZ1Sz7WOkVq1sS5tJoZ5XUj+SRb0SCiCThd/pLAvSmwed1rgIewdR/yDuHSS9+ehKuRJPoFYM5ctKteQ0z58Yin1pQCR/EsXuYcXoRSn5oRnToX1ypqzOCp5zKn7BYDVSDC7QvUCmXLUvUBTGrSKyys93CiMYtr+ySiK0sX1zd24VTCZmpTJqUiwMB2H/Q5Dn7nNnTgJbd46gW9bZuKhM/jb2mcUT9VAW+y1r5SNa6rqyLQKWwMro/oNYFMxGgs1ffH+zi1Ruz10IhI7Xi12JK98pkC4suoAVIIhjJGjG1I2WxmjWUfmWxrcUfhtTPO/OySSwFHZM3WipMEidao2lpDWT+K0Yrmzd8r2D7LARbl7zWnB91i40zy72lhr6Z/UqZV2cHLIk/xHL6fBSHE0vWCVY0YQterKWtXJZ2IvonGV8hUt30ABCjbN+RgXcbs4WFWv1Lmzz5eCyp34/Ij/ZtqweV3L1RKxfYCRcUmA6bLWcqqzipQk4SLHUX7B9TOVcUXwp8v49k2doemZeBIRHG2gwA2u9lJXtr7Bw4GjSq2Ar7KZA2+AS9YpVPOuAFrZRQ6XqRu940GhuUAaZvsGKdsk5ey/i8WfCD2hRLqlp75p8KuiOVQfQWG5PZd4IbuhUYfpCWHF86LC5hhXd2kAEmiqsdhpbY3SG2LfGTrWmWmCjYdIerwSRne2waNhKMOfbHqSIEn+HUUEPzNkngsZKfIHtprItkG7y0k5+5HnORKcvFRg/1gtYjEoXnyVjLrK7iyWab+tBz/ONY4yLtmWFsgyTBuI5/puxNe/N0fYYZnBE4C5sYQji+PWROHYP0KFj57WwweDtWVzvxV4S3MSN/g3d6P9bKK8DBGbeokFk4RD5648/X+PMuwEQi28KxOX4MAlxswj5JxDJl6KIdrOtJg2Q372tpf7+Drwq2FwfAZC7r7ABzd5vP74XH+GVfbWOXwvE+ioSdxuj8uiMvmAvfxO/H+PC7e5vX0cFd9x4VTAbPb1Dpb+lD7my6ru9fbmwXZzCcjrvmnb7nPtN997dx/+yK/3cdAZgSpM/iMd+/0NsncLKi23JO/ctwyVXDvv13rZOzbcA9/8GUEsHCPcAZm4ABgAAohAAAFBLAQIUABQACAAIAJx0aD/WN725GQAAABcAAAAWAAAAAAAAAAAAAAAAAAAAAABnZW9nZWJyYV9qYXZhc2NyaXB0LmpzUEsBAhQAFAAIAAgAnHRoP/cAZm4ABgAAohAAAAwAAAAAAAAAAAAAAAAAXQAAAGdlb2dlYnJhLnhtbFBLBQYAAAAAAgACAH4AAACXBgAAAAA=" framePossible = "false" showResetIcon = "true" showAnimationButton = "true" enableRightClick = "false" errorDialogsActive = "true" enableLabelDrags = "false" showMenuBar = "false" showToolBar = "false" showToolBarHelp = "false" showAlgebraInput = "false" useBrowserForJS = "true" allowRescaling = "true" /> | <ggb_applet width="634" height="463" version="4.0" ggbBase64="UEsDBBQACAAIAJx0aD8AAAAAAAAAAAAAAAAWAAAAZ2VvZ2VicmFfamF2YXNjcmlwdC5qc0srzUsuyczPU0hPT/LP88zLLNHQVKiuBQBQSwcI1je9uRkAAAAXAAAAUEsDBBQACAAIAJx0aD8AAAAAAAAAAAAAAAAMAAAAZ2VvZ2VicmEueG1sxVjrbts2FP7dPsWBUAwJmtiibrY7u8VaoFiBtCuQbii2rgAl0TYbmVRFypde/u199k57kh2SkiynaZauw2ok4e3wXL5z4YmnD7arAtasUlyKmUcGvgdMZDLnYjHzaj0/HXsP7t+eLphcsLSiMJfViuqZFxlKnuMkz/NxnmanGQ2T0ygYZ6eThJHTbD4nScYiGqfEA9gqfk/IZ3TFVEkzdp4t2YqeyYxqK3ipdXlvONxsNoNW1EBWi+FikQ62KvcA1RRq5jWTe8ju4NImtOSB75Phy6dnjv0pF0pTkTEPjAk1v3/71nTDRS43sOG5Xs68JEYzlowvlmhT4uNiaIhKBKRkmeZrpvBqb2lt1qvSs2RUmPNbbgZFZ44HOV/znFUzzx+E/iiZBN3fkQey4kzohpY0Moctt+mas41ja2YOZQ+0lEVKDUf48AECP/DhxAzEDQEOSeKOfLfnh24I3BC5IXY0kbseOdLI0USOJgo9WHPF04LNvDktFCLIxbxC73VrpXcFs/o0G3vryQnapPg7JA4NpA5y3Pf9E/OLOJ9ELdY9I0lPqq7qLxTaikzC6OYig68yNGxlBleZGcSfMTO5Rqiz+yZ2krgnE0XZH/v7icTwOjMvS3TrrxOYRP+LidNhmyrTJjtALQ1tEz2arZTJl3AC8cSEPYEYcyMZYZTHQCY4jALAbAASQxTjkowhMeMIwhEeRBDCGAwdCcEmRzzGP9HIMksgRmZmd4Q5CQQFRRCHQGxORYCZBDYvMUeDECniGGK8ZMSTwLAIE4gSXIVjiFBHk5IjgoQhXsQ1ig8gJBCay2QEQQKJ4Ucik+rJ2KiOLANIfEiIYYhZjRntshnpxxAaa5IGLi7KWh9AlK3ydqpl2fkCqbEe7aueq08HRfHWtKApK/CdODeeBFjTwmSEFTSXQkPrxMDtLSpaLnmmzpnWeEvBG7qmZ1Sz7WOkVq1sS5tJoZ5XUj+SRb0SCiCThd/pLAvSmwed1rgIewdR/yDuHSS9+ehKuRJPoFYM5ctKteQ0z58Yin1pQCR/EsXuYcXoRSn5oRnToX1ypqzOCp5zKn7BYDVSDC7QvUCmXLUvUBTGrSKyys93CiMYtr+ySiK0sX1zd24VTCZmpTJqUiwMB2H/Q5Dn7nNnTgJbd46gW9bZuKhM/jb2mcUT9VAW+y1r5SNa6rqyLQKWwMro/oNYFMxGgs1ffH+zi1Ruz10IhI7Xi12JK98pkC4suoAVIIhjJGjG1I2WxmjWUfmWxrcUfhtTPO/OySSwFHZM3WipMEidao2lpDWT+K0Yrmzd8r2D7LARbl7zWnB91i40zy72lhr6Z/UqZV2cHLIk/xHL6fBSHE0vWCVY0YQterKWtXJZ2IvonGV8hUt30ABCjbN+RgXcbs4WFWv1Lmzz5eCyp34/Ij/ZtqweV3L1RKxfYCRcUmA6bLWcqqzipQk4SLHUX7B9TOVcUXwp8v49k2doemZeBIRHG2gwA2u9lJXtr7Bw4GjSq2Ar7KZA2+AS9YpVPOuAFrZRQ6XqRu940GhuUAaZvsGKdsk5ey/i8WfCD2hRLqlp75p8KuiOVQfQWG5PZd4IbuhUYfpCWHF86LC5hhXd2kAEmiqsdhpbY3SG2LfGTrWmWmCjYdIerwSRne2waNhKMOfbHqSIEn+HUUEPzNkngsZKfIHtprItkG7y0k5+5HnORKcvFRg/1gtYjEoXnyVjLrK7iyWab+tBz/ONY4yLtmWFsgyTBuI5/puxNe/N0fYYZnBE4C5sYQji+PWROHYP0KFj57WwweDtWVzvxV4S3MSN/g3d6P9bKK8DBGbeokFk4RD5648/X+PMuwEQi28KxOX4MAlxswj5JxDJl6KIdrOtJg2Q372tpf7+Drwq2FwfAZC7r7ABzd5vP74XH+GVfbWOXwvE+ioSdxuj8uiMvmAvfxO/H+PC7e5vX0cFd9x4VTAbPb1Dpb+lD7my6ru9fbmwXZzCcjrvmnb7nPtN997dx/+yK/3cdAZgSpM/iMd+/0NsncLKi23JO/ctwyVXDvv13rZOzbcA9/8GUEsHCPcAZm4ABgAAohAAAFBLAQIUABQACAAIAJx0aD/WN725GQAAABcAAAAWAAAAAAAAAAAAAAAAAAAAAABnZW9nZWJyYV9qYXZhc2NyaXB0LmpzUEsBAhQAFAAIAAgAnHRoP/cAZm4ABgAAohAAAAwAAAAAAAAAAAAAAAAAXQAAAGdlb2dlYnJhLnhtbFBLBQYAAAAAAgACAH4AAACXBgAAAAA=" framePossible = "false" showResetIcon = "true" showAnimationButton = "true" enableRightClick = "false" errorDialogsActive = "true" enableLabelDrags = "false" showMenuBar = "false" showToolBar = "false" showToolBarHelp = "false" showAlgebraInput = "false" useBrowserForJS = "true" allowRescaling = "true" /> | ||
+ | |||
+ | = The awkward integral when <math>k=-1</math> = | ||
+ | |||
+ | We shall make use of Euler’s limit to evaluate the awkward integral when <math>k=-1</math>. | ||
+ | |||
+ | Let <math> x = \left( 1+\frac{y}{n}\right)^n.</math> | ||
+ | |||
+ | Rearranging this we have <math>x^{\frac{1}{n}} = 1 + \frac{y}{n}</math> and then <math>y = n \left(x^{\frac{1}{n}} -1\right).</math> | ||
+ | |||
+ | Let <math>m=\frac{1}{n}</math> so <math>y = \frac{1}{m} \left(x^{m}-1\right)=\frac{x^m-1}{m}.</math> | ||
+ | |||
+ | Now, set <math>m=k+1</math>, to give <math>y =\frac{x^{k+1}-1}{k+1}</math> and we are coming somewhat closer to the right hand side of the formula for our integral. | ||
+ | |||
+ | If <math>k\rightarrow -1</math>, then <math>m\rightarrow 0</math> and <math>n=\frac{1}{m}\rightarrow \infty</math>. But if <math>n\rightarrow \infty</math> then we may apply to and conclude that <math>x=e^y</math>. Taking the inverse gives <math>y=\ln(x)</math>. Hence <math>\lim_{k\rightarrow -1} \frac{x^{k+1}-1}{k+1} = \ln(x)</math> and so is indeed valid for all <math>k</math>, including the special case <math>k=-1</math>, when we define the value of in this case by the limit <math>\lim_{k\rightarrow -1}</math>. This limit happens to give the right answer to the integral for <math>k=-1</math>. | ||
+ | |||
+ | We can see this limit in the applet below. Mve the slider so that <math>k \rightarrow -1</math>. | ||
+ | |||
+ | <ggb_applet width="634" height="463" version="4.0" ggbBase64="UEsDBBQACAAIAPBzaD8AAAAAAAAAAAAAAAAWAAAAZ2VvZ2VicmFfamF2YXNjcmlwdC5qc0srzUsuyczPU0hPT/LP88zLLNHQVKiuBQBQSwcI1je9uRkAAAAXAAAAUEsDBBQACAAIAPBzaD8AAAAAAAAAAAAAAAAMAAAAZ2VvZ2VicmEueG1szVh7j9M4EP8bPsUoQicQ2zbOs+Va0IGEDmnhkJY7odMJyU3c1DSJQ+y0XR7f/cZ2kqa7sCwsOm61W7/G8/yNZ7rzR/sihy2rJRflwiFj1wFWJiLlZbZwGrUaTZ1HD2/PMyYytqwprERdULVwAk3J04XDYhLT6YqMZilZjoJwSUeUpNORN40jL8WZF64cgL3kD0rxghZMVjRhZ8maFfRUJFQZwWulqgeTyW63G3eixqLOJlm2HO9l6gCqWcqF004eILujSzvfkHuuSyavn59a9iNeSkXLhDmgTWj4w9u35jtepmIHO56q9cKJQjRjzXi2RpsiFxcTTVShQyqWKL5lEq8OlsZmVVSOIaOlPr9lZ5D35jiQ8i1PWb1w3LHvxtHM6z9jB0TNWalaWtLKnHTc5lvOdpatnhmJgQNKiHxJNUf4+BE813PhRA/EDh4OUWSPXLvn+nbw7BDYIbQ0gb0eWNLA0gSWJvAd2HLJlzlbOCuaS/QgL1c1Rq9fS3WeM6NPu3GwnpygTZK/R2Jfu9S6HPdd90T/oZ9Pgs7XAyPJQKqqm28U2omM/OD6Ir0bGep3Mr3PmemFXzAzukKotfs6dpJwIBNFmV/zd0mif5WZFyXa9c0ERsF/YuJ80qXKvM0OkGtN26JHsULqfPFnEM407AmEmBtRjCgPgcxwiD3AbAASQhDikkwh0mMMfowHAfgwBU1HfDDJEU7xI4gNswhCZKZ3Y8xJICgogNAHYnIqAMwkMHmJOer5SBGGEOIlLZ54moUfQRDhyp9CgDrqlIwJEvp4Edco3gOfgK8vkxi8CCLNjwQ61aOpVh1ZehC5EBHNELMaM9pmM9JPwdfWRK27eFk16shFSZF2UyWqPhZIje/R4dWz79PRo3hrntMly7FOnOlIAmxprjPCCFqJUkEXRM/uZTWt1jyRZ0wpvCXhLd3SU6rY/ilSy062oU1EKV/WQj0ReVOUEiARudvrLHIymHu91rjwBwfB8CAcHESDefxZuQJPoJEM5YtaduQ0TZ9pisPTgJ78o8zPH9eMbirBj82YT0zJmbMmyXnKafkXglVL0X6BvgLp56qrQIEfdoqIOj07l4hg2P/NaqH9GI5nRz8OnNsjbxaP3eEPBlwmVCdfMDu+NMVL7VHoHl8iVjTb9hGie9Ybn9U6sVvD9eKZfCzyw5Yx/wmtVFOb3gHfxlob9VuZ5cxAxCQ2FuZksxT7M4sN3/J6dV7hyrUKLDPjdsCnwQtDJGjHpR0Njdasp3INjWso3A5sPO3PycwzFGZc2tFQIXqtaq2lpDOTuJ0YLs2D5jpHaWOgr8t8U3J12i0UTzYHSzX9i6ZYsh5AxyzJD2I5n1wA2HzD6pLlLZ4xko1opE3PAdRTlvACl/agdQjVwfoTFbC7Kctq1umdm67MusucukOoXto2rJ7WonhWbl8hEi4oMJ90Ws5lUvNKAw6WWAM27ICplEuKJSQd3tMJiKYnulSge5R2DaZmo9aiNo0Xvig46rzLWYFtFigDrrIpWM2T3tEb08GhUk2n9zi2mmsvg1i+xafuQnAOUcTjL8APaF6tqWbX5lNOz1l95BrD7blIW8Etncx1wwgFxwo4QrgXdI85ivyWEl9BhS0zxqI8tMxWs/YVwQZEN+R7XTNjPTvH1Dfd/IrvBx5FJ/H3CAp6ZM0hDxS+0BtsQ6VpjVSblmbyO09TVvbq0hLhY4KAj1Sl7XU1QivGLLb7uxU6wLwIg9i3odFB2lc1itN8WifjtwXcxPHu/h4sYP/m7gbuA7kHE2hntj4dh3fVlAYSzoHN1bEcpMIwmN2Lczmc7jXD6X6vT692C9ZJRVrX/PKuEerXO/APdlTJh/2bD5v75NMn8wl37KFz2UWag3PM7ue6iMtT+oq9tnsHNJrmQWKyrvpe0RQLt20a+/v45a5WL3XdAY18/KYVXyqQ+pvq9EKZg/f2m+5XfZ61/s4sFO8esDj6dkBm3+ft4etiPP9/AqR3FSBH5Nsh6f1kJ/1gSJJx4B2BL+ggeV0QwsJZt05eWxTmQsPxM868hLj1zZ15g+y+XCCuWx++hlXyNbBOhk2CacTb/yk9/BdQSwcINd+bIxoGAADwEgAAUEsBAhQAFAAIAAgA8HNoP9Y3vbkZAAAAFwAAABYAAAAAAAAAAAAAAAAAAAAAAGdlb2dlYnJhX2phdmFzY3JpcHQuanNQSwECFAAUAAgACADwc2g/Nd+bIxoGAADwEgAADAAAAAAAAAAAAAAAAABdAAAAZ2VvZ2VicmEueG1sUEsFBgAAAAACAAIAfgAAALEGAAAAAA==" framePossible = "false" showResetIcon = "true" showAnimationButton = "true" enableRightClick = "false" errorDialogsActive = "true" enableLabelDrags = "false" showMenuBar = "false" showToolBar = "false" showToolBarHelp = "false" showAlgebraInput = "false" useBrowserForJS = "true" allowRescaling = "true" /> | ||
+ | |||
+ | = The mathematical moral = | ||
+ | |||
+ | They key to all of this is to choose the constant of integration in so that the resulting formula also holds in the awkward limiting case, <math>k=-1</math>. Like many tricks in mathematics it is a technique which occasionally has useful applications elsewhere. | ||
+ | |||
+ | There is a deep consistency to elementary mathematics, sometimes we just have to search to find it. | ||
+ | |||
+ | = References = | ||
+ | |||
+ | [1] C. H. Edwards. The Historical Development of the Calculus. Springer-Verlag, 1979. | ||
+ | |||
+ | [2] D. R. Stoutemyer. Crimes and misdemeanors in the computer algebra trade. Notices of the American Mathematical Society, 38(7):778-785, September 1991. |
Aktuelle Version vom 11. November 2011, 00:04 Uhr
We all expect a certain tidiness to mathematics, and for many this provides a gratifying aesthetic pleasure. However, there are some wrinkles. One of these arises in the integral of . Calculus provides us with the formula However, this equation is only correct for . If we try to take the right hand side is meaningless because we have a zero on the denominator of the fraction, i.e. . But in this case a separate argument gives the answer
Our expectation is that these two formulae should be reconciled. Indeed, if we let approach in the first we should end up with the second, but that fails to happen. For each , is undefined.
This vingette explores this issue.
Inhaltsverzeichnis |
Main article
Exceptional cases such as this are particularly problematic when we attempt to automate knowledge, that is when we seek to devise a computer algebra system which will symbolically manipulate expressions. If we want to integrate with respect to we certainly have two choices
- Ignore the special case .
- Ask the user if .
Both of these approaches are problematic. is only one point out of an infinity of other values for . However, in mathematics such exceptional cases occur rather more often than you might expect if you just stuck a pin in the number line, or even if you chose an integer at random. Asking the user can also be problematic as they many not yet know the value of , or the integral might occur deep within some other function and the user may have no idea what the question is asking! Automatically using the formula might result in an error creeping into a much longer calculation. Mistakes like this are very hard to track down and can be serious. These issues are addressed more fully by [2].
What we would really like is a formula that works for all . Here we are going to consider the comment in [2] that really, at a symbolic level, the answer to might better be written as
First consider the (easy) case . If we fix and differentiate then we obtain So is an anti-derivative of . Really we have chosen a slightly different constant of integration.
Whereas, for we originally had the form , the new formula gives . To understand what is going on we need to look at this expression for close to and take a limit, i.e. If we can make sense of this limit, then we will define the value of the expression at by the value of the limit, thus removing the exceptional case in a meaningful way. This sort of thing is a standard mathematician’s move.
Euler’s formula for
To consider the special case we shall need to make use of a special formula discovered by Euler.
If we consider the Taylor series for , on the left hand side we have
On the right, we can expand out the product using the Binomial Theorem
Consider the coefficients of . This equals As each term here containing converges to giving or exactly the Taylor series for .
You might think this argument has made somewhat cavalier use of limits. The following argument, adapted from [1], is close to Euler’s original. He was even more bold in his use of limits!
Euler unhesitatingly accepts the existence of both infinitely small and infinitely large numbers, and uses them to such effect that the modern reader’s own hesitation must be tinged with envy. [1]
Thinking about near , and its gradient, we see that, for infinitely small Let be any given number then is infinitely large. So Written in terms of limits, the infinitely large would give us our formula.
In this applet move the slider so that .
The awkward integral when
We shall make use of Euler’s limit to evaluate the awkward integral when .
Let
Rearranging this we have and then
Let so
Now, set , to give and we are coming somewhat closer to the right hand side of the formula for our integral.
If , then and . But if then we may apply to and conclude that . Taking the inverse gives . Hence and so is indeed valid for all , including the special case , when we define the value of in this case by the limit . This limit happens to give the right answer to the integral for .
We can see this limit in the applet below. Mve the slider so that .
The mathematical moral
They key to all of this is to choose the constant of integration in so that the resulting formula also holds in the awkward limiting case, . Like many tricks in mathematics it is a technique which occasionally has useful applications elsewhere.
There is a deep consistency to elementary mathematics, sometimes we just have to search to find it.
References
[1] C. H. Edwards. The Historical Development of the Calculus. Springer-Verlag, 1979.
[2] D. R. Stoutemyer. Crimes and misdemeanors in the computer algebra trade. Notices of the American Mathematical Society, 38(7):778-785, September 1991.