第一类椭圆积分计算机算法实现

鲁阳

（赤峰学院 教务处，内蒙古 赤峰 024000）

第一类椭圆积分计算机算法实现

鲁阳

（赤峰学院 教务处，内蒙古 赤峰 024000）

在某些科研工作中，常要用到大量数值计算，其中椭圆积分实际应用面不大，但计算机复杂，而且C语言数据函数集中没有收录，在网上也找不到此类算法，本文提供了一个C语言函数程序，语句简单易读，只需稍懂一点C语言，就可编出相应的主程序，通过调用该函数得到计算结果.

算法；第一类完全椭圆积分；计算机

1 函数定义

分.

2 形参说明

k：双精度实型变量.0≤k≤1.

f：双精度实型变量.该积分中的φ.

3 函数程序tyjf.c

该函数简单易读，已在TURBO C 2.0上运行通过.

〔1〕杨庆扬，等.数值分析.华是理工大学出版社，1986.

〔2〕[美]W.H.普雷斯，等.数值方法大全，科学计算的艺术.半州大学出版社，1991.

〔3〕徐士良.FORTRAN 常用算法程序集.清华大学出版社，1991.

TP391

A

1673-260X（2010）05-0021-02实际工作中我们常常要用到椭圆积分，通常可以通过查表等手段计算，但在某些科研工作中，比如单摆实验中的一些数据，就要遇到大量的此类问题.那么，一个能用计算机计算的方法就显得很有必要，而把该算法用C语言函数的形式实现，处理具体问题时调用此函数则方便宜行.下面我们用变步长辛卜生求积法实现此算法.