code 34390;
    procedure CHLSOL2(A, N, B); value N; integer N; array A, B;
    begin integer I;

        for I:= 1 step 1 until N do 
        B[I]:= (B[I] - TAMVEC(1, I - 1, I, A, B)) / A[I,I];
        for I:= N step - 1 until 1 do 
        B[I]:= (B[I] - MATVEC(I + 1, N, I, A, B)) / A[I,I]
    end CHLSOL2;
        eop