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);
