Crescimento logarítmico de uma população com ou sem vacas

Sabeis que:

dlogistic<-function(K, rd, N0, t)

{

N <- c(N0, numeric (t))

for(i in 1:t)

N [i+1] <- N [i]+rd*N[i]*(1-N[i]/K)

N

}

t<-30

rd=1.5

K=100

n<- dlogistic(K, rd, 1, t)

plot(0:t, n, type="l")

growth <-n[2:201]/n[1:200]

plot (n[1:200], growth)

dn<-n[2:201]-n[1:200]

plot(n[1:200],dn)

dá:

E esta função dá CAOS não COWS! (segundo o H, se nunca vos esquecerdes desta piada nunca mais terei problemas com estas distribuições!)

E isto é Ecologia Teórica!