蒙特卡洛求积分c语言

蒙特卡洛求积分是一种基于随机抽样的数值积分方法，它利用统计学原理来估计函数的定积分。在C语言中，可以这样实现：

1. **基本步骤**：
- 定义函数域和积分区间。
- 创建大量随机数（通常选择在0到1之间的均匀分布）。
- 计算每个随机点上函数值乘以区间宽度。
- 将所有点的结果相加并除以总的点数，得到近似积分值。
- 为了提高精度，可以增加样本点的数量。

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

// 假设我们要计算f(x) = x^2在[0, 1]上的积分
double func(double x) {
    return x * x;
}

int main() {
    int n_points = 1000000; // 样本点数量
    double sum = 0.0;
    srand(time(0)); // 设置随机种子

    for (int i = 0; i < n_points; ++i) {
        double u = (double) rand() / RAND_MAX; // 随机数在0到1之间
        double v = (double) rand() / RAND_MAX;
        if (u <= v) { // 只需考虑一维情况，这里简化条件
            sum += func(u);
        }
    }

    double result = n_points * sum / (double)n_points; // 平均值即为估算值
    printf("蒙特卡洛积分结果：%.6f\n", result);
    return 0;
}