我用了如下程序将三对角  aa bb cc  进行了LU分解，
得到相应的下三角和上三角矩阵，
但不知道下面该如何利用这两个矩阵 求出 x 
请英雄指点一二。


subroutine tr-mat
    dd(1)=bb(1)                       !aa,bb,cc依次为三对角矩阵三条对角线元素，mm为单位下三角阵L次对角元素，dd为上三角阵U对角线元素,U非对角元素为cc。
    do i=2,n                          !列循环
       mm(i)=aa(i)/dd(i-1)            !计算乘数
       dd(i)=bb(i)-mm(i)*cc(i-1)      !实施变换
    end do