Hcit=0.245653
Hjme=[1 0.952811 1.453925 0.742619 0.275628]
Phi=[8 0 8 0 1]
epsilon=0.508847
Zcit=Hjme/Hcit-epsilon*Phi
Zcit=Zcit([2:end])
Zjme=Hjme/Hcit+epsilon*Phi
syms s
Zcits=vpa(poly2sym(Zcit,s))
Zjmes=vpa(poly2sym(Zjme,s))
c1=maple('convert',Zcits/Zjmes,'confrac',s,'simple')
pretty(vpa(c1,5))