Let Φ = f + ig. Then fx = 3x2 - 3y2 (so that f(x,y) = x3 - 3xy2 + h(y) for some function h) and fy = -6xy, i.e., -6xy + h'(y) = -6xy. We deduce h'(y) = 0, forcing h(y) to be constant. If we assume Φ0,0) = 0, then f(0,0) = 0, implying h(y) = 0. Hence f(x,y) = x3 - 3xy2.
Since Φ is analytic, the Cauchy - Riemann equations are satisfied, so
(1) gx = -fy = 6xy
(2) gy = fx = 3x2 - 3y2
Integrating equation (1) yields g(x,y) = 3x2y + k(y) for some function k. Then (2) gives 3x2 + k'(y) = 3x2 - 3y2, or, k'(y) = - 3y2. Therefore k(y) = -y3 + C for some constant C. The condition Φ0,0) = 0 implies g(0,0) = 0, whence C = 0. Consequently, g(x,y) = 3x2y - y3, and
Φ(z) = (x3 - 3xy2) + (3x2y - y3)i = x3 + 3x(iy)2 + 3x2(iy) + (iy)3 = (x + iy)3 = z3
as desired.
Note: All complex potentials associated to v are of the form z3 + C where C is an arbitrary constant. Prescribing a value of Φ at a point then uniquely determines Φ. This is what I've done here, assuming Φ(0,0) = 0.