<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