Beispiel - SAM - ETH Zürich

Fallsc A (0) > c B (0) ➢
2. Reaktion kontrolliert Langzeitdynamik
Numerische
Mathemtik
Numerisches Experiment, MATLAB,t 0 = 0,T = 1,k 1 = 10 4 ,k 2 = 10 3 ,k 3 = 10,k 4 = 1
MATLAB-CODE : Explizite Integration steifer chemischer Reaktionsgleichungen
fun = @(t,y) ([-k1*y(1)*y(2) + k2*y(3) - k3*y(1)*y(3) + k4*y(4);
-k1*y(1)*y(2) + k2*y(3);
k1*y(1)*y(2) - k2*y(3) - k3*y(1)*y(3) + k4*y(4);
k3*y(1)*y(3) - k4*y(4)]);
tspan = [0 1];
y0 = [1;1;10;0];
options = odeset(’reltol’,0.1,’abstol’,0.001,’stats’,’on’);
[t,y] = ode45(fun,[0 1],y0,options);
R. Hiptmair
rev 35327,
24. Juni
2011
3.5
p. 378