数值积分C语言算法

 

一个计算Pi的程序
用的是数值积分算法

1. #include <math.h>
2. #include <stdio.h>
3. double f(double x);
4. main()
5. {
6. int     i=0;
7. int     a=0; 
8. int     b=0;
9. int     n=0;
10. double  x=0;
11. double  s=0;
12. double  h=0;
13. printf("/nplease input a=");
14. scanf("%d",&a);
15. printf("/nplease input b=");
16. scanf("%d",&b);
17. printf("/nplease input a number N=");
18. scanf("%d",&n);
19. h=(b-a)/(2.0*n);
20. s=f(a)-f(b);
21. for(x=a+h,i=1;i<(2*n+1);i++,x=x+h)
22. {
23.     if((i%2)==1) 
24.     { 
25.     s=s+f(x)*4;
26.     }
27.     else
28.     {
29.     s=s+f(x)*2;
30.     }
31. }
32. s=s*(b-a)/(6.0*n);
33. printf("%5.51f",4*s);
34. }
36. double f(double x)
37. {
38. double value=0;
39. value= 1/(1+x*x);
40. return  value;
41. }