造父變星型振子

造父變星型振子 ， 1955年天文學家W. S 。 Krogdahl在天體物理學學報上發表 《 天體的搏動是一種極限環現象 》 。 Krogdahl用下列非線性常微分方程組模擬造父變星的半徑隨時間的變化 ： 造父變星的半徑變化用下式表達 ： 此微分方程組無解析解 ， 但可以用龍格 － 庫塔法定性求得數值圖解 。 設參數params: = a = .1 * ( 1/3 ), b = .5 ， c = 10 ； 初始條件 ： init: = x ( 0 ) = 1.1, y ( 0 ) = 0