单位圆的面积为π,因此可以通过求单位圆面积的近似值来求π的近似值

一个学妹问的一道c语言题目

标题单位圆的面积为π,因此可以通过求单位圆面积的近似值来求π的近似值。考虑将单位圆沿x轴分成n份，从而求得n个梯形的面积和作为单位圆的一个近似值；
再将单位圆分成2n份，从而求得2n个梯形的面积和作为单位圆的另一近似值。此过程可以一直持续下去，将单位圆分成4n份、8n份等等，当二次求得的近似值很接近时，可认为得到单位圆面积较精确的近似值。将下列程序补充完整，利用上述方法求π的近似值

#include <math.h>#include <iostream>using namespace std;int main(){const  double epsilon = 1E-4;//这里我特意将极限写成负四次方，不然执行时间太长int n  = 2;double dist_s = 1.0;//默认1double pi = 0.0;for (n; epsilon < dist_s; n ++ ){   double x1 = 2.0/n;//高    cout << "x1 = " << x1 << endl;    double x2 = 2.0/(2*n);//高    cout << fixed<<"x2 = " << x2 << endl;    double s1 = 0.0;//x/n高度梯形的面积    double s2 = 0.0;//x/2n高度梯形的面积    cout << "n = " << n << endl;        for (double i = -1;i < 1;i = i + x1){            double h1 = 2*sqrt(1-pow(i,2));//上底            double h2 = 2*sqrt(1-pow(i + x1,2));//下底            s1 += (h1+h2)*x1/2.0;        }        cout << "s1 = " << s1 << endl;        for (double i = -1;i < 1;i = i + x2){            double h1 = 2*sqrt(1-pow(i,2));//上底            double h2 = 2*sqrt(1-pow(i + x2,2));//下底            s2 += (h1+h2)*x2/2.0;        }        cout << "s2 = " << s2 << endl;    dist_s = fabs(s1 - s2);    cout <<"dist =" << dist_s << endl;}cout << "pi = " << dist_s << endl;}

codeblock中运行