Model

\[P_{(t)} = P_{(0)}e^{rt}\]

Malthusian Growth Model

library(ggplot2)
n = 25
P = rep(0,n)
P[1] = 1000
N = seq(from = 1, to = n, by = 1)
e = exp(1)
r = 0.2
nmax = n - 1
for (t in 1:nmax) P[t+1] <- P[t]*e^r
popsim <- data.frame(N,P)
ggplot(data=popsim, aes(x=N, y=P)) +
  geom_line()

NA
LS0tDQp0aXRsZTogIk1hbHRodXNpYW4gR3Jvd3RoIE5CIg0KYXV0aG9yOiAiUm9uYWxkbyBCYWx0YXIiDQpkYXRlOiAnYHIgU3lzLkRhdGUoKWAnDQpvdXRwdXQ6DQogIGh0bWxfbm90ZWJvb2s6IGRlZmF1bHQNCiAgaHRtbF9kb2N1bWVudDogZGVmYXVsdA0KLS0tDQoNCg0KIyBNb2RlbA0KDQoNCiQkUF97KHQpfSA9IFBfeygwKX1lXntydH0kJA0KDQojIE1hbHRodXNpYW4gR3Jvd3RoIE1vZGVsDQpgYGB7ciBjb2RlMSwgZWNobz1UUlVFfQ0KbGlicmFyeShnZ3Bsb3QyKQ0KDQpuID0gMjUNClAgPSByZXAoMCxuKQ0KUFsxXSA9IDEwMDANCk4gPSBzZXEoZnJvbSA9IDEsIHRvID0gbiwgYnkgPSAxKQ0KZSA9IGV4cCgxKQ0KciA9IDAuMg0KDQpubWF4ID0gbiAtIDENCg0KZm9yICh0IGluIDE6bm1heCkgUFt0KzFdIDwtIFBbdF0qZV5yDQoNCnBvcHNpbSA8LSBkYXRhLmZyYW1lKE4sUCkNCg0KZ2dwbG90KGRhdGE9cG9wc2ltLCBhZXMoeD1OLCB5PVApKSArDQogIGdlb21fbGluZSgpDQogIA0KYGBgDQo=