# This example does a chemical kinetics simulation, for the reaction: # # A + 2B <==> C --> D # # The output displays the concentration of the species `C' # as a function of time. # # You may run this example by doing: # # ode < chem.ode | graph -T X -C # # or alternatively, to get a real-time plot, # # ode < chem.ode | graph -T X -C -x 0 10 -y 0 0.03 # # To improve the shape of the plotted curve, you may # wish to spline it, by doing e.g. # # ode < chem.ode | spline | graph -T X -C # # Alternatively, you could remove the `every 10' clause below. # The three rate constants are: # kf : A + 2B --> C # kb : C --> A + 2B # kd : C --> D a' = kb*c - kf*a*b^2 b' = kb*c - kf*a*b^2 c' = kf*a*b^2 - kb*c - kd*c d' = kd*c c = 0 d = 0 a = 0.1 b = 1 kf = 1 kb = 1 kd = 1 print t,c every 10 step 0,10
Generated by dwww version 1.15 on Wed Jun 26 01:22:12 CEST 2024.