*
* Kiso Enshu 2002, Sample program #5
************************************************************************
      program plotauto
*
*     for drawing an autocorrelation function
*
      implicit none
*
      real        alambda
      integer     nmax, maxl
      parameter ( alambda = 0.2 )
      parameter ( nmax = 500, maxl=15 )
*
      real*4      acor( nmax+1 )
      real*4      x( nmax+1 )
      real*4      t
      integer     i, iws
*
      do i = 1, nmax+1
         x( i ) = float( i-1 )
      enddo
      acor( 1 ) = 1.
      do i = 2, nmax+1
         t = float( i )
         acor( i ) = exp( -t*alambda )
      enddo
*
*     Dennou
*
      write( *,* ) ' WORKSTATION ID (I)  ?; '
      call sgpwsn
      read( *,* ) iws
      call gropn( iws )
      call grfrm

      call usspnt( maxl, x, acor )
      call uspfit
      call grstrf

      call ussttl( 'Lag', ' ', 'Cor. Coe.', ' ' )
      call usdaxs

      call uuslnt( 1 )
      call uuslni( 9 )
      call uulin( maxl, x, acor )

      call grcls

      stop
      end