def Eigen_Vectors2(A,lambdas=[]): if (not lambdas): lambdas=A.Eigen_Values() v0=Vector([ lambdas[0]-A[0][0], A[1][0] ]) if (abs(lambdas[0]-A[0][0])==0.0): v0=Vector([1.0,0.0]) v0=v0.Normalize() v1=v0.Transverse2() return Matrix([ [ v0[0],v1[0] ], [ v0[1],v1[1] ], ])