anchor.diff:  !	-nx_isothermal: Calculates the WIMP density in the Spergel-Press scheme
   ! 	-Etrans_sp: calculates the WIMP transported energy (eps_x) given the WIMP temperature (Tx)
   !	-Tx_integral: to be used in newtons_meth
  -!	-newtons_meth: solves Tx_integral=0 which defines Tx 
  +!	-newtons_meth: solves Tx_integral=0 which defines Tx
   
   ! All units are cgs except tab_r and tab_dr
   ! I apologize for the long function calls.
  
positive.code: r = tab_r*Rsun ! cm
  phi = -tab_vesc**2/2.d0 ! erg/g
  mxg = mdm*1.782662d-24
  
  ! WIMP number density in isothermal approximation
  nx_isothermal = exp(-mxg*phi/kB/T_x) 
  n_0 = Nwimps/trapz(r, 4.d0*pi*r**2.d0*nx_isothermal, nlines) ! Normalize so that integral(nx) = Nwimps
  nx_isothermal = n_0*nx_isothermal
  
  
negative.code: module spergelpressmod
  use capmod
  implicit none
  
  double precision, parameter :: kB=1.38064852d-16, mnucg=1.6726219e-24
  
  contains
  