*表題 非回転2次元浅水重力波の位相速度, 固有値を計算する
*
*履歴 2002/03/15 Hiroshi Taniguchi

      integer n
      real c
      real l
      real omega
      real k
      real pi

      parameter(pi=3.1415927, l=3.0)

      write(*,*) 'n=?'
      read(*,*) n
      write(*,*) 'k=?'
      read(*,*) k

      c     = sqrt( 1.   + ( (n*pi)**2 )/( (l*k)**2 ) )
      omega = sqrt( k**2 + ( (n*pi)**2 )/(     l**2 ) )

      write(*,*) 'c =±',c
      write(*,*) 'ω=±',omega

      stop
      end