Я пытаюсь создать программу c для интеграции sin(x)/sqrt(x) между 0 и бесконечностью. Я использую правило трапеции, отрезая конечные точки, когда функция стремится к бесконечности.
Однако общая сумма возврата слишком высока, и я не уверен, почему. Вот код:
#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);
}
Какие-нибудь мысли?
scanf()
, поэтому может случиться так, что ваш код вызывает неопределенное поведение. - person Iharob Al Asimi   schedule 17.03.2015