x1[1,3)	x2[1,3)	x3[1,3)	x4[1,3)	x5[1,3)	x6[1,3)	p_d/(4*pi*epsilon)*3*z/r^5*sqrt(x^2+y^2)
