c----------------------------------------------------------------------
	subroutine VHART(u)
c finds the Hartree potential from the density, and calculates the
c Hartree energy
c
c author: Kieron Burke, 8.14.99
c     in: dens(0:n)
c    out: vh(0:n)
c  calls: RNTGRL,RNGTR2
	include 'basics.inc'
	include 'dens.inc'
	include 'pot.inc'
	dimension f(0:m),g(0:m),rad(0:m),r2(0:m),dum(0:m)