The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable ...The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable and fast convergent computation of the bivariate normal integral with simple recursion. When tested computation on a random sample of one million parameter triplets with uniformly distributed components and using double precision arithmetic, the maximum absolute error was 3.45 × 10<sup>-</sup><sup>16</sup>. In additional testing, focusing on cases with correlation coefficients close to one in absolute value, when the computation may be very sensitive to small rounding errors, the accuracy was retained. In rare potentially critical cases, a simple adjustment to the computation procedure was performed—one potentially critical computation was replaced with two equivalent non-critical ones. All new series are suitable for vector and high-precision computation, assuming they are supplemented with appropriate efficient and accurate computation of the arctangent and standard normal cumulative distribution functions. They are implemented by the R package Phi2rho, available on CRAN. Its functions allow vector arguments and are ready to work with the Rmpfr package, which enables the use of arbitrary precision instead of double precision numbers. A special test with up to 1024-bit precision computation is also presented.展开更多
For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime...For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.展开更多
In this paper, the authors show some monotonicity and concavity of the classical psi function, by which several known results are improved and some new asymptotically sharp estimates are obtained for this function. In...In this paper, the authors show some monotonicity and concavity of the classical psi function, by which several known results are improved and some new asymptotically sharp estimates are obtained for this function. In addition, applying the new results to the psi function, the authors improve the well-known lower and upper bounds for the approximate evaluation of Euler＇s constant γ.展开更多
In this paper,we prove the Srivastava-Pint'er's addition theorems(see Applied Mathematic Lett.17(2004),375-380) by applying the another methods.We also provide some analoges of these addition theorems and dedu...In this paper,we prove the Srivastava-Pint'er's addition theorems(see Applied Mathematic Lett.17(2004),375-380) by applying the another methods.We also provide some analoges of these addition theorems and deduce the corresponding special cases.展开更多
As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the fi...As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the first time we propose the definitions of the fractional sum and fractional difference on non-uniform lattices by two different ways.The analogue of Euler’s Beta formula,Cauchy’Beta formula on non-uniform lattices are established,and some fundamental theorems of fractional calculas,the solution of the generalized Abel equation on non-uniform lattices are obtained etc.展开更多
Planck’s radiation law provides an equation for the intensity of the electromagnetic radiation from a physical body as a function of frequency and temperature. The frequency that corresponds to the maximum intensity ...Planck’s radiation law provides an equation for the intensity of the electromagnetic radiation from a physical body as a function of frequency and temperature. The frequency that corresponds to the maximum intensity is a function of temperature. At a specific temperature, for the frequencies correspond to much less than the maximum intensity, an equation was derived in the form of the Lambert <em>W</em> function. Numerical calculations validate the equation. A new form of solution for the Euler’s transcendental equation was derived in the form of the Lambert <em>W</em> function with logarithmic argument. Numerical solutions to the Euler’s equation were determined iteratively and iterative convergences were investigated. Numerical coincidences with physical constants were explored.展开更多
文摘The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable and fast convergent computation of the bivariate normal integral with simple recursion. When tested computation on a random sample of one million parameter triplets with uniformly distributed components and using double precision arithmetic, the maximum absolute error was 3.45 × 10<sup>-</sup><sup>16</sup>. In additional testing, focusing on cases with correlation coefficients close to one in absolute value, when the computation may be very sensitive to small rounding errors, the accuracy was retained. In rare potentially critical cases, a simple adjustment to the computation procedure was performed—one potentially critical computation was replaced with two equivalent non-critical ones. All new series are suitable for vector and high-precision computation, assuming they are supplemented with appropriate efficient and accurate computation of the arctangent and standard normal cumulative distribution functions. They are implemented by the R package Phi2rho, available on CRAN. Its functions allow vector arguments and are ready to work with the Rmpfr package, which enables the use of arbitrary precision instead of double precision numbers. A special test with up to 1024-bit precision computation is also presented.
基金Foundation item： Supported by the National Natural Science Foundation of China（10671056）
文摘For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.
基金Supported by the National Natural Science Foundation of China(11171307)
文摘In this paper, the authors show some monotonicity and concavity of the classical psi function, by which several known results are improved and some new asymptotically sharp estimates are obtained for this function. In addition, applying the new results to the psi function, the authors improve the well-known lower and upper bounds for the approximate evaluation of Euler＇s constant γ.
基金Supported by the PCSIRT of Education of China(IRT0621)Supported by the Innovation Program of Shanghai Municipal Education Committee of China(08ZZ24)Supported by the Henan Innovation Project for University Prominent Research Talents of China(2007KYCX0021)
文摘In this paper,we prove the Srivastava-Pint'er's addition theorems(see Applied Mathematic Lett.17(2004),375-380) by applying the another methods.We also provide some analoges of these addition theorems and deduce the corresponding special cases.
基金Supported by the National Natural Science Foundation Fujian province of China(2016J01032).
文摘As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the first time we propose the definitions of the fractional sum and fractional difference on non-uniform lattices by two different ways.The analogue of Euler’s Beta formula,Cauchy’Beta formula on non-uniform lattices are established,and some fundamental theorems of fractional calculas,the solution of the generalized Abel equation on non-uniform lattices are obtained etc.
文摘Planck’s radiation law provides an equation for the intensity of the electromagnetic radiation from a physical body as a function of frequency and temperature. The frequency that corresponds to the maximum intensity is a function of temperature. At a specific temperature, for the frequencies correspond to much less than the maximum intensity, an equation was derived in the form of the Lambert <em>W</em> function. Numerical calculations validate the equation. A new form of solution for the Euler’s transcendental equation was derived in the form of the Lambert <em>W</em> function with logarithmic argument. Numerical solutions to the Euler’s equation were determined iteratively and iterative convergences were investigated. Numerical coincidences with physical constants were explored.
