Untitled

V=100; %[L]
portataV= 2; %[L/min]
Ca0= 1; %[mol/L]
k=1; %[1/min]
t0=0; %[min]
t1=3; %[min]
t=[0:0.1:3];
dCdt= @(t,C) [-C(1); %dCa/dt
    2*C(1)-C(2)];    %dCb/dt
Cb0=0;
C0=[Ca0, Cb0];
tspan= [0 3],
[t,C]= ode45(dCdt, tspan, C0);
Ca=C(:,1); %estraggo la prima colonna della matrice 
Cb=C(:,2); %estraggo la seconda colonna della matrice
%serve per mettere Ca e Cb sullo stesso grafico
figure;
hold on;
plot(t,Ca, '-r');
plot(t,Cb,'-b');
xlabel('tempo');
ylabel('concentrazione');
grid on;
Cb_finale= C(end,2);
Cb_media= trapz(t,Cb)/3; %divide l'intervallo in 3 perhce guardando la regola degli infinitesimi il trapezio si divide in rettangoli e qui soo 3 minuti quindi diviso 3 
fprintf('Cb media= %.4f mol/L\n', Cb_media)