#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")