Newmark-β法で1次元の波動方程式

「Newmark-β法で1次元の波動方程式を読む」 (github.io)

Schweickart, James, Marschner の "Animating Elastic Rods with Sound" によると、Backward Euler法で弾性ロッドの音のシミュレーションを行うと、ダンピングによって受け入れられないほどミュートされた音が出るそうです。そこで論文で紹介されていたNewmark-β法から立てた画像の方程式について離散化と実装を行いました。

画像の一番上の式は、ばねとダンパの項を加えた波動方程式です。