dsolve(diff(y(x),x) = ((y(x)-Si(x))^2+sin(x))/x,y(x))