x1[1,3)	x2[1,3)	x3[1,3)	x4[3,5)	x5[1,2)	q*Ef/(m*(omega_0^2-omega^2))
