基于SymPy的指数增长模型

Question

我试图用SymPy实现指数增长模型。

n = Function('N')
t = symbols('t')
r = 1.7

eq = Eq(n(t).diff(t), r*n(t))
sol = dsolve(eq, ics={n(0): 10})
plot(sol.rhs, t, title='Exponential growth', xlabel='Time', ylabel='Population')
t = np.linspace(0,50,5)

这是我的输出：

​

​

为什么人口似乎从0开始？我怎样才能把这个从10开始修复呢？

Answer 1

它从10开始，但是指数增长很快，所以你的地块有一个1e9垂直标度，这使得10看起来很小。这就是为什么指数增长的事物通常是以对数标度绘制的。尝试减少x限制或用yscale="log"绘图以获得更好的视图。