DSolve[Derivative[1][y][x] == (x^2/2 + (x^2*Cos[2*y[x]])/2 - Sin[2*y[x]]/2)/x,y[x],x]