home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / sci / math / symbolic / 2297 < prev    next >
Encoding:
Text File  |  1992-09-01  |  2.9 KB  |  89 lines
  1. Newsgroups: sci.math.symbolic
  2. Path: sparky!uunet!haven.umd.edu!darwin.sura.net!convex!news.utdallas.edu!aicklen
  3. From: aicklen@utdallas.edu (Greg Aicklen)
  4. Subject: Re: Mma Solve[] oddity (a much shorter version)
  5. Message-ID: <1992Sep1.131440.15465@utdallas.edu>
  6. Sender: usenet@utdallas.edu
  7. Nntp-Posting-Host: leakey.utdallas.edu
  8. Organization: Univ. of Texas at Dallas
  9. References: <28914@option.GBA.NYU.EDU> <1992Aug23.215353.3637@u.washington.edu> <PARKES.92Aug26173241@ptah.itp.ethz.ch>
  10. Date: Tue, 1 Sep 1992 13:14:40 GMT
  11. Lines: 76
  12.  
  13. In article <PARKES.92Aug26173241@ptah.itp.ethz.ch> parkes@itp.ethz.ch (Andrew Parkes) writes:
  14. >
  15. >David Petry writes:
  16. >> In article <28914@option.GBA.NYU.EDU> jcao@option.GBA.NYU.EDU (Jingbin Cao) writes:
  17. >> >The following log documents some weired Solve[] oddities. They are:
  18. >> >
  19. >> >    1. Solve[] cannot find roots in some cases, at least not within
  20. >> >       reasonable period of time. But if you divide the equations
  21. >> >       by something, Solve[] finds some roots!
  22.  
  23.     ... deletions ...
  24.  
  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.     ... deletions ...
  31.  
  32. >> David Petry
  33. >> 
  34. >Here's an even simpler version of the above 
  35. >
  36. >eq  = (a x^2 + b x + c)/(d1 + d2 x) 
  37. >answer = Solve[ eq == 0 , x]
  38. >
  39. >Mma 1.2 gets the answer immediately of course.
  40. >
  41. >But I find that our  Mma 2.0  gets  gets nowhere. 
  42. >I killed the job after 40mins when it got to 30MB.
  43.  
  44.     ... deletions ...
  45.  
  46. >In case this is version dependent, our system is:-
  47. >$Version =  SPARC 2.0 (July 8, 1991)
  48. >$VersionNumber = 2
  49. >
  50. >Andrew Parkes,              | E-mail: parkes@itp.ethz.ch 
  51. >Theoretische Physik,        | 
  52. >ETH-Honggerberg,            |                           
  53. >CH-8093 Zurich, Switzerland.|                           
  54.  
  55. With Mma 2.1 on a SPARC, I get quick solutions to both problems given above.
  56. For Andrew Parkes "even simpler version" I get
  57.  
  58. /tmp_mnt/home/tulia/aicklen % math
  59. Mathematica 2.1 for SPARC
  60. Copyright 1988-92 Wolfram Research, Inc.
  61.  -- X11 windows graphics initialized --
  62.  
  63. In[1]:= eq  = (a x^2 + b x + c)/(d1 + d2 x)
  64.  
  65.                      2
  66.         c + b x + a x
  67. Out[1]= --------------
  68.           d1 + d2 x
  69.  
  70. In[2]:= answer = Solve[ eq == 0 , x]
  71.  
  72.                             2                               2
  73.                  b    Sqrt[b  - 4 a c]           b    Sqrt[b  - 4 a c]
  74.                -(-) + ----------------         -(-) - ----------------
  75.                  a           a                   a           a
  76. Out[2]= {{x -> -----------------------}, {x -> -----------------------}}
  77.                           2                               2
  78.  
  79.  
  80. The original problem posed by David Petry is also solved readily by Mma 2.1,
  81. but the answer is way too long to post.
  82.  
  83. This leads to the question ... how can I find out what the changes were from
  84. version 2.0 to 2.1?
  85.  
  86. Greg Aicklen
  87. University of Texas at Dallas
  88. Erik Jonsson School of Engineering and Computer Science
  89.