dsolve((x^2+1)*diff(y(x),x)+x*sin(y(x))*cos(y(x))-x*(x^2+1)*cos(y(x))^2 = 0,y(x))