Ejemplo de un modelo de reactor lote. En Matlab

function lotecomplejo
[t,x]=ode45(@lote,[0 2],[16 0 0 0])
subplot(4,1,1); plot(t,x(:,1))
subplot(4,1,2); plot(t,x(:,2))
subplot(4,1,3); plot(t,x(:,3))
subplot(4,1,4); plot(t,x(:,4))


function dx=lote(t,x)

% definicion de parametros

k1=5;
k2=4;
k3=3;
k4=3;

dx=zeros(4,1);

% ecuaciones de diseño (balances de masa), x1 = Ca, x2 = Cb, x3 = Cc, x4 = Cd

dx(1) = - k1*x(1);
dx(2) = k1*x(1) - k2*x(2) + k3*x(3) - k4*x(2);
dx(3) = k2*x(2) - k3*x(3);
dx(4) = k4*x(2);