一个计算Pi的程序
用的是数值积分算法
用的是数值积分算法
- #include <math.h>
- #include <stdio.h>
- double f(double x);
- main()
- {
- int i=0;
- int a=0;
- int b=0;
- int n=0;
- double x=0;
- double s=0;
- double h=0;
- printf("/nplease input a=");
- scanf("%d",&a);
- printf("/nplease input b=");
- scanf("%d",&b);
- printf("/nplease input a number N=");
- scanf("%d",&n);
- h=(b-a)/(2.0*n);
- s=f(a)-f(b);
- for(x=a+h,i=1;i<(2*n+1);i++,x=x+h)
- {
- if((i%2)==1)
- {
- s=s+f(x)*4;
- }
- else
- {
- s=s+f(x)*2;
- }
- }
- s=s*(b-a)/(6.0*n);
- printf("%5.51f",4*s);
- }
- double f(double x)
- {
- double value=0;
- value= 1/(1+x*x);
- return value;
- }