我想连续绘制几个雅可比多项式的(复数)根。这些多项式的参数和阶数随尺寸因子n的变化而变化。到目前为止,我已经有了计算根并在复平面上绘制它们的代码,但我希望每个绘图都覆盖上一个绘图,以便查看是否会形成根的模式。你知道吗
我写的代码现在生成了这些图,但它们都是独立的。你知道吗
import sympy
from sympy import *
import numpy
x = sympy.symbols('x')
def show_roots(n,a,b,c):
for k in range (1,n+1):
p = sympy.jacobi(a*k,-(a+b)*k,(b+c)*k,x)
coeff = sympy.poly(p).coeffs()
roots = numpy.roots(coeff)
X = [x.real for x in roots]
Y = [x.imag for x in roots]
plt.scatter(X,Y, color='red')
plt.show()
plt.pause(0.5)
目前没有回答
相关问题 更多 >
编程相关推荐