とある講義の課題で作った代用電荷法で2次元ラプラス方程式を解くプログラムのメイン部分だけを取り出したものです.
importしているライブラリは
- math
- matplotlib
- numpy
- scipy
です.
2次元ラプラス方程式に対して,単位円上に拘束点を取った場合に解く計算の本体です. 実行後に標準入力から与えられる変数Nを拘束点と特異点の数として等間隔に設置します.
2次元ラプラス方程式を(-1,-1),(1,-1),(1,1),(-1,1)の四角形上に拘束点を取った場合に解く計算の本体です. 実行後に標準入力から与えられる変数Nを拘束点と特異点の数として等間隔に設置します.
境界条件をこの中で定義しています.
数値解析に用いるためのグリーン関数の定義モジュールです. 他のモジュールからimportしていない1次元,3次元の場合のグリーン関数も定義しています.
- 村島定行,「代用電荷法とその応用」,1983,森北出版