一个计算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;
}