MINICERN/mathlib/gen/f/figi.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:34  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE FIGI(NM,N,T,D,E,E2,IERR)
  11       INTEGER I,N,NM,IERR
  12       REAL T(NM,3),D(N),E(N),E2(N)
  13       IERR = 0
  14       DO 100 I = 1, N
  15          IF (I .EQ. 1) GO TO 90
  16          E2(I) = T(I,1) * T(I-1,3)
  17          IF (E2(I)) 1000, 60, 80
  18    60    IF (T(I,1) .EQ. 0.0 .AND. T(I-1,3) .EQ. 0.0) GO TO 80
  19          IERR = -(3 * N + I)
  20    80    E(I) = SQRT(E2(I))
  21    90    D(I) = T(I,2)
  22   100 CONTINUE
  23       GO TO 1001
  24  1000 IERR = N + I
  25  1001 RETURN
  26       END