function [Tm,mu_l,rho_l] = n_properties(n,nw,Ts,Mw)
n0=0.34979;
Y0=-6288460;
Yinf=418.07;
b=8.929364;
g=0.0690406;
Tm=Yinf-(Yinf-Y0)*exp(-b*(n-n0)^g);
Tl=(Ts+2*Tm)/3;

n0=-2.293981;
Yinf0=57.8516;
b=2.476409;
g=0.0112117;
dY0= -602.688 + 77866.8/Tl + 198.006*log(Tl) - 4.18077*10^(-5)*Tl^2 - 2.49477*10^6/Tl^2;
dYinf= 0.0290196 - 241.023/Tl + 0.0440959*log(Tl) - 1.84891*10^(-7)*Tl^2 + 56561.7/Tl^2;
mu_l=exp(Yinf0+dYinf*(nw-n0)-dY0*exp(-b*(nw-n0)^g))*10^-3;

n0=-1.388524;
Yinf0=0;
b=5.519846;
g=0.0570632;
dY0= 8592.3 - 85.7292/Tl + 0.280284*log(Tl) - 4.48451*10^(-4)*Tl^2;
dYinf= 12.7924 + 0.0150627/Tl - 1.30794*10^(-5)*log(Tl) + 1.59611*10^(-8)*Tl^2;
Y=Yinf0 + dYinf*(n+n0) - dY0*exp(-b*(n+n0)^g);
rho_l=1000*Mw/Y;
end