三对角线方程组追赶法

function x=chasing_method(A,d)
    %{
        追赶法，适用于三对角线矩阵
        param A：三对角线矩阵
        param d：该求解的列向量
        return x：解向量
    %}
    rows=size(A,1);

    [x,y,b,c]=deal(zeros(rows,1));

    %先确定b,c,y的值
    b(1)=A(1,1);
    y(1)=d(1)/b(1);

    for ii=2:rows
        c(ii)=A(ii-1,ii)/b(ii-1);
        b(ii)=A(ii,ii)-A(ii,ii-1)*c(ii);

        y(ii)=(d(ii)-A(ii,ii-1)*y(ii-1))/b(ii);
    end
    
    % %确定y的值
    % y(1)=d(1)/b(1);
    % for ii=2:rows
    %     y(ii)=(d(ii)-A(ii,ii-1)*y(ii-1))/b(ii);
    % end

    %确定x的值
    x(rows)=y(rows);
    for ii=rows-1:-1:1
        x(ii)=y(ii)-c(ii+1)*x(ii+1);
    end
end