用matlab写出追赶法,不能用内置函数

追赶法是一种解三对角线方程组的方法，也称为托马斯算法。以下是用MATLAB实现追赶法的示例代码：

function x = tridiagonal_solver(a, b, c, d)
    n = length(d);
    x = zeros(n, 1);
    
    c_hat = zeros(n, 1);
    d_hat = zeros(n, 1);
    
    c_hat(1) = c(1) / b(1);
    d_hat(1) = d(1) / b(1);
    
    for i = 2:n
        denominator = b(i) - a(i) * c_hat(i-1);
        c_hat(i) = c(i) / denominator;
        d_hat(i) = (d(i) - a(i) * d_hat(i-1)) / denominator;
    end
    
    x(n) = d_hat(n);
    for i = n-1:-1:1
        x(i) = d_hat(i) - c_hat(i) * x(i+1);
    end
end

% 示例使用：
a = [0; -1; -1];  % 下对角线
b = [3; 4; 4];    % 主对角线
c = [-1; -1; 0];  % 上对角线
d = [7; 14; 13];  % 右侧向量

x = tridiagonal_solver(a, b, c, d);
disp(x);

在上述代码中，`tridiagonal_solver` 函数接受四个输入参数 `a`、 `b`、 `c` 和 `d`，分别代表下对角线、主对角线、上对角线和右侧向量的系数。函数中使用循环计算出副对角线和待求解的中间结果，然后再反向求解出最终的解向量 `x`。

请确保输入的参数 `a`、 `b`、 `c` 和 `d` 的长度相同，并且 `b` 中的元素不为零，否则可能会引发错误。