T=real(3404.255((1-a)/107.99^2)^(1/4) (1+3/4 t)^(1/4) ) ; t=0.385 ; a=0.306