c----------------------------------------------------------------------
	function DLHOOKE(n,d,r)
c constructs the log of the reduced density for a given solution to 
c Hooke's atom.
c
c author: Kieron Burke, 8.14.99
c     in: n - number of gridpoints ?
c         d - expansion coefficients (calculated with HOOKE)
c         r - spatial point 
c    out: return value of function: log of reduced density at point r
c  calls: DERFX
	implicit real*8(a-h,o-z)
	parameter(nm=150)
	dimension d(0:2*nm),p(0:2*nm+1),q(0:2*nm)
c