Introduction
The standard objective of a substitution involving a quadratic radical is rationalization. We replace the radical by a rational function of a new parameter and then hope that the resulting rational integral is manageable.
For several important classes of integrals, however, the Unified Substitution Method (USM) does something substantially stronger.
It does not merely produce a rational function. After the radical and the Jacobian interact, the transformed integrand becomes a finite Laurent polynomial:
\[L(p)=\sum_{k=m}^{N}c_kp^k,\]
where \(p\) is one of the USM parameters \(t\), \(r\), or \(s\).
Consequently,
\[\int L(p)\,dp = \sum_{\substack{k=m\\k\neq-1}}^{N} \frac{c_k}{k+1}p^{k+1} +c_{-1}\ln|p|+C.\]
There is no partial-fraction decomposition, no Hermite reduction, and no recursive trigonometric integration. The calculation is reduced to integrating powers of a single parameter.
This “Laurent collapse” can be viewed as the systematic principle already visible in Examples 4 and 6 of the public arXiv preprint A Unified Substitution Method for Integration (see arXiv:2505.03754v3). In those examples, a quotient radical collapses to a Laurent polynomial in \(t\), while a circular-radical integral collapses to one in \(r\).
The parametrizations behind the collapse
Let
\[X=x+b,\qquad a>0.\]
The five families considered below arise primarily from the difference-radical and circular-radical transformations. A closely related Laurent mechanism also appears under Transform 3 for the sum radical.
The difference-radical parameter
On either exterior component \(X>a\) or \(X<-a\), let
\[R=\sqrt{X^2-a^2},\qquad \varepsilon=\operatorname{sgn}(X),\]
and set
\[t=\frac{X-\varepsilon R}{a}.\]
Then
\[X=\frac{a}{2}\left(t+t^{-1}\right),\]
\[R=-\varepsilon\frac{a}{2}\left(t-t^{-1}\right),\]
and
\[\,dx=\frac{a(t^2-1)}{2t^2}\,dt.\]
Thus \(0<t<1\) on the component \(X>a\), whereas \(-1<t<0\) on the component \(X<-a\). The formulas remain Laurent in \(t\), with the component-dependent sign carried by \(\varepsilon\).
The circular-radical parameter
For \(|X|<a\), set
\[r=\frac{X}{a+\sqrt{a^2-X^2}}.\]
Equivalently,
\[X=\frac{2ar}{1+r^2},\]
\[\sqrt{a^2-X^2} =a\frac{1-r^2}{1+r^2},\]
and
\[\,dx=2a\frac{1-r^2}{(1+r^2)^2}\,dr.\]
If \(X=a\sin\theta\) with the principal choice \(\theta=\sin^{-1}(X/a)\in(-\pi/2,\pi/2)\), then \(r=\tan(\theta/2)\), the classical tangent half-angle parameter.
The arXiv preprint explains how the corresponding USM transformations recover Euler's first and second substitutions, up to explicit sign, scaling, or reciprocal reparametrizations. Thus, the comparison below is not between unrelated methods: USM reorganizes the classical substitutions in a way that makes the cancellations easier to detect.
Transform 3: Laurent polynomials for the sum radical
The Laurent-polynomial phenomenon is not limited to the five families associated with Transforms 2 and 5. Transform 3 produces the same kind of algebraic simplification for integrals containing
\[\sqrt{X^2+a^2}.\]
Define
\[s=\frac{X+\sqrt{X^2+a^2}}{a}>0.\]
Then
\[X=\frac a2\left(s-s^{-1}\right),\]
\[\sqrt{X^2+a^2} =\frac a2\left(s+s^{-1}\right),\]
and
\[\,dx=\frac a2\left(1+s^{-2}\right)\,ds.\]
Each of the three fundamental quantities \(X\), \(\sqrt{X^2+a^2}\), and \(\,dx\) is therefore already a finite Laurent polynomial in \(s\).
It follows immediately that an expression of the form
\[Q\left(X,\sqrt{X^2+a^2}\right)\,dx,\]
where \(Q\) is a polynomial in its two arguments, becomes a finite Laurent polynomial in \(s\).
The same remains true for many integrands containing integer powers of
\[X+\sqrt{X^2+a^2}=as,\]
because multiplication or division by such a power merely introduces a Laurent monomial in \(s\).
For example, Example 3 of the arXiv preprint considers
\[\int_0^\infty \frac{\,dx}{\left(x+\sqrt{1+x^2}\right)^2}.\]
With
\[s=x+\sqrt{1+x^2},\]
we have
\[\,dx=\frac12(1+s^{-2})\,ds\]
and
\[\left(x+\sqrt{1+x^2}\right)^{-2}=s^{-2}.\]
Thus,
\[\int_0^\infty \frac{\,dx}{\left(x+\sqrt{1+x^2}\right)^2} = \frac12\int_1^\infty\left(s^{-2}+s^{-4}\right)\,ds = \frac23.\]
This is another Laurent collapse: the transformed integrand contains only two negative powers of \(s\).
Traditional comparison
The classical hyperbolic substitution
\[x=\sinh u\]
gives
\[x+\sqrt{1+x^2} =\sinh u+\cosh u =e^u.\]
The integral then becomes
\[\int_0^\infty e^{-2u}\cosh u\,du.\]
This is already manageable; a standard next step is to expand \(\cosh u\) into exponentials. Setting \(s=e^u\) then produces the same Laurent integral obtained directly by Transform 3.
Transform 3 therefore packages two classical steps,
\[x=\sinh u, \qquad s=e^u,\]
into the single algebraic parameter
\[s=x+\sqrt{1+x^2}.\]
Its advantage is not that hyperbolic substitution fails, but that the final Laurent structure is exposed from the beginning.
The five Laurent-polynomial families
The cancellation mechanism applies systematically to the following classes. Each statement is understood on a connected real component on which its integrand and substitution are defined.
For every polynomial \(P\), on either component \(X>a\) or \(X<-a\),
\[\int P(x)\sqrt{X^2-a^2}\,dx\]
becomes a finite Laurent polynomial in \(t\).
Indeed,
\[P(x) = P\left(\frac{a}{2}(t+t^{-1})-b\right)\]
is itself a Laurent polynomial, while
\[\sqrt{X^2-a^2}\,dx = -\varepsilon\frac{a^2}{4}t^{-3}(t^2-1)^2\,dt,\]
where
\[\varepsilon=\operatorname{sgn}(X)\]
is constant on the selected component.
Polynomial times a quotient radical
For every polynomial \(P\), on either component \(X>a\) or \(X<-a\),
\[\int P(x)\sqrt{\frac{X-a}{X+a}}\,dx\]
also becomes a finite Laurent polynomial in \(t\).
The decisive cancellation is
\[\sqrt{\frac{X-a}{X+a}}\,dx = -\frac a2(1-t^{-1})^2\,dt.\]
The linear factor \(1+t\) contributed by the denominator of the radical cancels a matching factor in the Jacobian.
Reciprocal powers divided by a circular radical
For every integer \(n\geq1\), on a connected component of \(0<|X|<a\),
\[\int\frac{\,dx}{X^n\sqrt{a^2-X^2}}\]
becomes
\[2^{1-n}a^{-n} \int r^{-n}(1+r^2)^{n-1}\,dr.\]
Because \(n-1\geq0\), the factor \((1+r^2)^{n-1}\) expands finitely.
A circular radical divided by reciprocal powers
For every integer \(n\geq3\), on a connected component of \(0<|X|<a\),
\[\int\frac{\sqrt{a^2-X^2}}{X^n}\,dx\]
becomes
\[2^{1-n}a^{2-n} \int r^{-n}(1-r^2)^2(1+r^2)^{n-3}\,dr.\]
The restriction \(n\geq3\) is exactly what prevents a residual denominator involving \(1+r^2\).
Reciprocal powers multiplied by a root ratio
For every integer \(n\geq2\), on a connected component of \(0<|X|<a\) and with principal square roots,
\[\int \frac1{X^n} \sqrt{\frac{a\pm X}{a\mp X}}\,dx\]
becomes
\[2^{1-n}a^{1-n} \int r^{-n}(1\pm r)^2(1+r^2)^{n-2}\,dr.\]Once again, the exponent \(n-2\) is nonnegative, so the transformed expression is a finite Laurent polynomial.
These formulas show that the restrictions on \(n\) are structural rather than arbitrary: they are precisely the thresholds at which all remaining powers of \(1+r^2\) move into the numerator.
Five examples: USM versus the traditional routes
Example 1: A polynomial times a difference radical
Consider
\[I_1=\int x^2\sqrt{x^2-1}\,dx, \qquad x>1.\]
Set
\[t=x-\sqrt{x^2-1}.\]
Then
\[x=\frac12(t+t^{-1}),\]
and direct substitution gives
\[I_1 = \int\left( -\frac{t^3}{16} +\frac1{8t} -\frac1{16t^5} \right)\,dt.\]
The transformed integrand contains only three Laurent monomials. Therefore,
\[I_1 = -\frac{t^4}{64} +\frac18\ln t +\frac{t^{-4}}{64} +C.\]
Using
\[t^{-1}=x+\sqrt{x^2-1},\]
we obtain
\[\boxed{ I_1= \frac18\left[ x\sqrt{x^2-1}(2x^2-1) -\ln\left(x+\sqrt{x^2-1}\right) \right]+C }.\]
Traditional comparison
The hyperbolic substitution \(x=\cosh u\) produces
\[\int \cosh^2u\,\sinh^2u\,du,\]
which is commonly handled with multiple-angle identities before back-substitution.
The trigonometric substitution \(x=\sec\theta\) produces an integral involving \(\sec^3\theta\tan^2\theta\), normally handled through identities and reduction formulas.
Both methods are valid. The USM advantage is that the complete algebraic structure is exposed immediately: three powers of \(t\), integrated term by term.
Example 2: The quotient-radical cancellation
Consider the integral appearing as Example 4 in the arXiv preprint:
\[I_2=\int\sqrt{\frac{x+1}{x+3}}\,dx, \qquad x>-1.\]
Here
\[X=x+2, \qquad a=1,\]
so the USM parameter is
\[t=x+2-\sqrt{x^2+4x+3}.\]
The radical and Jacobian combine to give
\[I_2 = -\frac12\int\left(1-2t^{-1}+t^{-2}\right)\,dt.\]
Thus,
\[I_2 = -\frac t2+\ln t+\frac1{2t}+C.\]
Since
\[\frac1{2t}-\frac t2 =\sqrt{x^2+4x+3},\]
the result is
\[\boxed{ I_2= \sqrt{x^2+4x+3} +\ln\left(x+2-\sqrt{x^2+4x+3}\right)+C }.\]
The important feature is not merely that the integral is solvable. It is that the transformed integrand is already a Laurent polynomial.
Traditional comparison
The natural direct rationalization is
\[u=\sqrt{\frac{x+1}{x+3}}.\]
Solving for \(x\) and differentiating gives
\[I_2=\int\frac{4u^2}{(u^2-1)^2}\,du.\]
The radical has disappeared, but the price is a repeated quadratic denominator. A standard continuation uses partial fractions or introduces a second hyperbolic substitution.
The USM parameter packages both steps into one transformation:
\[-\frac12\left(1-2t^{-1}+t^{-2}\right).\]
Moreover, multiplying the original integrand by any polynomial \(P(x)\) preserves the Laurent form, because
\[x=\frac12(t+t^{-1})-2.\]
The direct \(u\)-substitution, by contrast, generally increases the degree of the rational denominator as the polynomial factor becomes more complicated.
This is the first seed of the general Laurent-polynomial phenomenon: Example 4 demonstrates the cancellation in one particular integral, while the corresponding family shows that it persists for every polynomial multiplier.
Example 3: Reciprocal powers and a circular radical
Now consider Example 6 from the arXiv preprint:
\[I_3=\int\frac{\,dx}{x^3\sqrt{4-x^2}}, \qquad 0<x<2.\]
Set
\[r=\frac{2-\sqrt{4-x^2}}{x} =\frac{x}{2+\sqrt{4-x^2}}.\]
Transform 5 gives
\[I_3 = \frac1{32}\int\left(r^{-3}+2r^{-1}+r\right)\,dr.\]
Therefore,
\[I_3 = -\frac1{64r^2} +\frac1{16}\ln r +\frac{r^2}{64} +C.\]
After combining the reciprocal powers,
\[\boxed{ I_3= \frac1{16}\left[ \ln\left(\frac{2-\sqrt{4-x^2}}{x}\right) -\frac{2\sqrt{4-x^2}}{x^2} \right]+C }.\]
Traditional comparison
The usual circular substitution
\[x=2\sin\theta\]
produces
\[I_3=\frac18\int\csc^3\theta\,d\theta.\]
That integral is standard and is commonly evaluated by a reduction formula or an integration-by-parts derivation, followed by the reconstruction of both \(\csc\theta\) and \(\cot\theta\) in terms of \(x\).
The USM route bypasses the reduction formula entirely. The radical and Jacobian cancel before integration, leaving exactly three Laurent monomials.
This is the second worked example in the arXiv preprint from which the general Laurent classification naturally emerges.
Example 4: A circular radical over a high reciprocal power
Consider
\[I_4=\int\frac{\sqrt{1-x^2}}{x^5}\,dx, \qquad 0<x<1.\]
Use
\[r=\frac{x}{1+\sqrt{1-x^2}}.\]
This is the fourth family with \(a=1\) and \(n=5\). The general formula gives\[I_4 = \frac1{16}\int r^{-5}(1-r^2)^2(1+r^2)^2\,dr.\]
But
\[(1-r^2)^2(1+r^2)^2=(1-r^4)^2,\]
so
\[I_4 = \frac1{16}\int\left(r^{-5}-2r^{-1}+r^3\right)\,dr.\]
Hence
\[I_4 = \frac{r^4-r^{-4}}{64} -\frac18\ln r+C.\]
Returning to \(x\),
\[\boxed{ I_4= \frac18\left[ \ln\left(\frac{1+\sqrt{1-x^2}}{x}\right) -\frac{(2-x^2)\sqrt{1-x^2}}{x^4} \right]+C }.\]
Traditional comparison
With \(x=\sin\theta\), one obtains
\[I_4 = \int\frac{\cos^2\theta}{\sin^5\theta}\,d\theta = \int\left(\csc^5\theta-\csc^3\theta\right)\,d\theta.\]
A standard evaluation uses reduction formulas for two odd powers of \(\csc\theta\). It is entirely feasible, but the method treats the integral as a trigonometric-reduction problem.
USM instead detects the algebraic factorization
\[r^{-5}-2r^{-1}+r^3\]
automatically. The same pattern persists for every \(n\geq3\); increasing \(n\) merely enlarges a finite binomial expansion.
Example 5: Reciprocal powers times a root ratio
Consider
\[I_5= \int\frac1{x^3}\sqrt{\frac{1+x}{1-x}}\,dx, \qquad 0<x<1.\]
Again set
\[r=\frac{x}{1+\sqrt{1-x^2}}.\]
Because
\[\sqrt{\frac{1+x}{1-x}} =\frac{1+r}{1-r},\]
the fifth Laurent family gives
\[I_5 = \frac14\int r^{-3}(1+r)^2(1+r^2)\,dr.\]
Expanding,
\[I_5 = \frac14\int\left(r^{-3}+2r^{-2}+2r^{-1}+2+r\right)\,dr.\]
Thus,
\[I_5 = -\frac1{8r^2} -\frac1{2r} +\frac12\ln r +\frac r2 +\frac{r^2}{8} +C.\]
Combining reciprocal pairs gives
\[\boxed{ I_5= -\frac{\sqrt{1-x^2}}{2x^2} -\frac{\sqrt{1-x^2}}{x} +\frac12\ln\left(\frac{x}{1+\sqrt{1-x^2}}\right)+C }.\]
The conjugate family
\[\frac1{x^n}\sqrt{\frac{1-x}{1+x}}\]
is handled identically, with \((1+r)^2\) replaced by \((1-r)^2\).
Traditional comparison
The substitution \(x=\sin\theta\) gives
\[\sqrt{\frac{1+\sin\theta}{1-\sin\theta}} =\frac{1+\sin\theta}{\cos\theta}\]
on the relevant component. Consequently,
\[I_5 = \int\left(\csc^3\theta+\csc^2\theta\right)\,d\theta.\]
A standard evaluation combines the odd-cosecant reduction formula with a separate elementary integration of \(\csc^2\theta\).
USM reaches the same trigonometric half-angle parameter without first converting the problem into trigonometric notation. More importantly, it reveals immediately that the entire class for \(n\geq2\) is Laurent-polynomial after substitution.
Compressing \(t^n\pm t^{-n}\): the binomial difference and sum formulas
Let \(n\in\mathbb N\) and \(y\in(-\infty,-1]\cup[1,\infty)\). Termwise integration often produces reciprocal pairs such as
\[t^n-t^{-n} \qquad\text{and}\qquad t^n+t^{-n}.\]
Let
\[w=\sqrt{y^2-1}.\]
For \(y\geq1\), one has \(t=y-w\) and \(t^{-1}=y+w\); for \(y\leq-1\), these two reciprocal roots exchange roles. Treating the two roots as an unordered pair avoids unnecessary expression swell.
Indeed,
\[(y-w)(y+w)=1,\]
so
\[\{y-w,y+w\}=\{t,t^{-1}\}.\]
Define
\[D_n(y)=(y-w)^n-(y+w)^n.\]
By the binomial theorem, all even powers of \(w\) cancel, leaving only the odd terms:
\[D_n(y) = -2\sum_{j=0}^{\lfloor(n-1)/2\rfloor} \binom{n}{2j+1} y^{n-2j-1}w^{2j+1}.\]
Since
\[w^{2j+1}=(y^2-1)^jw,\]
this may also be written as
\[D_n(y) = -2\sqrt{y^2-1} \sum_{j=0}^{\lfloor(n-1)/2\rfloor} \binom{n}{2j+1} y^{n-2j-1}(y^2-1)^j.\]
For the USM parameter chosen as the small root on both exterior components,
\[t= \begin{cases} y-\sqrt{y^2-1},&y\geq1,\\[4pt] y+\sqrt{y^2-1},&y\leq-1, \end{cases}\]
the result is
\[t^n-t^{-n}=\sigma D_n(y),\]
where
\[\sigma= \begin{cases} 1,&y\geq1,\\ -1,&y\leq-1. \end{cases}\]
The factor \(\sigma\) is essential: the roles of \(y-w\) and \(y+w\) are interchanged on the negative exterior component.
The first few identities are
\[D_1(y)=-2\sqrt{y^2-1},\]
\[D_2(y)=-4y\sqrt{y^2-1},\]
\[D_3(y)=-2(4y^2-1)\sqrt{y^2-1},\]
and
\[D_4(y)=-8y(2y^2-1)\sqrt{y^2-1}.\]
There is an equally useful companion formula for the sum. Define
\[S_n(y)=(y-w)^n+(y+w)^n.\]
The odd powers now cancel, giving
\[S_n(y) = 2\sum_{j=0}^{\lfloor n/2\rfloor} \binom{n}{2j} y^{n-2j}w^{2j},\]
or
\[S_n(y) = 2\sum_{j=0}^{\lfloor n/2\rfloor} \binom{n}{2j} y^{n-2j}(y^2-1)^j.\]
Unlike the difference, the sum is unchanged when the two reciprocal roots are interchanged. Therefore,
\[t^n+t^{-n}=S_n(y)\]
on both exterior components, with no additional sign factor. Equivalently, in terms of the Chebyshev polynomials \(T_n\) and \(U_{n-1}\),
\[S_n(y)=2T_n(y),\qquad D_n(y)=-2\sqrt{y^2-1}\,U_{n-1}(y).\]
For example,
\[t^2+t^{-2}=4y^2-2,\]
and
\[t^4+t^{-4}=16y^4-16y^2+2.\]
These formulas shorten the final back-substitution substantially. In Example 1, the Laurent primitive contains
\[\frac{t^{-4}-t^4}{64}.\]
For \(x>1\), the difference formula gives
\[t^4-t^{-4} =-8x(2x^2-1)\sqrt{x^2-1}.\]
Hence
\[\frac{t^{-4}-t^4}{64} = \frac18x(2x^2-1)\sqrt{x^2-1},\]
which recovers the algebraic part of the antiderivative in one step.
The arXiv preprint explicitly presents the binomial-difference identity as a device for simplifying terms of the form \(t^n-t^{-n}\) during back-substitution (see arXiv:2505.03754v3).
Why Laurent polynomials are better than general rational functions
The difference is computationally meaningful.
A general rational function
\[\frac{A(p)}{B(p)}\]
may require factorization of \(B\), polynomial division, repeated-factor decomposition, irreducible quadratic terms, or Hermite reduction.
A Laurent polynomial
\[\sum c_kp^k\]
requires none of these operations. Its integration complexity is visible at inspection:
• every exponent except \(-1\) produces another power;
• the exponent \(-1\) produces a logarithm;
• no new algebraic denominator is introduced;
• the number of terms is controlled by finite polynomial or binomial expansions.
There is also an advantage during back-substitution. Positive and negative powers often occur in pairs such as
\[t^m-t^{-m} \qquad\text{or}\qquad t^m+t^{-m}.\]
The binomial difference and sum formulas convert these pairs directly into expressions involving \(y\) and \(\sqrt{y^2-1}\). This avoids separately expanding two large reciprocal powers and then simplifying the result.
That is why the final answers in the examples above can be written using only the original radical and one logarithm, despite passing through several Laurent powers.
What USM is—and is not—claiming
The Laurent-polynomial property is a precise structural statement. For the stated families and domain restrictions, the transformed integrand has a finite Laurent expansion.
It does \emph{not} imply that every individual integral will have a shorter final antiderivative than every classical derivation. For small exponents, a familiar trigonometric substitution may be equally short. Nor does it constitute a universal theorem about computer-algebra performance.
The real advantage is systematic predictability:
1. the appropriate parameter is selected from the radical geometry;
2. the branch and real component are stated explicitly;
3. cancellation occurs before integration;
4. the resulting integral is termwise;
5. reciprocal powers can be recombined before back-substitution;
6. the procedure scales uniformly with \(P(x)\) or \(n\).
Traditional substitutions usually teach these examples as separate categories: Euler substitutions for difference radicals, circular substitutions for \(\sqrt{a^2-x^2}\), hyperbolic substitutions for \(\sqrt{x^2+a^2}\), half-angle identities for root ratios, and reduction formulas for powers of secant or cosecant.
USM shows that they are manifestations of a common algebraic mechanism.
Final perspective
Examples 4 and 6 of the arXiv preprint display two striking simplifications:
\[\sqrt{\frac{x+1}{x+3}}\,dx \quad\longrightarrow\quad -\frac12\left(1-2t^{-1}+t^{-2}\right)\,dt,\]
and
\[\frac{\,dx}{x^3\sqrt{4-x^2}} \quad\longrightarrow\quad \frac1{32}\left(r^{-3}+2r^{-1}+r\right)\,dr.\]
Transform 3 supplies a parallel sum-radical example:
\[\frac{1}{\left(x+\sqrt{1+x^2}\right)^2}\,dx \quad\longrightarrow\quad \frac12\left(s^{-2}+s^{-4}\right)\,ds.\]
The five systematic families explain the deeper pattern behind the first two examples, while Transform 3 shows that the Laurent principle extends naturally beyond those families.
The best substitution does not merely remove the radical. It makes the radical, the powers of the centered variable, and the Jacobian simplify one another before the integration begins.
Once the integration has been performed, the binomial difference and sum formulas complete the process by compressing reciprocal pairs such as
\[t^n\pm t^{-n}\]
directly into the original variable and radical.
For these classes, USM moves the problem beyond rationalization. What initially appears to be a radical integral becomes finite arithmetic with powers—and, at worst, one logarithm.
