home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #19 / NN_1992_19.iso / spool / sci / math / symbolic / 2316 < prev    next >
Encoding:
Internet Message Format  |  1992-09-03  |  1.7 KB
  1. Path: sparky!uunet!haven.umd.edu!darwin.sura.net!wupost!usc!elroy.jpl.nasa.gov!news.claremont.edu!ucivax!news.service.uci.edu!beckman.com!dn66!a_rubin
  2. Newsgroups: sci.math.symbolic
  3. Subject: Re: ode solver, different results from different packages
  4. Message-ID: <a_rubin.715544670@dn66>
  5. From: a_rubin@dsg4.dse.beckman.com (Arthur Rubin)
  6. Date: 3 Sep 92 18:24:30 GMT
  7. References: <1992Sep3.102037.20539@news.Hawaii.Edu>
  8. Nntp-Posting-Host: dn66.dse.beckman.com
  9. Lines: 39
  10.  
  11. In <1992Sep3.102037.20539@news.Hawaii.Edu> qhuang@wiliki.eng.hawaii.edu (Qiming Huang) writes:
  12.  
  13. >Hi, I am trying to learn these different symbolic math packages
  14. >such as mathmatica, mapleV, and macsyma. I tested one example
  15. >listed on the macsyma manual for ordinary differential equations
  16. >which writes as
  17.  
  18. >    x*dy/dx +a*x*y^2+2*y+b*x=0
  19.  
  20. >the solution given by macsyma is
  21. >                          
  22. >                          (2 a sqrt(b) x)/sqrt(-a)
  23. >    y=((2 %c a b x - 2 %c sqrt(-a) sqrt(b)) %e
  24.  
  25. >                                          (2 a sqrt(b) x)/sqrt(-a)
  26. >+ sqrt(-a) a sqrt(b) x -a)/(2 %c sqrt(-a) a sqrt(b) x %e
  27. >   
  28. >   2
  29. >+ a x)
  30.  
  31.  
  32. Of course, the correct answer is:
  33.  
  34. if a=0, y = -b x /3 + %c / x^2
  35.  
  36. if a!=0, y = -1/(a x) + l/a produces the equation
  37.  
  38. dl/dx + l^2 + a b = 0,
  39.  
  40. which Mathematica also cannot solve, but "clearly" has solutions
  41.  
  42. (a b < 0) sqrt(-a b) tanh ( sqrt(-a b) (x + %c) )   (including c=+/- infinity)
  43. (a b = 0) 1/(x + %c)
  44. (a b > 0) sqrt(a b) tan ( sqrt(a b) (x + %c) )      
  45. --
  46. Arthur L. Rubin: a_rubin@dsg4.dse.beckman.com (work) Beckman Instruments/Brea
  47. 216-5888@mcimail.com 70707.453@compuserve.com arthur@pnet01.cts.com (personal)
  48. My opinions are my own, and do not represent those of my employer.
  49. My interaction with our news system is unstable; if you want to be sure I see a post, mail it.
  50.