CC****************************************************************
      PROGRAM smb 
CC  Calcula el mapa standard para varias c.i. en y 
CC****************************************************************
      IMPLICIT REAL*8(a-h,o-z)
      OPEN(21,FILE='smb.dat')
      READ(21,*)xk,n1
      READ(21,*)N,icada
      READ(21,*)x0,y0
      CLOSE(21)
      OPEN(31,FILE='smb.out')
      dospi=8.d0*DATAN(1.d0)
      xk1=xk/dospi
      dy=1.d0/n1
      WRITE(31,*)x0,y0
      DO 5 j=1,n1
      y=y0+(j-1)*dy;x=x0
      DO 4 i=1,N
      y=y+xk1*DSIN(dospi*x)
      x=x+y
       x=DMOD(x,1.d0)
       IF(x.lt.0.d0) x=x+1.d0
       y=DMOD(y,1.d0) 
       IF(y.LT.0.d0)y=y+1.d0 
      IF(MOD(i,icada).eq.0) WRITE(31,*)x,y
    4 CONTINUE
    5 CONTINUE
      END

