■計算を続行■ とにかく、プログラムを実行してみて下さい。 するとファイルSOLUTION.BEMが出力されます。 このデータをベースにして u(0,y)を表計算ソフトで作図すると下図が得られます。 図から解析解としてΔP=0.6983 がえられます。 この場合の厳密解は ΔP=0.6499 です。 わりと荒っぽい要素分割ですが まーまーな数値解が得られていますね。 なにせ u(x) の最大値がdoublet の中にありますから計算には気をつかいます。
