I think levap pretty much covers it (on how to find a potential function), but the only thing I wanted to add is another way to determine if a vector field is conservative.
Another way to check if $F$ is conservative is to calculate curl$F$ = curl$\nabla f$ = $\bigg(\dfrac{\partial F_z}{\partial y} - \dfrac{\partial F_y}{\partial z}\bigg)i + \bigg(\dfrac{\partial F_x}{\partial z} - \dfrac{\partial F_z}{\partial x}\bigg)j + \bigg(\dfrac{\partial F_y}{\partial x} - \dfrac{\partial F_x}{\partial y}\bigg)k = 0$. If the curl is not $0$, then $F$ cannot be conservative.
But if the curl does end up being $0$, you also have to check that $F$ is defined over all of $\mathbb{R}^3$ (or, for all $(x, y, z) \in \mathbb{R}^3$, $F(x, y, z)$ exists), and that all its component functions have continuous partials. In your case, since all the component functions are products and sums of continuous functions that are defined on the entire domain of $\mathbb{R}^3$, you don't have to worry about this condition, but I would just state it on the side.
Indeed,
$\bigg(\dfrac{\partial F_z}{\partial y} - \dfrac{\partial F_y}{\partial z}\bigg) = cos(yz) - yzsin(yz) - cos(yz) + yzsin(yz) = 0;$
$\bigg(\dfrac{\partial F_x}{\partial z} - \dfrac{\partial F_z}{\partial x}\bigg) = 0 - 0 = 0;$ and
$\bigg(\dfrac{\partial F_y}{\partial x} - \dfrac{\partial F_x}{\partial y}\bigg) = 1 - 1 = 0$.
Then $F$ is indeed conservative.
I would also refer to this previous post: Question on conservative fields, which gives detailed information on equivalent definitions of conservative. The previous solutions use property 3, there exists $f$ st $\nabla f = F$.