SUBROUTINE GRUNG2
      DOUBLE PRECISION T1,TD,ZO,EMA,PI,W,WC,X,Y
      DIMENSION X(50),Y(50)
      COMMON /GRUNGE/ I,TD,ZO,EMA,PI,W,WC
      COMMON /FIT/ X,Y
      T1=DCOS(W*TD*1D-6)
      X(I)=(T1-1.)*W*W
      Y(I)=T1+2.*W*W/(WC*WC)
      RETURN
      END