Write $f = u+iv$, with real-valued functions $u,\, v$. The Cauchy-Riemann equations are
$$\frac{\partial u}{\partial x} = \frac{\partial v}{\partial y}; \qquad \frac{\partial u}{\partial y} = - \frac{\partial v}{\partial x}.$$
So with the given $u(x,y) = x^3y - xy^3$, the first yields
$$\frac{\partial u}{\partial x} = 3x^2y - y^3 = \frac{\partial v}{\partial y},$$
and thus $v(x,y) = \frac{3}{2}x^2y^2 - \frac{1}{4}y^4 + \varphi(x)$ with an unknown function $\varphi$. The second CR equation then yields
$$-\frac{\partial u}{\partial y} = 3xy^2 - x^3 = 3xy^2 + \varphi'(x),$$
so we obtain $\varphi(x) = -\frac14 x^4 + c$ with an unknown (real) constant $c$.
Assembling,
$$v(x,y) = - \frac14x^4 + \frac32 x^2y^2 - \frac14y^4 + c = -\frac14(x^4 - 6x^2y^2 + y^4) + c,$$
so
$$f(x+iy) = x^3y - xy^3 -\frac{i}{4}(x^4 - 6x^2y^2 + y^4) + ic = \frac{1}{4i}(x+iy)^4 + ic.$$