θ=-sgn((Subscript[z, 0]))a cos(((Sqrt[(Subscript[x, 0]^2)+(Subscript[y, 0]^2)])/(\[LeftBracketingBar](Subscript[v, 0])\[RightBracketingBar]))).