# Exercici 12.14

kill(all)$

Designem
V  el volum de l'aula
x(t) , quantitat de CO2 a l'aula
C(t), concentració de CO2 a l'aula. C(t) = x(t)/V
C0 = concentracio de CO2 de l'aire ambient = 4*10^(-4)
q, cabal de bombeix de l'aire per a la renovació

Quantitat de CO2 que entra en un temps dt : q*C0*dt
Quantitat de CO2 que surt en un temps dt : q*C(t)*dt
Pèrdua de CO2 en dt: - dx

V = V:384;
C0 = C0 : 4*10^(-4);

EDO:  ('diff(x,t))  + q*(x/V - C0)= 0;

ode2(EDO, x, t);

expand(%);

define (x(t) , c*%e^(-(q*t)/384)+96/625);

Eq1 : x(0) = 0.0012*V;
solve ([Eq1] , [c]);

c : 192/625;

Eq2: x(10) = 0.0006*V;

solve (Eq2 , q), numer;

q = q: 53.2337;

define (C(t) , x(t)/V);

wxplot2d([C(t)], [t,0,10], [y, 0, 0.0015] , [ylabel, ""] , [legend, "C(t)"] , [style, [lines, 2]] ,
 [gnuplot_preamble, "set grid; set zeroaxis;"])$