#PlotScaledRoots returns a plot of the roots of the (scaled) m-pt Strodt poly
#Note: If the polynomials were not squarefree, this program would loop forever
#Note: There is no particular reason for the color and size of the points
from sage.rings.polynomial.complex_roots import complex_roots
def PlotScaledRootsP(m,n):
L=len(complex_roots(P(m,n,n*x),skip_squarefree=True))
l=[]
for j in range(L):
l.append((real(complex_roots(P(m,n,n*x),skip_squarefree=True)[j][0]).center(), imaginary(complex_roots(P(m,n,n*x),skip_squarefree=True)[j][0]).center()))
return point(l,rgbcolor=hue(1),size=30)