read `common.map`: rhB:=oms*rhg+Arh*rhg^2: qB:=qg+q0+rhg*Aq+Cq*rhg^2: pB:=pg+p0+rhg*Ap+Cp*rhg^2: uB:=ug+u0+rhg*Au+Cp*rhg^2: for a7 to 4 do for b7 to 4 do jac[a7,b7]:=diff(xB[a7],xg[b7]): od:od: gg:=matrix(4,4): gt:=matrix(4,4): assume(K>0);assume(oms4>0); for a11 to 4 do for b11 to 4 do gt[a11,b11]:=subs(rh=rhg*oms,gBondi[a11,b11]): od:od: for a9 to 4 do for b9 to 4 do gg[a9,b9]:=sum(sum(gt[a10,b10]*jac[a10,a9]*jac[b10,b9],a10=1..4),b10=1..4): gg[a9,b9]:=convert(series('leadterm'(gg[a9,b9]),rhg),polynom): od:od: FA:=array(1..4,[0,rhg*fg/2*(sK1/2-J/2/sK1),I*rhg*fg/2*(sK1/2+J/2/sK1),0]): FA_Bondi:=array(1..4): for a11 to 4 do FA_Bondi[a11]:=sum(FA[a12]*jac[a11,a12],a12=1..4): FA_Bondi[a11]:=convert(simplify(series('leadterm'(FA_Bondi[a11]),rhg)),polynom): FA_Bondi[a11]:=simplify(subs(sK1^2=K+1,FA_Bondi[a11])): od: t1:=simplify((FA_Bondi[2]-I*FA_Bondi[3])/oms/rhg*sqrt(2)/fo): t1:=subs(sK1=sqrt(K+1),t1): t2:=simplify(subs(oms=oms0,rethAX,t1)): t3:=simplify(subs(rA,t2)): eidelta:=simplify(subs(qg=(zeta+zetab)/2,pg=(zetab-zeta)*I/2,t3));