Absorbed some phases in the definition of X and Y to get more symmetric reduced stability matrix. Fixed error of missing lower diagonal ph eneriges. Added more printout