home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / sci / math / symbolic / 2268 < prev    next >
Encoding:
Text File  |  1992-08-26  |  2.2 KB  |  65 lines
  1. Newsgroups: sci.math.symbolic
  2. Path: sparky!uunet!europa.asd.contel.com!darwin.sura.net!Sirius.dfn.de!chx400!bernina!bernina!parkes
  3. From: parkes@itp.ethz.ch (Andrew Parkes)
  4. Subject: Re: Mma Solve[] oddity (a much shorter version)
  5. In-Reply-To: petry@runners.math.washington.edu's message of Sun, 23 Aug 1992 21:53:53 GMT
  6. Message-ID: <PARKES.92Aug26173241@ptah.itp.ethz.ch>
  7. Sender: news@bernina.ethz.ch (USENET News System)
  8. Organization: Institute for Theoretical Physics, ETH, Zuerich, Switzerland
  9. References: <28914@option.GBA.NYU.EDU> <1992Aug23.215353.3637@u.washington.edu>
  10. Date: Wed, 26 Aug 1992 16:32:41 GMT
  11. Lines: 52
  12.  
  13.  
  14. David Petry writes:
  15. > In article <28914@option.GBA.NYU.EDU> jcao@option.GBA.NYU.EDU (Jingbin Cao) writes:
  16. > >The following log documents some weired Solve[] oddities. They are:
  17. > >
  18. > >    1. Solve[] cannot find roots in some cases, at least not within
  19. > >       reasonable period of time. But if you divide the equations
  20. > >       by something, Solve[] finds some roots!
  21. > 
  22. > [long transcript deleted]
  23. > 
  24. > Here's a simple example of what Jingbin Cao is getting at (I think):
  25. > 
  26. > eqs  =  {  a*x*y + b*x + c*y ==0,  d*x*y + e*x + f*y + g == 0}
  27. > 
  28. > Solve[ eqs, {x,y}]
  29. > 
  30. > Mathematica 2.  seems to be unable to solve the above equations, but if we 
  31. > divide out the "a" and "d", so that
  32. > 
  33. > eqs = { x*y + b*y + c*x == 0,  x*y + e*x + f*y + g == 0}
  34. > 
  35. > then Mathematica gets the answer right away.
  36. > 
  37. > A frustrating deficiency!
  38. > 
  39. > 
  40. > David Petry
  41. > 
  42. Here's an even simpler version of the above 
  43.  
  44. eq  = (a x^2 + b x + c)/(d1 + d2 x) 
  45. answer = Solve[ eq == 0 , x]
  46.  
  47. Mma 1.2 gets the answer immediately of course.
  48.  
  49. But I find that our  Mma 2.0  gets  gets nowhere. 
  50. I killed the job after 40mins when it got to 30MB.
  51.  
  52. If I give d1 a numerical value then the answer is again immediate.  I
  53. haven't looked whether giving some options to Solve[] makes things
  54. work correctly (as I don't think changing options should be necessary
  55. to solve a quadratic equation!).
  56.  
  57. In case this is version dependent, our system is:-
  58. $Version =  SPARC 2.0 (July 8, 1991)
  59. $VersionNumber = 2
  60.  
  61. Andrew Parkes,              | E-mail: parkes@itp.ethz.ch 
  62. Theoretische Physik,        | 
  63. ETH-Honggerberg,            |                           
  64. CH-8093 Zurich, Switzerland.|                           
  65.