c - 与c的数值积分

我必须对文本文件中的一组数据点进行数值积分。
我的数据点看起来像

0.5   0.479425539
1     0.841470985
1.5   0.997494987
2     0.909297427
2.5   0.598472144
3     0.141120008
3.5   -0.350783228
4     -0.756802495
4.5   -0.977530118
5     -0.958924275

我的尝试是

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

double trapezoidalRule (double size, double *x, double *y)
{
    double sum = 0.0,increment;
    int k;
    for (k=1; k<size; k++)
    {
        increment = 0.5 * (x[k]-x[k-1]) * (y[k]+y[k-1]);
        sum += increment;
    }
    return sum;
    _getch();
}
int main ( int argc, char * argv[])
{
    char*  fileName = argc > 1 ? argv[1] : "C:\\Users\\g\\Desktop\\test.txt";
    FILE*  inputFile = fopen (fileName, "r");
    int  k;
    double size,*x, *y;
    double integral;
    if ( inputFile ==NULL)
    {
        fprintf (stderr, "Open failed for %s\n", fileName);
        exit(666);
    }
    fscanf (inputFile, "%d", &size);
    printf (" Number of points: %d\n", size);

    x = (double *) calloc (size, sizeof(*x));
    y = (double *) calloc (size, sizeof(*y));

    for (k=0; k< size; k++)
        fscanf (inputFile, "%lg%lg" , x+k, y+k);
    integral = trapezoidalRule (size, x, y);
    printf ("Integral:", "\n", integral);
    printf ("\n");
    //printf ("check: ","\n", x[size-1], log(x[size-1]) );
    _getch();
}

但我没有得到所需的输出。。。我想不出有什么不对。。。它可以计算积分值，但它不能。。。另外，点数是错误的…它只是取第一个数，不算点数。。。请帮忙。。。

最佳答案

我认为你离解决方案不远。这个公式至少看起来不错。
也许你在代码中犯的最大错误是你在数据中丢失了要读取的点数。所以你的代码可能读“0.5”作为点数。然后，k上的循环就是k=0（然后k=1>0.5），这可能就是为什么你只有一个点。为了让它工作，我做了以下更改：
在数据文件的开头添加点数。
更改int size的大小类型
打印积分值printf ("Integral: %lg \n", integral);
这对我很有用。
（大编辑，因为它已经被作为C而不是C++被延迟）

关于c - 与c的数值积分，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/14603874/