C++求积分代码

基本上是靠定义来求的，而且代码也简明扼要，精度为1/1000；

如下：

1. #include"iostream"
2. #include"cstdlib"
3. #include"cmath"
4. #define N 1000
5. using namespace std;
6. 
   
7. double f(double x)
8. {
9. return x*x*x;
10. }
11. double  intergrate(double(*fp)(double), double a, double b)
12. {
13. fp = f;
14. int N = (b - a) * N;
15. double dx = (double)1/N;
16. double s = 0;
17. for (int i = 1; i <= N; i++)
18. {
19. s += f(a + dx)*dx;
20. a += dx;
21. }
22. return s;
23. }
24. 
    
25. int main()
26. {
27. cout << "f(x)=x^3" << endl;
28. cout << "f(x)在(0,1)的积分：" << intergrate(f, 0, 1);
29. system("pause");
30. return 0;
31. }

想要更改精度的话可以直接修改N的值