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

function modelo5=modelo5(C) % C=[V CA CB]
    F1=2.5; %m3/s
    A=2; %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",[VB,CA0,CB0],t) (octave)
res=ode45("modelo5",t,[V0,CA0,CB0]) (matlab)