code: def profileMain(self, M, z):
          
          c = self.c_M_z(M, z)
          r200 = self.r200_M(M)
          rho0 = self.rho0_c(c)
          Rs = r200/c
          return r200, rho0, c, Rs
docs: returns all needed parameter (in comoving units modulo h) to draw the profile of the main halo
          r200 in co-moving Mpc/h
          rho_s in  h^2/Mpc^3 (co-moving)
          Rs in Mpc/h co-moving
          c unit less