用复合梯形求积公式和复合辛普生求积公式计算定积分(+小实验 7) (大实验见 p317,实验 6 数值积分(1)(2) 【小实验 7】 求函数在指定点处的数值导数:𝑓𝑓(𝑥𝑥) = � 𝑥𝑥 𝑥𝑥2 𝑥𝑥3 1 2𝑥𝑥 3𝑥𝑥2 0 2 6𝑥𝑥 �, 𝑥𝑥 = 1,2,3
大实验 #include<iostream> using namespace std; double fun(double x) { return 4 / (1 + x * x); } double Simpson(double a, double b, int n) { double h = (b - a) / n; double s1 = fun(a + h / 2); double s2 = 0; for (int i = 1; i < n; i++) { s1 += fun(a + i * h + h / 2); s2 += fun(a + i * h); } return h * (fun(a) + 4 * s1 + 2 * s2 + fun(b)) / 6; } int main() { int n = 10; //区间等分数 double a = 0, b = 1; //上下限 cout << "被积函数 f = 4 / (1 + x * x)" << endl; cout << "区间等分数为 10" << endl; cout << "积分区间为(0,1)" << endl; cout << "积分值为:" << Simpson(a, b, n) << endl; return 0; } 【小实验 7】
通过这次实验,我更加深入地了解了复合梯形求积公式和复合辛普生求积公式,以及如何运用他们进行解题,并熟练的运用MATLAB进行解题,锻炼了我的思维能力,使我更加灵活熟练度使用MATLAB |
计算方法实验:用拉格朗日插值法计算插值问题
最新推荐文章于 2024-09-16 17:46:36 发布