18 lines
445 B
Plaintext
Executable File
18 lines
445 B
Plaintext
Executable File
#Iniciando Modelo SI para plantas vs babosas que se las comen
|
|
N<-1000
|
|
P<-seq(1:1000)
|
|
B<-seq(1:1000)
|
|
P[1]<-1000
|
|
B[1]<-100
|
|
x<-0.166
|
|
y1<-0.05
|
|
y2<-0.03
|
|
h<-1
|
|
for(i in 1:N){P[i+1]<-P[i]+h*(x*P[i]-y1*P[i]*B[i])
|
|
B[i+1]<-B[i]+h*(y2*P[i]*B[i])}
|
|
m<-max(S,I)
|
|
plot(P,type="l", ylim=c(0,m),xlab="Tiempo",ylab="Plantas")
|
|
lines(B,col="red")
|
|
text(8, 12, "Plantas que quedan",col="black")
|
|
text(8, 10, "Plantas ingeridas por la babosa",col="red")
|