追赶法的简单MATLAB代码

%追赶法
function [x]=zhuiganfa（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（