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