delta rho = 1-mu|x|^r