A method for extending the accuracy of certain mathematical tables. (Q1513630)

Eine Tafel (working-table) enthalte die bis auf 7 Decimalstellen berechneten Werte einer Function \(u\) von \(x\) für Werte von \(x\) von 0,00 bis 0,50, deren Intervalle gleich 0,01 sind; eine andere (checking table) enthalte die bis auf 11 Stellen berechneten Werte derselben Function für Werte, deren Intervalle grösser, nämlich 0,1 sind. Es soll eine neue Tafel hergestellt werden, welche soviel wie möglich die Vorzüge beider vereinigt, die grössere Genauigkeit der zweiten mit den kleineren Intervallen der ersten. Dabei dient die working-table als Basis zur Berechnung der ersten und zweiten Differenzen für die neue Tafel; die Werte dieser werden gefunden durch successive Addition der gefundenen Differenzen und controllirt mittels der checking-table. Ist \(u=f(x)=e^x\), \(\frac1{\sqrt{2\pi}}e^{-\frac12x^2}\), \(\sin\frac12\pi x\), so wird jeder Wert aus einem oder zwei vorhergehenden Werten berechnet. Das ist nicht immer möglich. Aber es giebt immer eine andere Formel für die Berechnung, wenn (vorausgesetzt, dass das Intervall \(h\) klein genug ist) der Differentialquotient der Function durch die Function oder durch die Function und das Argument ausgedrückt werden kann. Es sei \(u=f(x)\) und \(v=\frac{du}{dx}= \varphi(x,u)\), wo \(\varphi(x,u)\) eine Function bedeutet, deren Wert für einen gegebenen Wert von \(x\) berechnet werden kann, wenn \(u\) für diesen Wert bekannt ist. Setzt man \(u_0=f(x_0)\), \(u_{\pm n}=f(x_0\pm nh)\), wo \(x_0\) ein in der Tafel vorkommender Wert ist, so ist \[ \begin{aligned} u_1 &= f(x_0+h) = u_0 + hf'(x_0) + \frac{h^2}{2!}f''(x_0) + \frac{h^3}{3!}f'''(x_0) +\cdots,\\ hv_0 &= hf'(x_0),\quad hv_{-1} = hf'(x_0-h) = hf'(x_0) - h^2f''(x_0) + \frac{h^3}{2!}f'''(x_0) -\cdots,\\ hv_{-2} &= hf'(x_0) - 2h^2f''(x_0) + \frac{4h^3}{2!}f'''(x_0) -\cdots,\end{aligned} \] daher \[ \varDelta hv_{-1} = h^2f''(x_0) - \frac12h^3f'''(x_0) + \frac16h^4f^{\text{IV}}(x_0) -\cdots, \] \[ \varDelta^2 hv_{-2} = h^3f'''(x_0) - h^4f^{\text{IV}}(x_0) +\cdots,\quad \varDelta^3 hv_{-32} = h^4f^{\text{IV}}(x_0) -\cdots; \] demnach \[ u_1 = u_0 + hv_0 + \frac12\varDelta hv_{-1} + \frac5{12}\varDelta^2 hv_{-3} + \frac{251}{720}\varDelta^4 hv_{-4} + \frac{95}{288}\varDelta^5 hv_{-5} +\cdots. \] \(hv_0\), \(\varDelta hv_{-1}\), \(\varDelta hv_{-2}\), ... sind bekannt, wenn \(u_0\), \(u_{-1}\), \(u_{-2}\), ... bekannt sind; demnach kann jeder Wert von \(u\) aus dem ihm vorhergehenden Werte nach dieser Formel berechnet werden. II. Für seine Methode der Centraldifferenzen führt Verf. drei Operatoren \(\delta\), \(\sigma\) und \(\mu\) ein durch die Formeln: \[ \delta f(x) = f(x+\frac12h) - f(x-\frac12h),\quad \mu f(x) = \frac12\{f(x+\frac12h) + f(x-\frac12h)\}, \] \[ \delta\sigma f(x) = f(x); \] dann ergiebt sich \[ \begin{aligned} \sigma u_{-\frac12} &= \cdots+ u_{-2} + u_{-1},\,\sigma u_\frac12 = \cdots+ u_{-2} + u_{-1} + u_0,\,\sigma u_\frac32 = \cdots+ u_{-2} + u_{-1} + u_0 + u_1,\\ .\quad.\quad&.\quad.\quad.\quad.\quad.\quad.\quad.\quad.\quad.\quad.\\ \mu\sigma u_{-1} &= \cdots+ u_{-2} + \frac12u_{-1},\,\mu\sigma u_0 = \cdots+ u_{-2} + u_{-1} + \frac12u_0,\\ \mu\sigma u_1 &= \cdots+ u_{-2} + u_{-1} + u_0 + \frac12i_1.\end{aligned} \] Die Operation \(\sigma\) involvirt die Einführung einer willkürlichen Constante; ist diese passend gewählt, so ist \[ \sigma\delta u_0 = u_0,\quad \sigma\delta^2u_\frac12 = \delta u_\frac12,\quad \sigma\delta^3u_0 = \delta^2u_0,\dots, \] daher \(\sigma=\delta^{-1}\). Für \(u\) ergeben sich die Formeln: \[ u = \mu\left(\sigma - \frac12\delta + \frac{11}{720}\delta^3 - \frac{191}{60480}\delta^5 + \frac{2497}{3628800}\delta^7 -\cdots\right)hv,\tag{1} \] \[ u = \left(\sigma + \frac1{24}\delta - \frac{17}{5760}\delta^3 + \frac{367}{967680}\delta^5 - \frac{27859}{464486400}\delta^7 -\cdots\right)hv.\tag{2} \] Man berechnet \(hv_0\), \(hv_1\), ... oder \(hv_\frac12\), \(hv_\frac32\), ... und ihre Differenzen, dann \(\sigma hv_\frac12\) oder \(\sigma hv_0\) mittels des durch die checking-table gegebenen Wertes von \(u_0\) nach der Formel \[ \sigma hv_\frac12 = u_0 + \frac12hv_0 + \frac12\mu\delta hv_0 - \frac{11}{720}\mu\delta^3hv_0 +\cdots \] oder \[ \sigma hv_0 = u_0 - \frac1{24}\delta hv_0 + \frac{17}{5760}\delta^3hv_0 -\cdots, \] in ähnlicher Weise \(\sigma hv_{n+\frac12}\), \(\sigma hv_{2n+\frac12}\), ... oder \(\sigma hv_n\), \(\sigma hv_{2n}\), ... Wenn die Bedingungen \[ \begin{aligned} \sigma hv_{n+\frac12} - \sigma hv_\frac12 &= hv_1 + hv_2 +\cdots+ hv_n,\\ \sigma hv_{2n+\frac12} - \sigma hv_0 &= hv_{n+1} + hv_{n+2} +\cdots+ hv_{2n}\end{aligned} \] oder \[ \begin{aligned} \sigma hv_n - \sigma hv_0 &= hv_\frac12 + hv_\frac32 +\cdots+ hv_{n-\frac12},\\ \sigma hv_{2n} - \sigma hv_n &= hv_{n+\frac12} + hv_\frac32 +\cdots+ hv_{2n-\frac12}\end{aligned} \] erfüllt sind, so wird \(\sigma hv\) und dann \(u\) nach einer der obigen Formeln berechnet. III. Bezeichnet \(Du\) die Ableitung von \(u\), so ergiebt sich \[ u = \left(\sigma^2 + \frac1{12} - \frac1{240}\delta^2 + \frac{31}{60480}\delta^4 - \frac{289}{3628800}\delta^6 +\cdots\right)h^2D^2u; \] mit dieser Formel kann die Tabelle von \(u\) weiter ausgedehnt werden, wenn \(u\) einer Gleichung von der Form \[ d^2u/dx^2 = \psi(x,u) \] genügt. Man hat zunächst die Tafel der Werte von \(h^2w = h^2d^2u/dx^2\) nach dieser Gleichung unter Benutzung der aus der working-table bekannten Werte von \(u\), dann die Tafel der Werte von \(\sigma h^2w\), \(\sigma^2h^2w\) und dann \(u\) nach der Formel \[ u = \left(\sigma^2 + \frac1{12} - \frac1{240}\delta^2 + \frac{31}{60480}\delta^4 - \frac{289}{3628800}\delta^6 +\cdots\right)h^2w \] zu berechnen.