魯陽(yáng)
(赤峰學(xué)院 教務(wù)處,內(nèi)蒙古 赤峰 024000)
第一類(lèi)橢圓積分計(jì)算機(jī)算法實(shí)現(xiàn)
魯陽(yáng)
(赤峰學(xué)院 教務(wù)處,內(nèi)蒙古 赤峰 024000)
在某些科研工作中,常要用到大量數(shù)值計(jì)算,其中橢圓積分實(shí)際應(yīng)用面不大,但計(jì)算機(jī)復(fù)雜,而且C語(yǔ)言數(shù)據(jù)函數(shù)集中沒(méi)有收錄,在網(wǎng)上也找不到此類(lèi)算法,本文提供了一個(gè)C語(yǔ)言函數(shù)程序,語(yǔ)句簡(jiǎn)單易讀,只需稍懂一點(diǎn)C語(yǔ)言,就可編出相應(yīng)的主程序,通過(guò)調(diào)用該函數(shù)得到計(jì)算結(jié)果.
算法;第一類(lèi)完全橢圓積分;計(jì)算機(jī)
分.
k:雙精度實(shí)型變量.0≤k≤1.
f:雙精度實(shí)型變量.該積分中的φ.
該函數(shù)簡(jiǎn)單易讀,已在TURBO C 2.0上運(yùn)行通過(guò).
〔1〕楊慶揚(yáng),等.數(shù)值分析.華是理工大學(xué)出版社,1986.
〔2〕[美]W.H.普雷斯,等.數(shù)值方法大全,科學(xué)計(jì)算的藝術(shù).半州大學(xué)出版社,1991.
〔3〕徐士良.FORTRAN 常用算法程序集.清華大學(xué)出版社,1991.
TP391
A
1673-260X(2010)05-0021-02實(shí)際工作中我們常常要用到橢圓積分,通??梢酝ㄟ^(guò)查表等手段計(jì)算,但在某些科研工作中,比如單擺實(shí)驗(yàn)中的一些數(shù)據(jù),就要遇到大量的此類(lèi)問(wèn)題.那么,一個(gè)能用計(jì)算機(jī)計(jì)算的方法就顯得很有必要,而把該算法用C語(yǔ)言函數(shù)的形式實(shí)現(xiàn),處理具體問(wèn)題時(shí)調(diào)用此函數(shù)則方便宜行.下面我們用變步長(zhǎng)辛卜生求積法實(shí)現(xiàn)此算法.