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