''import sympy as sp
# 斐波那契数列
# 定义变量k、
k = sp.var('k', positive=True, integer=True)# 新建符号变量时可以指定其定义域
a = sp.Matrix([[0,1], [1,1]])
val = a.eigenvals() # 求特征值
vec = a.eigenvects() # 求特征向量
P,D = a.diagonalize() # 把a相似对角化
ak = P @(D**k)@(P.inv())
F = ak@sp.Matrix([1,1])
s = sp.simplify(F[0]) # 表达式
print(s)
sm = []
for i in range(20):
sm.append(s.subs(k, i).n()) #符号表达式S中的符号变量OLD替换为新的值NEW。 # .n()得到计算结果，否则得到计算表达式
print(sm)
