在C語言中,進行數值積分的方法有很多種,其中一些常見的方法包括辛普森法則(Simpson’s Rule)、梯形法則(Trapezoidal Rule)和布倫特-薩爾明法則(Brent-Salamin Algorithm)等。
以下是使用辛普森法則進行數值積分的一個示例代碼:
#include <stdio.h>
double f(double x) {
return x * x;
}
double simpson(double a, double b, int n) {
double h = (b - a) / n;
double sum = f(a) + f(b);
for (int i = 1; i < n; i += 2) {
sum += 4 * f(a + i * h);
}
for (int i = 2; i < n - 1; i += 2) {
sum += 2 * f(a + i * h);
}
return h / 3 * sum;
}
int main() {
double a = 0, b = 1;
int n = 100; // 積分點的個數
double result = simpson(a, b, n);
printf("Integral of f(x) from %f to %f is %f\n", a, b, result);
return 0;
}
在上述代碼中,f(x)
是被積函數,simpson()
函數使用辛普森法則計算定積分的值。main()
函數中指定了積分的區間 [a, b]
和積分點的個數 n
,然后調用 simpson()
函數計算積分值并輸出結果。
請注意,這只是一個示例代碼,實際應用中可能需要根據具體問題進行調整和優化。同時,數值積分方法的選擇也會對結果精度和計算效率產生影響。