c        1         2         3         4         5         6         7
c23456789012345678901234567890123456789012345678901234567890123456789012
c        1         2         3         4         5         6         7

c----------------------------------------------------------------------
c TESTDENS tests a density.
c Kieron Burke, 5.12.99
	
	implicit real*8(a-h,o-z)
	parameter(m=1000000)
	dimension f(0:m)
	pi4=16d0*datan(1d0)

	write(6,*) 'TESTING expdens'
	write(6,*) ' n=no of points'
	do n=20,80,20
	  do i=1,n
	    y=dreal(i)/dreal(n)
	    r=1d0/y-1d0
	    f(i)=pi4*r*r*expdens(r)
	  enddo
	  ans=rntgrl(n,1d0,f)-1d0
	  write(6,*) ' n=',n,' ans=',ans
	enddo
	end
c----------------------------------------------------------------------