拡散の係数kの値を自由に与える 昨日の記事の係数(k1,..,k4)を各点で与える行列を作る Nx<-30 Ny<-30 Nt<-30 U<-tempU<-matrix(0,Nx,Ny) U[10:20,10:20]<-1 k1<-matrix(0.2,Nx-1,Ny) k2<-matrix(0.2,Nx,Ny-1) k3<-k4<-matrix(0.1,Nx-1,Ny-1) k1[1:ceiling(Nx…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。