魯陽
(赤峰學院 教務處,內蒙古 赤峰 024000)
第一類橢圓積分計算機算法實現
魯陽
(赤峰學院 教務處,內蒙古 赤峰 024000)
在某些科研工作中,常要用到大量數值計算,其中橢圓積分實際應用面不大,但計算機復雜,而且C語言數據函數集中沒有收錄,在網上也找不到此類算法,本文提供了一個C語言函數程序,語句簡單易讀,只需稍懂一點C語言,就可編出相應的主程序,通過調用該函數得到計算結果.
算法;第一類完全橢圓積分;計算機

分.

k:雙精度實型變量.0≤k≤1.
f:雙精度實型變量.該積分中的φ.



該函數簡單易讀,已在TURBO C 2.0上運行通過.

〔1〕楊慶揚,等.數值分析.華是理工大學出版社,1986.
〔2〕[美]W.H.普雷斯,等.數值方法大全,科學計算的藝術.半州大學出版社,1991.
〔3〕徐士良.FORTRAN 常用算法程序集.清華大學出版社,1991.
TP391
A
1673-260X(2010)05-0021-02實際工作中我們常常要用到橢圓積分,通常可以通過查表等手段計算,但在某些科研工作中,比如單擺實驗中的一些數據,就要遇到大量的此類問題.那么,一個能用計算機計算的方法就顯得很有必要,而把該算法用C語言函數的形式實現,處理具體問題時調用此函數則方便宜行.下面我們用變步長辛卜生求積法實現此算法.