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