计算方法实验：用拉格朗日插值法计算插值问题

实验内容 用复合梯形求积公式和复合辛普生求积公式计算定积分（+小实验 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