使用辛普森规则计算定积分的Java实现

使用辛普森规则计算定积分的Java实现

在数值计算中，辛普森规则是一种常用的方法，用于估算函数的定积分。它基于将函数曲线划分为若干个小区间，并在每个小区间内采用二次多项式来逼近原函数。这篇文章将介绍如何使用Java语言实现辛普森规则来计算定积分，并提供相应的源代码。

首先，我们需要了解辛普森规则的原理。对于一个连续函数 f(x)，我们可以将其在区间 [a, b] 上进行积分。辛普森规则的基本思想是，在每个小区间 [x_i, x_{i+1}] 内，我们使用二次多项式来逼近原函数。对于每个小区间，我们可以使用以下公式来计算定积分的近似值：

∫[x_i, x_{i+1}] f(x) dx ≈ (h/3) * [f(x_i) + 4f(x_{i+1/2}) + f(x_{i+1})]

其中，h = (b - a) / n 是子区间的宽度，n 是划分区间的数量。使用这个公式，我们可以将整个区间 [a, b] 分成 n 个小区间，并对每个小区间应用辛普森规则，最后将结果累加得到定积分的近似值。

现在，让我们来编写Java代码实现这个算法。首先，我们需要定义一个函数，用于计算 f(x) 的值。假设我们要计算 f(x) = x^2 在区间 [0, 1] 上的积分，可以使用以下代码：

public static double function(double x) {
   
    return x * x;
}

接下来，我们需要编写一个辛普森规则的函数，用于计算定积分的近似值。代码如下：