Skip to article frontmatterSkip to article content

Интегрирование рациональных функций

Higher School of Economics

До сих пор мы использовали либо какие-то простые наблюдения, либо какие-то трюки, чтобы интегрировать форму. Разумеется, наукой это назвать нельзя. Здесь мы систематически разработаем технику интегрирования для очень важного класса форм.

Элементарные свойства полиномов

Напомним, что полиномом P(x)P(x) от переменной xx над R\mathbb{R} называется выражение вида

αnxn+αn1xn1++α1x+α0, \alpha_nx^n + \alpha_{n-1}x^{n-1} + \cdots + \alpha_1x + \alpha_0,

где все αiR\alpha_i \in \mathbb{R} и можно положить, что αn0\alpha_n \ne 0 и в таком случае говорят, что полином P(x)P(x) имеет степень nn и пишут deg(P(x))=n.\mathrm{deg}(P(x))= n.

Множество всех полиномов от переменной xx над R\mathbb{R} обозначают так: R[x].\mathbb{R}[x].

Proof

Пусть

A(x)=αnxn+αn1xn1++α1x+α0,B(x)=βkxk+βk1xk1++β1x+β0,\begin{align*} A(x) &= \alpha_nx^n + \alpha_{n-1}x^{n-1} + \cdots + \alpha_1x + \alpha_0, \\ B(x) &= \beta_kx^k + \beta_{k-1}x^{k-1} + \cdots + \beta_1x + \beta_0, \end{align*}

где n=deg(A(x))n = \mathrm{deg}(A(x)), k=deg(B(x))k = \mathrm{deg}(B(x)), так что αn,βk0\alpha_n, \beta_k \ne 0.

Применим индукцию по nn.

(1) Если n=0n =0 (т.е. A(x)=α0A(x) = \alpha_0) и 0=deg(A(x))<deg(B(x))0=\mathrm{deg}(A(x)) < \mathrm{deg}(B(x)), то положим Q(x):=0Q(x):=0 и R(x):=A(x)R(x) : = A(x), т.е. имеем

α0=0B(x)+α0. \alpha_0 = 0 \cdot B(x) + \alpha_0.

(2) Если deg(A(x))=deg(B(x))=0\mathrm{deg}(A(x)) = \mathrm{deg}(B(x)) = 0, т.е. A(x)=α0,B(x)=β0A(x) = \alpha_0, B(x) = \beta_0, то положим R(x):=0R(x) :=0 и R(x):=α0β0R(x): = \frac{\alpha_0}{\beta_0}, т.е.

α0=α0β0β0+0. \alpha_0 = \frac{\alpha_0}{\beta_0} \beta_0 +0.

(3) Пусть теперь n>0n>0. Если deg(A(x))<deg(B(x))\mathrm{deg}(A(x)) < \mathrm{deg}(B(x)), то положим Q(x):=0Q(x) : = 0, а R(x):=A(x).R(x):=A(x).

Итак, пусть теперь теорема доказана в случае n>0n>0 и deg(A(x))deg(B(x))\mathrm{deg}(A(x)) \ge \mathrm{deg}(B(x)). Тогда можем записать

A(x)=αnxn+αn1xn1++α1x+α0=αnβkxnk(βkxk)+αn1xn1++α1x+α0=αnβkxnk(B(x)βk1xk1β0)+αn1xn1++α1x+α0=αnβkxnkB(x)+A~(x),\begin{align*} A(x) &= \alpha_nx^n + \alpha_{n-1}x^{n-1} + \cdots + \alpha_1x + \alpha_0 \\ &= \frac{\alpha_n}{\beta_k} x^{n-k} (\beta_k x^k) + \alpha_{n-1}x^{n-1} + \cdots + \alpha_1x + \alpha_0 \\ &= \frac{\alpha_n}{\beta_k} x^{n-k} \Bigl( B(x) - \beta_{k-1}x^{k-1} - \cdots - \beta_0 \Bigr) + \alpha_{n-1}x^{n-1} + \cdots + \alpha_1x + \alpha_0 \\ &= \frac{\alpha_n}{\beta_k}x^{n-k} B(x) + \widetilde{A}(x), \end{align*}

где deg(A~(x))=n1.\mathrm{deg}(\widetilde{A}(x)) = n-1. Но тогда по предположению индукции мы можем найти такие Q~(x)\widetilde{Q}(x) и R(x){R}(x), что

A~(x)=Q~(x)B(x)+R(x), \widetilde{A}(x) = \widetilde{Q}(x)B(x) + {R}(x),

и deg((x))<deg(B(x))\mathrm{deg}((x))< \mathrm{deg}(B(x)). Тогда получаем

A(x)=αnβkxnkB(x)+A~(x)=αnβkxnkB(x)+Q~(x)B(x)+R(x)=Q(x)B(x)+R(x),\begin{align*} A(x) &= \frac{\alpha_n}{\beta_k}x^{n-k} B(x) + \widetilde{A}(x) \\ &= \frac{\alpha_n}{\beta_k}x^{n-k} B(x) + \widetilde{Q}(x)B(x) + {R}(x)\\ &= Q(x) B(x) + R(x), \end{align*}

где Q(x):=αnβkxnk+Q~(x)Q(x): = \frac{\alpha_n}{\beta_k}x^{n-k} + \widetilde{Q}(x), чем доказательство существования Q(x)Q(x) и R(x)R(x) закончено.

(4) Докажем теперь единственность. Предположим, что

A(x)=Q1(x)B(x)+R1(x)=Q2(x)B(x)+R2(x), A(x) = Q_1(x) B(x) + R_1(x) = Q_2(x) B(x) + R_2(x),

где deg(R1(x)),deg(R2(x))<deg(B(x))\mathrm{deg}(R_1(x)), \mathrm{deg}(R_2(x)) < \mathrm{deg}(B(x)). Тогда получаем

(Q1(x)Q2(x))B(x)=R2(x)R1(x), \Bigl(Q_1(x) - Q_2(x)\Bigr) B(x) = R_2(x) - R_1(x),

следовательно,

deg((Q1(x)Q2(x))B(x))=deg(R2(x)R1(x)), \mathrm{deg}\left((Q_1(x) - Q_2(x)\Bigr) B(x) \right) = \mathrm{deg}(R_2(x) -R_1(x)),

но

deg((Q1(x)Q2(x))B(x))=deg(Q1(x)Q2(x))+deg(B(x)), \mathrm{deg}\left((Q_1(x) - Q_2(x)\Bigr) B(x) \right) = \mathrm{deg}(Q_1(x) - Q_2(x)) + \mathrm{deg}(B(x)),

и так как deg(R2(x)R1(x))<deg(B(x))\mathrm{deg}(R_2(x) - R_1(x)) < \mathrm{deg}(B(x)), то последнее равенство возможно лишь в случае, когда Q1(x)=Q2(x)Q_1(x) = Q_2(x), т.е. Q1(x)=Q2(x)Q_1(x) = Q_2(x), и следовательно, R1(x)=R2(x)R_1(x) = R_2(x), что и требовалось доказать.

Мы считаем известными (или мы просто их принимаем) следующие два факта из алгебры:

Тогда мы получаем очень важное для нас следствие:

Proof

Согласно теореме Безу, нужно просто доказать, что любой полином P(x)R[x]P(x) \in \mathbb{R}[x] делится либо на линейный полином вида (xa)(x-a), либо на квадратный (x2+bx+c)(x^2 + bx + c). Но по теореме Безу, если P(x)P(x) имеет корень aa, то P(x)P(x) делится на (xa)(x-a), потому все числа a1,,apa_1,\ldots, a_p — это просто все корни уравнения P(x)=0P(x) = 0.

Далее, рассмотрим теперь P(x)P(x) как полином в множестве C[x]\mathbb{C}[x], это возможно, потому что RC\mathbb{R} \subset \mathbb{C}. Тогда по основной теореме алгебры существует комплексное число ζ=α+β1\zeta = \alpha + \beta \sqrt{-1}, что P(z)=0P(z) =0. Тогда по теореме Безу P(x)P(x) делится на (xαβ1)(x- \alpha - \beta \sqrt{-1}).

Пусть P(x)=αnxn++α1x+α0P(x) = \alpha_nx^n + \cdots + \alpha_1 x + \alpha_0, тогда если P(ζ)=0P(\zeta) =0, где ζC\zeta \in \mathbb{C}, то получаем

P(ζ)=αnζn++α1ζ+α0=αnζn++α1ζ+α0=αnxn++α1x+α0=P(ζ)=0.\begin{align*} P(\overline{\zeta}) &= \alpha_n {\overline{\zeta}}^n + \cdots + \alpha_1 \overline{\zeta} + \alpha_0 \\ &= \overline{\alpha_n \zeta^n} + \cdots + \overline{\alpha_1 \zeta} + \overline{\alpha_0} \\ &= \overline{\alpha_nx^n + \cdots + \alpha_1 x + \alpha_0} \\ &= \overline{P(\zeta)} = 0. \end{align*}

Это означает, что ζ\overline{\zeta} тоже является корнем уравнения P(x)=0P(x) = 0, тогда по теореме Безу P(x)P(x) делится на xζ=x(αβ1)x - \overline{\zeta} = x - (\alpha -\beta \sqrt{-1}).

Таким образом, полином P(x)P(x) делится на произведение (xζ)(xζ)(x - \zeta)(x - \overline{\zeta}).

Находим

(xζ)(xζ)=(xαβ1)(xα+β1)=x22αx+α2+β2, (x - \zeta)(x - \overline{\zeta}) = (x - \alpha - \beta \sqrt{-1})(x - \alpha + \beta \sqrt{-1}) = x^2 - 2 \alpha x + \alpha^2 + \beta^2,

т.е. полином P(x)P(x) делится на полином второй степени. Это завершает доказательство.

Рациональная функция и её разложение

Мы умеем интегрировать формы вида P(x)dxP(x) \mathrm{d}x, где P(x)P(x) — полином, сейчас мы хотим развить теорию интегрирования форм вида P(x)Q(x)dx\frac{P(x)}{Q(x)}\mathrm{d}x, где Q(x)Q(x) тоже полином и всегда подразумевается, что Q(x)0Q(x) \ne 0.

Выражения вида P(x)Q(x)\frac{P(x)}{Q(x)} называются дробями, и как и в случае числовых дробей, целесообразно рассматривать так называемые несократимые дроби. Формализацией этого желания является следующее определение.

Множество всех рациональных функций образует поле \footnote{Доказательство этого факта очевидно и при этом банальное, нужно лишь проверить выполнение аксиом i.1 — i.9 определения поля (см. Определение Definition 1).}, которое так и называется полем рациональных функций и обозначается так: R(x).\mathbb{R}(x).

Так как имеется процесс деления полиномов друг на друга, то введём следующее понятие:

Используя теперь теорему Theorem 4, мы можем разложить знаменатель на более простые множители. Оказывается, это влечёт разложение и самой дроби. Введём следующее понятие.

Рассмотрим дробь P(x)Q(x)\frac{P(x)}{Q(x)}, согласно теореме Theorem 4, мы можем записать знаменатель в виде

Q(x)=(xa1)k1(xap)kp(x2+b1x+c2)m1(x2+bqx+cq)mq, Q(x) = (x - a_1)^{k_1} \cdots (x-a_p)^{k_p} (x^2 + b_1x + c_2)^{m_1}\cdots (x^2 + b_qx + c_q)^{m_q},

где k1++kp+2m1++2mq=deg(P(x))k_1 + \cdots + k_p + 2m_1 + \cdots + 2m_q = \mathrm{deg}(P(x)) и все ki,mjZ0k_i, m_j \in \mathbb{Z}_{\ge 0}, и более того, согласно теореме Безу, все aia_i — это всё корни уравнения Q(x)=0.Q(x) =0.

(1) Пусть хотя бы один kik_i больше нуля, обозначим его просто через kk, тогда можно записать Q(x)=(xa)kQ~(x)Q(x) = (x-a)^k \widetilde{Q}(x), где aa — это соответствующее число из чисел aia_i. Тогда aa не является корнем уравнения Q~(x)=0.\widetilde{Q}(x) =0.

Допустим теперь, что существует такое число AA и такой полином P~(x)\widetilde{P}(x), что

P(x)Q(x)=A(xa)k+P~(x)(xa)k1Q~(x). \frac{P(x)}{Q(x)} = \frac{A}{(x-a)^k} + \frac{\widetilde{P}(x)}{(x-a)^{k-1}\widetilde{Q}(x)}.

Для доказательства этого равенства достаточно подобрать эти неизвестные A,P~(x)A, \widetilde{P}(x) так, чтобы выполнялось равенство

P(x)AQ~(x)=(xa)P~(x). P(x) - A \widetilde{Q}(x) = (x-a)\widetilde{P}(x).

Так как AA это число, то оно не должно зависеть от xx, поэтому положим в этом равенстве x=ax = a, и тогда мы получаем, что

P(a)AQ~(a)=0 P(a) - A \widetilde{Q}(a) = 0

откуда A=P(a)Q~(a)A = \frac{P(a)}{\widetilde{Q}(a)}. Это выражение корректно, так как aa был выбран так, чтобы aa — корень уравнения Q(x)=0Q(x) =0, но не корень уравнения Q~(x)=0.\widetilde{Q}(x) = 0.

Далее, полином P~(x)\widetilde{P}(x) можно теперь определить следующим образом:

P~(x):=P(x)AQ~(x)xa. \widetilde{P}(x): = \frac{P(x) - A \widetilde{Q}(x)}{x-a}.

(2) Пусть теперь Q(x)Q(x) содержит хотя бы один сомножитель вида (x2+bx+c)m(x^2 + bx + c)^m, тогда запишем Q(x)=(x2+bx+c)mQ^(x)Q(x)= (x^2 + bx + c)^m \widehat{Q}(x), где уже Q^(x)\widehat{Q}(x) не делится на x2+bx+cx^2 + bx + c. Тогда подберём числа B,CB,C и полином P^(x)\widehat{P}(x) так, чтобы

P(x)Q(x)=Bx+C(x2+bx+c)m+P^(x)(x2+bx+c)m1Q^(x). \frac{P(x)}{Q(x)} = \frac{Bx + C}{(x^2 + bx + c)^m} + \frac{\widehat{P}(x)}{(x^2 + bx + c)^{m-1} \widehat{Q}(x)}.

Это то же самое, что подобрать эти же неизвестные, чтобы выполнялось равенство

P(x)(Bx+C)Q^(x)=(x2+bx+c)P^(x). P(x) - (Bx +C)\widehat{Q}(x) = (x^2 + bx + c) \widehat{P}(x).

Поступим следующим образом. Разделим полиномы P(x)P(x), Q^(x)\widehat{Q}(x) на x2+bx+cx^2 + bx + c с остатком;

P(x)=F(x)(x2+bx+c)+αx+β,Q^(x)=H(x)(x2+bx+c)+γx+δ.\begin{align*} P(x) &= F(x)(x^2 + bx + c) + \alpha x + \beta, \\ \widehat{Q}(x) &= H(x)(x^2 + bx + c) + \gamma x + \delta. \end{align*}

Тогда, подставляя в предыдущее равенство, получаем

F(x)(x2+bx+c)+αx+β(Bx+C)(H(x)(x2+bx+c)+γx+δ)=(x2+bx+c)P^(x). F(x)(x^2 + bx + c) + \alpha x + \beta - (Bx +C) \bigl( H(x)(x^2 + bx + c) + \gamma x + \delta\bigr) = (x^2 + bx + c) \widehat{P}(x).

Потребуем теперь, чтобы полином

R(x)=αx+β(Bx+C)(γx+δ)=0 R(x) = \alpha x + \beta - (Bx + C) (\gamma x +\delta ) =0

делился на x2+bx+cx^2 + bx + c без остатка.

Итак, имеем

R(x)=αx+β(Bx+C)(γx+δ)=γBx2+(αδBγC)x+(βδC).\begin{align*} R(x) &= \alpha x + \beta - (Bx + C) (\gamma x +\delta ) \\ &= - \gamma B x^2 + (\alpha - \delta B - \gamma C)x + (\beta - \delta C). \end{align*}

Разделив теперь R(x)R(x) на x2+bx+cx^2 + bx +c на x2+bx+cx^2 + bx +c, мы получим в остатке следующее выражение:

((bγδ)BγC+α)x+cγBδC+β.\Bigl( (b \gamma - \delta)B - \gamma C + \alpha \Bigr)x + c \gamma B - \delta C + \beta.

Тогда мы получаем систему (относительно неизвестных B,CB,C) линейных уравнений

{(bγδ)BγC=αcγBδC=β. \begin{cases} (b \gamma - \delta)B - \gamma C =- \alpha \\ c \gamma B - \delta C =- \beta. \end{cases}

Определитель этой системы имеет вид

Δ=bγδγcγδ=δ2bγδ+cγ2 \Delta = \begin{vmatrix} b \gamma - \delta & \gamma \\ c\gamma & -\delta \end{vmatrix} = \delta^2 - b\gamma \delta + c \gamma^2

Пусть γ0\gamma \ne 0, тогда

Δ=γ2((δγ)2+b(δγ)+c), \Delta = \gamma^2 \left( \left( -\frac{\delta}{\gamma} \right)^2 + b \left(-\frac{\delta}{\gamma}\right) + c \right),

но это есть значение полинома x2+bx+cx^2 + bx +c в точке x=δγx = -\frac{\delta}{\gamma} и, следовательно Δ0\Delta \ne 0, ибо мы предположили, что x2+bx+cx^2 + bx +c не имеет корней. Таким образом, система имеет решение, и числа с необходимым требованием существуют.

Если же γ=0\gamma =0, то Δ=δ2\Delta = \delta^2, но так как Q^(x)=H(x)(x2+bx+c)+γx+δ\widehat{Q}(x) = H(x)(x^2 + bx + c) + \gamma x + \delta, то δ0\delta \ne 0 ибо Q^(x)\widehat{Q}(x) на x2+bx+cx^2 + bx +c не делится.

Итак, в любом случае, решение системы существует, а значит, можно подобрать так B,CB,C, чтобы полином P(x)(Bx+C)Q^(x)P(x) - (Bx + C)\widehat{Q}(x) делится на x2+bx+cx^2 + bx + c без остатка. В таком случае полином

P^(x):=P(x)(Bx+C)Q^(xx2+bx+c. \widehat{P}(x): = \frac{P(x) - (Bx + C)\widehat{Q}(x}{x^2 + bx + c}.

Таким образом, доказательство теоремы сводится к повторному применению случаев (1) и (2), которые обеспечивают возможность последовательного выделения простых дробей из данной правильной дроби, вплоть до её исчерпывания. Это доказывает теорему.

Proof

С помощью замены y=xαy =x-\alpha интегрирование формы ω=A(xα)k\omega = \frac{A}{(x-\alpha)^k} сводится к интегрированию формы Ayk\frac{A}{y^k}, что и даёт нужный результат, так как dy=d(xα)=dx.\mathrm{d}y = \mathrm{d}(x-\alpha) = \mathrm{d}x.

Чтобы научится интегрировать форму Ax+B(x2+ax+b)ndx\dfrac{Ax + B}{(x^2 + ax + b)^n} \mathrm{d}x, нам нужно сначала научится интегрировать форму dx(x2+α2)n.\dfrac{\mathrm{d}x}{(x^2 + \alpha^2)^n}.

(1) Так как (arctan(y))=1y2+1(\arctan(y))' = \frac{1}{y^2 + 1}, то

ω1=dxx2+α2= dxα2((x2α2)+1)=1α2αd(xα)(x2α2)+1=1αd(xα)(xα)2+1=1αarctan(xα)+C.\begin{align*} \int \omega_1 &= \int \frac{\mathrm{d}x}{x^2 + \alpha^2} = \ \frac{\mathrm{d}x}{\alpha^2\cdot\left (\left( \frac{x^2}{\alpha^2} \right) + 1\right)}\\ &= \frac{1}{\alpha^2} \int \frac{\alpha \cdot\mathrm{d}\left(\frac{x}{\alpha}\right)}{\left( \frac{x^2}{\alpha^2} \right) + 1} = \frac{1}{\alpha} \int \frac{ \mathrm{d}\left(\frac{x}{\alpha}\right)}{\left( \frac{x}{\alpha} \right)^2 + 1} = \frac{1}{\alpha} \cdot \arctan\left( \frac{x}{\alpha}\right) + C. \end{align*}

(2) Пусть теперь n1n \ge 1, будем интегрировать ωn\omega_n по частям, т.е. воспользуемся правилом

udv=uvvdu. \int u \mathrm{d}v = uv - \int v \mathrm{d}u.

Положили u=1(x2+α2)nu = \frac{1}{(x^2 + \alpha^2)^n}, v=xv = x, находим

du=(1(x2+α2)n)dx=2nx(x2+α2)n+1dx,dv=dx. \mathrm{d}u = \left( \frac{1}{(x^2 + \alpha^2)^n} \right)'\mathrm{d}x = - \frac{2nx}{(x^2 + \alpha^2)^{n+1}}\mathrm{d}x, \qquad \mathrm{d}v = \mathrm{d}x.

Тогда

ωn=dx(x2+α2)n=x(x2+α2)n+2nx2(x2+α2)n+1dx=x(x2+α2)n+2n(x2+α2)α2(x2+α2)n+1dx=x(x2+α2)n+2n((x2+α2)(x2+α2)n+1dxα2dx(x2+α2)n+1)=x(x2+α2)n+2n(dx(x2+α2)nα2dx(x2+α2)n+1)=x(x2+α2)n+2nωn2nα2ωn+1,\begin{align*} \int \omega_n &= \int \frac{\mathrm{d}x}{(x^2 + \alpha^2)^n} = \frac{x}{(x^2 + \alpha^2)^n} + 2n \cdot \int \frac{x^2}{(x^2 + \alpha^2)^{n+1}}\mathrm{d}x \\ &= \frac{x}{(x^2 + \alpha^2)^n} + 2n \cdot \int \frac{(x^2+\alpha^2) - \alpha^2}{(x^2 + \alpha^2)^{n+1}}\mathrm{d}x \\ &= \frac{x}{(x^2 + \alpha^2)^n} + 2n \cdot \left( \int \frac{(x^2+\alpha^2)}{(x^2 + \alpha^2)^{n+1}}\mathrm{d}x - \alpha^2 \int \frac{\mathrm{d}x}{(x^2 + \alpha^2)^{n+1}}\right) \\ &= \frac{x}{(x^2 + \alpha^2)^n} + 2n \cdot \left( \int \frac{\mathrm{d}x}{(x^2 + \alpha^2)^{n}} - \alpha^2 \int \frac{\mathrm{d}x}{(x^2 + \alpha^2)^{n+1}}\right) \\ &=\frac{x}{(x^2 + \alpha^2)^n} + 2n\cdot \int \omega_n - 2n\alpha^2 \int \omega_{n+1}, \end{align*}

т.е. мы получили рекуррентное соотношение

ωn=x(x2+α2)n+2nωn2nα2ωn+1, \int \omega_n = \frac{x}{(x^2 + \alpha^2)^n} + 2n\cdot \int \omega_n - 2n\alpha^2 \int \omega_{n+1},

из которого и следует требуемое.

Выделим в выражении x2+ax+bx^2 + ax + b полный квадрат

x2+ax+b=(x+a2)2+(ba24), x^2 + ax + b = \left(x+ \frac{a}{2} \right)^2 + \left(b - \frac{a^2}{4} \right),

так как по условию x2+ax+b=0x^2 + ax + b =0 не имеет корней, то a24b<0a^2 - 4b <0, тогда положим

c2:=ba24,c=+ba24 c^2: = b- \frac{a^2}{4}, \qquad c = + \sqrt{b- \frac{a^2}{4}}

тогда сделаем замену

y:=x+a2, y:= x+ \frac{a}{2},

находим

dy=(x+a2)dx=dx,x2+ax+b=(x+a2)2+(ba24)=y2+c2,Ax+B=Ay+(BAa2).\begin{align*} \mathrm{d}y &= \left( x+ \frac{a}{2}\right)' \mathrm{d}x = \mathrm{d}x,\\ x^2 + ax + b &= \left(x+ \frac{a}{2} \right)^2 + \left(b - \frac{a^2}{4} \right) = y^2 + c^2, \\ Ax + B &= & Ay + \left(B - \frac{Aa}{2} \right). \end{align*}

Рассмотрим два случая.

(1) n=1n = 1, тогда получаем

Ax+Bx2+ax+bdx=Ax+(BAa2)y2+c2dy=A22ydyy2+c2+(BAa2)dyy2+c2=A2ln(y2+c2)+1c(BAa2)arctan(yc)+C,\begin{align*} \int \frac{Ax + B}{x^2 + ax + b}\mathrm{d}x = \int \frac{Ax + \left( B - \frac{Aa}{2} \right)}{y^2 + c^2}\mathrm{d}y &= \frac{A}{2} \int \frac{2y\mathrm{d}y}{y^2 + c^2} + \left( B - \frac{Aa}{2} \right) \int \frac{\mathrm{d}y}{y^2 + c^2} \\ &= \frac{A}{2} \ln(y^2 + c^2) + \frac{1}{c} \cdot \left( B - \frac{Aa}{2} \right)\arctan\left(\frac{y}{c}\right) + C, \end{align*}

или, возвращаясь к xx и подставляя вместо cc его значение:

Ax+Bx2+ax+b)dx=A2ln(x2+ax+b)+2BAa4ba2arctan(2x+a4ba2)+C. \int \frac{Ax + B}{x^2 + ax + b)}\mathrm{d}x = \frac{A}{2}\ln(x^2 + ax + b) + \frac{2B-Aa}{\sqrt{4b-a^2}}\arctan\left( \frac{2x+a}{\sqrt{4b-a^2}} \right) + C.

(2) Пусть n>1n >1, делая ту же замену, получаем

Ax+B(x2+ax+b)ndx=Ax+(BAa2)(y2+c2)ndy=A22ydy(y2+c2)n+(BAa2)dy(y2+c2)n. \int \frac{Ax + B}{(x^2 + ax + b)^n}\mathrm{d}x = \int \frac{Ax + \left( B - \frac{Aa}{2} \right)}{(y^2 + c^2)^n}\mathrm{d}y = \frac{A}{2} \int \frac{2y\mathrm{d}y}{(y^2 + c^2)^n} + \left( B - \frac{Aa}{2} \right) \int \frac{\mathrm{d}y}{(y^2 + c^2)^n}.

Видим, что второй интеграл это интеграл от формы ω который найден в лемме , первый же интеграл легко берётся с помощью замены t:=y2+c2t:=y^2 + c^2 , тогда dt=(y2+c2)dy=2ydy\mathrm{d}t = (y^2 + c^2)'\mathrm{d}y = 2y\mathrm{d}y, следовательно ydy=12dty\mathrm{d}y = \frac{1}{2}\mathrm{d}t, и мы получаем

2ydy(y2+c2)n=dttn=1n11tn1+C. \int \frac{2y\mathrm{d}y}{(y^2 + c^2)^n} = \int \frac{\mathrm{d}t}{t^n} = - \frac{1}{n-1}\cdot \frac{1}{t^{n-1}} +C.

Тем самым лемма доказана.

Proof

Действительно, если дробь неправильная, то поделив P(x)=B(x)Q(x)+R(x)P(x) = B(x) Q(x) + R(x) на Q(x)Q(x) мы получаем

P(x)Q(x)dx=B(x)dx+R(x)Q(x)dx. \int \frac{P(x)}{Q(x)}\mathrm{d}x = \int B(x)\mathrm{d}x + \int \frac{R(x)}{Q(x)}\mathrm{d}x.

Первый интеграл находится легко, так как это интеграл от полинома, а во втором интеграле присутствует уже правильная дробь. Но согласно теореме Theorem 5 каждая правильная дробь представляется в виде суммы простых, тогда из лемм }, Lemma 3 вытекает утверждение теоремы.