File: HHUK.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

      SUBROUTINE HHUK(A,IASEP,NR,SIG,BETA)
C 
C        PURPOSE
C           TO CALCULATE SIGMA,BETA AND THE VECTOR U(K) 
C 
      DIMENSION A(2)
C 
    9 SIG=0.
      JA = 1
      DO 10 I=1,NR
      SIG = SIG + A(JA)* A(JA)
   10 JA = JA + IASEP 
      SIG=SQRT(SIG) 
C 
C     GENERATE CORRECT SIGN OF PARAMETER SIG
   11 H=A(1)
      IF(H)12,13,13 
   12 SIG=-SIG
   13 BETA=H+SIG
      A(1) = BETA 
      BETA=1./(SIG*BETA)
  100 RETURN
      END 


---