home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #26 / NN_1992_26.iso / spool / sci / math / symbolic / 2921 < prev    next >
Encoding:
Internet Message Format  |  1992-11-08  |  1.2 KB

  1. Path: sparky!uunet!charon.amdahl.com!pacbell.com!ames!sun-barr!rutgers!psuvax1!psuvm!sml108
  2. From: SML108@psuvm.psu.edu
  3. Newsgroups: sci.math.symbolic
  4. Subject: Mathematica Question
  5. Message-ID: <92311.165716SML108@psuvm.psu.edu>
  6. Date: 6 Nov 92 21:57:16 GMT
  7. Organization: Penn State University
  8. Lines: 20
  9.  
  10. Hi, I am posting a question for a friend, please send any replies by
  11. email so I don't clutter up the net...  Here goes...
  12.  
  13. We're trying to solve the following kinetics rate law:
  14. dx/dt = k*(ao - x)*(bo - x)*(x + xo):(1/2)
  15. where k,ao,bo, and xo are constants and x is a function of time t.
  16.  
  17. We have integrated the inverse of the right hand side of the equation with
  18. respect to x.  The resulting expression, which I shall call y(x), is then
  19. equal to time t.  Ultimately, we want to plot x as a function of time t, but
  20. we have been unable to manage this so far.  The Plot[Evaluate] function does
  21. not seem to work in this case.  There seems to be a imaginary part to the
  22. expression y(x).  Any suggestions on how to extract the real part would be
  23. helpful.
  24.  
  25. Once x is found as a function of t, we plan to differentiate x twice with
  26. respect to t in order to find where d2x/dt2 = 0, so any solution for x
  27. must be compatible with this next step.
  28.  
  29. Thanks in advance...
  30.