File: EX1COR.FT of Tape: Various/ETH/eth11-3

(Source file text) 

---

      DIMENSION ICOS(32),ISIN(32),ICOS2(32),ISIN2(32)
      DIMENSION ISTORE(64),COR(32)
      COMPLEX*8 ERRMES(2,3)

      DATA ERRMES/'LESS THA','N EIGHT ','EXCEEDS ','F.MAXN  ',
     1'NOT A PO','WER OF 2'/
      DATA PI/3.14159/,T/0./
      DT=PI/16.

      DO 1 I=1,32
      ICOS(I)=COS(T)*1000.
      ISIN(I)=SIN(T)*1000.
      ICOS2(I)=ICOS(I)
      ISIN2(I)=ISIN(I)
    1 T=T+DT

      CALL CORREL(IERROR,32,ISIN,ICOS,ISTORE,ISCALE)
      IF(IERROR.NE.0) GO TO 20
   10 TYPE 1000,ISIN
      FAC=2.**IABS(ISCALE)
      IF(ISCALE.LT.0) FAC=1./FAC
      DO 2 I=1,32
    2 COR(I)=FAC*ISIN(I)
      TYPE 1001,COR

      CALL CORREL(IERROR,32,ISIN2,ISIN2,ISTORE,ISCALE)
      IF(IERROR.NE.0) GO TO 20
      FAC=2.**IABS(ISCALE)
      IF(ISCALE.LT.0)FAC=1./FAC
      CSIN0=FAC*ISIN2(1)

      CALL CORREL(IERROR,32,ICOS2,ICOS2,ISTORE,ISCALE)
      IF(IERROR.NE.0) GO TO 20
      FAC=2.**IABS(ISCALE)
      IF(ISCALE.LT.0)FAC=1./FAC
      CCOS0=FAC*ICOS2(1)

      FAC=SQRT(CSIN0*CCOS0)
      DO 3 I=1,32
    3 COR(I)=COR(I)/FAC
      TYPE 1002,CSIN0,CCOS0,COR
      STOP

   20 TYPE 2000,(ERRMES(I,IERROR),I=1,2)
      STOP

 1000 FORMAT(///' UNSCALED RESULTS OF CORRELATION OF SINE AND COSINE'/,
     1 /,(4I17))
 1001 FORMAT(///' SCALED RESULTS OF CORRELATION OF SINE AND COSINE'//,
     1 (1P4E17.5))
 1002 FORMAT(///' VALUE OF AUTO-CORRELLATION OF SIN*1000 WITH ZERO SHIFT
     1 = ',1PE10.4,//,' VALUE OF AUTO-CORRELLATION OF COS*1000 WITH ZERO
     2 SHIFT = ',1PE10.4,//,' CROSS-CORRELATION COEFFICIENTS FOR SIN*100
     30 BY COS*1000 ARE AS FOLLOWS:',//,(1P4E17.5))
 2000 FORMAT(///' ****ERROR**** ARRAY LENGTH ',4A4)
      END

---