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
videotutoriales de matlab, octave, numpy, python, POA, algoritmia, simulacion de procesos en ingenieria quimica, y algo de ingenieria quimica
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
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);