PROPERTIES OF ARCTAN - UFL MAEuhk/PROPERTIES-OF-ARCTAN.pdf · PROPERTIES OF ARCTAN(Z) We know from...
-
Upload
phungtuong -
Category
Documents
-
view
225 -
download
8
Transcript of PROPERTIES OF ARCTAN - UFL MAEuhk/PROPERTIES-OF-ARCTAN.pdf · PROPERTIES OF ARCTAN(Z) We know from...
PROPERTIES OF ARCTAN(Z) We know from elementary calculus that the function z=tan(θ) has an inverse θ=arctan(z). In differentiating z once we have-
∫+
=+=z
zdzzequivalentitsorddz
02
2
1)arctan(])tan(1[ θθ
On setting the upper limit to 1/N with N<1 we find the infinite series expansion for arctan given by-
∑+
−=∫ ∑ −=∞
= ++
=
∞
=
+
0 1212/1
0 0
1
)12(1)1()1()/1arctan(
n nnnN
z n
n
NndzzN
or the equivalent-
dtt
tNn
NN
t
mm
n n
n∫
++∑
+
−=
=
−
= +
/1
02
21
0 12 1)12()1()/1arctan(
This series will converge quite rapidly when N>>1. Thus-
+−+−= ....
)239(71
)239(51
)239(311
2391)239/1arctan( 642
However for N=1, the series just equals that of Gregory which is known to be notoriously slowly convergent-
4
.....91
71
51
311)1arctan( π
=−+−+−=
If one takes the first hundred terms(m=100) in the Gregory series, the integral remainder will still be-
percentsomeordtt
t
t3/1...0024999.0
1
1
02
200=∫
+=
In general the larger N becomes the more rapidly the infinite series for arctan(z) will converge. Thus the series for (π/8) =arctan{ 1/[1+sqrt(2)]} reads -
−
++
+−
+= ...
)21(1
)21(11
)21(1
8 42π
which converges somewhat faster than the Gregory series. Lets examine some of the other analytical characteristics of arctan(z). Its plot for z real looks like this-
We see that arctan(z) varies linearly with z for small z starting with value zero and becomes non-linear in its variation with increasing z, eventually approaching Pi/2 as Pi/2-1/z as z approaches infinity. The function has odd symmetry since arctan(-z)=-arctan(z). Its derivative is just 1/(1+z^2) and hence represents a special case of the Witch of Agnesi ( this curve was studied by the Italian mathematician Maria Agnesi 1718-1799 and received its name due to a mistranslation of the Italian word versiero for curve by an English translator who mixed it up with the Italian word for witch). Using the multiple angle formula for tangent , one also has-
)](tan[arctan)](tan[arctan1)](tan[arctan)](tan[arctan)]arctan()(tan[arctanyxyxyx
−+
=+
or the equivalent -
)arctan()arctan(])1()(arctan[ yx
xyyx
+=−+
On setting x=z and y=∞ we find –
)1arctan()arctan(2 z
z +=π
so that, for example, arctan(2)=π/2-arctan(0.5)=π/2-0.46364..= 1.1071…If x=1 and y=-1/3 one obtains the well known identity-
)31arctan()
21arctan(
4+=
π
and x=1/7, y=-1/8 produces-
)571arctan()
81arctan()
71arctan( +=
Consider next the complex number z=x+iy. Writing this out in polar form yields-
)]arctan(exp[22xyiyxiyx +=+
so that-
+
+−=
22ln)/arctan(
yx
iyxixy
This result relates the arctan to the logarithm function so that-
42
1ln πii=
+
Looking at the near linear relation between arctan(z) and z for z<<1 suggests that arctan(1/N)=m*arctan(1/(m*N) +small correction of the order 1/N^3 for large N. This is indeed the case. By looking at the imaginary part of-
[ ] )ln()ln()()(ln 221121 21 iNpiNpiNiN pp +++=++
one finds-
))34(
1arctan()21arctan(2)1arctan( 2 +
−=NNNN
,
)11827
8arctan()31arctan(3)1arctan( 24 −+
−=NN
NNN
and-
))15160256(
180arctan()41arctan(4)1arctan( 24
2
−−
−−=
NNNN
NN
We next solve an integral in terms of arctan to get-
Therefore one finds-
)117arctan(
72)]
73arctan()
75[arctan(
72
43
1
02 =−=∫
++=t ttdt
It is also possible to manipulate the original integral form for artctan(z) into a variety of different range integrals. Consider the substitutions t=u/N and Nt=tanh(v). These produce the integrals-
−
+
−=∫
−++
=∫++ 22
2
22
2 4
2arctan4
2
)4
4()2
(
1
bac
bat
bacabac
abt
dtacbtat
dt
∫
+−
+=∫
+=∫
+=
∞
=== 02
22
1
022
/1
02
)]1
1()[cosh(111)1arctan(
vu
N
tN
v
dvNN
uNduNdt
tN
Expanding the term in the denominator of the last integral leads to an alternate series for arctan(1/N). In compact form, it reads-
∑++
+=
∞
=0 2
2
2 )1()!12(!4
1)1arctan(
n n
n
Nnn
NN
N
and produces the identity-
...!92!4
!72!3
!52!2
!32!11
2
42322212+++++=
π
Also using the variable substitution u=w/sqrt(w^2+1) yields the symmetric form-
∫
+−+++
=+∞=
−∞=
w
wN
ww
dwNNN
]1
1)1[()1()1(2)/1arctan(
2222
so that-
∫ ∫+
−−
=∫++
=∞
=
∞
=
∞+
∞− 0 02221)(cosh()
21)(cosh(
21)5.0( v v v
dv
v
dv
ww
dwπ
It seems that this last integral in w can form the starting point for an AGM approach for finding precise values of π. It can also be expanded as the series-
...97531
!47531
!3531!2
31!1
1!0
2+
⋅⋅⋅⋅+
⋅⋅⋅+
⋅⋅+
⋅+=
π
which shows an interesting pattern but is unfortunately only slowly convergent. A much more rapidly convergent series is found for larger N. Indeed, we have in general that-
∑+⋅⋅⋅+
=∑ ∫++
=
++−
∫++
=
∞
= +
∞
=
∞
=++
∞
0 120 02/3212
220
2/322
)]12...(531[)1(!2
)1()1(
)1)(1(11
1)1(1
)/1arctan(
n n
n
n w nn nNNn
wdw
NN
wNwdw
NNN
This yields at N=239 the result-
065...54527410483821495928074723864541840760020.)1239()!22(
2)!1)(!(239)239/1arctan(0 12
12
=
=∑++
+=
∞
= +
+
n n
n
nnn
which is accurate to 43 places when adding up just the first nine terms in the infinite series. By telescoping the arctan(1/N) series terms by two, one finds the even faster convergent form-
++
+++∑
+Γ+
+Γ=
∞
= + )31616()58()34(
)212()1(
)21(!4
)/1arctan( 2
2
0 222 nnnNn
nN
nnNN
n n
n
Also using our earlier discussed four term arctan formula for π ( see-NUMERICAL EVALUATION OF PI BY A FOUR TERM ARCTAN FORMULA) we have that-
795028...462643383235897932383.14159265
)71825(51456
)57122(13384
)3250(9120
)1445(3648
)!2()!1(!4
41
1
=
+++∑
−=
∞
= nnnnn
n
nnnπ
with each additional term taken in this series improving the accuracy of π by about 3 places. Note the summation procedure requires no taking of roots and simply involves summation, multiplication, and division of integers. Arctan(z) also relates to the hypergeometric series. Matching term by term of the infinite series for F(a,b,c,x) with the the first infinite series expansion for arctan given earlier, one has-
∫+−
=−==
1
0 22
)1)(1(2);2/3;1,2/1()arctan(
t tzt
dtzzzFz
Also it follows that the second order differential equation-
0)53()1(2 2
2=−−+− w
dzdwz
dzwdzz
has-
)arctan(1 zz
w =
as a solution. Finally, one integration of arctan yields-
)1ln(21)arctan()arctan( 2zzzz +∫ −=
which is easily verified by differentiating both sides. March 2009