prac5

La función dCa se expresa en forma de vector. Cada componente del vector (z en la función)
está dado por la ecuación diferencial para cada reacción.

function z = dCa (t, C)
k1=0.2;
k2=0.1;
z = [-k1*C(1); k1*C(1)-k2*C(2)^2; k2*C(2)^2];
Esta función se ingresa como parámetro de la función ODE45. Observe el ejemplo:
C0=[1 0 0];
[t C]= ode45('dCa',[0 100],C0);
Ca = C(:,1);
Cb = C(:,2);Cc = C(:,3);
plot (t,Ca,t,Cb,t,Cc)

1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0

0

10

20

30

40

50

60

70

80

90

100

Esta función implementa el método de Runge Kutta (similar al ODE45)
function [t,z]=edoRK(f,a,b,z0,n)
dx=(b-a)/n;
t=linspace(a,b,n+1);
z(:,1)=z0;
for i=1:n
k1=dx*feval(f,t(i),z(:,i));
k2=dx*feval(f,t(i)+(dx/2),z(:,i)+(k1/2));k3=dx*feval(f,t(i)+(dx/2),z(:,i)+(k2/2));
k4=dx*feval(f,t(i)+dx,z(:,i)+k3);
z(:,i+1)=z(:,i)+(1/6)*(k1+(2*k2)+(2*k3)+k4);
end
z=z';
Observe el ejemplo:
C0=[1 0 0];
n = 100;
[t C]= edoRK('dCa', 0, 100,C0, n);
Ca = C(:,1);
Cb = C(:,2);
Cc = C(:,3);
plot (t,Ca,t,Cb,t,Cc)
Se va a preparar acetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100 gal de una solución decomposición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidad es de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor se calienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2 gpm (17.4lb/min). La solución se extrae a la misma velocidad volumétrica. A este nivel, los datos de velocidad son:Se va a prepararacetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100 gal de una solución de composición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidad es de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor se calienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2 gpm (17.4lb/min). La soluciónse extrae a la misma velocidad volumétrica. A este nivel, los datos de velocidad son:Se va a preparar acetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100 gal de una solución de composición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidad es de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor secalienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2 gpm (17.4lb/min). La solución se extrae a la misma velocidad volumétrica. A este nivel, los datos de velocidad son:Se va a preparar acetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100 gal de una solución de composición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidades de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor se calienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2 gpm (17.4lb/min). La solución se extrae a la misma velocidad volumétrica. A este nivel, los datos de velocidad son:Se va a preparar acetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100gal de una solución de composición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidad es de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor se calienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2 gpm (17.4lb/min). La solución se extrae a la misma velocidad volumétrica. A este nivel, los datos develocidad son:Se va a preparar acetato de etilo en un reactor de 110 gal. El reactor contiene originalmente 100 gal de una solución de composición igual a 20% en peso de etanol y 35% en peso de ácido acético. Su densidad es de 8.7 Ib/gal; suponga que este valor permanece constante para todas las composiciones. El reactor se calienta y se mantiene a 100°C. El etanol puro se adiciona a una velocidad de 2...