写一个用pthread实现梯形积分法的实例

好的，下面是一个使用pthread库实现梯形积分法的示例代码：

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>

#define NUM_THREADS 4
#define N 1000000

double a = 0.0, b = 1.0, h = (b - a) / N;
double result = 0.0;
pthread_mutex_t mutex;

void *worker(void *arg) {
    long id = (long)arg;
    double local_result = 0.0;

    for (long i = id; i < N; i += NUM_THREADS) {
        double x = a + i * h;
        double fx = 1 / (1 + x * x);
        double x1 = a + (i + 1) * h;
        double fx1 = 1 / (1 + x1 * x1);
        local_result += (fx + fx1) * h / 2;
    }

    pthread_mutex_lock(&mutex);
    result += local_result;
    pthread_mutex_unlock(&mutex);

    pthread_exit(NULL);
}

int main() {
    pthread_t threads[NUM_THREADS];
    pthread_mutex_init(&mutex, NULL);

    for (long i = 0; i < NUM_THREADS; i++) {
        pthread_create(&threads[i], NULL, worker, (void*)i);
    }

    for (long i = 0; i < NUM_THREADS; i++) {
        pthread_join(threads[i], NULL);
    }

    printf("The result is: %f\n", result);

    pthread_mutex_destroy(&mutex);
    pthread_exit(NULL);
}

该程序将区间 [0, 1] 分成 N 个小区间，每个线程计算一部分的积分值，并将结果累加到全局变量 result 中。为了避免多个线程同时修改 result 变量导致竞争条件，使用互斥锁来保护该变量。最后，主线程输出积分结果。

希望这个示例能够帮助你理解如何使用 pthread 库实现梯形积分法。