Question

2) Chemical-kinetics problems often lead to sets of coupled,firsorder ODEs. For example, consider the reaction network Assuming a first-order reaction-rate expression for each transformation, material balances for each species lead to the following set of ODEs dt dB dt dC dt Solve the set of ODE given above, with k,-5, k2-2, k,-1, A(0) concentrations versus time graph 1, B(0)-C-(0)-0. Plot and label the

I need to solve this using matlab.

Answer 1

plz rate:

Script 1:

function dx=ala_func(t,x,Para)
a=Para(1);
b=Para(2);
c=Para(3);
dx=[-a*x(1)+b*x(2);a*x(1)-b*x(2)-c*x(2);c*x(2)];
end

Script2:

Para=[5 2 1]';
x0=[1 0 0]';
dt=0.01;
tspan=0:dt:30;
[t,x]=ode45(@(t,x) ala_func(t,x,Para),tspan,x0);
f=x(:,2);
h=x(:,1);
y=x(:,3);
plot(t,h)

hold on
plot(t,f)

hold on
plot(t,y)
legend(' Conc of A vs t',' Conc of B vs t',' Conc of C vs t')

1.2 Conc of A vs Conc ofB vs t Conc of C vs t 0.8 0.6 0.4 0.2 아 0.2 10 15 20 25 30