MathGroup Archive 1993

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Eigensystem Misbehavior

  • To: mathgroup at
  • Subject: Re: Eigensystem Misbehavior
  • From: Mark E. Kotanchek <mek at>
  • Date: Fri, 8 Oct 93 14:12:40 -0400


I've figured out the problem with Mma's Eigensystem routine! It returns a list  
of the eigenvectors--which in Mma is equivalent to a matrix. However, this list  
of lists is the TRANSPOSE of V in the definition: 

   However, since the eigenvalues are defined as solving R.V =
   V.A, this should imply that A = V'.R.V where the apostrophe
   denotes the Hermitian transpose. (Since V'.V = V.V' = I where I
   is the identity matrix.)

Thus, to get back to the diagonal matrix of eigenvalues, I should use 

   Conjugate[V].R.Transpose[V] //Chop // MatrixForm

rather than,

      Conjugate[Transpose[V]].R.V //Chop // MatrixForm

IMHO, the documentation which could have helped me out on this one is a little  
on the sparse side. However, at the moment it appears that the eigensystem  
routine is working -- although, Jesus Rivero may disagree with me.

Thanks for your time and patience,

Mark Kotanchek
Guidance & Control Dept - 363 ASB
Applied Research Lab/Penn State
P.O. Box 30
State College, PA 16804

e-mail:	mek at (NeXTmail)
TEL:	(814)863-0682
FAX:	(814)863-7843

  • Prev by Date: Re: Several Functions on the same graph
  • Next by Date: Re: Producing expr == 0 from Eliminate
  • Previous by thread: Eigensystem Misbehavior
  • Next by thread: Integration failure