- 拡散方程式をDOP853(ドルマンプリンス-ルンゲクッタ8次)で計算(Excel VBA)
du_i/dt = d (u_i+1 + u_i-1 - 2*u_i) / h^2
- maximaで微分方程式の解
dx/dt=-kx,t=0にてx(t)=aの解
maxima[return]
使い方は、? solve
?? desolve
終了はquit();
--- let's try. ---
'diff(x(t),t)=-k*x(t);
atvalue(x(t),t=0,a);
desolve(%o1,x(t);
x(t)=a %exp(-kt)と解けました
- 独り言
ay''+by'+cy=αexp(。x)・・・①
a,b,c,α,βは定数
|