According to Wolfram,
$$\frac{dy}{dx}-y^2=\frac{b}{x}+a$$
The equation above has the solution;
Where $a$ and $b$ are constants. $_1F_1(a;b;x)$ and $U(a,b,x)$ are the Kummer confluent hypergeometric function and confluent hypergeometric function of the second kind, respectively.
I don't know how Wolfram got this solution, can anyone here help me please?
Thank you guys

$$\frac{dy}{dx}=y^2+\frac{b}{x}+a$$ This is a Riccati ODE. The usual change of function is : $$y(x)=-\frac{u'(x)}{u(x)}\quad;\quad y'=-\frac{u''}{u}+\frac{(u')^2}{u^2}$$ $$y'=\left(-\frac{u'}{u}\right)^2+\frac{b}{x}+a=-\frac{u''}{u}+\frac{(u')^2}{u^2}$$ $$\frac{b}{x}+a=-\frac{u''}{u}$$ $$u''-\left(\frac{b}{x}+a\right)u=0$$ This is the general confluent hypergeometric differential equation :
https://mathworld.wolfram.com/GeneralConfluentHypergeometricDifferentialEquation.html
From Eq.$(1)$ in this page determine $A$, $f(x)$ and $h(x)$. This leads to :
Then computing $y(x)=-\frac{u'(x)}{u(x)}$ will lead to the solution.