<Matlab-3:追赶法（Doolittle分解）工具箱

function x=chase (d,e,f,b) 
% ---------------------------------------------------------------
%the method of chaase*******************************
%  d, e, f,分别是是方程组的下对角，主对角，上对角的三条
%  l，r， f分别是Doolittle分解的主和次对角线
%    /1          0    ...             ...        0  0     \
%    |l2               1  ...          ...       0  0     | 
%  L=|0         l3              ...1                      |
%    |                                                    |
%    |                                                    |
%    \                                   ln            1  /
%------------------------------------------------------------------
r(1)=e(1);y(1)=b(1);
for i=2:length(e)
l(i-1)=d(i-1)/r(i-1);r(i)=e(i)-l(i-1)*f(i-1);
y(i)=b(i)-l(i-1)*y(i-1);
end
x(length(e))=y(length(e))/r(length(e))
for j=length(e)-1:-1:1
x(j)=(y(j)-f(j)*x(j+1))/r(j)
end