dM = slash(rho*4;3)*pi*(r_c^3-(r_c-dr)^3)

dr= dM/(4*rho*pi*r_c^2)

dr=(1,5*10^11)^2*(C_soleil/(4*rho*L*a^2))*dt