y'(x)=\frac{(y(x)-\text{Si}(x))^2+\sin (x)}{x}