x1[1,3)	x2[1,3)	x3[1,3)	x4[1,3)	x5[1,3)	x6[1,3)	x7[1,3)	x8[1,3)	mom*H/(kb*T)+(mom*alpha)/(epsilon*c^2*kb*T)*M
