谈Stirling式的改良.PDFVIP

談Stirling 公式的改良 蔡永裕 壹. 緣 起 從學生時代接觸到 Stirling 公式後, 就 n! = e−nnn √2πn + h (1) 被它深深吸引–利用一個簡單的計算式, 三招 . −n n √2πn; (2) = e n 兩式下就能求得原來必須反覆計算才能求解 n = 1, 2, 3, . . . 的數據。 但除非 n 階乘的 n 真的夠大, 大到 註: 遠遠超乎我們的實用範圍, 否則其精確度始 終令人不太敢恭維! 例如相對誤差只要求萬 e = 2.71828 18284 59045 23536 分之一, n 竟然就要高於 800, (此時 n! 的值 02874 71352 66249 77572 是幾乎是 10 的 2000次方)。 在現代, 個人電 47093 . . . 腦如此普及的情況下想算得 n! 精確一點的 π = 3.14159 26535 89793 23846 數據, 難道還是這麼可望而不可及嗎? 在一 26433 83279 50288 41971 次偶然的機緣中找到了其改良式, 原已準備 發表, 卻因為搬家而造成了論文資料遺失 (包 69399... h 1 含文件失蹤、 磁片發霉、 電腦故障等三衰齊 (0 ) n! 12n 至)。 因為筆者個性隨緣, 便也隨匆匆歲 月而 使用 (2) 式在所謂 「大的 n 值」 時, 會 淡忘此事。 最近失蹤已久的文件忽重現江湖, 得到 n! 的漸近相等值。 但是 n 值究竟要大到 在賢妻的督促下, 遂配合電腦重新撰寫本文, 多少, 代入式子中才會有令人滿意的結果呢? 順便為本刊 66 期蔡聰明先生大作 「談 Stir- 100? 500? 1000? . . ., 經過驗算後, 雖然這 ling 公式」 做一補充。 些數據已大到幾乎超越一般工程運算所能使 用到的範圍, 卻總有幾分美中不足之憾。 也就 貳. Stirling’s formula 是說– 54 談 Stirling 公式的改良 55 1. 「n 值如果不夠大」 (例如小於 70), 則以 如果這個假設可行, E 值如何求得呢? 一般工程用計算機直接將答案按出即可, 就算再大些, 使用 Microsoft 的視窗軟 體工具 「小算盤」 等類似工具, 亦可算到