x1[1,5)	x2[1,5)	x3[1,5)	x4[1,5)	-rho_c_0*q*A_vec/m
