x1[1,3)	x2[1,3)	x3[1,3)	x4[1,3)	1/2*m*(omega^2+omega_0^2)*1/2*x^2
