#!/usr/bin/python
import numpy as np
def main():
	a=np.array([[1,2],[3,4]])
	eigenvalues,eigenvectors=np.linalg.eig(a)
	print(eigenvalues)
	print (eigenvectors)
	print (eigenvectors[:,0]) #Ta idioduanusmata tis protis idiotimis
	print (eigenvectors[:,1]) #Tis deyterhs idiotimis
	#Verify
	d=eigenvectors[:,0]*eigenvalues[0]
	e=a@eigenvectors[:,0]
	print (d,e)
	print (d==e)
	print (np.allclose(d,e))
if __name__=="__main__":
	main()
