Difference between revisions of "ArcCosc"

From TORI
Jump to: navigation, search
 
m (Text replacement - "\$([^\$]+)\$" to "\\(\1\\)")
 
Line 1: Line 1:
−
[[File:acoscmapT300.png|600px|thumb|[[complex map]] of $u+\mathrm i v=\mathrm{acosc}(x+\mathrm i y)$]]
 +
[[File:acoscmapT300.png|600px|thumb|[[complex map]] of \(u+\mathrm i v=\mathrm{acosc}(x+\mathrm i y)\)]]
 
[[File:acoscplotT.png|600px|thumb| [[explicit plot]] of acosc and its asymptitics]]
  
[[File:acoscplotT.png|600px|thumb| [[explicit plot]] of acosc and its asymptitics]]
 
'''ArcCosc''' denoted aslo as '''acosc''' is [[inverse function]] of [[Cosc]] defined with
  
'''ArcCosc''' denoted aslo as '''acosc''' is [[inverse function]] of [[Cosc]] defined with
−
: $ \displaystyle \mathrm{cosc}(z)=\frac{\cos(z)}{z}$
 +
: \( \displaystyle \mathrm{cosc}(z)=\frac{\cos(z)}{z}\)
   
 
ArcCosc appears in the calculation of mode of wave channeled between absorbing walls. However, it may have also other applications. Tterefore, the robust [[C++]] implementation is supplied in the description of the complex map (click on the map at right).
  
ArcCosc appears in the calculation of mode of wave channeled between absorbing walls. However, it may have also other applications. Tterefore, the robust [[C++]] implementation is supplied in the description of the complex map (click on the map at right).
Line 8: Line 8:
 
==History and ethimology==
  
==History and ethimology==
 
The previous version of this article uses name [[ArcCip]] instead of '''ArcCosc'''. Then it happened that the inverse function of [[sinc]] has similar application (for antisymmetric mode);
  
The previous version of this article uses name [[ArcCip]] instead of '''ArcCosc'''. Then it happened that the inverse function of [[sinc]] has similar application (for antisymmetric mode);
−
: $\displaystyle \mathrm{sinc}(z)=\frac{\sin(z)}{z}$
 +
: \(\displaystyle \mathrm{sinc}(z)=\frac{\sin(z)}{z}\)
 
Function [[sinc]] already has well established name
  
Function [[sinc]] already has well established name
 
<ref name="wolfram">
  
<ref name="wolfram">
Line 22: Line 22:
   
 
The first real extremum of function cosc is denoted with name [[Sazae]]; its value is denoted with name [[Tarao]]:
  
The first real extremum of function cosc is denoted with name [[Sazae]]; its value is denoted with name [[Tarao]]:
−
: $\mathrm{cosc}(\mathrm{Sazae})=\mathrm{Tarao}~$ , $~ ~ \mathrm{cosc}'(\mathrm{Sazae})=0$
 +
: \(\mathrm{cosc}(\mathrm{Sazae})=\mathrm{Tarao}~\) , \(~ ~ \mathrm{cosc}'(\mathrm{Sazae})=0\)
−
: $\mathrm{Sazae} \approx ~ 2.798386045783887$
 +
: \(\mathrm{Sazae} \approx ~ 2.798386045783887\)
−
: $\mathrm{Tarao} \approx\! -0.33650841691839534$
 +
: \(\mathrm{Tarao} \approx\! -0.33650841691839534\)
−
: $\mathrm{acosc}(\mathrm{Tarao})=\mathrm{Sazae}$
 +
: \(\mathrm{acosc}(\mathrm{Tarao})=\mathrm{Sazae}\)
   
−
In the figure with the [[explicit plot]] of acosc, it is compared with its asymptotic denoted with name $\mathrm{Left}_3$.
 +
In the figure with the [[explicit plot]] of acosc, it is compared with its asymptotic denoted with name \(\mathrm{Left}_3\).
−
: $\displaystyle \mathrm{Left}_3(z)=\mathrm{Sazae}
+
: \(\displaystyle \mathrm{Left}_3(z)=\mathrm{Sazae}
 
- \Bigg(
  
- \Bigg(
 
\frac{2(z\!-\!\mathrm{Tarao})}
  
\frac{2(z\!-\!\mathrm{Tarao})}
 
{- \mathrm{Tarao} } \Bigg)^{\!\!1/2}
  
{- \mathrm{Tarao} } \Bigg)^{\!\!1/2}
 
- \frac{2 (z\! -\!\mathrm{Tarao})}
  
- \frac{2 (z\! -\!\mathrm{Tarao})}
−
{3 ~\rm Sazae~Tarao} $
 +
{3 ~\rm Sazae~Tarao} \)
−
The subscript 3 indicates that the only 3 terms taken into account in the expansion. Function $\mathrm{Left}_3(z)$ approximates $\mathrm{acosc}(z)$ in the left hand side of the figure;
 +
The subscript 3 indicates that the only 3 terms taken into account in the expansion. Function \(\mathrm{Left}_3(z)\) approximates \(\mathrm{acosc}(z)\) in the left hand side of the figure;
−
: $\displaystyle \mathrm{acosc}(z)= \mathrm{Left}_3(z) + O(z -\mathrm{Tarao})^{3/2}$
 +
: \(\displaystyle \mathrm{acosc}(z)= \mathrm{Left}_3(z) + O(z -\mathrm{Tarao})^{3/2}\)
   
 
==ArcCosc of real argument==
  
==ArcCosc of real argument==
 
Along the real axis, ArcCosc is defined for argument not smaller than
 
Along the real axis, ArcCosc is defined for argument not smaller than
−
: $\mathrm{Tarao} \approx -0.33650841691839534$
 +
: \(\mathrm{Tarao} \approx -0.33650841691839534\)
   
−
AcrCosc is monotonously decreasing function, it decreases from value Sazae at Tarao till zero at $+\infty$.
 +
AcrCosc is monotonously decreasing function, it decreases from value Sazae at Tarao till zero at \(+\infty\).
   
−
: $\mathrm{acosc}(\mathrm{Tarao}) = \mathrm{Sazae} \approx 2.798386045783887$
 +
: \(\mathrm{acosc}(\mathrm{Tarao}) = \mathrm{Sazae} \approx 2.798386045783887\)
   
 
At Tarao, ArcCosc has singularity (branch point), its derivative becomes infinite.
  
At Tarao, ArcCosc has singularity (branch point), its derivative becomes infinite.
   
 
Specific value ArcCosc has at zero,
  
Specific value ArcCosc has at zero,
−
: $\mathrm{acosc}(0)= \pi/2 \approx 1.5707963267948966$
 +
: \(\mathrm{acosc}(0)= \pi/2 \approx 1.5707963267948966\)
   
 
The expansion of ArcCosc at zero has form
  
The expansion of ArcCosc at zero has form
−
: $\displaystyle \mathrm{acosc}(z)= \frac{\pi}{2} - \frac{\pi}{2} z
 +
: \(\displaystyle \mathrm{acosc}(z)= \frac{\pi}{2} - \frac{\pi}{2} z
 
+ \frac{\pi}{2} z^2
 
+ \frac{\pi}{2} z^2
 
- \frac{\pi ~ (24\!+\!\pi^2)}{48} z^3
  
- \frac{\pi ~ (24\!+\!\pi^2)}{48} z^3
 
+ \frac{\pi ~ (6\!+\!\pi^2)}{12} z^4
  
+ \frac{\pi ~ (6\!+\!\pi^2)}{12} z^4
−
+...$
 +
+...\)
 
[[Mathematica]] built-in operation [[InverseSeries]] allows straightforward calculation of many terms of this
 
[[Mathematica]] built-in operation [[InverseSeries]] allows straightforward calculation of many terms of this
−
expansion. The radius of convergence of this series, id est, $- \mathrm{Tarao}$, is determined by the distance to the closest singularity.
 +
expansion. The radius of convergence of this series, id est, \(- \mathrm{Tarao}\), is determined by the distance to the closest singularity.
   
 
At infinity, ArcCosc decays becomming almost inversly proportional to its arument.
 
At infinity, ArcCosc decays becomming almost inversly proportional to its arument.
Line 65: Line 65:
 
==More expansions of ArcCosc==
  
==More expansions of ArcCosc==
 
In vicinity of Tarao, ArcCosc can be approximated with series
  
In vicinity of Tarao, ArcCosc can be approximated with series
−
: $\!\!\!\!\!\!\!\!\!\! \displaystyle \mathrm{acosc}(\mathrm{Tarao}\!+\!t)=\mathrm{Sazae}
 +
: \(\!\!\!\!\!\!\!\!\!\! \displaystyle \mathrm{acosc}(\mathrm{Tarao}\!+\!t)=\mathrm{Sazae}
 
-\Bigg(\frac{- 2 ~ t}{\mathrm{Tarao}}\Bigg)^{\! 1/2}
  
-\Bigg(\frac{- 2 ~ t}{\mathrm{Tarao}}\Bigg)^{\! 1/2}
 
+ \frac{2 ~ t}{\mathrm{Sazae} ~ \mathrm{Tarao}}
 
+ \frac{2 ~ t}{\mathrm{Sazae} ~ \mathrm{Tarao}}
 
- \frac{3~\mathrm{Sazae}^2-4}{9~ \mathrm{Sazae}^2} \Bigg( \frac{- ~ t}{2~ \mathrm{Tarao}} \Bigg)^{3/2} +
  
- \frac{3~\mathrm{Sazae}^2-4}{9~ \mathrm{Sazae}^2} \Bigg( \frac{- ~ t}{2~ \mathrm{Tarao}} \Bigg)^{3/2} +
−
..$
 +
..\)
−
The truncated series with tree terms gives approximation called $\mathrm{Left}_3$ shown at the figure; it approximates ArcCosc at the left hand side of the figure.
 +
The truncated series with tree terms gives approximation called \(\mathrm{Left}_3\) shown at the figure; it approximates ArcCosc at the left hand side of the figure.
   
 
Similar expansion can be written for the point named i Wakame; [[Wakame]] is real solution of equation
  
Similar expansion can be written for the point named i Wakame; [[Wakame]] is real solution of equation
−
: $\mathrm{cosc}'(\mathrm{i} ~ \mathrm{Wakame})=0$
 +
: \(\mathrm{cosc}'(\mathrm{i} ~ \mathrm{Wakame})=0\)
−
: $\mathrm{Wakame} \approx 1.50887956153832$
 +
: \(\mathrm{Wakame} \approx 1.50887956153832\)
−
: $\mathrm{cosc}(\mathrm{i} ~ \mathrm{Wakame})= - \mathrm{i} ~ \mathrm{Fune} \approx - \mathrm{i} ~ 1.199678640257734$
 +
: \(\mathrm{cosc}(\mathrm{i} ~ \mathrm{Wakame})= - \mathrm{i} ~ \mathrm{Fune} \approx - \mathrm{i} ~ 1.199678640257734\)
 
For calculations and evaluation of the coefficients of this expansion, it is convenient to deal with
  
For calculations and evaluation of the coefficients of this expansion, it is convenient to deal with
 
function [[ArcCohc]], which is inverse function of
 
function [[ArcCohc]], which is inverse function of
−
:$ \displaystyle \mathrm{cohc}(z)=\frac{ \cosh(z) }{z}= - \mathrm{i}~ \mathrm{cosc}(\mathrm{i} z)$
 +
:\( \displaystyle \mathrm{cohc}(z)=\frac{ \cosh(z) }{z}= - \mathrm{i}~ \mathrm{cosc}(\mathrm{i} z)\)
   
 
A dozen terms of these expansions can be calculated with [[Mathematica]].
 
A dozen terms of these expansions can be calculated with [[Mathematica]].
−
The radius $r$ of convergence of these series is determined by the distance to the nearby singularity,
 +
The radius \(r\) of convergence of these series is determined by the distance to the nearby singularity,
−
$r=\sqrt{\mathrm{Tarao}^2+\mathrm{Wakame}^2} \approx 1.5459480734762723$
 +
\(r=\sqrt{\mathrm{Tarao}^2+\mathrm{Wakame}^2} \approx 1.5459480734762723\)
   
 
At Infinity, ArcCosc can be expanded as follows
  
At Infinity, ArcCosc can be expanded as follows
−
: $\displaystyle \mathrm{acosc}(z) = \frac{1}{z}-\frac{1}{2
 +
: \(\displaystyle \mathrm{acosc}(z) = \frac{1}{z}-\frac{1}{2
 
z^3}+\frac{13}{24
  
z^3}+\frac{13}{24
 
z^5}-\frac{541}{720
  
z^5}-\frac{541}{720
 
z^7}+\frac{9509}{8064
  
z^7}+\frac{9509}{8064
−
z^9} + ..$
 +
z^9} + ..\)
−
This expansion can be used for the precise approcimation of ArcCosc while $|z|\gg 1$, $\Re(z)\!>\!0 $.
 +
This expansion can be used for the precise approcimation of ArcCosc while \(|z|\gg 1\), \(\Re(z)\!>\!0 \).
   
 
==Numerical implementation==
  
==Numerical implementation==
Line 99: Line 99:
   
 
The bransh points are Tarao,
 
The bransh points are Tarao,
−
$\mathrm{i}\,\mathrm{Wakame}$
+
\(\mathrm{i}\,\mathrm{Wakame}\)
 
and
  
and
−
$-\mathrm{i}\,\mathrm{Wakame}$, where
+
\(-\mathrm{i}\,\mathrm{Wakame}\), where
   
−
: $ \mathrm{Wakame}=\mathrm{cohc}(\mathrm{Fune})$
 +
: \( \mathrm{Wakame}=\mathrm{cohc}(\mathrm{Fune})\)
 
Name Fune is used for the real solution of equation
  
Name Fune is used for the real solution of equation
−
: $ \mathrm{cohc}'( \mathrm {Fune} )=0$
 +
: \( \mathrm{cohc}'( \mathrm {Fune} )=0\)
 
where
 
where
−
: $ \displaystyle \mathrm{cohc}(z)= \frac{\cosh(z)}{z}$
 +
: \( \displaystyle \mathrm{cohc}(z)= \frac{\cosh(z)}{z}\)
   
 
Constants Wakame and Fine have the following values:
  
Constants Wakame and Fine have the following values:
−
: $ \mathrm{Wakame} \approx 1.50887956153832$
 +
: \( \mathrm{Wakame} \approx 1.50887956153832\)
−
: $ \mathrm{Fune} \approx 1.199678640257734$
 +
: \( \mathrm{Fune} \approx 1.199678640257734\)
 
in such a way,
  
in such a way,
−
: $ \mathrm{Fune} = \mathrm{acohc}(\mathrm{Wakame})$
 +
: \( \mathrm{Fune} = \mathrm{acohc}(\mathrm{Wakame})\)
 
where acohc=[[ArcCohc]] is inverse function of [[cohc]].
  
where acohc=[[ArcCohc]] is inverse function of [[cohc]].
   
−
For the rough approximation of $\mathrm{acosc}(z)$, first, one of the four expansions mentioned above is used:
 +
For the rough approximation of \(\mathrm{acosc}(z)\), first, one of the four expansions mentioned above is used:
−
at Tarao, at $\pm \mathrm i \, \mathrm{Wakame}$ and at $+\infty~$.
 +
at Tarao, at \(\pm \mathrm i \, \mathrm{Wakame}\) and at \(+\infty~\).
−
Then the approximation $s_0$ is adjusted with the [[Newton algorithm]],
 +
Then the approximation \(s_0\) is adjusted with the [[Newton algorithm]],
−
: $\displaystyle s_{n+1}=s_n+ \frac{z-\mathrm{cosc}(s_n)}{ \mathrm{cosc}'(s_n)}$
 +
: \(\displaystyle s_{n+1}=s_n+ \frac{z-\mathrm{cosc}(s_n)}{ \mathrm{cosc}'(s_n)}\)
−
While there exist efficient implementations of $\mathrm{cosc}$ and $\mathrm{cosc}'$, it is faster to make few additional iterations than to increase number of terms in the truncated series of the primary approximation; the goal of the primary approximation is to indicate the correct branch of the inverse function. Then, the Newtom method converges quickly: from 4 to 7 iterations are applied giving at least 14 correct decimal digits in the result.
 +
While there exist efficient implementations of \(\mathrm{cosc}\) and \(\mathrm{cosc}'\), it is faster to make few additional iterations than to increase number of terms in the truncated series of the primary approximation; the goal of the primary approximation is to indicate the correct branch of the inverse function. Then, the Newtom method converges quickly: from 4 to 7 iterations are applied giving at least 14 correct decimal digits in the result.
   
 
==Constants==
  
==Constants==
Line 126: Line 126:
 
Four mathematical constants are defined above. In order to avoid confusions with single–letter notations, for
 
Four mathematical constants are defined above. In order to avoid confusions with single–letter notations, for
 
these constants the following names are suggested: [[Sazae]], [[Tarao]], [[Fune]], [[Wakame]].
  
these constants the following names are suggested: [[Sazae]], [[Tarao]], [[Fune]], [[Wakame]].
−
: $\mathrm{Sazae} \approx ~ 2.798386045783887$
 +
: \(\mathrm{Sazae} \approx ~ 2.798386045783887\)
−
: $\mathrm{Tarao} \approx\! -0.33650841691839534$
 +
: \(\mathrm{Tarao} \approx\! -0.33650841691839534\)
−
: $ \mathrm{Fune} \approx 1.199678640257734$
 +
: \( \mathrm{Fune} \approx 1.199678640257734\)
−
: $ \mathrm{Wakame} \approx 1.50887956153832$
 +
: \( \mathrm{Wakame} \approx 1.50887956153832\)
 
with the following properties:
  
with the following properties:
−
: $ \mathrm{cosc}'(\mathrm{Sazae}) =0$
 +
: \( \mathrm{cosc}'(\mathrm{Sazae}) =0\)
−
: $ \mathrm{cosc}(\mathrm{Sazae}) =\mathrm{Tarao}$
 +
: \( \mathrm{cosc}(\mathrm{Sazae}) =\mathrm{Tarao}\)
−
: $ \mathrm{acosc}(\mathrm{Tarao}) =\mathrm{Sazae}$
 +
: \( \mathrm{acosc}(\mathrm{Tarao}) =\mathrm{Sazae}\)
−
: $ \mathrm{cohc}'(\mathrm{Fune}) =0$
 +
: \( \mathrm{cohc}'(\mathrm{Fune}) =0\)
−
: $ \mathrm{cohc}(\mathrm{Fune}) =\mathrm{Wakame}$
 +
: \( \mathrm{cohc}(\mathrm{Fune}) =\mathrm{Wakame}\)
−
: $ \mathrm{acohc}(\mathrm{Wakame}) =\mathrm{Fune}$
 +
: \( \mathrm{acohc}(\mathrm{Wakame}) =\mathrm{Fune}\)
 
Knowledge of these constants simplifies evaluation of functions [[ArcCosc]] and [[ArcCohc]].
  
Knowledge of these constants simplifies evaluation of functions [[ArcCosc]] and [[ArcCohc]].
   

Latest revision as of 18:27, 30 July 2019

complex map of \(u+\mathrm i v=\mathrm{acosc}(x+\mathrm i y)\)
explicit plot of acosc and its asymptitics

ArcCosc denoted aslo as acosc is inverse function of Cosc defined with

\( \displaystyle \mathrm{cosc}(z)=\frac{\cos(z)}{z}\)

ArcCosc appears in the calculation of mode of wave channeled between absorbing walls. However, it may have also other applications. Tterefore, the robust C++ implementation is supplied in the description of the complex map (click on the map at right).

History and ethimology

The previous version of this article uses name ArcCip instead of ArcCosc. Then it happened that the inverse function of sinc has similar application (for antisymmetric mode);

\(\displaystyle \mathrm{sinc}(z)=\frac{\sin(z)}{z}\)

Function sinc already has well established name [1][2]. For this analogy, the name of function Cosc is used. As usually, the name of the inverse function is created adding prefic Arc or a.

For this reason the function is called ArcCosc.

Notations

In formulas, the old name ArcCip and the short name acosc also may be used.

The first real extremum of function cosc is denoted with name Sazae; its value is denoted with name Tarao:

\(\mathrm{cosc}(\mathrm{Sazae})=\mathrm{Tarao}~\) , \(~ ~ \mathrm{cosc}'(\mathrm{Sazae})=0\)
\(\mathrm{Sazae} \approx ~ 2.798386045783887\)
\(\mathrm{Tarao} \approx\! -0.33650841691839534\)
\(\mathrm{acosc}(\mathrm{Tarao})=\mathrm{Sazae}\)

In the figure with the explicit plot of acosc, it is compared with its asymptotic denoted with name \(\mathrm{Left}_3\).

\(\displaystyle \mathrm{Left}_3(z)=\mathrm{Sazae} - \Bigg( \frac{2(z\!-\!\mathrm{Tarao})} {- \mathrm{Tarao} } \Bigg)^{\!\!1/2} - \frac{2 (z\! -\!\mathrm{Tarao})} {3 ~\rm Sazae~Tarao} \)

The subscript 3 indicates that the only 3 terms taken into account in the expansion. Function \(\mathrm{Left}_3(z)\) approximates \(\mathrm{acosc}(z)\) in the left hand side of the figure;

\(\displaystyle \mathrm{acosc}(z)= \mathrm{Left}_3(z) + O(z -\mathrm{Tarao})^{3/2}\)

ArcCosc of real argument

Along the real axis, ArcCosc is defined for argument not smaller than

\(\mathrm{Tarao} \approx -0.33650841691839534\)

AcrCosc is monotonously decreasing function, it decreases from value Sazae at Tarao till zero at \(+\infty\).

\(\mathrm{acosc}(\mathrm{Tarao}) = \mathrm{Sazae} \approx 2.798386045783887\)

At Tarao, ArcCosc has singularity (branch point), its derivative becomes infinite.

Specific value ArcCosc has at zero,

\(\mathrm{acosc}(0)= \pi/2 \approx 1.5707963267948966\)

The expansion of ArcCosc at zero has form

\(\displaystyle \mathrm{acosc}(z)= \frac{\pi}{2} - \frac{\pi}{2} z + \frac{\pi}{2} z^2 - \frac{\pi ~ (24\!+\!\pi^2)}{48} z^3 + \frac{\pi ~ (6\!+\!\pi^2)}{12} z^4 +...\)

Mathematica built-in operation InverseSeries allows straightforward calculation of many terms of this expansion. The radius of convergence of this series, id est, \(- \mathrm{Tarao}\), is determined by the distance to the closest singularity.

At infinity, ArcCosc decays becomming almost inversly proportional to its arument.

Few asymptotics of ArcCosc are shown in the figure.

More expansions of ArcCosc

In vicinity of Tarao, ArcCosc can be approximated with series

\(\!\!\!\!\!\!\!\!\!\! \displaystyle \mathrm{acosc}(\mathrm{Tarao}\!+\!t)=\mathrm{Sazae} -\Bigg(\frac{- 2 ~ t}{\mathrm{Tarao}}\Bigg)^{\! 1/2} + \frac{2 ~ t}{\mathrm{Sazae} ~ \mathrm{Tarao}} - \frac{3~\mathrm{Sazae}^2-4}{9~ \mathrm{Sazae}^2} \Bigg( \frac{- ~ t}{2~ \mathrm{Tarao}} \Bigg)^{3/2} + ..\)

The truncated series with tree terms gives approximation called \(\mathrm{Left}_3\) shown at the figure; it approximates ArcCosc at the left hand side of the figure.

Similar expansion can be written for the point named i Wakame; Wakame is real solution of equation

\(\mathrm{cosc}'(\mathrm{i} ~ \mathrm{Wakame})=0\)
\(\mathrm{Wakame} \approx 1.50887956153832\)
\(\mathrm{cosc}(\mathrm{i} ~ \mathrm{Wakame})= - \mathrm{i} ~ \mathrm{Fune} \approx - \mathrm{i} ~ 1.199678640257734\)

For calculations and evaluation of the coefficients of this expansion, it is convenient to deal with function ArcCohc, which is inverse function of

\( \displaystyle \mathrm{cohc}(z)=\frac{ \cosh(z) }{z}= - \mathrm{i}~ \mathrm{cosc}(\mathrm{i} z)\)

A dozen terms of these expansions can be calculated with Mathematica. The radius \(r\) of convergence of these series is determined by the distance to the nearby singularity, \(r=\sqrt{\mathrm{Tarao}^2+\mathrm{Wakame}^2} \approx 1.5459480734762723\)

At Infinity, ArcCosc can be expanded as follows

\(\displaystyle \mathrm{acosc}(z) = \frac{1}{z}-\frac{1}{2 z^3}+\frac{13}{24 z^5}-\frac{541}{720 z^7}+\frac{9509}{8064 z^9} + ..\)

This expansion can be used for the precise approcimation of ArcCosc while \(|z|\gg 1\), \(\Re(z)\!>\!0 \).

Numerical implementation

The C++ implementation with complex double precision is available in the descriptions of the figures (Just click at the picture).

The expansions of the function at its branch points is used to get the primary approximation, then the value is adjusted with the Newton algorithm.

The bransh points are Tarao, \(\mathrm{i}\,\mathrm{Wakame}\) and \(-\mathrm{i}\,\mathrm{Wakame}\), where

\( \mathrm{Wakame}=\mathrm{cohc}(\mathrm{Fune})\)

Name Fune is used for the real solution of equation

\( \mathrm{cohc}'( \mathrm {Fune} )=0\)

where

\( \displaystyle \mathrm{cohc}(z)= \frac{\cosh(z)}{z}\)

Constants Wakame and Fine have the following values:

\( \mathrm{Wakame} \approx 1.50887956153832\)
\( \mathrm{Fune} \approx 1.199678640257734\)

in such a way,

\( \mathrm{Fune} = \mathrm{acohc}(\mathrm{Wakame})\)

where acohc=ArcCohc is inverse function of cohc.

For the rough approximation of \(\mathrm{acosc}(z)\), first, one of the four expansions mentioned above is used: at Tarao, at \(\pm \mathrm i \, \mathrm{Wakame}\) and at \(+\infty~\). Then the approximation \(s_0\) is adjusted with the Newton algorithm,

\(\displaystyle s_{n+1}=s_n+ \frac{z-\mathrm{cosc}(s_n)}{ \mathrm{cosc}'(s_n)}\)

While there exist efficient implementations of \(\mathrm{cosc}\) and \(\mathrm{cosc}'\), it is faster to make few additional iterations than to increase number of terms in the truncated series of the primary approximation; the goal of the primary approximation is to indicate the correct branch of the inverse function. Then, the Newtom method converges quickly: from 4 to 7 iterations are applied giving at least 14 correct decimal digits in the result.

Constants

Four mathematical constants are defined above. In order to avoid confusions with single–letter notations, for these constants the following names are suggested: Sazae, Tarao, Fune, Wakame.

\(\mathrm{Sazae} \approx ~ 2.798386045783887\)
\(\mathrm{Tarao} \approx\! -0.33650841691839534\)
\( \mathrm{Fune} \approx 1.199678640257734\)
\( \mathrm{Wakame} \approx 1.50887956153832\)

with the following properties:

\( \mathrm{cosc}'(\mathrm{Sazae}) =0\)
\( \mathrm{cosc}(\mathrm{Sazae}) =\mathrm{Tarao}\)
\( \mathrm{acosc}(\mathrm{Tarao}) =\mathrm{Sazae}\)
\( \mathrm{cohc}'(\mathrm{Fune}) =0\)
\( \mathrm{cohc}(\mathrm{Fune}) =\mathrm{Wakame}\)
\( \mathrm{acohc}(\mathrm{Wakame}) =\mathrm{Fune}\)

Knowledge of these constants simplifies evaluation of functions ArcCosc and ArcCohc.

Conclusion

The knowledhe of the asymptotic behavior and the high efficiency of the evaluation algorithm allows to deal with ArcCosc as it it would be elementary function. It has sense to include such functions as Tetration, ArcTetration, SuperFactorial, AbelFactorial, and ArcCosc to the set of Special functions.

See also

cos, ArcCos, acosq, Sinc, complex map, inverse function

References

  1. ↑ http://mathworld.wolfram.com/SincFunction.html
  2. ↑ http://en.wikipedia.org/wiki/Sinc_function
Retrieved from "https://mizugadro.mydns.jp/t/index.php?title=ArcCosc&oldid=28040"