subroutine foo(x)
      double precision r1mach
      x = sqrt(r1mach(3))
      end