“…Changing to cylindrical coordinates x = r cos θ, y = r sin θ, z = z, the system (2) writes as (16)ṙ = εg 1 (θ, r, z), θ = 1 + ε r cos θ(e + h|z| + f r cos θ) − sin θ(a + d|z|+ (b − g)r cos θ) − cr sin 2 θ , z = ir cos θ + ε(j + m|z| + kr cos θ + lr sin θ), where g 1 (θ, r, z) = cos θ a + d|z| + br cos θ + sin θ e + h|z|+ (c + f )r cos θ + gr sin 2 θ, and taking θ as the new independent variable system (2) becomes (17) r ′ = εg 1 (θ, r, z), z ′ = ir cos θ + εg 2 (θ, r, z), where the prime denotes derivative with respect to θ, and g 2 (θ, r, z) = j + m|z| − f ir cos 3 θ + lr sin θ − i cos 2 θ e+ h|z| + (g − b)r sin θ + cos θ kr + di|z| sin θ+i sin θ(a + cr sin θ) . The unperturbed system is(18)…”