This will contain code to plot the density profile rho = rho(r) when initial rho is 10.