function u=PiT(Tb,rhob,omega,mu0)

R=8.314;
u=R*Tb*(4.44-0.78*omega+log(Tb))*rhob/mu0;