/*函數：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