变步长数值积分算法+动态演示

数值算法课程：变步长数值积分算法+动态演示，C语言+flash开发

/*函数：f（x）=sin（x）/x	积分区间：（ab）=（06）*/
#include
#include

double f（double x）
{
   return sin（x）/x;
}

int main（）
{
	double abehT0T1; //a积分起点、b积分终点、e精度
	printf（“请输入积分区间[ab]：“）;
	scanf（“%lf%lf“&a&b）;
	printf（“请输入精度e：“）;
	scanf（“%lf“&e）;


	//解决输入的积分区间端点中存在0的情况
	if（a==0）
		a+=0.0000000001;
	if（b==0）
		b+=0.0000000001;
	//解决输入的积分区间端点中存在0的情况


	h=b-a;
	T0=h*（f（a）+f（b））/2;
	T1=T0/2+h/2*f（a+h/2）;
	int n=1;
	while（fabs（T1-T0）>=e）
	{
		h=h/2;
		n=2*n;
		T0=T1;
		T1=0;
		for（int k=1;k<=n;k++）
			T1=T1+f（a+（k-1/2）*h）;
		T1=T0/2+h*T1/2;
	}

	printf（“积分结果为：%lf\n“T1）;
	return 0;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        217  2012-04-21 14:58  变步长数值积分算法+动态演示\动态演示\使用说明.txt

     文件      16012  2011-06-01 20:49  变步长数值积分算法+动态演示\动态演示\变步长数值积分算法.swf

     文件        718  2011-05-31 16:35  变步长数值积分算法+动态演示\变步长数值积分算法.cpp

     目录          0  2012-04-21 14:57  变步长数值积分算法+动态演示\动态演示

     目录          0  2012-04-21 14:57  变步长数值积分算法+动态演示

----------- ---------  ---------- -----  ----

                16947                    5