function spraydryer=solver(t,x)

global P Gdry A B C rhoL Kpw nAvp h MWwater MWair dHev CpL CpG rhoG muG dp mp0 vg

mp=x(1);
Gvap=x(2);
Tp=x(3);
Tg=x(4);
vs=x(5);
z=x(6);

Pw=MWair/MWwater*Gvap*P/Gdry;
P0=exp(A-B/(Tp+C))/760;                                                     % mmHg--->atm
Sp=3.14*((mp/rhoL*6/3.14)^(1/3))^2;

spraydryer(1)=Kpw*Sp*(Pw-P0);
spraydryer(2)=-Kpw*Sp*(Pw-P0)*nAvp;
spraydryer(3)=(h*Sp*(Tg-Tp)+Kpw*Sp*(Pw-P0)*dHev)/mp/CpL;
spraydryer(4)=-h*Sp*(Tg-Tp)*nAvp/(Gvap+Gdry)/CpG;
spraydryer(5)=((1-rhoG/rhoL)*9.81-3*vs*muG*3.14*dp/mp0);
spraydryer(6)=(vs+vg);

spraydryer=spraydryer';
end