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