崩壊過程

R

KABIRA2011-05-13

  • この記事の2.数学的な扱い方に忠実なモデル
  • 原子崩壊の過程
    • 崩壊していない原子数 N
    • t = 0 の時の原子数 N_0
    • 崩壊時間の確率分布 p(t) = e^{-\frac{t}{\tau}}
  • 一つ一つの原子 X[i] が崩壊する時刻を乱数 (rexp) で与える
  • 残存している原子数 sum(X > t[i]) の経時変化をプロットする
  • 赤い線は期待値
N0<-100
tau<-10
X<-rexp(N0,rate=1/tau)
#hist(X,breaks=N0)

t<-seq(0,100,by=0.1)
N<-c()
for(i in 1:length(t)){
	N[i]<-sum(X>t[i])
	}
	
plot(N0*exp(-t/tau),col=2)
par(new=TRUE)
plot(N)