基于GSL,使用梯形法对正弦函数积分

首先需要引入GSL库，然后写出计算正弦函数积分的代码：

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

double integrand (double x, void * params) {
  return sin(x);
}

int main (void) {
  gsl_integration_workspace * w 
    = gsl_integration_workspace_alloc (1000);

  double result, error;

  gsl_function F;
  F.function = &integrand;
  F.params = NULL;

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

  printf ("result          = % .18f\n", result);
  printf ("estimated error = % .18f\n", error);

  gsl_integration_workspace_free (w);

  return 0;
}

在这个代码中，我们定义了一个integrand函数来计算被积函数sin(x)，并将其作为参数传递给gsl_function结构体的function成员。然后我们使用gsl_integration_qags函数来计算积分，其中第一个参数是gsl_function结构体，第二个和第三个参数是积分区间的下限和上限，第四个参数是积分误差上限，第五个参数是积分区间的细分数，第六个参数是gsl_integration_workspace结构体指针，用于存储计算过程中的中间结果，最后两个参数分别是积分结果和误差。