// Reactor semibatch con reacción A --> B
k = 1 // min-1

function y = semibatch(t,x,v)
  V = x(1); Ca = x(2)
  dVdt = v
  dVCadt = -k*V*Ca + v*Cain
  dCadt = (dVCadt - Ca*dVdt)/V
  y = [dVdt dCadt]
endfunction

v = input('v entrada = ')
Cain = input('Cain = ')
Vo = 1 // L
Cao = 1 // concentración inicial, mol/L
xo = [Vo Cao]'
t=0:0.1:10
x = ode(xo,0,t,semibatch)
clf
subplot(2,1,1)
plot2d(t,x(1,:))
xlabel('tiempo (min)')
ylabel('volumen (L)')

subplot(2,1,2)
plot2d(t,x(2,:))
xlabel('tiempo (min)')
ylabel('Ca (mol/L)')
