As $M=2x\sin(y) +1$, $N= 4x^2\cos(y) + 3x\cot(y)+10 \sin(y)\cos(y)$ we have
$$
M_y-N_x= 2x\cos(y)-[8x\cos(y)+3\cot(y)]=-6x \cos(y) -3\frac{\cos(y)}{\sin(y)}
\implies
\frac{N_x-M_y}{M}= 3\cot(y).
$$
Therefore we are looking for an integrating factor $\mu(y)$, satisfying
$$
\frac{\mu'}{\mu}= 3\cot(y)\implies \ln (\mu)=3\int \cot(y)\,dy=3\ln (\sin(y))
$$
(we take constant equal $0$). Then $\mu =\sin^3(y)$.
After multiplication we get
$$
\bigl[2x\sin^4(y) +\sin^3(y)\bigr]\,dx +
\bigl[4x^2\sin^3(y)\cos(y) + 3x\sin^2(y)\cos(y)+10 \sin^4(y)\cos(y)\bigr]\,dy=0\,.
$$
Looking for $H(x,y)$ satisfying
\begin{align}
&H_x = 2x\sin^4(y) +\sin^3(y)\,,\label{eq-1-1}\\
&H_y=4x^2\sin^3(y)\cos(y) + 3x\sin^2(y)\cos(y)+10 \sin^4(y)\cos(y)\,.
\label{eq-1-2}
\end{align}
(\ref{eq-1-1}) implies
$$H(x,y)= \int \bigl[2x\sin^4(y) +\sin^3(y)\bigr]\,dx=
x^2\sin^4(y)+ x\sin^3(y) +h(y).$$
Then
$$
H_y= 4x^2\sin^3(y)\cos(y)+ 3x\sin^2\cos(y) +h'(y)
$$
and comparing with (\ref{eq-1-2}) we see that
$$
h'(y)=10\sin^4(y)\cos(y)\implies h(y)=\int 10\sin^4(y)\cos(y)\,dy= 2\sin^5(y);
$$
we pick-up constant $0$.
Finally
$$
H(x,y)=
x^2\sin^4(y)+ x\sin^3(y) +2\sin^5(y) =C
$$
is a solution to the original problem.