Numerically simulating stellar pulsation using a "one zone" model

[patrykh18]

Homework Statement

Write a code that will simulate stellar pulsation using a "one zone" model, and plot various graphs

Relevant Equations

The relevant equations are given in the image attached

The question is taken from An Introduction to Modern Astrophysics by Carroll and Ostlie. I did manage to do the entire question and plot the relevant graphs but I just want to to investigate a bit more. For example I want to look at how the graph would like in the case of the Sun. I don't know what the right value for pressure I would need to input in order to get the right graph, with the right amplitude and period. The value for P in part d of the question seems arbitrary to me.

 

[Nate810]

Any help would be appreciated. Question: Consider a gaseous sphere in hydrostatic equilibrium. Assume that the gas is ideal, with no magnetic field and that the equation of state is given by P = Kρ^γ, where γ is a constant.(a) Derive an expression for the density as a function of radius.(b) Assume that the density at the center of the sphere is ρ_c and that the sound speed is given by c_s = (γP/ρ)^1/2. Obtain an expression for the frequency of radial pulsations of the sphere.(c) Show that for small amplitudes, the frequency is independent of the amplitude.(d) Assuming that P = P_0 cos (ωt), plot the behavior of the density as a function of r for two different values of ω.