code 34424;
  procedure SOLTRI(SUB, DIAG, SUPER, N, B);
  value N; integer N; array SUB, DIAG, SUPER, B;
  begin integer I;
     real R;
     R:= B[1]:= B[1] / DIAG[1];
     for I:= 2 step 1 until N do 
     R:= B[I]:= (B[I] - SUB[I - 1] * R) / DIAG[I];
     for I:= N - 1 step -1 until 1 do 
     R:= B[I] := B[I] - SUPER[I] * R
  end SOLTRI;
        eop