レヴィ飛行

• 浮動と言えばこちら
• 数セミの話題その３
• 動物の餌探索に関して、酔歩とその一種であるレヴィ飛行についての記事もある「べき則で食べる」
• Wiki記事はこちら

Ns<-100
XY<-matrix(0,Ns,2)
# 標準コーシー分布のパラメタ
a<-0
b<-1
Ls<-rcauchy(Ns-1,a,b)
ts<-runif(Ns-1)*2*pi
for(i in 2:Ns){
 XY[i,]<-XY[i-1,]+c(Ls[i]*cos(ts[i]),Ls[i]*sin(ts[i]))
}
plot(XY,type="b")