“…In the following, we expand functions g(r, z) and f (r, z) as power series of a small expansion parameter ε. In order to get the first-order solution, substituting function g(r, z) truncated as g(r, z) = εg 1 (r, z) and f (r, z) as f (r, z) = 1 + ε 2 f 2 (r, z) [32] in Eqs. (7) and (8) yields…”