2013-04-27 153 views

回答

2

我不完全理解你的問題。 LSODE解決了dx/dt = f(x,t)[1]形式的任何泛型ODE,就像odeint中的方法一樣。然而,LSODE檢查問題的剛度並選擇正確的方案(預測校正器用於非硬性和僵硬問題的後向差異)。 odeint不提供這種自動選擇,而是依賴於用戶爲問題選擇正確的方案。就我所知,這兩種實現都沒有使用可能的現有分析解決方案。

[1] https://computation.llnl.gov/casc/odepack/odepack_home.html

0

大多數ODE求解器確實意味着某種超過原LSODE限制(在Fortran語言實現)。我知道的最接近的選擇是來自GNU的Octave。