Hola. Combinando estas ideas que han propuesto sale directo:
Podemos ver que
\( u(x, y) = \sin(\text{Re}(z^3)) \cosh(\text{Im}(z^3)) \)
desarrollando el seno y el coseno hiperbólico a partir de la última expresión que tienes llegas a una expresión de la forma \( g(z)=f(z)+f(\bar z) \) para una función holomorfa \( f \)
Sacando la constante se puede llegar a la expresión
\( u(x,y) = \frac{1}{2}(f(z)+f(\overline z)) \) donde \( f(z)=\sen(z^3) \).
Pero \( f \) conmuta con la conjugación, así que
\( u(x,y) = \frac{1}{2}(f(z)+\overline{f(z)}) = \text{Re}(f(z)) \).
Lo que prueba que \( u(x,y) \) es armónica.
Para hallar su conjugado armónico, expandes el cubo en \( \sen((x+iy)^3) \), separas seno de la suma (de real + imaginaria) y tomas su parte imaginaria.
A partir de una función armónica \( u(x,y) \) se puede reconstruir una función holomorfa (módulo constantes) con \( u \) como parte real calculando
\( 2\, u \left( \dfrac{z}{2}, \, \dfrac{z}{2i}\right) - u(0,0). \)
Esto se puede justificar usando series de potencia.