（その１）に掲げた，区間［ａ−ｘ，ａ＋ｘ］に対して，

　　２ｘ／ａ＜∫ｄｔ／ｔ＜ｘ（１／（ａ＋ｘ）＋１／（ａ＋ｘ））

の左辺は矩形公式，右辺は台形公式に拠っていることはおわかり頂けるであろう．念のため説明するが，ｙ＝１／ｔのグラフで，ｔ＝［ａ−ｘ，ａ，ａ＋ｘ］となる３点をとればよい．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

　（その１）では

　　２ｘ／ａ＜∫ｄｔ／ｔ＜ｘ（１／（ａ＋ｘ）＋１／（ａ＋ｘ））

を２回適用し，

［１］ａ＝５／４，ｘ＝１／４→２／５＜∫(1,3/2)ｄｔ／ｔ＜５／１２

［２］ａ＝７／４，ｘ＝１／４→２／７＜∫(3/2,2)ｄｔ／ｔ＜７／２４

より

　　２４／３５＜ｌｏｇ２＝∫(1,2)ｄｔ／ｔ＜１７／２４

　　０．６８５＜ｌｏｇ２＜０．７０８

と絞り込んだが，１回適用では無理なのだろうか？

　　ａ−ｘ＝１，ａ＋ｘ＝２→ａ＝３／２，ｘ＝１／２

　　１／３＜ｌｏｇ２＝∫(1,2)ｄｔ／ｔ＜３／２

　　０．３３３＜ｌｏｇ２＜１．５

となり，近似の精度を上げることはできない．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

　（その２）の台形則，シンプソン則，中点則をそのまま適用してみると，

［１］台形則

　　∫(a,b)ｆ（ｘ）ｄｘ≒（ｂ−ａ）／２・｛ｆ（ａ）＋ｆ（ｂ）｝

＝１／２・（１／１＋１／２）＝３／４＝０．７５

［２］シンプソン則（１／３則）

　　∫(a,b)ｆ（ｘ）ｄｘ≒（ｂ−ａ）／６・｛ｆ（ａ）＋４ｆ（（ａ＋ｂ）／２）＋ｆ（ｂ）｝

　　＝１／６・（１／１＋４・２／３＋１／２）＝２５／３６＝０．６９４

［３］中点則

　　∫(a,b)ｆ（ｘ）ｄｘ≒（ｂ−ａ）ｆ（（ａ＋ｂ）／２）

　　＝２／３＝０．６

　さらに

［４］シンプソン則（３／８則）

　　∫(a,b)ｆ（ｘ）ｄｘ≒３ｈ／８・｛ｆ（ａ）＋３ｆ（ａ＋ｈ）＋３ｆ（ａ＋２ｈ）ｆ（ｂ）｝　　（ｈ＝（ｂ−ａ）／３）

　　＝１／８（１／１＋３・３／４＋３・３／５＋１／２）＝１１１／１６０＝０．６９４

［５］ブール則

　　∫(a,b)ｆ（ｘ）ｄｘ≒２ｈ／４５・｛７ｆ（ａ）＋３２ｆ（ａ＋ｈ）＋１２ｆ（ａ＋２ｈ）＋３２ｆ（ａ＋３ｈ）＋７ｆ（ｂ）｝　　（ｈ＝（ｂ−ａ）／４）

　　＝１／９０（７＋３２・４／５＋１２・２／３＋３２・４／７＋７／２）＝４３６７／６３００＝０．６９３

　シンプソン則（１／３則）の大健闘ということであろう．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝