Python：将分析解决方案绘制为三阶

Question

我目前正在阅读关于国粹的评论文章，我想尝试重现他们的一些情节并玩弄免费参数。我使用的等式是：

和

我试图获得这个情节：

然而，在玩完初始条件之后，我能得到的最合适的是：

我的方法是简单地将等式中的和扩展为w，直到n = 3（但我不确定这是否正确），我的代码是：

import numpy as np
import matplotlib.pyplot as plt

w0 = -0.667
omega0 = 0.6

a = np.linspace(0,1,1000)
omega = omega0*(a**(-3.0*w0))/(omega0*(a**(-3.0*w0)) + 1.0 - omega0)
w = w0
for n in range(1,4):
    w += ((((-1.0)**n)*w0*(1.0 - w0**2.0))/(1.0 - (n + 1.0)*w0 + 2.0*n*(n + 1.0)*(w0**2.0)))*((omega/(1.0 - omega))**n)

plt.plot(a,w)
plt.grid(True)
plt.axis([0,1,-0.8,-0.6])
plt.xlabel('a')
plt.ylabel('w')
plt.show()

编辑：我的问题是为什么我的代码产生的图形与文章中的图形略有不同（注意文章末尾的渐变略有反转）？我真的不知道我正在做的是解决第1，第2和第3顺序的分析解决方案（如论文标题中所述）

非常感谢任何帮助。