t=0:1:30; %s
V0=1.5; %m
CA0=2.5; %mol/m3
CB0=0; %mol/m3

function modelo5=modelo5(C) %C=[V CA CB]
  F1=2.5; %m3/s
  A=1; %m2
  F2=0.5*(C(1)/A); %m3/s
  CA1=1.5; %mol/m3
  CB1=1; %mol/m3

  modelo5(1)=F1-F2;
  modelo5(2)=(F1*CA1-F2*C(2)-C(2)*(F1-F2))/C(1);
  modelo5(3)=(F1*CB1-F2*C(3)-C(3)*(F1-F2))/C(1);
endfunction

res=lsode("modelo5",[V0,CA0,CB0],t)