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