時間について法で差分近似した方程式は
(6) |
これに対する弱形式は以下の通り。
FreeFem++ で対応するコードは
heatB.edp 内の problem |
problem heat(u,v,init=n)= int2d(Th)(u*v)-int2d(Th)(uold*v) +int2d(Th)(tau*(dx(u)*dx(v)+dy(u)*dy(v))) -int2d(Th)(tau*f*v)-int1d(Th,2,3)(tau*g2*v) +on(1,4,u=g1);(ラベル , が を, ラベル , が を構成している。) |
このコードが弱形式 (7) の実現であることは読み取ろう。