Bernoulli's Equation follows the form dy/dx + yP(x) = (yn)Q(x).
Here, x(dy/dx) + y = 1/y2 can go to dy/dx + y(1/x) = (y-2)(1/x), which
obeys the form above.
Make the transformation v = y-(-2)+1 or y3 and y = v1/3.
Then dy/dx = (1/3)v-2/3(dv/dx).
Rewrite dy/dx + y(1/x) = (y-2)(1/x) as (1/3)v-2/3(dv/dx) + v1/3(1/x) = (v1/3)-2(1/x).
Multiply through this last by 3v2/3 to gain dv/dx + 3v(1/x) = 3/x.
Note that dv/dx + 3v(1/x) = 3(1/x) is a linear differential equation of the first order with
integrating factor e∫(3/x)dx.
Then build the general solution or "primitive" as ve3ln|x| equals ∫(3/x) • e3ln|x|dx.
Rewrite this last as y3(x3) = ∫[(3/x) • x3]dx or x3y3 = x3 + C or y3 = 1 + Cx-3.
Finally, express y3 = 1 + Cx-3 as y = (1 + Cx-3)1/3.