krebs Paul Cochrane Simulation of a simplification of the Krebs cycle of biochemical reactions. Adapted for xmds from "Mathematica computer programs for physical chemistry", William H. Cropper, Springer Verlag (1998) This is a cyclic reaction scheme with the following reactions: A + R -> B + S B -> C R + C -> S + D D -> A with this reaction included which produces R at a constant rate from a reactant P whose concentration is large enough to be nearly constant. P -> R The overall reaction is R -> S With equations: d[A]_dt = k4[D] - k1[A][R] d[B]_dt = k1[A][R] - k2[B] d[C]_dt = k2[B] - k3[C][R] d[D]_dt = k3[C][R] - k4[D] d[R]_dt = k5 - k1[A][R] - k3[C][R] t yes yes yes yes main 1 main double A B C D R RK4IP 250 100000 1000 main Aout Bout Cout Dout Rout