Input Data

P := RandomGoodPoly(30):
N := 150:
a := [0.4, 0.1, 0.2, 0.3]:

EllipseProperties(P, alpha = a, scaling = "eval"):
AnimateConvergence(P, N, alpha = a, scaling = "eval");

Output Data

Under eigenvalue rescaling, the limiting ellipse will have normal form equation:
0.77 x^2 + -0.93 xy + 0.30 y^2 = 1
which has:
semi-major axis length = 8.1588
semi-minor axis length = 0.9733
eccentricity = 0.99286
with the semi-major axis tilted at an angle of 121.55 degrees counter-clockwise from the x-axis.