
      double precision one
      parameter (one=1.0D0)
      
