program testp
	implicit double precision (a-h, o-z)
	external f

c lower limit of integration
        a=0.0d+00
c upper limit of integration
        b=1.0d+00
        do 10 l=1, 3000

	call dqk15(f, a, b, result, abserr, resabs, resasc)

	call dqk21 (f, a, b, result, abserr, resabs, resasc)

	call dqk31 (f, a, b, result, abserr, resabs, resasc)

	call dqk41 (f, a, b, result, abserr, resabs, resasc)

	call dqk51 (f, a, b, result, abserr, resabs, resasc)

	call dqk61 (f, a, b, result, abserr, resabs, resasc)

   10   continue

        stop
 	end