\[ y'(x)=\frac {(y(x)-\text {Si}(x))^2+\sin (x)}{x} \] ✓ Mathematica : cpu = 0.0807745 (sec), leaf count = 17
DSolve[Derivative[1][y][x] == (Sin[x] + (-SinIntegral[x] + y[x])^2)/x,y[x],x]
\[\left \{\left \{y(x)\to \text {Si}(x)+\frac {1}{-\log (x)+c_1}\right \}\right \}\] ✓ Maple : cpu = 0.069 (sec), leaf count = 15
dsolve(diff(y(x),x) = ((y(x)-Si(x))^2+sin(x))/x,y(x))
\[y \left (x \right ) = \operatorname {Si}\left (x \right )+\frac {1}{-\ln \left (x \right )+c_{1}}\]