ita https://blog.hatena.ne.jp/ita/ ita’s diary https://ita.hatenadiary.jp/ 理系ネタ J これを数値的に解いてみる。http://d.hatena.ne.jp/ita/20071218/p1 長さLの一次元系があり、両端が時間Tだけずれてつながっている。差分法で単純なオイラー法, f[x] +=dt*dfdt(f[])を考える。 拡散方程式 方程式は。差分だと #define XNUM 1000 #define dt 0.001 #define dx (L/XNUM) double f[XNUM]; f[x] += dt* (2*f[x]-f[x-1]-f[x+1])/(dx*dx);普通はt=0の初期値を設定して時間発展させる。しかし両端の時間がずれていると未来の値が必要になり計… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fita.hatenadiary.jp%2Fentry%2F20071227%2Fp1" title=" タイムトンネルの数値計算 - ita’s diary" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/i/ita/20071229/20071229125140.png Hatena Blog https://hatena.blog 2007-12-27 00:00:00 rich https://ita.hatenadiary.jp/entry/20071227/p1 1.0 100%