The given differential equation can be written as
dydx=yx[xcos(yx)+ysin(yx)][ysin(yx)−xcos(yx)]
This being a homogeneous differential equation, we can solve the equation by putting
y=vx ⇒dydx=v+xdvdx
Now putting this, then the equation changes to
v+xdvdx=v[cos(v)+vsin(v)][vsin(v)−cos(v)]
xdvdx=v[cos(v)+vsin(v)][vsin(v)−cos(v)]−v
xdvdx=v[cos(v)+vsin(v)−vsin(v)+cos(v)]vsin(v)−cos(v)
xdvdx=2vcos(v)vsin(v)−cos(v)
⇒12∫[tan(v)−1v]dv−∫dxx=C
⇒12[ln|sec(v)|−ln|v|]−ln|x|=C
⇒−ln|vcos(v)|−2ln|x|=C
⇒ln|x2yxcos(yx)|=−C=ln(C1)
⇒xycos(yx)=C1
Here both C and C1 are constants.