This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / f / trbak1.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1996/04/01 15:02:38  mclareni
6 * Mathlib gen
7 *
8 *
9 #include "gen/pilot.h"
10       SUBROUTINE TRBAK1(NM,N,A,E,M,Z)
11       INTEGER I,J,K,L,M,N,NM
12       REAL A(NM,N),E(N),Z(NM,M)
13       REAL H,S
14       IF (N .EQ. 1) GO TO 200
15       DO 140 I = 2, N
16          L = I - 1
17          H = E(I) * A(I,L)
18          IF (H .EQ. 0.0) GO TO 140
19          DO 130 J = 1, M
20             S = 0.0
21             DO 110 K = 1, L
22   110       S = S + A(I,K) * Z(K,J)
23             S = S / H
24             DO 120 K = 1, L
25   120       Z(K,J) = Z(K,J) + S * A(I,K)
26   130    CONTINUE
27   140 CONTINUE
28   200 RETURN
29       END