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

分.

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



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

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