def main():
"""
Generates a matrix with given eigenvalues
and launches the animation.
"""
print 'Running the power method...'
n = input('Give the dimension : ')
# m = input('Convergence factor ? ')
m = 1.2
L = np.random.uniform(0,2*np.pi,n)
j = complex(0,1)
K = np.exp(j*L)
for i in range(1,len(K)): K[i] = K[i]/m
A = generate_matrix(K)
p = Power_Method(n,eps=1.0e-2,A=A);
animate(p)
def main():
"""
Prompts the user for a dimension
and the number of iterations.
"""
print 'Running the power method...'
n = input('Give the dimension : ')
m = input('Convergence factor ? ')
L = np.random.uniform(0,2*np.pi,n)
j = complex(0,1)
K = np.exp(j*L)
for i in range(1,len(K)): K[i] = K[i]/m
A = generate_matrix(K)
p = Power_Method(n,A=A);
for i in range(1000):
p.next()
print p
if p.accurate(): break
