追赶法求解三对角方程的MATLAB程序代码

function [k]=zhuigan（abcd）
%对角线下方的元素，个数比A少一个
% %对角线元素
%对角线上方的元素，个数比A少一个
%d为方程常数项
%用追赶法解三对角矩阵方程
r=size（a）;
m=r（2）;
r=size（b）;
n=r（2）;
if size（a）~=size（c）|m~=n-1|size（b）~=size（d）
    error（‘变量不匹配检查变量输入情况!‘）;
end
%%
%LU分解
u（1）=b（1）;
for i=2:n
    l（i-1）=a（i-1）/u（i-1）;
    u（i）=b（i）-l（