用c语言用追赶法解三对角方程

好的，以下是用C语言实现追赶法解三对角方程的代码：

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

void tridiagonal_solve(int n, double *a, double *b, double *c, double *f, double *x) {
    // 下三角矩阵L和上三角矩阵U
    double *L = (double*)malloc(n * sizeof(double));
    double *U = (double*)malloc(n * sizeof(double));

    // 前向追赶法求解Ly = f
    L[0] = a[0];
    x[0] = f[0] / L[0];

    for (int i = 1; i < n; i++) {
        L[i] = a[i] - b[i] * c[i-1] / L[i-1];
        x[i] = (f[i] - b[i] * x[i-1]) / L[i];
    }

    // 后向追赶法求解Ux = y
    U[n-1] = c[n-1];
    x[n-1] = x[n-1] / U[n-1];

    for (int i = n-2; i >= 0; i--) {
        U[i] = c[i] - b[i] * a[i+1] / U[i+1];
        x[i] = (x[i] - b[i] * x[i+1]) / U[i];
    }

    // 释放动态内存
    free(L);
    free(U);
}

int main() {
    int n = 5;  // 方程组的未知数个数
    double a[] = {2, 3, 4, 5, 6};  // 主对角线上的元素
    double b[] = {1, 1, 1, 1};  // 次对角线上的元素
    double c[] = {1, 1, 1, 1};  // 副对角线上的元素
    double f[] = {10, 20, 30, 40, 50};  // 方程组的右端向量
    double x[n];  // 方程组的解向量

    // 求解三对角线性方程组
    tridiagonal_solve(n, a, b, c, f, x);

    // 输出解向量
    for (int i = 0; i < n; i++) {
        printf("x[%d] = %lf\n", i, x[i]);
    }

    return 0;
}

以上代码中，`tridiagonal_solve()`函数实现了追赶法求解三对角线性方程组的过程。`n`为未知数个数，`a[]`、`b[]`、`c[]`分别为主对角线、次对角线和副对角线上的元素，`f[]`为方程组的右端向量，`x[]`为解向量。

在`main()`函数中，我们定义了一个五个未知数的三对角线性方程组，并通过调用`tridiagonal_solve()`函数求解出了方程组的解向量。最后，我们输出解向量的值。

希望以上代码对您有所帮助！