This is my work on Chain Integration and Integration by Parts (higher-order integrals). Also included is a differential analysis on cycloids.
Chain Integration Formula and Examples
edit
See Fresnel integral .
We could approximate the tangent and secant integrals
T
(
x
)
:=
∫
0
x
tan
(
t
2
)
d
t
{\displaystyle T(x):=\int _{0}^{x}\tan {\big (}t^{2}{\big )}~{\text{d}}t}
and
Sc
(
x
)
:=
∫
0
x
sec
(
t
2
)
d
t
{\displaystyle {\text{Sc}}(x):=\int _{0}^{x}\sec {\big (}t^{2}{\big )}~{\text{d}}t}
by using the Cauchy Principal Value and integration by parts :
T
(
x
)
=
ln
|
sec
x
2
|
2
x
+
∫
0
x
ln
|
sec
t
2
|
d
t
2
t
2
{\displaystyle T(x)={\frac {\ln |\sec x^{2}|}{2x}}+\int _{0}^{x}{\frac {\ln |\sec t^{2}|~{\text{d}}t}{2t^{2}}}}
and
Sc
(
x
)
=
ln
|
tan
x
2
+
sec
x
2
|
2
x
+
∫
0
x
ln
|
tan
t
2
+
sec
t
2
|
d
t
2
t
2
{\displaystyle {\text{Sc}}(x)={\frac {\ln |\tan x^{2}+\sec x^{2}|}{2x}}+\int _{0}^{x}{\frac {\ln |\tan t^{2}+\sec t^{2}|~{\text{d}}t}{2t^{2}}}}
with the square Chain Integration formula
∫
0
x
f
(
t
2
)
d
t
=
F
(
x
2
)
2
x
+
∫
0
x
F
(
t
2
)
d
t
2
x
2
{\displaystyle \int _{0}^{x}f{\big (}t^{2}{\big )}~{\text{d}}t={\frac {F{\big (}x^{2}{\big )}}{2x}}+{\frac {\int _{0}^{x}F{\big (}t^{2}{\big )}~{\text{d}}t}{2x^{2}}}}
where
F
(
x
)
=
∫
0
x
f
(
t
)
d
t
{\displaystyle F(x)=\int _{0}^{x}f(t)~{\text{d}}t}
.
In fact, integrals through CPV are defined when subtracting balanced pole functions - pure powers of
1
/
x
{\displaystyle 1/x}
- from original poles, results in functions with remaining possible singularities of size
O
(
n
−
α
)
{\displaystyle O(n^{-\alpha })}
with
α
<
1
{\displaystyle \alpha <1}
strict).
Notice that both
tan
x
2
{\displaystyle \tan x^{2}}
and
sec
x
2
{\displaystyle \sec x^{2}}
have simple nonzero isolated poles, limiting to scalar multiples of
[
∓
1
/
x
]
x
=
0
{\displaystyle [\mp 1/x]_{x=0}}
. Then
∫
1
x
∓
d
t
t
=
∓
ln
x
{\displaystyle \int _{1}^{x}\mp {\frac {{\text{d}}t}{t}}=\mp \ln x}
and
∫
0
x
∓
ln
t
d
t
=
∓
x
(
ln
x
−
1
)
{\displaystyle \int _{0}^{x}\mp \ln t~{\text{d}}t=\mp x(\ln x-1)}
with
lim
x
→
0
∓
[
x
(
ln
x
−
1
)
]
=
0
{\displaystyle \lim _{x\rightarrow 0}\mp {\big [}x(\ln x-1){\big ]}=0}
.
Therefore,
∫
0
x
F
(
t
2
)
d
t
2
x
2
{\displaystyle {\frac {\int _{0}^{x}F{\big (}t^{2}{\big )}~{\text{d}}t}{2x^{2}}}}
is bounded for
F
(
x
)
=
ln
|
sec
x
|
=
∫
0
x
tan
t
d
t
{\displaystyle F(x)=\ln |\sec x|=\int _{0}^{x}\tan t~{\text{d}}t}
and
F
(
x
)
=
ln
|
tan
x
+
sec
x
|
=
∫
0
x
sec
t
d
t
{\displaystyle F(x)=\ln |\tan x+\sec x|=\int _{0}^{x}\sec t~{\text{d}}t}
; as we have the related bounded integral
∫
0
1
∓
ln
x
d
x
=
∓
[
x
(
ln
x
−
1
)
]
0
1
=
±
1
{\displaystyle \int _{0}^{1}\mp \ln x~{\text{d}}x=\mp [x(\ln x-1)]_{0}^{1}=\pm 1}
and we can do scalar multiple comparisons . So by CPV , the above integrals are defined except at isolated poles. Graphs of these integrals for
x
≥
0
{\displaystyle x\geq 0}
are found below:
Graphs of
T
(
x
)
{\displaystyle T(x)}
and
Sc
(
x
)
{\displaystyle {\text{Sc}}(x)}
Integral of Tan x^2
Integral of Sec x^2
[0,15] × [-5, 5]
[0,15] × [-5, 5]
Finally, we have, approximately:
lim inf
x
→
∞
∫
0
x
tan
(
t
2
)
d
t
=
lim
x
→
∞
∫
0
x
ln
|
sec
t
2
|
d
t
2
t
2
≈
0.535914345467
{\displaystyle \liminf _{x\rightarrow \infty }\int _{0}^{x}\tan {\big (}t^{2}{\big )}~{\text{d}}t=\lim _{x\rightarrow \infty }\int _{0}^{x}{\frac {\ln |\sec t^{2}|~{\text{d}}t}{2t^{2}}}\approx 0.535914345467}
lim inflection
x
→
∞
∫
0
x
sec
(
t
2
)
d
t
=
lim
x
→
∞
∫
0
x
ln
|
tan
t
2
+
sec
t
2
|
d
t
2
t
2
≈
0.836834795146
{\displaystyle {\underset {x\rightarrow \infty }{\text{lim inflection}}}\int _{0}^{x}\sec {\big (}t^{2}{\big )}~{\text{d}}t=\lim _{x\rightarrow \infty }\int _{0}^{x}{\frac {\ln |\tan t^{2}+\sec t^{2}|~{\text{d}}t}{2t^{2}}}\approx 0.836834795146}
compared to
π
/
8
≈
0.62665706866
{\displaystyle {\sqrt {\pi /8}}\approx 0.62665706866}
for both
S
(
x
)
{\displaystyle S(x)}
and
C
(
x
)
{\displaystyle C(x)}
.
Modulated Integrals
edit
We could also approximate the cotangent and cosecant integrals
Ct
(
x
)
:=
∫
0
x
cot
(
t
2
)
d
t
{\displaystyle {\text{Ct}}(x):=\int _{0}^{x}\cot {\big (}t^{2}{\big )}~{\text{d}}t}
and
Cc
(
x
)
:=
∫
0
x
csc
(
t
2
)
d
t
{\displaystyle {\text{Cc}}(x):=\int _{0}^{x}\csc {\big (}t^{2}{\big )}~{\text{d}}t}
by using an integration by parts, but need to isolate a pole of order 2 at
x
=
0
{\displaystyle x=0}
for each function. We do so by subtracting
1
/
x
2
{\displaystyle 1/x^{2}}
from each function to yield bounded functions
f
(
x
)
{\displaystyle f(x)}
at
x
=
0
{\displaystyle x=0}
(in fact, with
lim
x
→
0
f
(
x
)
=
0
{\displaystyle \lim _{x\rightarrow 0}f(x)=0}
for both functions!), applying the same treatment
∫
0
x
f
(
t
2
)
d
t
=
F
(
x
2
)
2
x
+
∫
0
x
F
(
t
2
)
d
t
2
x
2
{\displaystyle \int _{0}^{x}f{\big (}t^{2}{\big )}~{\text{d}}t={\frac {F{\big (}x^{2}{\big )}}{2x}}+{\frac {\int _{0}^{x}F{\big (}t^{2}{\big )}~{\text{d}}t}{2x^{2}}}}
as previously (it works similarly at all other poles), and then adding the antiderivative
−
1
/
x
{\displaystyle -1/x}
of
1
/
x
2
{\displaystyle 1/x^{2}}
back.
Since both integrals have right limit
−
∞
{\displaystyle -\infty }
as
x
{\displaystyle x}
approaches
0
{\displaystyle 0}
, we instead add constants so that critical points/inflection points of the cotangent/cosecant integrals, respectively, approach
y
=
0
{\displaystyle y=0}
as
x
→
∞
{\displaystyle x\rightarrow \infty }
(canonicalization). We thus have:
Ct
(
x
)
:≈
1.29414
−
(
ln
2
−
ln
(
|
cot
x
2
+
csc
x
2
|
)
−
2
ln
(
|
x
|
)
2
x
2
+
1
x
+
∫
0
x
ln
2
−
ln
(
|
cot
t
2
+
csc
t
2
|
)
−
2
ln
(
|
t
|
)
2
t
2
d
t
)
{\displaystyle {\text{Ct}}(x):\approx 1.29414-\left({\frac {\ln 2-\ln \left(\left|\cot x^{2}+\csc x^{2}\right|\right)-2\ln \left(\left|x\right|\right)}{2x^{2}}}+{\frac {1}{x}}+\int _{0}^{x}{\frac {\ln 2-\ln \left(\left|\cot t^{2}+\csc t^{2}\right|\right)-2\ln \left(\left|t\right|\right)}{2t^{2}}}dt\right)}
Cc
(
x
)
:≈
ln
2
−
ln
(
|
cot
x
2
+
csc
x
2
|
)
−
2
ln
(
|
x
|
)
2
x
−
1
x
+
∫
0
x
ln
2
−
ln
(
|
cot
t
2
+
csc
t
2
|
)
−
2
ln
(
|
t
|
)
2
t
2
d
t
+
0.5360919
{\displaystyle {\text{Cc}}(x):\approx {\frac {\ln 2-\ln \left(\left|\cot x^{2}+\csc x^{2}\right|\right)-2\ln \left(\left|x\right|\right)}{2x}}-{\frac {1}{x}}+\int _{0}^{x}{\frac {\ln 2-\ln \left(\left|\cot t^{2}+\csc t^{2}\right|\right)-2\ln \left(\left|t\right|\right)}{2t^{2}}}dt+0.5360919}
With graphs seen below:
Graphs of
Ct
(
x
)
{\displaystyle {\text{Ct}}(x)}
and
Cc
(
x
)
{\displaystyle {\text{Cc}}(x)}
Integral of Cot x^2
Integral of csc x^2
[0,15] × [-5, 5]
[0,15] × [-5, 5]
Higher Order (Stacked) Integrals
edit
A general formula for the second-order real integral (second antiderivative ) is
∬
0
x
f
(
t
)
d
t
2
=
x
∫
0
x
f
(
t
)
d
t
−
∫
0
x
t
f
(
t
)
d
t
{\displaystyle \iint _{0}^{x}f(t)~{\text{d}}t^{2}=x\int _{0}^{x}f(t)~{\text{d}}t-\int _{0}^{x}tf(t)~{\text{d}}t}
Proof
d
d
x
(
x
∫
0
x
f
(
t
)
d
t
−
∫
0
x
t
f
(
t
)
d
t
)
=
x
f
(
x
)
+
∫
0
x
f
(
t
)
d
t
−
x
f
(
x
)
=
∫
0
x
f
(
t
)
d
t
;
{\displaystyle {\frac {\text{d}}{{\text{d}}x}}\left(x\int _{0}^{x}f(t)~{\text{d}}t-\int _{0}^{x}tf(t)~{\text{d}}t\right)=xf(x)+\int _{0}^{x}f(t)~{\text{d}}t-xf(x)=\int _{0}^{x}f(t)~{\text{d}}t;}
by Product Rule and Fundamental Theorem of Calculus I .
A general formula for the third-order real integral (third antiderivative ) is
∭
0
x
f
(
t
)
d
t
3
=
1
2
(
x
2
∫
0
x
f
(
t
)
d
t
−
2
x
∫
0
x
t
f
(
t
)
d
t
+
∫
0
x
t
2
f
(
t
)
d
t
)
{\displaystyle \iiint _{0}^{x}f(t)~{\text{d}}t^{3}={\frac {1}{2}}\left(x^{2}\int _{0}^{x}f(t)~{\text{d}}t-2x\int _{0}^{x}tf(t)~{\text{d}}t+\int _{0}^{x}t^{2}f(t)~{\text{d}}t\right)}
Proof
d
d
x
(
1
2
(
x
2
∫
0
x
f
(
t
)
d
t
−
2
x
∫
0
x
t
f
(
t
)
d
t
+
∫
0
x
t
2
f
(
t
)
d
t
)
)
=
1
2
(
x
2
f
(
x
)
+
2
x
∫
0
x
f
(
t
)
d
t
−
2
x
2
f
(
x
)
−
2
∫
0
x
t
f
(
t
)
d
t
+
x
2
f
(
x
)
)
=
1
2
(
2
x
∫
0
x
f
(
t
)
d
t
−
2
∫
0
x
t
f
(
t
)
d
t
)
=
∬
0
x
f
(
t
)
d
t
2
;
{\displaystyle {\frac {\text{d}}{{\text{d}}x}}\left({\frac {1}{2}}\left(x^{2}\int _{0}^{x}f(t)~{\text{d}}t-2x\int _{0}^{x}tf(t)~{\text{d}}t+\int _{0}^{x}t^{2}f(t)~{\text{d}}t\right)\right)={\frac {1}{2}}\left(x^{2}f(x)+2x\int _{0}^{x}f(t)~{\text{d}}t-2x^{2}f(x)-2\int _{0}^{x}tf(t)~{\text{d}}t+x^{2}f(x)\right)={\frac {1}{2}}\left(2x\int _{0}^{x}f(t)~{\text{d}}t-2\int _{0}^{x}tf(t)~{\text{d}}t\right)=\iint _{0}^{x}f(t)~{\text{d}}t^{2};}
by Product Rule and Fundamental Theorem of Calculus I .
A general formula for the
n
{\displaystyle n}
th-order real integral (
n
{\displaystyle n}
th antiderivative ) is
∫
⋯
∫
0
x
n
times
f
(
t
)
d
t
n
=
1
(
n
−
1
)
!
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
[
x
n
−
1
−
i
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle {\underset {n~{\text{times}}}{\int \cdots \int _{0}^{x}}}f(t)~{\text{d}}t^{n}={\frac {1}{(n-1)!}}\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}\left[x^{n-1-i}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
Proof
Assume by induction , that
∫
⋯
∫
0
x
n
−
1
times
f
(
t
)
d
t
n
−
1
=
1
(
n
−
2
)
!
∑
i
=
1
n
−
1
(
−
1
)
i
(
n
−
2
i
)
[
x
n
−
2
−
i
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle {\underset {n-1~{\text{times}}}{\int \cdots \int _{0}^{x}}}f(t)~{\text{d}}t^{n-1}={\frac {1}{(n-2)!}}\sum _{i=1}^{n-1}(-1)^{i}{\binom {n-2}{i}}\left[x^{n-2-i}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
Then,
d
d
x
(
1
(
n
−
1
)
!
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
[
x
n
−
i
−
1
∫
0
x
t
i
f
(
t
)
d
t
]
)
=
1
(
n
−
1
)
!
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
[
x
n
−
i
−
1
x
i
f
(
x
)
+
(
n
−
i
−
1
)
x
n
−
i
−
2
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle {\frac {\text{d}}{{\text{d}}x}}\left({\frac {1}{(n-1)!}}\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}\left[x^{n-i-1}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]\right)={\frac {1}{(n-1)!}}\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}\left[x^{n-i-1}x^{i}f(x)+(n-i-1)x^{n-i-2}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
=
1
(
n
−
1
)
!
[
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
x
n
−
1
f
(
x
)
+
∑
i
=
1
n
(
−
1
)
n
(
n
−
1
i
)
∫
0
x
t
i
f
(
t
)
d
t
]
=
1
(
n
−
1
)
!
[
x
n
−
1
f
(
x
)
(
1
+
(
−
1
)
)
n
−
1
+
∑
i
=
1
n
(
−
1
)
n
(
n
−
1
)
!
i
!
(
n
−
i
−
1
)
!
(
n
−
i
−
1
)
x
n
−
i
−
2
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle ={\frac {1}{(n-1)!}}\left[\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}x^{n-1}f(x)+\sum _{i=1}^{n}(-1)^{n}{\binom {n-1}{i}}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]={\frac {1}{(n-1)!}}\left[x^{n-1}f(x){\big (}1+(-1){\big )}^{n-1}+\sum _{i=1}^{n}(-1)^{n}{\frac {(n-1)!}{i!(n-i-1)!}}(n-i-1)x^{n-i-2}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
=
1
(
n
−
1
)
!
[
x
n
−
1
f
(
x
)
⋅
0
n
−
1
+
∑
i
=
1
n
(
−
1
)
n
(
n
−
1
)
!
i
!
(
n
−
i
−
2
)
!
x
n
−
i
−
2
∫
0
x
t
i
f
(
t
)
d
t
]
=
1
(
n
−
1
)
!
[
(
n
−
1
)
∑
i
=
1
n
(
−
1
)
n
(
n
−
2
)
!
i
!
(
n
−
i
−
2
)
!
x
n
−
i
−
2
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle ={\frac {1}{(n-1)!}}\left[x^{n-1}f(x)\cdot 0^{n-1}+\sum _{i=1}^{n}(-1)^{n}{\frac {(n-1)!}{i!(n-i-2)!}}x^{n-i-2}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]={\frac {1}{(n-1)!}}\left[(n-1)\sum _{i=1}^{n}(-1)^{n}{\frac {(n-2)!}{i!(n-i-2)!}}x^{n-i-2}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
=
1
(
n
−
2
)
!
∑
i
=
1
n
−
1
(
−
1
)
i
(
n
−
2
i
)
[
x
n
−
2
−
i
∫
0
x
t
i
f
(
t
)
d
t
]
=
∫
⋯
∫
0
x
n
−
1
times
f
(
t
)
d
t
n
−
1
{\displaystyle ={\frac {1}{(n-2)!}}\sum _{i=1}^{n-1}(-1)^{i}{\binom {n-2}{i}}\left[x^{n-2-i}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]={\underset {n-1~{\text{times}}}{\int \cdots \int _{0}^{x}}}f(t)~{\text{d}}t^{n-1}}
by Product Rule , Binomial formula , combination formulas , and Fundamental Theorem of Calculus I , which implies
∫
⋯
∫
0
x
n
times
f
(
t
)
d
t
n
=
1
(
n
−
1
)
!
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
[
x
n
−
1
−
i
∫
0
x
t
i
f
(
t
)
d
t
]
{\displaystyle {\underset {n~{\text{times}}}{\int \cdots \int _{0}^{x}}}f(t)~{\text{d}}t^{n}={\frac {1}{(n-1)!}}\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}\left[x^{n-1-i}\int _{0}^{x}t^{i}f(t)~{\text{d}}t\right]}
and we are done. Intuitively, terms in the derivative of the form
x
n
−
1
f
(
x
)
{\displaystyle x^{n-1}f(x)}
come from differentiating each integral (but not the power in front of it, thereby obtaining the right side
f
g
′
{\displaystyle fg'}
in
(
f
g
)
′
=
f
′
g
+
f
g
′
{\displaystyle (fg)'=f'g+fg'}
contributing from each integral term), and cancel out by means of
0
=
0
n
−
1
=
(
1
+
(
−
1
)
)
n
−
1
=
∑
i
=
1
n
(
−
1
)
i
(
n
−
1
i
)
{\displaystyle 0=0^{n-1}={\big (}1+(-1){\big )}^{n-1}=\sum _{i=1}^{n}(-1)^{i}{\binom {n-1}{i}}}
; but terms of the form
x
n
−
i
−
1
∫
0
x
t
i
f
(
t
)
d
t
{\displaystyle x^{n-i-1}\int _{0}^{x}t^{i}f(t)~{\text{d}}t}
are few, with only one contribution from each integral term.
So to compute higher-order integrals, no other nonelementary integrals need to be considered except for those possibly equal to
∫
0
x
t
n
f
(
t
)
d
t
{\displaystyle \int _{0}^{x}t^{n}f(t)~{\text{d}}t}
for
f
∈
R
,
n
∈
N
{\displaystyle f\in {\mathcal {R}},~n\in \mathbb {N} }
(Riemann-integrable functions). To compute from another bound
c
≠
0
{\displaystyle c\neq 0}
, it is sufficient for integrals only of the form
∫
0
x
(
t
−
c
)
n
f
(
t
)
d
t
{\displaystyle \int _{0}^{x}(t-c)^{n}f(t)~{\text{d}}t}
to be considered.
Below, the first five antiderivatives of
f
(
x
)
=
e
x
2
+
1
/
x
2
{\displaystyle f(x)=e^{x^{2}+1/x^{2}}}
are computed and graphed, using this method. Since
f
{\displaystyle f}
is undefined at
x
=
0
{\displaystyle x=0}
(in fact,
x
=
0
{\displaystyle x=0}
is an essential singularity , and every antiderivative diverges as
x
→
0
{\displaystyle x\rightarrow 0}
), we start at
x
=
1
{\displaystyle x=1}
instead (a critical minimum point for
f
{\displaystyle f}
, so most suitable):
e
x
2
+
1
/
x
2
{\displaystyle e^{x^{2}+1/x^{2}}}
∫
1
x
e
t
2
+
1
/
t
2
d
t
{\displaystyle \int _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t}
∬
1
x
e
t
2
+
1
/
t
2
d
t
2
=
(
x
−
1
)
∫
1
x
e
t
2
+
1
/
t
2
d
t
−
∫
1
x
(
t
−
1
)
e
t
2
+
1
/
t
2
d
t
{\displaystyle \iint _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t^{2}=(x-1)\int _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t-\int _{1}^{x}(t-1)e^{t^{2}+1/t^{2}}~{\text{d}}t}
∭
1
x
e
t
2
+
1
/
t
2
d
t
3
=
1
2
[
(
x
−
1
)
2
∫
1
x
e
t
2
+
1
/
t
2
d
t
−
2
(
x
−
1
)
∫
1
x
(
t
−
1
)
e
t
2
+
1
/
t
2
d
t
+
∫
1
x
(
t
−
1
)
2
e
t
2
+
1
/
t
2
d
t
]
{\displaystyle \iiint _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t^{3}={\frac {1}{2}}\left[(x-1)^{2}\int _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t-2(x-1)\int _{1}^{x}(t-1)e^{t^{2}+1/t^{2}}~{\text{d}}t+\int _{1}^{x}(t-1)^{2}e^{t^{2}+1/t^{2}}~{\text{d}}t\right]}
∫
⋯
∫
1
x
4
times
e
t
2
+
1
/
t
2
d
t
4
=
1
6
[
(
x
−
1
)
3
∫
1
x
e
t
2
+
1
/
t
2
d
t
−
3
(
x
−
1
)
2
∫
1
x
(
t
−
1
)
e
t
2
+
1
/
t
2
d
t
+
3
(
x
−
1
)
∫
1
x
(
t
−
1
)
2
e
t
2
+
1
/
t
2
d
t
−
∫
1
x
(
t
−
1
)
3
e
t
2
+
1
/
t
2
d
t
]
{\displaystyle {\underset {4~{\text{times}}}{\int \cdots \int _{1}^{x}}}e^{t^{2}+1/t^{2}}~{\text{d}}t^{4}={\frac {1}{6}}\left[(x-1)^{3}\int _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t-3(x-1)^{2}\int _{1}^{x}(t-1)e^{t^{2}+1/t^{2}}~{\text{d}}t+3(x-1)\int _{1}^{x}(t-1)^{2}e^{t^{2}+1/t^{2}}~{\text{d}}t-\int _{1}^{x}(t-1)^{3}e^{t^{2}+1/t^{2}}~{\text{d}}t\right]}
∫
⋯
∫
1
x
5
times
e
t
2
+
1
/
t
2
d
t
5
=
1
24
[
(
x
−
1
)
4
∫
1
x
e
t
2
+
1
/
t
2
d
t
−
4
(
x
−
1
)
3
∫
1
x
(
t
−
1
)
e
t
2
+
1
/
t
2
d
t
+
6
(
x
−
1
)
2
∫
1
x
(
t
−
1
)
2
e
t
2
+
1
/
t
2
d
t
−
4
(
x
−
1
)
∫
1
x
(
t
−
1
)
3
e
t
2
+
1
/
t
2
d
t
+
∫
1
x
(
t
−
1
)
4
e
t
2
+
1
/
t
2
d
t
]
{\displaystyle {\underset {5~{\text{times}}}{\int \cdots \int _{1}^{x}}}e^{t^{2}+1/t^{2}}~{\text{d}}t^{5}={\frac {1}{24}}\left[(x-1)^{4}\int _{1}^{x}e^{t^{2}+1/t^{2}}~{\text{d}}t-4(x-1)^{3}\int _{1}^{x}(t-1)e^{t^{2}+1/t^{2}}~{\text{d}}t+6(x-1)^{2}\int _{1}^{x}(t-1)^{2}e^{t^{2}+1/t^{2}}~{\text{d}}t-4(x-1)\int _{1}^{x}(t-1)^{3}e^{t^{2}+1/t^{2}}~{\text{d}}t+\int _{1}^{x}(t-1)^{4}e^{t^{2}+1/t^{2}}~{\text{d}}t\right]}
Note that the general chain integral formula cannot be used since (1)
g
(
x
)
=
x
2
+
1
/
x
2
{\displaystyle g(x)=x^{2}+1/x^{2}}
is not strictly monotone, (2) no balanced essential singularity can be isolated from
e
x
2
+
1
/
x
2
{\displaystyle e^{x^{2}+1/x^{2}}}
, not even
e
1
/
x
2
{\displaystyle e^{1/x^{2}}}
, and (3) even the balanced essential singularity
e
1
/
x
2
{\displaystyle e^{1/x^{2}}}
is not elementary-integrable.
The work-energy formula is intrinsic to any twice continuously differentiable increasing function:
y
′
(
x
|
y
=
b
)
2
−
y
′
(
x
|
y
=
a
)
2
2
=
∫
y
=
a
y
=
b
y
″
(
x
|
y
)
d
y
{\displaystyle {\frac {y'(x|_{y=b})^{2}-y'(x|_{y=a})^{2}}{2}}=\int _{y=a}^{y=b}y''(x|_{y})~{\text{d}}y}
This is by Chain Rule , since
d
f
′
(
f
−
1
(
x
)
)
2
2
d
x
=
f
″
(
f
−
1
(
x
)
)
{\displaystyle {\frac {{\text{d}}f'(f^{-1}(x))^{2}}{2~{\text{d}}x}}=f''(f^{-1}(x))}
for any twice continuously differentiable increasing function
f
{\displaystyle f}
.
Using the work-energy formula
KE
=
1
2
m
(
v
f
2
−
v
0
2
)
=
∫
x
1
x
2
F
(
x
)
d
x
{\displaystyle {\text{KE}}={\frac {1}{2}}m(v_{f}^{2}-v_{0}^{2})=\int _{x_{1}}^{x_{2}}F(x)~{\text{d}}x}
with
m
=
1
{\displaystyle m=1}
, initial position
x
1
=
r
{\displaystyle x_{1}=r}
, and initial velocity
v
0
(
x
1
)
=
0
{\displaystyle v_{0}(x_{1})=0}
; we have
(
d
y
d
x
)
2
=
2
r
y
−
1
=
v
2
=
2
∫
r
y
F
(
x
)
d
x
⟹
d
d
y
[
2
r
y
−
1
]
=
d
d
y
[
2
∫
r
y
F
(
x
)
d
x
]
⟺
−
2
r
y
2
=
2
F
(
y
)
⟹
F
(
y
)
=
−
r
y
2
{\displaystyle \left({\frac {{\text{d}}y}{{\text{d}}x}}\right)^{2}={\frac {2r}{y}}-1=v^{2}=2\int _{r}^{y}F(x)~{\text{d}}x\Longrightarrow {\frac {\text{d}}{{\text{d}}y}}\left[{\frac {2r}{y}}-1\right]={\frac {\text{d}}{{\text{d}}y}}\left[2\int _{r}^{y}F(x)~{\text{d}}x\right]\Longleftrightarrow -{\frac {2r}{y^{2}}}=2F(y)\Longrightarrow F(y)=-{\frac {r}{y^{2}}}}
by the Fundamental Theorem of Calculus I . So in fact the cycloid is the solution
(
x
,
y
(
x
)
)
{\displaystyle {\big (}x,y(x){\big )}}
to Newton's Gravitational Law (where time is measured in
x
{\displaystyle x}
) if a particle is bounded in a heavy point object's gravitational field (negative net energy), with zero angular momentum; if we set
r
=
G
M
{\displaystyle r=GM}
,
y
=
R
{\displaystyle y=R}
, to obtain
F
(
R
)
=
−
G
M
m
R
2
⟹
a
(
R
)
=
−
G
M
R
2
{\displaystyle F(R)=-{\frac {GMm}{R^{2}}}\Longrightarrow a(R)=-{\frac {GM}{R^{2}}}}
.
Other zero-angular-momentum solutions for identical mass include
y
(
x
)
=
x
2
3
9
r
2
3
{\displaystyle y(x)=x^{\frac {2}{3}}{\sqrt[{3}]{\frac {9r}{2}}}}
for zero net energy, and
(
x
,
y
(
x
)
)
=
(
r
(
sinh
t
−
t
)
,
r
(
cosh
t
−
1
)
)
{\displaystyle {\big (}x,y(x){\big )}={\big (}r(\sinh t-t),r(\cosh t-1){\big )}}
for positive net energy. For the first case, indeed
(
d
y
d
x
)
2
=
[
d
d
x
(
x
2
3
9
r
2
3
)
]
2
=
[
2
3
9
r
2
x
3
]
2
=
4
9
81
r
2
4
x
2
3
=
2
r
9
3
81
⋅
4
2
3
⋅
r
3
=
2
r
9
r
2
3
x
2
3
=
2
r
y
{\displaystyle \left({\frac {{\text{d}}y}{{\text{d}}x}}\right)^{2}=\left[{\frac {\text{d}}{{\text{d}}x}}\left(x^{\frac {2}{3}}{\sqrt[{3}]{\frac {9r}{2}}}\right)\right]^{2}=\left[{\frac {2}{3}}{\sqrt[{3}]{\frac {9r}{2x}}}\right]^{2}={\frac {4}{9}}{\sqrt[{3}]{\frac {81r^{2}}{4x^{2}}}}={\frac {2r}{\sqrt[{3}]{{\frac {9^{3}}{81}}\cdot {\frac {4}{2^{3}}}\cdot r}}}={\frac {2r}{{\sqrt[{3}]{\frac {9r}{2}}}x^{\frac {2}{3}}}}={\frac {2r}{y}}}
by the Power Rule , where the additional constant of
−
1
{\displaystyle -1}
is specific to the cycloid only (which is the total energy); and for the second case, indeed
(
d
y
d
x
)
2
=
[
d
d
t
(
r
(
cosh
t
−
1
)
)
d
d
t
(
r
(
sinh
t
−
t
)
)
]
2
=
[
r
sinh
t
r
(
cosh
t
−
1
)
]
2
=
[
sinh
t
cosh
t
−
1
]
2
=
sinh
2
t
cosh
2
t
−
2
cosh
t
+
1
=
cosh
2
t
−
1
cosh
2
t
−
2
cosh
t
+
1
{\displaystyle \left({\frac {{\text{d}}y}{{\text{d}}x}}\right)^{2}=\left[{\frac {{\frac {\text{d}}{{\text{d}}t}}{\big (}r(\cosh t-1){\big )}}{{\frac {\text{d}}{{\text{d}}t}}{\big (}r(\sinh t-t){\big )}}}\right]^{2}=\left[{\frac {r\sinh t}{r(\cosh t-1)}}\right]^{2}=\left[{\frac {\sinh t}{\cosh t-1}}\right]^{2}={\frac {\sinh ^{2}t}{\cosh ^{2}t-2\cosh t+1}}={\frac {\cosh ^{2}t-1}{\cosh ^{2}t-2\cosh t+1}}}
=
(
cosh
t
−
1
)
(
cosh
t
+
1
)
(
cosh
t
−
1
)
2
=
cosh
t
+
1
cosh
t
−
1
=
2
cosh
t
−
1
+
1
=
2
r
y
+
1
{\displaystyle ={\frac {(\cosh t-1)(\cosh t+1)}{(\cosh t-1)^{2}}}={\frac {\cosh t+1}{\cosh t-1}}={\frac {2}{\cosh t-1}}+1={\frac {2r}{y}}+1}
by the parametric derivative , and the additional constant of
+
1
{\displaystyle +1}
is specific to this second path only (which is the total energy).