こちらのつづき kの値を各位置について与える まず を各位置で与え、その から k を計算する Nx<-30 Ny<-30 Nt<-100 U<-tempU<-matrix(0,Nx,Ny) U[11:20,11:20]<-1 #初期分布 #これがl方向とt方向の拡散の定数 kl<-0.20 kt<-0.05 #線維走向の傾き theta<-mat…