ここを見ていたら無性にやりたくなったこちらのネタ。
madoka<- function(M,W,C,D,B,K,F,Time){
data<- matrix(0,Time,2)
data[1,1]<- M; data[1,2]<- W
for(t in 2:Time){
t.counter<- t
W<- W + F*W + (B-K)*M
M<- M + C - (D+B)*M
if((W<0)||(M<0)){
break
}
data[t,1]<- M; data[t,2]<- W
}
return(data[1:(t.counter-1),])
}
Data<- madoka(M=1000,W=1000,C=1000,D=0.1,B=0.1,K=0.3,F=0.4,Time=200)
png("battle.png")
matplot(Data,type="l",frame=FALSE,xlab="Time",ylab="Population")
title("Magic girls and Witches battle")
legend(1,max(Data)/2,c("magic girls","witches"),col=1:2,lty=1:2)
dev.off()