c----------------------------------------------------------------------
	subroutine DENSFILL(T,Vee,Vext,Eion,ecusp)
c fills the density in, which one determined by several parameters.
c
c author: Kieron Burke, 8.14.99
c     in: alpha -
c         num   -
c	  ncorr -
c	  np    - 
c	  next  -
c	  nmat  -
c	  cyrus.Z, if next=1,num=2,ncorr=1 -
c    out: dens(0:n)
c	  T - exact kinetic energy
c	  Vee - exact e-e repulsion
c	  Eion - the difference: E(N)-E(N-1)
c	  ecusp - err in e-e cusp if next=2,num=2,ncorr=1
c calls: DSPLINE,DSPLINT,HOODL,DLKAIS
	include 'basics.inc'
	include 'con.inc'
	include 'dens.inc'
	include 'sys.inc'
	dimension rc(m),rhol(m),rho2(m)