(* PT coefficients for Barbanis potential *)
nqx = 9
 
nqy = 1
 
omx = 1
 
Null
 
ncoef = 5
 
Null
en[1] = 19/2 + (3*omy)/2
en[2] = (3*(5 - 76*omy - 12*omy^2))/(-8*omy^2 + 32*omy^4)
en[3] = (3*(-55 + 760*omy - 2133*omy^2 - 16720*omy^3 + 48124*omy^4 + 47880*omy^5 - 84096*omy^6 - 23712*omy^7 - 1152*omy^8))/(32*omy^5*(-1 + omy^2)*(-1 + 4*omy^2)^3)
en[4] = (3*(11745 - 188100*omy + 574650*omy^2 + 6541928*omy^3 - 41292455*omy^4 - 85075540*omy^5 + 864661020*omy^6 + 220896432*omy^7 - 8013337392*omy^8 + 4181347936*omy^9 + 28776773696*omy^10 - 18298997888*omy^11 - 43589405184*omy^12 + 24213351936*omy^13 + 26388099072*omy^14 - 8338802688*omy^15 - 4705714176*omy^16 - 543522816*omy^17 - 14155776*omy^18))/(128*(-1 + omy)^2*omy^8*(1 + omy)^2*(-3 + 2*omy)*(-1 + 2*omy)^5*(1 + 2*omy)^5*(3 + 2*omy)*(-1 + 4*omy)*(1 + 4*omy))
en[5] = (3*(56212380 - 1028298240*omy + 3509407125*omy^2 + 49996569600*omy^3 - 396024922555*omy^4 - 637125397920*omy^5 + 13626317948539*omy^6 - 8796612400000*omy^7 - 257425688946245*omy^8 + 417181313629920*omy^9 + 2932871076227316*omy^10 - 6530037676319424*omy^11 - 19502319849589280*omy^12 + 53342111158145408*omy^13 + 69000151348526464*omy^14 - 229259037947236224*omy^15 - 130024016132747520*omy^16 + 546917119847751168*omy^17 + 123298615299318784*omy^18 - 747366001447481344*omy^19 - 45322435132129280*omy^20 + 581633770387488768*omy^21 - 1231545913573376*omy^22 - 247493130438672384*omy^23 - 334555651768320*omy^24 + 51445921037549568*omy^25 + 3246076628828160*omy^26 - 3675551982157824*omy^27 - 792424486010880*omy^28 - 53654241411072*omy^29 - 913217421312*omy^30))/(2048*(-2 + omy)*(-1 + omy)^3*omy^11*(1 + omy)^3*(2 + omy)*(-3 + 2*omy)^2*(-1 + 2*omy)^7*(1 + 2*omy)^7*(3 + 2*omy)^2*(-1 + 4*omy)^2*(1 + 4*omy)^2)