program 11
    real(kind=n)::aa,bb,cc,dd,mm
    aa(1)=0,aa(2)=1,aa(3)=2,bb(1)=3,bb(2)=4,bb(3)=5,cc(1)=1,cc(2)=1,cc(3)=0
    real n=3
    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)   !实施变换
    write(*,*) mm(1) 
    end do    

end