PROGRAM SCKN3
      COMMON T,T1,T2,E1(4)
      T1=0.50025
      T2=2.0
      T=.0499975
      E1(1)=1.0
      E1(2)=-1.0
      E1(3)=-1.0
      E1(4)=-1.0
      I=10000
      N3=14*I
      IF(N3)39,39,31                                                    SCK   46
   31 DO 38 I=1,N3,1                                                    SCK   47
   38 CALL PA(E1)                                                       SCK   48
   39 CONTINUE                                                          SCK   49
      CALL POUT(N3,N2,N2,E1(1),E1(2),E1(3),E1(4))                       SCK   50
      END
      SUBROUTINE PA(E)
      COMMON T,T1,T2
      DIMENSION E(4)
      J=0
    1 E(1)=(E(1)+E(2)+E(3)-E(4))*T
      E(2)=(E(1)+E(2)-E(3)+E(4))*T
      E(3)=(E(1)-E(2)+E(3)+E(4))*T
      E(4)=(-E(1)+E(2)+E(3)+E(4))/T2
      J=J+1
      IF(J-6)1,2,2
    2 CONTINUE
      RETURN
      END
      SUBROUTINE POUT(N,J,K,X1,X2,X3,X4)
         PRINT *,N,J,K
         PRINT *,X1,X2,X3,X4
         RETURN
      END