	implicit none
	real s
	external f

	call qtrap(f, 0.0, 1.0, s)
	write (*,*) s
	call qsimp(f, 0.0, 1.0, s)
	write (*,*) s
	end

	function f(x)
	implicit none
	real f,x
	f=x**2+exp(x**2)
	end
