double precision function f(x)
      double precision x
      f = x*(x*x - 2.) - 5.
      return
      end

.