double precision function f(x)
	implicit double precision (a-h, o-z)

        f = 1.0d+00 + x * x + 1.0d+00 / (1.0d+00 + 1.0d+02 * x * x)

        return
	end