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（