3 $\theta$ 法のコードを作る

弱形式は次のようになる。

$$\frac{1}{\Delta t}\left(u^{n+1}-u^{n},v\right) +\langle u^{n+\theta},v\rangle-(f,v)-[g_2,v]=0,\quad u^{n+\theta}:=\theta u^{n+1}+(1-\theta)u^n.$$ (8)

すなわち

$$\left(u^{n+1},v\right)-\left(u^{n},v\right) +\Delta t\theta\langl... ...gle +\Delta t(1-\theta)\langle u^{n},v\rangle -\Delta t(f,v)-\Delta t[g_2,v]=0.$$ (9)

次のコードがこの弱形式を実現したものであることは理解できるであろう。

2024/6/18 実習が進むまで秘密