A new formula with derivatives for numerical integration was presented. Based on this formula and the Richardson extrapolafion process, a numerical integration method was established. It can converge faster than the R...A new formula with derivatives for numerical integration was presented. Based on this formula and the Richardson extrapolafion process, a numerical integration method was established. It can converge faster than the Romberg's. With the same accuracy, the computation of the new numerical integration with derivatives is only half of that of Romberg's numerical integration.展开更多
The classical Archimedean approximation of π uses the semiperimeter or area of regular polygons inscribed in or circumscribed about a unit circle in R~2 and it is well-known that by using linear combinations of these...The classical Archimedean approximation of π uses the semiperimeter or area of regular polygons inscribed in or circumscribed about a unit circle in R~2 and it is well-known that by using linear combinations of these basic estimates,modern extrapolation techniques can greatly speed up the approximation process.Similarly,when n vertices are randomly selected on the circle,the semiperimeter and area of the corresponding random inscribed and circumscribing polygons are known to converge to π almost surely as n→∞,and by further applying extrapolation processes,faster convergence rates can also be achieved through similar linear combinations of the semiperimeter and area of these random polygons.In this paper,we further develop nonlinear extrapolation methods for approximating π through certain nonlinear functions of the semiperimeter and area of such polygons.We focus on two types of extrapolation estimates of the forms χ_n=S_n~αA_n~β and Y_n(p)=(αS_n~p+βA_n~p)~(1/p) where α+β=1,p≠0,and Sn and An respectively represents the semiperimeter and area of a random n-gon inscribed in the unit circle in R~2,and Xn may be viewed as the limit of Y_n(p) when p→0.By deriving probabilistic asymptotic expansions with carefully controlled error estimates for Xn and Y_n(p),we show that the choice α=4/3,β=-1/3 minimizes the approximation error in both cases,and their distributions are also asymptotically normal.展开更多
文摘A new formula with derivatives for numerical integration was presented. Based on this formula and the Richardson extrapolafion process, a numerical integration method was established. It can converge faster than the Romberg's. With the same accuracy, the computation of the new numerical integration with derivatives is only half of that of Romberg's numerical integration.
基金supported in part by the National Natural Science Foundation of China (No.12131003)。
文摘The classical Archimedean approximation of π uses the semiperimeter or area of regular polygons inscribed in or circumscribed about a unit circle in R~2 and it is well-known that by using linear combinations of these basic estimates,modern extrapolation techniques can greatly speed up the approximation process.Similarly,when n vertices are randomly selected on the circle,the semiperimeter and area of the corresponding random inscribed and circumscribing polygons are known to converge to π almost surely as n→∞,and by further applying extrapolation processes,faster convergence rates can also be achieved through similar linear combinations of the semiperimeter and area of these random polygons.In this paper,we further develop nonlinear extrapolation methods for approximating π through certain nonlinear functions of the semiperimeter and area of such polygons.We focus on two types of extrapolation estimates of the forms χ_n=S_n~αA_n~β and Y_n(p)=(αS_n~p+βA_n~p)~(1/p) where α+β=1,p≠0,and Sn and An respectively represents the semiperimeter and area of a random n-gon inscribed in the unit circle in R~2,and Xn may be viewed as the limit of Y_n(p) when p→0.By deriving probabilistic asymptotic expansions with carefully controlled error estimates for Xn and Y_n(p),we show that the choice α=4/3,β=-1/3 minimizes the approximation error in both cases,and their distributions are also asymptotically normal.