我正在尝试创建一个 C 程序来对 0 和无穷大之间的 sin(x)/sqrt(x) 进行积分。我使用梯形规则,当函数趋于无穷大时切断端点。
但是返回的总数太高,我不知道为什么。代码如下:
#include<math.h>
#include<stdio.h>
double func(double u)
{ double a;
a = ((sin(u))/(sqrt(u)));
return a;}
void main()
{
int i, N;
double sum, u, a, b, h, Fa, Fb, F;
printf("Enter value of N\n");
scanf("%d" ,&N);
a=0.01;
b=1000;
h=(b-a)/(N-1);
sum=0;
F=func(a);
u=a;
for(i=0; i<N; i++)
{
sum=sum+F;
u=u+h;
F=fabs(func(u));
}
Fa=func(a);
Fb=func(b);
sum=sum-(0.5*Fa)-(0.5*Fb);
sum=sum*h;
printf("I: %lf\n", sum);
}
有什么想法吗?
最佳答案
只需删除 F=fabs(func(u));
行中的 fabs
即可。
并且您应该在末尾使用 int main(void)
和 return 0;
而不是 void main()
。
关于c - 如何用C实现无限极限梯形积分?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/29087126/