$$x \left(x - 1\right) \frac{d}{d x} y{\left(x \right)} + y{\left(x \right)} = x^{2} \left(2 x - 1\right)$$

Разделим обе части ур-ния на множитель при производной y':
$$x \left(x - 1\right)$$
Получим уравнение:
$$\frac{x \left(x - 1\right) \frac{d}{d x} y{\left(x \right)} + y{\left(x \right)}}{x \left(x - 1\right)} = \frac{x \left(2 x - 1\right)}{x - 1}$$
Это дифф. уравнение имеет вид:

y' + P(x)y = Q(x)

где
$$P{\left(x \right)} = \frac{1}{x \left(x - 1\right)}$$
и
$$Q{\left(x \right)} = \frac{x \left(2 x - 1\right)}{x - 1}$$
и называется линейным однородным
дифф. уравнением 1го порядка:
Решим сначала надо соответствующее линейное однородное ур-ние

y' + P(x)y = 0

с разделяющимися переменными
Данное ур-ние решается следущими шагами:

Из y' + P(x)y = 0 получаем

$$\frac{dy}{y} = - P{\left(x \right)} dx$$, при y не равным 0
$$\int \frac{1}{y}\, dy = - \int P{\left(x \right)}\, dx$$
$$\log{\left(\left|{y}\right| \right)} = - \int P{\left(x \right)}\, dx$$
Или,
$$\left|{y}\right| = e^{- \int P{\left(x \right)}\, dx}$$
Поэтому,
$$y_{1} = e^{- \int P{\left(x \right)}\, dx}$$
$$y_{2} = - e^{- \int P{\left(x \right)}\, dx}$$
Из выражения видно, что надо найти интеграл:
$$\int P{\left(x \right)}\, dx$$
Т.к.
$$P{\left(x \right)} = \frac{1}{x \left(x - 1\right)}$$, то
$$\int P{\left(x \right)}\, dx$$ =
= $$\int \frac{1}{x \left(x - 1\right)}\, dx = \left(- \log{\left(x \right)} + \log{\left(x - 1 \right)}\right) + Const$$
Подробное решение интеграла
Зн., решение однородного линейного ур-ния:
$$y_{1} = \frac{x e^{C_{1}}}{x - 1}$$
$$y_{2} = - \frac{x e^{C_{2}}}{x - 1}$$
что соотв. решению
с любой константой C, не равной нулю:
$$y = \frac{C x}{x - 1}$$
Мы нашли решение соотв. однородного ур-ния
Теперь надо решить наше неоднородное уравнение

y' + P(x)y = Q(x)

Используем метод вариации произвольной постоянной
Теперь, считаем, что C - это функция от x

$$y = \frac{x C{\left(x \right)}}{x - 1}$$
И подставим в исходное уравнение.
Воспользовавшись правилами
- дифференцирования произведения;
- производной сложной функции,
находим, что
$$\frac{d}{d x} C{\left(x \right)} = Q{\left(x \right)} e^{\int P{\left(x \right)}\, dx}$$
Подставим Q(x) и P(x) в это уравнение.
Получим простейшее дифф. ур-ние для C(x):
$$\frac{d}{d x} C{\left(x \right)} = 2 x - 1$$
Зн., C(x) =
$$\int \left(2 x - 1\right)\, dx = \left(x^{2} - x\right) + Const$$
Подробное решение интеграла
подставим C(x) в
$$y = \frac{x C{\left(x \right)}}{x - 1}$$
и получим окончательный ответ для y(x):
$$\frac{x \left(x^{2} - x + Const\right)}{x - 1}$$

$$y{\left(x \right)} = \frac{x \left(C_{1} + x^{2} - x\right)}{x - 1}$$