使用C语言编写一个求积分的代码
时间: 2024-04-30 15:20:37 浏览: 139
下面是一个使用C语言求解定积分的示例代码:
```c
#include <stdio.h>
#include <math.h>
double f(double x) {
return sin(x); // 定义被积函数 f(x)
}
double integral(double a, double b, int n) {
double h = (b - a) / n;
double sum = 0.0;
for (int i = 0; i < n; i++) {
double x1 = a + i * h;
double x2 = x1 + h;
double y1 = f(x1);
double y2 = f(x2);
sum += (y1 + y2) / 2.0 * h;
}
return sum;
}
int main() {
double a = 0.0, b = 3.1415926; // 积分区间 [a, b]
int n = 1000; // 分割数
double result = integral(a, b, n);
printf("The integral of sin(x) from %f to %f is %f\n", a, b, result);
return 0;
}
```
该代码使用梯形法计算定积分,即将积分区间分割成若干个小区间,每个小区间近似看作一个梯形,然后计算所有梯形面积之和,最终得到定积分的近似值。
阅读全文