資源簡介
圓周率的計算C程序(數值法)梯形積分公式,能運行?。?!

代碼片段和文件信息
#include
#include
#include
#define?N????20
#define?Max??10
#define??a????0
#define??b????1
#define?eps???1e-6
#define?pi???3.14159
double?f(double?x)???????????????????????????????????????/*計算f(x)在各點值的函數*/
{double?y;
y=4/(1+x*x);
return(y);
}
double?trapz(double?a1double?b1long?int?n)??????????????????/*梯形公式求小區間上的積分函數*/
{int?i;
??double?sumh=(b1-a1)/n;
??for(i=1;i ??sum+=f(a1+i*h);
??sum+=(f(a1)+f(b1))/2;
??return(h*sum);
}
void?main()????????????????????????????????????????????/*Romberg積分程序主函數*/
{int?ik;
?long?int?n=Nm=0;
?double?T[Max+1][2];
?T[0][1]=trapz(abn);
?n*=2;
?for(m=1;m ??{
???for(i=0;i ??????{T[i][0]=T[i][1];}
???T[0][1]=trapz(abn);
???n*=2;
???for(i=1;i<=m;i++)
????T[i][1]=T[i-1][1]+(T[i-1][1]-T[i-1][0])/(pow(22*m)-1);
??if(fabs(T[m-1][1]-T[m][1]) ????{
?????printf(“pi的估計值為:%lf\n“T[m][1]/2);
?????break;
?if(?k?!=?1?)?exit(0);?}
????}
}
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1056??2009-05-24?14:44??Text1.c
-----------?---------??----------?-----??----
?????????????????1056????????????????????1
- 上一篇:RawCap.exe
- 下一篇:ViewUrl 上網信息檢查工具
評論
共有 條評論