x1[1,5)	x2[1,5)	x3[1,5)	x4[1,5)	x5[1,5)	n_rho*p_d^2*Ef/(3*kb*T)
