The function ln ⡠x can always be computed from its ascending power series after preliminary scaling. Suppose first 1 / 10 ⤠x ⤠10 . Then we take square roots repeatedly until | y | is sufficiently small, where
After computing ln â¡ ( 1 + y ) from (4.6.1)
For other values of x set x = 10m ⢠ξ , where 1 / 10 ⤠ξ ⤠10 and m â ⤠. Then
ExponentialsLet x have any real value. First, rescale via
4.45.4 m = â x ln â¡ 10 + 1 2 â , y = x â m ⢠ln â¡ 10 .Then
and since | y | ⤠1 2 ⢠ln ⡠10 = 1.15 ⢠⦠, ey can be computed straightforwardly from (4.2.19).
Trigonometric FunctionsLet x have any real value. We first compute ξ = x / Ï , followed by
Then
4.45.7 sin â¡ x = ( â 1 ) m ⢠sin ⡠θ , cos â¡ x = ( â 1 ) m ⢠cos ⡠θ ,and since | θ | ⤠1 2 â¢ Ï = 1.57 ⢠⦠, sin ⡠θ and cos ⡠θ can be computed straightforwardly from (4.19.1) and (4.19.2).
The other trigonometric functions can be found from the definitions (4.14.4)â(4.14.7).
Inverse Trigonometric FunctionsThe function arctan â¡ x can always be computed from its ascending power series after preliminary transformations to reduce the size of x . From (4.24.15) with u = v = ( ( 1 + x2 ) 1 / 2 â 1 ) / x , we have
Beginning with x 0 = x , generate the sequence
4.45.9 x n = x n â 1 1 + ( 1 + x n â 1 2 ) 1 / 2 , n = 1 , 2 , 3 , ⦠,until x n is sufficiently small. We then compute arctan â¡ x n from (4.24.3), followed by
Another method, when x is large, is to sum
compare (4.24.4).
As an example, take x = 9.47376 . Then
4.45.12 x 1 = 0.90000 ⢠⦠, x 2 = 0.38373 ⢠⦠, x 3 = 0.18528 ⢠⦠, x 4 = 0.09185 ⢠⦠.For the remaining inverse trigonometric functions, we may use the identities provided by the fourth row of Table 4.16.3. For example, arcsin â¡ x = arctan â¡ ( x ⢠( 1 â x2 ) â 1 / 2 ) .
Hyperbolic and Inverse Hyperbolic FunctionsThe hyperbolic functions can be computed directly from the definitions (4.28.1)â(4.28.7). The inverses arcsinh , arccosh , and arctanh can be computed from the logarithmic forms given in §4.37(iv), with real arguments. For arccsch , arcsech , and arccoth we have (4.37.7)â(4.37.9).
Other MethodsSee Luther (1995), Ziv (1991), Cody and Waite (1980), Rosenberg and McNamee (1976), Carlson (1972a). For interval-arithmetic algorithms, see Markov (1981). For Shift-and-Add and CORDIC algorithms, see Muller (1997), Merrheim (1994), Schelin (1983). For multiprecision methods, see Smith (1989), Brent (1976).
§4.45(ii) Complex VariablesFor ln ⡠z and ez
4.45.15 ln â¡ z = ln â¡ | z | + i ⢠ph â¡ z , â Ï â¤ ph â¡ z â¤ Ï , 4.45.16 ez = e â â¡ z ⢠( cos â¡ ( â â¡ z ) + i ⢠sin â¡ ( â â¡ z ) ) .See §1.9(i) for the precise relationship of ph â¡ z to the arctangent function.
The trigonometric functions may be computed from the definitions (4.14.1)â(4.14.7), and their inverses from the logarithmic forms in §4.23(iv), followed by (4.23.7)â(4.23.9). Similarly for the hyperbolic and inverse hyperbolic functions; compare (4.28.1)â(4.28.7), §4.37(iv), and (4.37.7)â(4.37.9).
For other methods see Miel (1981).
§4.45(iii) Lambert W -FunctionFor x â [ â 1 / e , â ) the principal branch W 0 â¡ ( x ) can be computed by solving the defining equation W ⢠eW = x numerically, for example, by Newtonâs rule (§3.8(ii)). Initial approximations are obtainable, for example, from the power series (4.13.6) (with t ⥠0 ) when x is close to â 1 / e , from the asymptotic expansion (4.13.10) when x is large, and by numerical integration of the differential equation (4.13.4) (§3.7) for other values of x .
Similarly for W â 1 â¡ ( x ) in the interval [ â 1 / e , 0 ) (with t ⤠0 in (4.13.6)).
See also Barry et al. (1995b) and Chapeau-Blondeau and Monir (2002).
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.5