add sample code to find Y_e(Y_lep,rho,T)