ϕr=jb1,r22b4a(1+2b(b2)+jb1,r23(2b4a)2)+O(1a5),