¿Cómo resolver ecuaciones diferenciales exactas? 2

En esta ocasión, vengo a mostrarles la solución del problema 7 de ecuaciones homogéneas del libro de Ecuaciones diferenciales con aplicaciones y notas históricas 2da edición, del autor George F. Simmons
El problema es el siguiente:
Determinar si la ecuación es exacta y resolverla
$$(sen(x)sen(y)-xe^{y})dy=(e^{y}+cos(x)cos(y))dx$$
Primero debemos igualar a cero para poder calcular las respectivas derivadas y tendremos:
$$-(e^{y}+cos(x)cos(y))dx+(sen(x)sen(y)-xe^{y})dy=0$$
Para que la ecuación sea exacta debemos tener en cuenta lo siguiente:
$$\frac{\partial^{2}f}{\partial y \partial x}=\frac{\partial^{2}f}{\partial x \partial y}$$
O como normalmente se muestra:
$$\frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}$$
Corresponde a la ecuación diferencial:
$$M(x,y)dx+N(x,y)dy=0$$
Luego podemos identificar las funciones $M=\frac{\partial f}{\partial x}$ y $N=\frac{\partial f}{\partial y}$
El problema es hallar la función $f$ que cumpla tal ecuación (comprobando primero si es una ecuación exacta para este caso).
Calculamos la derivada parcial de $M$ respecto a $y$, del problema que vamos a intentar solucionar:
$$\frac{\partial M}{\partial y}=\frac{\partial}{\partial y}(-(e^{y}+cos(x)cos(y)))=\frac{\partial}{\partial y}(-e^{y}-cos(x)cos(y))$$
$$=-e^{y}+cos(x)sen(y)$$
Ahora calculamos la derivada parcial de $N$ respecto a $x$
$$\frac{\partial N}{\partial x}=\frac{\partial}{\partial x}(sen(x)sen(y)-xe^{y})$$
$$=-e^{y}+cos(x)sen(y)$$
Luego podemos ver que ambos resultados son iguales
$$\frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}$$
Corresponde a una ecuación diferencial exacta, ahora pasamos a hallar nuestra función $f$ como sigue:
$$f=\int Mdx + g(y)$$
Reemplazamos $M$ e integramos:
$$f=\int (-e^{y}-cos(x)cos(y))dx+ g(y)$$
Aplicamos linealidad e integramos:
$$f=-\int e^{y}dx-\int cos(x)cos(y)dx+ g(y)$$
$$f=-xe^{y}- sen(x)cos(y)+ g(y)$$
Derivamos parcialmente la función $f$ respecto a $y$
$$\frac{\partial f}{\partial y}=\frac{\partial}{\partial y}(-xe^{y}- sen(x)cos(y)+ g(y))$$
$$\frac{\partial f}{\partial y}=-xe^{y}+sen(x)sen(y)+g'(y)$$
Recordamos que $\frac{\partial f}{\partial y}=N$ y reemplazamos para hallar la derivada de la función $g(y)$:
$$sen(x)sen(y)-xe^{y}=-xe^{y}+sen(x)sen(y)+g'(y)$$
$$g'(y)=0$$
Luego integrando respecto a $y$ tenemos:
$$\int g'(y)dy=\int 0 dy$$
$$g(y)=c$$
Luego la solución queda reemplazando el valor de $g(y)$ en $f$
$$-xe^{y}-sen(x)cos(y)+c=0$$
$$xe^{y}+sen(x)cos(y)=c$$