蒙特卡罗法求椭圆面积c语言代码

下面是使用蒙特卡罗方法求椭圆面积的 C 语言代码：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

int main()
{
    srand(time(NULL)); // 初始化随机数种子
    int n = 1000000; // 抛点总数
    int count = 0; // 落在椭圆内的点的个数
    double a = 2.0; // 椭圆长轴
    double b = 1.0; // 椭圆短轴
    double x, y;
    for (int i = 0; i < n; i++) {
        x = 2.0 * rand() / RAND_MAX - 1.0; // 在 [-1, 1] 区间内随机生成 x 坐标
        y = 2.0 * rand() / RAND_MAX - 1.0; // 在 [-1, 1] 区间内随机生成 y 坐标
        if (x * x / (a * a) + y * y / (b * b) <= 1.0) {
            count++; // 落在椭圆内
        }
    }
    double area = 4.0 * a * b * count / n; // 计算椭圆面积
    printf("The area of the ellipse is: %f\n", area);
    return 0;
}

该程序使用了随机抛点的方法，将落在椭圆内的点的个数与总点数的比例乘以矩形面积来估算椭圆的面积。由于随机抛点的方式具有随机性，因此得到的结果可能会有一定的误差。