function humps ( x )

      real humps
      real x

      humps = 1.0 / ( (x-0.3)**2 + 0.01 ) 
     &      + 1.0 / ( (x-0.9)**2 + 0.04 ) - 6.0

      return
      end