Question

我使用梯形法则计算0和无穷大之间函数的积分。我可以计算给定N值的积分值，现在我试图将N从2循环到给定值，但它不起作用。当N为2并重复而不是N的新值时，它会不断计算积分值。问题出在main（）中的for循环中我认为。

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <float.h>

double f(double x) {
 double a;
  a =1/((1+x)*pow(x,0.5));
  return a;}

double tra(double upper, double lower, int N) {
 double sum, step, integral,lowest;
  step=(upper-lower)/(N-1);
  lower=lower+step;

    if(lower==0) {
       lowest=DBL_EPSILON;}
    else {
    lowest=lower;}

    while(lower<upper) {
         sum=sum+f(lower);
         lower=lower+step;}

 integral=step*(sum+(f(upper)/2)+(f(lowest)/2));
 sum=0;
 return integral;}


main() {
int N;
double upper=DBL_EPSILON*2, lower=0, total=0;
for(N=2;N<20000;N+=100) {   /*Here im trying to loop N so that the        integral is calculated for increasing values of N*/

   while(upper<FLT_MAX) {
        total=total+tra(upper, lower, N);
        lower=upper;
        upper=upper*2;}
        printf("Integral is %.10f\n", total);
}
}

Answer 1

我建议你将变量初始化移到for循环中，如下所示：

int main(void) {
    int N;
    double upper, lower, total;
    for(N=2;N<20000;N+=100) {
        upper = DBL_EPSILON*2;
        lower = 0;
        total = 0;
        while(upper<FLT_MAX) {
            total=total+tra(upper, lower, N);
            lower=upper;
            upper=upper*2;
        }
        printf("Integral is %.10f\n", total);
    }
    return 0;
}

数值积分，使用C中的梯形法则

1 个答案: