KABIRA https://blog.hatena.ne.jp/KABIRA/ KABIRAの日記 https://kabira.hatenadiary.org/ シミュレーション 拡散 拡散をシミュレーションしていて こちらやこちらでNeumann条件になっていなかった 境界でゼロフラックスになるようにする 端からの流出が起きないようにする そのため流入の起こるマスを制限する 計算過程で流入の起こるマスの数も数えておく ソースの中のkが流入の起こるマスの数 計算しているのは N<-30 M<-30 T<-100 d<-1 #キョリの次元 dt<-1 #時間の次元 D<-0.12 #拡散係数 U<-array(0,c(N,M,T)) U[5:10,5:10,1]<-1 #初期条件 #拡散 for(t in 1:(T-1)){ for(i in 1:N){ for(j in 1:… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fkabira.hatenadiary.org%2Fentry%2F20110102%2F1293464224" title=" 全量が保存されるNeumann条件 - KABIRAの日記" 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.d.st-hatena.com/diary/KABIRA/2011-01-02.gif Hatena Blog https://hatena.blog 2011-01-02 00:37:04 rich https://kabira.hatenadiary.org/entry/20110102/1293464224 1.0 100%