- import numpy as np
- import matplotlib.pyplot as plt
- def rho(r, rho_c, R_star):
- return rho_c * ((r/R_star) - 1)**2
- def M(r, rho_c, R_star):
- integrand = lambda r_prime: 4*np.pi*r_prime**2 * rho(r_prime, rho_c, R_star)
- result, _ = np.trapz(integrand(np.linspace(0, r, 1000)), np.linspace(0, r, 1000))
- return result
- rho_c = 1.0 # Central density
- R_star = 1.0 # Radius of the star
- r_values = np.linspace(0, R_star, 100)
- M_values = [M(r, rho_c, R_star) for r in r_values]
- plt.plot(r_values, M_values)
- plt.xlabel('Radius (r)')
- plt.ylabel('Enclosed Mass (M(r))')
- plt.title('Enclosed Mass vs Radius')
- plt.grid(True)
- plt.show()
[text] ASTR
Viewer
*** This page was generated with the meta tag "noindex, nofollow". This happened because you selected this option before saving or the system detected it as spam. This means that this page will never get into the search engines and the search bot will not crawl it. There is nothing to worry about, you can still share it with anyone.
Editor
You can edit this paste and save as new: