domingo, 31 de enero de 2016

Resolver ecuaciones diferenciales con maple.

Tenemos la siguiente ecuación de variables separables a resolver.




solucion

dsolve({(diff(y(x), x))/y(x)^2 = -e^cos(x)*sin(x)});
print(`output redirected...`); # input placeholder
                               ln(e)       
                  |y(x) = -------------------
                                      cos(x) 
                         _C1 ln(e) - e
 

miércoles, 1 de julio de 2015

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

Método de McCabe Thiele para destilación fraccionada.

**** DESCARGAS ****