The Mathematica Journal
Feature Articles
New Products
New Publications
News Bulletins
Write Us
About the Journal
Staff and Contributors
Back Issues
Download this Issue

Speciation of Adenosine [Graphics:../Images/index_gr_100.gif]-triphosphate

Example 2.

ATP binds both [Graphics:../Images/index_gr_101.gif] and metal ions such as [Graphics:../Images/index_gr_102.gif] and presents a more complex example than the phosphoric acid example. Again, the system must be defined and makeInput[] is used to obtain the transpose of the stoichiometric number matrix [Graphics:../Images/index_gr_103.gif], the species vector, and the charge vector.


Again the parameters used by reactions[] must be defined. However, this time ATP is a biochemical reactant (i.e., it represents the sums of the molalities of the various ATP species) and as such is the entry for the parameter reactants. The vector pseudo specifies which of the species belong to the ATP pseudoisomer group. The values of [Graphics:../Images/index_gr_110.gif] for the species [Graphics:../Images/index_gr_111.gif] and [Graphics:../Images/index_gr_112.gif] have been constrained to be 7.0 and 3.0, respectively. In this case the transpose of the apparent stoichiometric number matrix [Graphics:../Images/index_gr_113.gif] as specified in ntprime should be set to {{1}}.


As before, ionicEquilib[] is used to obtain the molalities m of the species at equilibrium.


The mole fractions x within the ATP pseudoisomer group can be obtained by using the function moleFractions[]. All of the parameters have already been specified.


Converted by Mathematica    

[Article Index] [Prev Page][Next Page]