home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #16 / NN_1992_16.iso / spool / sci / math / symbolic / 2042 < prev    next >
Encoding:
Text File  |  1992-07-21  |  1.4 KB  |  39 lines

  1. Newsgroups: sci.math.symbolic
  2. Path: sparky!uunet!stanford.edu!EE.Stanford.EDU!siegman
  3. From: siegman@EE.Stanford.EDU (Anthony E. Siegman)
  4. Subject: Re: Mma and linear algebra
  5. Message-ID: <1992Jul22.025710.17724@EE.Stanford.EDU>
  6. Organization: Stanford University
  7. References: <SMITHW.92Jul21184319@sofya.hamblin.math.byu.edu>
  8. Date: Wed, 22 Jul 92 02:57:10 GMT
  9. Lines: 28
  10.  
  11. From Dr. William V. Smith):
  12.  
  13. >Someone asked about Mma and linear algebra, in particular, finding the
  14. >eigenvalues and eigenvectors of the matrix
  15. >
  16. >{{0,1,0},{0,0,1},{1,1,0}}
  17. >
  18. >I tried "Eigensystem" on this on my home NeXT, and it did fail.
  19. >Mathematica has problems with complex numbers.  I also just
  20. >tried it on a system at work (a MIPS) and it gave the algorithm
  21. >failed message. 
  22.  
  23. Try {{0.,1.,0.},{0.,0.,1.},{1,.1.,0.}}.  On a Mac SE/30 with mma 2.0
  24. both the eigenvalues and eigenvectors come out almost instantaneously.
  25. Without the dots the eigenvalues come out in a few seconds, and I gave
  26. up on the eigenvectors after a minute or so.
  27.  
  28. One thing I've learned from several helpful msgs from David Withoff at
  29. wri is that mma behaves _very_ differently in many situations
  30. depending upon whether you write integers (including especially zero)
  31. with or without a decimal point.  Even something like
  32.  
  33.     m = 1; rootTwo = Sqrt[2.]; f[x_] := Exp[ - m rootTwo x^2]
  34.  
  35. can behave differently for "m = 1" or "m = 1." despite the fact that
  36. you'd think the product "m rootTwo' would lose the pure number
  37. identification. 
  38.  
  39.