% Carico termico totale (W)
Qtot = 13000;
% Portata di vapore immessa nel locale (kg/s)
mv = 6/3600;

% Condizioni nel locale (temperatura °C, umidità relativa)
tA = 26;
PhiA = 0.55;

% Temperatura aria nelle condizioni di immissione (°C)
tI = 17;

% ===================
% Calore latente di evaporazione a 0°C (J/kg)
r0 = 2501e3;
% Calori specifici aria secca e vapore d'acqua (J/kg K)
cpa = 1006;
cpv = 1875;

% Formula per il calcolo della tensione di vapore 
Avs = 6.41542;
Bvs = 31.614894;
Cvs = 1.327603e-1;
Dvs = 1.5593343e-5;
pvs = @(t) exp(Avs+2.301585*t/(Bvs+Cvs*t+Dvs*t^2));

% Variabile temporanea utile per il calcolo
a = Qtot/mv;

% Fattore termico
R = 1-mv*r0/Qtot;
disp(['R=',num2str(R)])

% Umidità specifica presente nel locale (kgv/kga)
xA = 0.622*PhiA*pvs(tA)/(101325-PhiA*pvs(tA));

% Umidità specifica presente nelle condizioni di immissione,
% ottenuta dividendo il bilancio energetico del locale per il
% bilancio di massa acqua del locale e risolvendo per xI (kgv/kga)
xI = (cpa*(tA-tI)+cpv*xA*tA+(r0-a)*xA)/(r0-a+cpv*tI);
disp(['xI=',num2str(xI)])

% portata di massa di aria trattata (kg/s) da bilancio di massa acqua
ma = mv/(xA-xI);
disp(['ma=',num2str(ma)])