c----------------------------------------------------------------------
	function EDIFF(rho,rhop,v) 
c finds the dE=E[rhop]-E[rho], using v=dE/drho, to first order.
c 
c author: Kieron Burke, 8.14.99
c     in: rho(0:n),rhop(0:n),v(0:n)
c    out: ediff
c  calls: RNTGRL
	include 'basics.inc'
	dimension rho(0:m),rhop(0:m),v(0:m),dum(0:m)