(a) Use a numerical integration algorithm such as the Euler method to compute the density profile for the n=1.5 and n=3 polytropes. Be sure to correctly incorporate the boundary conditions in your integrations.
(b) Plot your results and compare them with the n=0, n=1, and n=5 analytic models determined in Problem 10.20.