c
c   fx.f  91/21/11 
c
c   Function to integrate in the integration example.
c
	double precision function f(x)
	double precision x

	f = 4.D0 / (1.D0 + x*x)

   	return 
	end
