 Solution

The middle component suggests a scalar potential with $-{\varphi }_y=2y{\cosh }^2(xz)$, so $\varphi =-y^2{\cosh }^2(xz)+h(x,z)$. The chain rule and $2\sinh t\cosh t=\sinh (2t)$ giveMatching the other two components requires $h_x=h_z=0$. ThereforeIts negative gradient is the given vector field in all three components; the additive constant has no effect on that gradient.