CC****************************************************************
      PROGRAM whis2 
CC****************************************************************
      IMPLICIT REAL*8(a-h,o-z)
      OPEN(21,FILE='whis2.dat')
      READ(21,*)xlambda,W
      READ(21,*)N,icada
      READ(21,*)s,t
      CLOSE(21)
      OPEN(31,FILE='whis2.out')
      dospi=8.d0*DATAN(1.d0)
      pi=dospi/2.d0
      delta0=1.d-9
      s2=s-delta0
      t2=t+delta0
      G=xlambda*DLOG(32.d0/W)
      DO 4 i=1,N
      xi=i
      s=s-DSIN(t)
      t=t-xlambda*DLOG(ABS(s))+G
      s2=s2-SIN(t2)
      t2=t2-xlambda*DLOG(ABS(s2))+G
       t=MOD(t,dospi)
       t2=MOD(t2,dospi)
       IF(t.lt.0.d0) t=t+dospi
       IF(t2.lt.0.d0) t2=t2+dospi
       IF(t.GT.pi)t=t-dospi
       IF(t2.GT.pi)t2=t2-dospi
       delta=dsqrt((s-s2)**2+(t-t2)**2)
       IF(mod(i,icada).eq.0) WRITE(31,*)i,delta
    4 CONTINUE
      END

