code 34261;
procedure PSTTFMMAT(A, N, V, B);
value N; integer N; array A, V, B;
begin integer I, I1, J;
    real H;

    I1:= N; V[N,N]:= 1;
    for I:= N - 1 step - 1 until 1 do 
    begin H:= B[I] * A[I,I1]; if H < 0 then 
        begin for J:= I1 step 1 until N do V[J,I]:= A[I,J] /
            H;
            for J:= I1 step 1 until N do 
            ELMCOL(I1, N, J, I, V, V, MATMAT(I1, N, I, J, A, V))
        end;
        for J:= I1 step 1 until N do V[I,J]:= V[J,I]:= 0;
        V[I,I]:= 1; I1:= I
    end 
end PSTTFMMAT;
        eop