RE: Simple question from a biologist
- To: mathgroup at smc.vnet.net
- Subject: [mg67180] RE: [mg67139] Simple question from a biologist
- From: "David Park" <djmp at earthlink.net>
- Date: Sun, 11 Jun 2006 02:18:37 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Although it hasn't gotten you into any trouble here, in general it is better
to use symbol names that start with small case letters because then you can
be sure they won't conflict with Mathematica defined symbols. So rewriting
and taking the first (and only) solution...
sol = Solve[{{{-x, 0, z}, {x, -y, 0}, {0, y, -z}} .
{a, b, c} == 0, a + b + c == 1}, {a, b, c}][[1]]
{a -> (y*z)/(x*y + x*z + y*z), b -> (x*z)/(x*y + x*z + y*z),
c -> (x*y)/(x*y + x*z + y*z)}
Then we could, for example, calculate all your rates in a list and check
that they are equal. We just write the expression for each rate and then
substitute the solution. You don't have to Solve again.
{a x, b y, c z} /. sol
Equal @@ %
{(x*y*z)/(x*y + x*z + y*z), (x*y*z)/(x*y + x*z + y*z),
(x*y*z)/(x*y + x*z + y*z)}
True
Or, perhaps more what you want, your could define the rate function as
follows...
rate[x_, y_, z_] = a x /. sol
(x*y*z)/(x*y + x*z + y*z)
Then you could, for example, plot the rate as a function of x and y for
fixed values of z.
Plot3D[rate[x, y, 1], {x, 0.00001, 2}, {y, 0.00001, 2}];
A couple of good books are 'The Beginner's Guide to Mathematica' by Jerry
Glynn & Theodore Gray and 'Mathematica Navigator: Graphics and Methods of
Applied Mathematics' by Heikki Ruskeepaa. But it is very worthwhile going
through Part I of The Mathematica Book. And using Help as much as possible.
David Park
djmp at earthlink.net
http://home.earthlink.net/~djmp/
From: tnad [mailto:terry_najdi at hotmail.com]
To: mathgroup at smc.vnet.net
I'm a bit new to this so please bear with me. I solved this eqation:
Sol = Solve[{{{-x, 0, z}, {x, -y, 0}, {0, y, -z}}.{A, B, C} == 0, A + B + C
== 1}, {A, B, C}]
and got the ouputs of A, B and C interms of x,y, and z each.
Now I want to express a term called "rate" where rate = Ax = By = Cz in
terms of x, y and Z only.
So I tried to do this:
Solve[rate == Ax , rate] /. Sol
but I cannot get the rate in terms of x,y and z.
Is there a better way to do this? Also if someone knows of a better
tuttorial (better than the built-in tuttorial) for mathematica, please let
me know.