2010-11-17から1日間の記事一覧

4行で書くブラウン運動のパス

R

Rだとサクっとかけるという話。実質1行。 N <- 10 T <- 100 path <- apply(matrix(rnorm(N*T),nrow=T,ncol=N),2,cumsum) matplot(path,type="l",ylab="Value",xlab="Time") 実行すると な感じで図が出力。path変数には各経路の軌跡がはいってる。 くそう、p…