In the Laurent expansion ζ(s, a) = 1 s− 1 + ∞ ∑ k=0 (−1)γk(a) k! (s− 1) , 0 < a ≤ 1, of the Riemann-Hurwitz zeta function, the coefficients γk(a) are known as Stieltjes, or generalized Euler, constants. [When a = 1, ζ(s, 1) = ζ(s) (the Riemann zeta function), and γk(1) = γk.] We present a new approach to high-precision approximation of γk(a). Plots of our results reveal much structure in the gr...

We show how ideas originating in the theory of dynamical systems inspire a new approach to numerical integration of functions. Any Lebesgue integral can be approximated by a sequence of integrals with respect to equidistributions, i.e. evenly weighted discrete probability measures concentrated on an equidistributed set. We prove that, in the case where the integrand is real analytic, suitable l...

In this work, a new family of time marching procedures based on Green’s function matrices is presented. The formulation is based on the development of new recurrence relationships, which employ time integral terms to treat initial condition values. These integral terms are numerically evaluated taking into account Newton-Cotes formulas. The Green’s matrices of the model are also numerically com...

The main objective of the study is to compare the Newton-Cotes methods such as the Trapezium rule, Simpson 1/3 rule and Simpson 3/8 rule to estimate the area under the Lorenz curve and Gini coefficient of income using polynomial function with degree 5. Comparing the Gini coefficients of income computed from the Polynomial function with degree 5 for the Trapezium, Simpson 1/3 and Simpson 3/8 met...

2. Write a Matlab script of no more than 10 lines to calculate the Newton-Cotes quadrature weights wj , j = 1, 2, . . . , n, n + 1; go to the end of this handout for more details about wj . For uniformity, we require that the script starts with the four lines n=2; h=1/n; nodes=(0:h:1)’; f=eye(n+1); where the j-th column of the matrix f are needed for the Lagrange interpolation in order to find ...

The main purpose of this paper is to examine the effectiveness of the Quarter-Sweep Gauss-Seidel (QSGS) method in solving the dense linear systems generated from the discretization of the linear Fredholm integral equations of the second kind. In addition, the applications of the various orders of closed Newton-Cotes quadrature discretization schemes will be investigated in order to form linear ...

and hosti 013.02.0 Abstract A numerical method based on an NM-set of general, hybrid of block-pulse function and Taylor series (HBT), is proposed to approximate the solution of nonlinear Volterra–Fredholm integral equations. The properties of HBT are first presented. Also, the operational matrix of integration together with Newton-Cotes nodes are utilized to reduce the computation of nonlinear ...

In many applications of ageand size-structured population models, there is an interest is obtaining good approximations of total population numbers rather than of their densities. Therefore, it is reasonable in such cases to solve numerically not the PDE model equations themselves, but rather their integral equivalents. For this purpose quadrature formulae are used in place of the integrals. Si...