double rho_ZHAO	(	double &	r,
		double	par[5]
	)

Definition at line 180 of file profiles.cc.

{
   //--- Returns rho(r) for a ZHAO (alpha,beta,gamma) [Msol/kpc^3]
   //  r         Distance from the centre of the profile [kpc]
   //  par[0]    rho_s: density normalisation [Msol/kpc^3]
   //  par[1]    r_s: scale radius [kpc]
   //  par[2]    alpha: transition slope
   //  par[3]    beta: outer slope
   //  par[4]    gamma: inner slope
   //
   //    => rho(r)= rho_s/[(r/r_s)^{gamma} [1+(r/r_s)^{alpha}]^{(beta-gamma)/alpha}]


   // Required to avoid divergences at r=0
   if (par[4] > 1.e-5) {
      if (r > par[1] * pow(par[0] / gDM_RHOSAT, 1. / par[4])) {
         return par[0] * pow(par[1] / r, par[4])
                / pow(1. + pow(r / par[1], par[2]), (par[3] - par[4]) / par[2]);
      } else
         return gDM_RHOSAT;
   } else
      return par[0] / pow(1. + pow(r / par[1], par[2]), par[3] / par[2]);
}