如何用c语言实现追赶法解三对角方程组

追赶法（也称为托马斯算法）是一种用于解三对角线性方程组的方法，可以用C语言实现。以下是一种基于追赶法的C语言代码示例：

#include <stdio.h>

void thomas(double* a, double* b, double* c, double* d, int n)
{
    int i;
    double m;

    // 前向消元
    for (i = 1; i < n; i++) {
        m = a[i] / b[i - 1];
        b[i] = b[i] - m * c[i - 1];
        d[i] = d[i] - m * d[i - 1];
    }

    // 回代求解
    d[n - 1] = d[n - 1] / b[n - 1];
    for (i = n - 2; i >= 0; i--) {
        d[i] = (d[i] - c[i] * d[i + 1]) / b[i];
    }
}

int main()
{
    double a[3] = {1, 2, 3};
    double b[3] = {4, 5, 6};
    double c[3] = {7, 8, 9};
    double d[3] = {10, 11, 12};

    thomas(a, b, c, d, 3);

    printf("Solution:\n");
    for (int i = 0; i < 3; i++) {
        printf("%f\n", d[i]);
    }

    return 0;
}

在这个示例中，我们定义了一个名为thomas的函数来实现追赶法。该函数需要传入四个数组：a、b、c和d，它们分别代表三对角线性方程组中的对角线、下对角线、上对角线和常数项。n参数表示方程组的大小。在函数中，我们首先进行前向消元来消除下对角线，然后进行回代求解来计算未知变量的值。最后，我们在主函数中调用thomas函数来解方程组，并输出解。