2016-07-06から1日間の記事一覧

離散的な拡散

# 部屋の数 n <- 1000 # 追跡時刻数 n.t <- 1000 # 砂の量を記録する行列 X <- matrix(0,n.t,n) # 初期時刻の砂の量(高さ)をところどころに与える # m箇所に砂がある m <- n/10 # 砂のある場所 X[1,10:50] <- 10 X[1,200:300] <- 100 X[1,500:800] <- 50 # …