c++ - C++中的数值积分

我需要集成一个函数(两个变量)。
我知道我可以通过使用Fubini定理来集成一个变量函数，然后使用诸如Rectangle方法或Trapezoidal规则之类的数值方法来做到这一点。

但是，在C++中是否有任何预建函数可以做到这一点？我需要对单位R2三角((0,0), (1,0), (0,1))进行集成。

最佳答案

您可以使用GNU Scientific Library，它支持许多“数值分析”功能，包括集成。

手册中非常简单的example of integration只是几行代码:

#include <stdio.h>
#include <math.h>
#include <gsl/gsl_integration.h>

double f (double x, void * params) {
  double alpha = *(double *) params;
  return log(alpha*x) / sqrt(x);
}

int
main (void)
{
  double result, error;
  double expected = -4.0;
  double alpha = 1.0;
  gsl_integration_workspace * w
    = gsl_integration_workspace_alloc (1000);

  gsl_function F;
  F.function = &f;
  F.params = &alpha;

  gsl_integration_qags (&F, 0, 1, 0, 1e-7, 1000,
                        w, &result, &error);

  printf ("result          = % .18f\n", result);
  printf ("exact result    = % .18f\n", expected);
  printf ("estimated error = % .18f\n", error);
  printf ("actual error    = % .18f\n", result - expected);
  printf ("intervals =  %d\n", w->size);

  gsl_integration_workspace_free (w);

  return 0;
}

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