yumulog

北海道の大学教員/情報科学研究者の日記

続・ラプラス方程式

だんだんわかってきたぞ。



発散するのは z の下端で勾配が急になるから。それを補正すると次に z の上端の勾配も急になるのでそれも補正。



そして、次は電位のピークによって発散してしまうようだ。格子点を縦横2倍(計4倍)にして計算すると何とかなるんじゃないかと思い、現在走らせています。



明日起きたらちゃんと結果が出てたらいいですね。おやすみなさい。