如何用python绘制地图的稳定流形和不稳定流形？

Question

我需要为Henon映射$x{n+1}=1-ax{n^2+y{n$和$y{n+1}=bx{n$绘制稳定和不稳定流形，其中$a$和$b$是参数。为此，我知道我可以沿着雅可比矩阵的特征向量的方向选择一行初始条件，在固定点位置进行计算。问题是我不知道如何做到这一点，我对python或任何编程语言都很陌生。我可以把初始条件画成一条直线，但我不知道它们在特征向量方向上会是什么样子

def HenonMap(a,b,x,y):
return 1 - a*x**2 + y, b*x


from pylab import *
import numpy.linalg as LA
a = 1.4
b = 0.3
nIterates = 100
xlabel('x(n)') # set x-axis label
ylabel('y(n)') # set y-axis label

nICs = 10000
xmin = -6
xmax = 0.6
xinc = (xmax - xmin) / nICs
ymin =  -0.18
ymax =   0.18
yinc = (ymax - ymin) / nICs
x = []
y = []

for i in xrange(0,nICs):
    x.append(xmin + (xmax - xmin) * xinc * i)
    y.append(ymin + (ymax - ymin) * yinc * i)
for n in range(0,nIterates):
    for i in range(0,nICs):
        x[i], y[i] = HenonMap(a,b,x[i],y[i])
plot(x,y, '-r')

show（）