こちらのつづき 3dでの拡散のシミュレーション ナナメ方向も移動する 近傍のうち、中心間キョリが の18マスについて移動がある 線維方向は考えていない library(rgl) Nx<-35 Ny<-35 Nz<-35 Nt<-30 U<-tempU<-array(0,c(Nx,Ny,Nz)) #トーラスを作る A<-array(…