library("abind")
par(mfrow=c(1,1))
j<-c(1:100)
n1<-NULL
a<-function(x){
b<-(1*(x/100))*(1-(x/100))
c<-b*100
}
for (i in j){
if (1==i){ n1<-as.vector(a(10))
}else{m<-a(n1[i-(i-1)])
n1<-abind(m,n1,along=)
}
}
n2<-n1[100:1]
par(bg="white",bty="l")
plot(j,n2,main="Crecimiento poblacional",type="l",lwd=2.5,lty=2,
xlab="Generaciones",ylab="Tamaño Poblacional",pch=21,bg="grey")
#se me olvidó decirte que el modelo tiene capacidad de carga y es de 100 individuos.