自然数の自然対数ln nを積分で近似したい

はじめに

この記事は 🐟️🍊みかん🍊🐟️様のこの記事 を参考にしています。

直接的には(多分)無理

この記事 にもある通り、例えば$\ln 5$を近似しようとして、
$${\displaystyle \frac{1}{5}}{\int }_0^4{x}^n(1-x{)}^{n}dx<{\int }_0^4{\displaystyle \frac{x^n(1-x{)}^n}{1+x}}<{\int }_0^4{x}^n(1-x{)}^{n}dx$$
と不等式を作ったとしても、$n=2,4,8$で順に、(大体)
$$\begin{array}{rl}0.9166& <\ln 5<20.5834\\ -6.6131& <\ln 5<387.721\\ -2275.5& <\ln 5<267854.5\\ \end{array}$$
みたいに全然うまくいかないのです。直接それだけを求めるというのは結構無理そうなので少しだけ変えた手法で近似しようと思います。

早速本題へ

とはいえ仕組みは結構単純なもので、${\displaystyle {\int }_0^1{\displaystyle \frac{x^n(1-x{)}^n}{1+x}}dx}$の分子は変えずに分母をいじるだけです。

例えば$\ln 3$を求めたいなら
$${\displaystyle \frac{1}{3}}{\int }_0^1{x}^n(1-x{)}^{n}dx<{\int }_0^1{\displaystyle \frac{x^n(1-x{)}^n}{2+x}}dx<{\displaystyle \frac{1}{2}}{\int }_0^1{x}^n(1-x{)}^{n}dx$$
とすればいいだけです。

懸念点

この式を見た瞬間に分かる人もいるとは思いますが、この積分、$\ln 2$と$\ln 3$が両方とも出るので$\ln 2$の近似が必要ですし、なんならその誤差をそのまま持ち越します。

何がともあれ実際にやってみよう

まず$\ln 2$の近似が必要なので
$${\displaystyle \frac{1}{2}}{\int }_0^1{x}^n(1-x{)}^{n}dx<{\int }_0^1{\displaystyle \frac{x^n(1-x{)}^n}{1+x}}dx<{\int }_0^1{x}^n(1-x{)}^{n}dx$$
で近似します。$\ln 3$等々を近似するときに必要になるので、ここは精度よく$n=8$にしましょう。(ちなみに計算結果は 参考の記事 にもあります。)
$${\displaystyle \frac{1}{2}}{\int }_0^1{x}^8(1-x{)}^{8}dx={\displaystyle \frac{8!8!}{2\times 17!}}={\displaystyle \frac{1}{437580}}(こっそりベータ関数の公式を使用しています。)$$であることと、
$$\begin{array}{rcl}& & {\int }_0^1{\displaystyle \frac{x^8(1-x{)}^8}{1+x}}dx={\int }_0^1\{x^{15}-9x^{14}+37x^{13}-93x^{12}+163x^{11}-219x^{10}+247x^9-255x^8\\ & & +256(x^7-x^6+x^5-x^4+x^3-x^2+x-1)+{\displaystyle \frac{256}{1+x}}\}dx\\ & & =[{\displaystyle \frac{x^{16}}{16}}-{\displaystyle \frac{3x^{15}}{5}}+{\displaystyle \frac{37x^{14}}{14}}-{\displaystyle \frac{93x^{13}}{13}}+{\displaystyle \frac{163x^{12}}{12}}-{\displaystyle \frac{219x^{11}}{11}}+{\displaystyle \frac{247x^{10}}{10}}-{\displaystyle \frac{85x^9}{3}}+32x^8\\ & & {-{\displaystyle \frac{256x^7}{7}}+{\displaystyle \frac{128x^6}{3}}-{\displaystyle \frac{256x^5}{5}}+64x^4-{\displaystyle \frac{256x^3}{3}}+128x^2-256x+256\ln (1+x)]}_0^1\\ & & =256\ln 2-{\displaystyle \frac{42629549}{240240}}\end{array}$$であることから、
$${\displaystyle \frac{1}{437580}}<256\ln 2-{\displaystyle \frac{42629549}{240240}}<{\displaystyle \frac{2}{437580}}$$
$${\displaystyle \frac{1}{256}}({\displaystyle \frac{1}{437580}}+{\displaystyle \frac{42629549}{240240}})<\ln 2<{\displaystyle \frac{1}{256}}({\displaystyle \frac{2}{437580}}+{\displaystyle \frac{42629549}{240240}})$$
$$(適当な小数表示0.6931471775<\ln 2<0.6931471864)$$
これを基に$\ln 3$を計算するのですが、$n=8$で計算すると
$${\displaystyle \frac{1}{1102449553920}}+{\displaystyle \frac{98165766053}{242106568704}}+\ln 2<\ln 3<{\displaystyle \frac{1}{734966369280}}+{\displaystyle \frac{98165766053}{242106568704}}+\ln 2$$
となり$\ln 2$の7桁精度に見合わないくらい高い精度で出てきてしまいます。
(適当な小数表示 $1.09861228562<\ln 3<1.09861229455$ ほぼ$\ln 2$の誤差しか残ってないことが分かる)

何故

ここでこの謎を解明するために$x^n(1-x{)}^n$を$m+x$で割った余りに着目したいと思います。

書き方が下手なんじゃ。まぁつまり何が言いたいかというと、$m$の値が大きくなるほど、積分した際の$\ln \left(m\right)$の係数が存外大きくなるということです。求める際にはその係数で割るので、思ったより精度良く出てしまうということだったんですね。

とりあえず一般化したもの。ベータ関数の公式使うと結構すっきりしますね。
ちなみに降順もいけます。

おまけ1:実践的にどうなの

ということでこちらの問題です。

$5^{\pi }$は整数か。のままだといろいろ大変ですのでこうしました。
この問題は言い換えると
$${\displaystyle \frac{\ln 156}{\ln 5}}<\pi <{\displaystyle \frac{\ln 157}{\ln 5}}$$
を示せばよいことになるので、早速近似の準備をしましょう。
$\ln 2$と$\ln 3$の近似はさっきのを使いまわします。計算量を出来るだけ削減したいならば$\ln 3$の近似は$n=7とかn=6$でも十分です。

$\ln 5$を求めよう

$\ln 5$は$\ln 4$の近似値が分かれば求められますが$\ln 4=2\ln 2$なのでもうわかってますね。${\displaystyle \frac{4!4!}{9!\times {20}^5}}\approx 5\times {10}^{-10}$なので$n=4$で十分。
予め${20}^4$で割ってあげると計算が楽です。
$${\displaystyle \frac{1}{{20}^4\cdot 5\cdot 630}}+2\ln 2<\ln 5+{\displaystyle \frac{1}{{20}^4}}({\displaystyle \frac{1}{8}}-{\displaystyle \frac{8}{7}}+{\displaystyle \frac{38}{6}}-{\displaystyle \frac{156}{5}}+{\displaystyle \frac{625}{4}}-{\displaystyle \frac{2500}{3}}+5000-40000)<{\displaystyle \frac{1}{{20}^4\cdot 4\cdot 630}}+2\ln 2$$
(小数表示$1.6094379060<\ln 5<1.6094379244$)

$\ln 156$を求めよう

$\ln 156=\ln (2^2\cdot 3\cdot 13)=2\ln 2+\ln 3+\ln 13$なので$\ln 13$の近似値が分かればよいです。
$\ln 13$の近似をするために$\ln 12$の値が必要ですが、$\ln 12=2\ln 2+\ln 3$なので分かっていますね。${\displaystyle \frac{3!3!}{7!\times (12\times 13{)}^4}}\approx {10}^{-11}$から$n=3$で十分。色々整理して、
$${\displaystyle \frac{1}{{156}^3}}({\displaystyle \frac{3646505}{12}}-{\displaystyle \frac{1}{12\cdot 140}})+2\ln 2+\ln 3<\ln 13<{\displaystyle \frac{1}{{156}^3}}({\displaystyle \frac{3646505}{12}}-{\displaystyle \frac{1}{13\cdot 140}})+2\ln 2+\ln 3$$
(小数表示$2.5649493482<\ln 13<2.5649493751$)
よって、
$$5.0498559888<\ln 156<5.0498560425$$
であることが分かる。

$\ln 157$を求めよう

といっても$\ln 156$が近似できているのでさっきの公式をそのまま使うだけです。ここまで$m$の値が大きければ$n=1$でも十分です。
$${\displaystyle \frac{1}{156\cdot 157}}({\displaystyle \frac{313}{2}}-{\displaystyle \frac{1!1!}{3!\cdot 156}})+\ln 156<\ln 157<{\displaystyle \frac{1}{156\cdot 157}}({\displaystyle \frac{313}{2}}-{\displaystyle \frac{1!1!}{3!\cdot 157}})+\ln 156$$
(小数表示:$5.0562457867<\ln 157<5.0562458407$)

あとは$\pi$と比較

上記から、$${\displaystyle \frac{\ln 156}{\ln 5}}<{\displaystyle \frac{5.0498560425}{1.6094379060}}<3.13765199$$
$${\displaystyle \frac{\ln 157}{\ln 5}}>{\displaystyle \frac{5.0562457867}{1.6094379244}}>3.14162212$$
であり、$\pi$は、 参考の記事 から、($n=2$の場合を引用して、)
$3.14159231<\pi <3.1415929$
(近似の方法の説明は上記の記事に丸投げします。すみません。)
であるので、
${\displaystyle \frac{\ln 156}{\ln 5}}<\pi <{\displaystyle \frac{\ln 157}{\ln 5}}$が示された。

まだ結構余裕があるので計算量減らせそうですね。

おまけ2

また国際信州学院大学です。
これも
${\log }_{10}7.14\le 1800{\log }_{10}2 -\lfloor 1800{\log }_{10}2\rfloor <{\log }_{10}7.15$
${\displaystyle \frac{\ln 714}{\ln 10}}-2\le 1800{\displaystyle \frac{\ln 2}{\ln 10}}-\lfloor 1800{\displaystyle \frac{\ln 2}{\ln 10}}\rfloor <{\displaystyle \frac{\ln 715}{\ln 10}}-2$
を示せれば良いです。
(おまけ1までの計算結果は使いまわします。)
$\ln 10$は$\ln 2+\ln 5$で表せるので、
$$2.3025850835<\ln 10<2.3025851108$$
となります。

$\ln 715$を求めよう

$715=5\times 11\times 13$より$\ln 11$が求まれば$\ln 715$が求まります。
$\ln 12\to \ln 11、n=3$で近似すると
$${\displaystyle \frac{1}{{132}^3}}({\displaystyle \frac{1}{12\times 140}}-{\displaystyle \frac{12007403}{60}})+\ln 12<\ln 11<{\displaystyle \frac{1}{{132}^3}}({\displaystyle \frac{1}{11\times 140}}-{\displaystyle \frac{12007403}{60}})+\ln 12$$
(小数表示$2.3978952636<\ln 11<2.3978952904$)

$\ln 715=\ln 5+\ln 11+\ln 13$より
$$6.5722825178<\ln 715<6.5722825899$$

$\ln 714$を求めよう

$\ln 715\to \ln 714$,$n=1$で近似すると
$${\displaystyle \frac{1}{714\cdot 715}}({\displaystyle \frac{1!1!}{3!\cdot 715}}-{\displaystyle \frac{1429}{2}})+\ln 715<\ln 714<{\displaystyle \frac{1}{714\cdot 715}}({\displaystyle \frac{1!1!}{3!\cdot 714}}-{\displaystyle \frac{1429}{2}})+\ln 715$$
(小数表示:$6.5708829379<\ln 714<6.5708830101$)

あとは

$${\displaystyle \frac{\ln 714}{\ln 10}}-2<{\displaystyle \frac{6.5708830101}{2.3025850835}}-2<0.85369824429$$
$${\displaystyle \frac{\ln 715}{\ln 10}}-2>{\displaystyle \frac{6.5722825178}{2.3025851108}}>0.85430600891$$
$$0.85398561295<1800{\displaystyle \frac{\ln 2}{\ln 10}}-\lfloor 1800{\displaystyle \frac{\ln 2}{\ln 10}}\rfloor <0.85399899469$$
であるから、
$${\displaystyle \frac{\ln 714}{\ln 10}}-2<0.85369824429<0.85398561295<1800{\displaystyle \frac{\ln 2}{\ln 10}}-\lfloor 1800{\displaystyle \frac{\ln 2}{\ln 10}}\rfloor <0.85399899469<0.85430600891<{\displaystyle \frac{\ln 715}{\ln 10}}-2$$
よって示された。

おわりに

これが出来ることによって$p,q$を有理数($q>0$)としたとき、$\ln p$とか${\log }_{q}p$とかの近似が高校範囲でできます。実用性までは考えてません。あとこれやってると昔やってたマージゲームを思い出します。若干違うけど。そして 🐟️🍊みかん🍊🐟️ 様には大変感謝です。

計算ミス等があれば教えてください。

参考

各学部の過去問 | 国際信州学院大学, 閲覧日2024年9月3日https://kokushin-u.jp/prospective-student/past-questions/

円周率を積分で近似計算する | 🐟️🍊みかん🍊🐟️, 閲覧日2024年9月3日 https://mathlog.info/articles/F5qZheDGn1M44uI6BFaD