# -*- coding: utf-8 -*-
from __future__ import division
import math
import jacobi
import random
import sys
#import psyco
#psyco.full()

if len(sys.argv)>1 : n=int(sys.argv[1])
else : n=10
M = [[0 for i in range(n)] for j in range(n)]
for i in range(n) :
	for j in range(i,n) :
		M[i][j]=random.random()
		M[j][i]=M[i][j]

E, V, sweeps = jacobi.jacobi(M) # wow min while-løkke virkede :)

print 'sweeps = ', sweeps
print 'E='
print E[len(E)-1]