a p r 4 beta p r alpha r r0 2 ur0 m r r1 ur1 y x1 t1 y x1 t1 k 1 d z uz q x2 t2 x2 x3 t2 x3 phi 0 r 0.25 0 r 2 1 1.3061225 r 1.125 r 1.125 1 r r0 0 1 r r1 0 1 z 0 0 f h t4