天天看點

龍格庫塔C 語言程式設計實作

最新一直在學一門很苦逼的課程《數值分析》,哎,高等數學沒學好現在後悔了呀,哎回來再惡補,現在說正事。

龍格庫塔很牛掰的名字,是兩個國外數學家的名字的合并,應該又是兩個大牛。

這個方法主要是用來解決微分方程的解,大體思路就是用差分代替微分。細節我就不說了,大家可以去百度,我這裡直接給出龍格庫塔家族的通式:

龍格庫塔C 語言程式設計實作

這個是龍格庫塔家族的通式,如果我們取累加和中的r為4則可以得到如下公式:

龍格庫塔C 語言程式設計實作

這就是四階龍格庫塔方程,據說這裡的精度很好,可以精确到小數點後6位。

方程組中的h為差分的步長。

下面給出c語言實作的代碼:

結果為:

龍格庫塔C 語言程式設計實作

繼續閱讀