The differential equation is:
[xsin2(xy)−y]dx+xdy=0
Rearranging:
xsin2(xy)dx+(xdy−ydx)=0
Since d(xy)=x2xdy−ydx, we have xdy−ydx=x2d(xy)
xsin2(xy)dx+x2d(xy)=0
Dividing by x2:
xsin2(xy)dx+d(xy)=0
Let v=xy:
xsin2(v)dx+dv=0
Dividing by sin2(v):
xdx+csc2(v)dv=0
Integrating:
∫xdx+∫csc2(v)dv=C
ln∣x∣−cot(v)=C
ln∣x∣−cot(xy)=C
Applying y=4π when x=1:
ln(1)−cot(1π/4)=C
0−1=C
C=−1
The particular solution is:
lnx−cot(xy)=−1
lnx=cot(xy)−1
x=ecot(xy)−1