Please excuse me for my stupidness - because omega ~ sqrt(k/m) , there should be 0.5*0.8^2*x^2+0.5*1.2^2*y^2 and 0.25*0.8^2*(x-y)^2 + 0.25*1.2^2*(x+y)^2 ; in these cases, results are OK. Regards Igor.