φ=sgn((Subscript[x, 0]))a cos(((Subscript[y, 0])/((Subscript[x, 0]^2)+(Subscript[y, 0]^2)))),