魔法少女の終わらない戦い

MikuHatsune2011-07-15

ここを見ていたら無性にやりたくなったこちらのネタ

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()
#Data<- madoka(M=1000,W=1000,C=1000,D=0.2,B=0.1,K=0.3,F=0.4,Time=120)