FitzHugh–Nagumo反応拡散モデルを解くためにAlternating direction implicit methodを実装することに取り組んでいます。私はblogでそれをPythonの実装例を発見したが、私はこの方法に誤りがあると思います - ここで紹介するステンシルで： 
それは fの反応用語を掛ける半分の時間ステップサイズすべきではありません？Pythonにおけるpdeの有限差分ソルバのための交互方向暗黙のメソッド
差分商によって差分商を取り付け、一方が元々ポーズタスクた式
U_t = D * (U_xx + U_yy) + f
されていない、両方の場合において
U_t = D/2 * U_xx + D/2 * U_yy + Δt*f
を取得します。
それはfためU_t、U_ppでD/(Δp^2)、p=x,yと1にあったように、係数は1/(Δt/2)でなければなりません。
数式は、差の商とそれが掛け合わされる次の段階とが混じっているようです。Δt/2。 fの要因はΔt/2でなければならないこと
そして、その次の式で1が実際にα_p=σ_p、p=x,yとして新しい定数を必要としません
http://math.stackexchange.com/でもっと答えが出るかもしれませんが、これは数学的な問題です。 – Daneel
またはhttp://scicomp.stackexchange.comには、大規模な科学計算の専門家がいます。 - また、分母 'D /（2 * dp^2）'の因子2であると思われ、リンクされたウィキペディアのページはこれを確認しているようです。 - そして、なぜdtがあるのでしょうか？他の言葉はすべて完全な差の商であり、無限小の要素はありません。 – LutzL