数值积分之Newton_Cotes闭合积分公式

本文主要是介绍数值积分之Newton_Cotes闭合积分公式，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

//Newton_Cotes闭合积分公式
#include <iostream>
#include <math.h>

using namespace std;

class trapezoidal
{
private:
 int n, k;
 double a, b, f, h, integral, sum;

public:
 double func(double v)
 {
  f = 24.6 / (v - 0.065) - 5.5 / pow(v, 2);
  return f;
 }
 void integration();
};

void main()
{
 trapezoidal trape;
 trape.integration();
}

void trapezoidal::integration()
{
 cout << "\n输入下限和上限：";
 cin >> a >> b;
 cout << "\n输入你要的段数：";
 cin >> n;
 h = (b - a) / n;
 sum = 0.0;
 for (k = 1; k < n; k++)
 {
  sum += func(a + k * h);
 }
 integral = (h / 2.0) * (func(a) + 2.0 * sum + func(b));
 cout << endl << "积分值 = " << integral << endl;
}

这篇关于数值积分之Newton_Cotes闭合积分公式的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---