x1[1,5)	x2[1,5)	x3[1,5)	sqrt(gamma*pr/rho)
