Page 256 - Thermodynamics of Biochemical Reactions
P. 256
256 Mathematica Solutions to Problems
This can be compared with the values calculated directly from In p:
p/.h->hh
6
I2.02478 10 , 208860., 22904.8, 2929.52, 495.647, 114.278, 32.5074, 10.5321, 3.97116, 1
1.29532, 1.09335, 1.02951, 1.00933, 1.00295, 1.00093, 1.0003}
The agreement is quite good as expected.
We can also plot the value of P calculated by the successive differentiation and integration to see whether it is the same as the
plot in part b.
Plot [Evaluate[Exp[Integrate[(nH/h) ,h]+23.716611, ~h,10A-2,10A-10},AxesLabel->I"h"r"g")l;
P
2X1O6 -
1. 5x106 -
1x106 -
500000 -
' h
0.002 0.004 0.006 0.008 0.01
(c) If we introduce the relation between h and pH, the expressions for p and nH will be converted to functions of pH.
10 -pH
B
10 6
1.99526 10 2.95121 10
1+ +
2 PH PH
10 10