4_34.cpp
上传用户:zipjojo
上传日期:2009-07-20
资源大小:70k
文件大小:1k
源码类别:

文章/文档

开发平台:

C/C++

  1. # include<stdio.h>
  2. # include<math.h>
  3. # include<iostream.h>
  4. double fx(double w);
  5. double integral(double s,double t,int m);
  6. void main()
  7. {
  8. int n;
  9. double a,b,I;
  10. printf("请输入积分区间(从a到b),并输入等分区间数n:n");
  11.     cin>>a>>b>>n;
  12.     I=integral(a,b,n);
  13. printf("I=%fn",I);
  14. }
  15. double fx(double w)
  16. {
  17. return (log(1.0+w)/(1.0+w*w));
  18. }
  19. double integral(double s,double t,int m)
  20. {
  21. double (*p)(double w);
  22. p=fx;
  23. double f=0.0,h,x,y1,y2,area;
  24. h=(t-s)/m;
  25. x=s;
  26. y1=(*p)(x);
  27. for(int i=1;i<=m;i++)
  28. {
  29. x+=h;
  30. y2=(*p)(x);
  31. area=(y1+y2)*h/2;
  32. y1=y2;
  33. f+=area;
  34. }
  35. return f;
  36. }