chenaoki
commented
5 years ago 下記のanisotropy_ratioは「心臓のフィジオーム」中のmuだと思うけど，p175の定義ではsigma_i = mu * sigma_eとなっているので定義が違うのでは？
anisotropy_ratio = sigma_l_i / sigma_t_i
rhs_vmem += anisotropy_ratio * i_ext_e
rhs_vmem *= 1 / (1 + anisotropy_ratio)
このモノドメインモデルをテストしたところ， タイムステップは通電刺激シミュレーション時で0.001ms, それ以外では0.005ms程に変更可能． タイムステップをそれ以上(0.01msとか)にすると計算が発散したことを確認．