function [nu,k]=nuBO(T,rhos,rhoc,xi)

rhosr=rhos/rhoc;
nu=exp(-0.42704*(xi*rhosr-1)+...
    2.089*(xi*rhosr-1).^2-0.42367*(xi*rhosr-1).^3)-1;%...

k=xi*(1+nu.^(-1)).*(-0.42704+4.1780*(xi*rhosr-1)-1.27101*(xi*rhosr-1).^2)/rhoc;