Some continuous programming problems in numerical analysis

Many problems in the design and implementation of computational schemes may be studied using the theory and methods of mathematical programming. One seeks to minimize bounds for the errors in the calculated results obtained from a given set of input data, exploiting analytical relations. We describe optimal quadrature rules and give an application to the evaluation of the sums of power series, belonging to an important class. We present results which are based on the theory of linear and semi-infinite programming. We also study the associated complexity issues and obtain simple qualitative results for the computational work required.     

A dynamical strategy for approximation methodsUne stratégie dynamique pour les méthodes d'approximation

The numerical result provided by an approximation method is affected by a global error, which consists of both a truncation error and a round-off error. Let us consider the converging sequence generated by successively dividing by two the step size used. If computations are performed until, in the convergence zone, the difference between two successive approximations is only due to round-off errors, then the global error on the result obtained is minimal. Furthermore its significant bits which are not affected by round-off errors are in common with the exact result, up to one. To cite this article: F. Jézéquel, C. R. Mecanique 334 (2006).     

Quadrature formulas for the generalized Riemann-Stieltjes integral

In this paper we propose a technique of approximation for the generalized Riemann-Stieltjes integral and we found an analogue for Newton-Cotes formulas in the case n = 2 and n = 3. *Beneficiary of a Socrates fellowship at the Department of Mathematics, University of Study of Cagliari, Via Ospedale, n. 72, Cagliari, 09124, Italy, in the period February – July 2002.     

On Domains in Which Harmonic Functions Satisfy Generalized Mean Value Properties

Potential Analysis - Assume that a bounded domain Ω??N (N ≥ 2) has the property that there exists a signed measure µ with compact support in Ω such that, for every...     

Revised Iterative Solution for Groundstate of Schrodinger Equation

A revised iterative method based on Green function defined by quadratures along a single trajectory is proposed to solve the low-lying quantum wave function for Schrodinger equation. Specially a new expression of the perturbed energy is obtained, which is much simpler than the traditional one. The method is applied to solve the unharmonic oscillator potential. The revised iteration procedure gives exactly the same result as those based on the single trajectory quadrature method. A comparison of the revised iteration method to the old one is made using the example of Stark effect. The obtained results are consistent to each other after making power expansion.     

The dynamical behavior of the discontinuous Galerkin method and related difference schemes

We study the dynamical behavior of the discontinuous Galerkin finite element method for initial value problems in ordinary differential equations. We make two different assumptions which guarantee that the continuous problem defines a dissipative dynamical system. We show that, under certain conditions, the discontinuous Galerkin approximation also defines a dissipative dynamical system and we study the approximation properties of the associated discrete dynamical system. We also study the behavior of difference schemes obtained by applying a quadrature formula to the integrals defining the discontinuous Galerkin approximation and construct two kinds of discrete finite element approximations that share the dissipativity properties of the original method.

     

变截面门式刚架结构分析的微分求积单元法

本文采用一种精确、简便的数值计算方法——微分求积单元法(DQEM)对变截面门式刚架结构进行了力学分析。首先建立了一般荷载作用下变截面构件的平衡微分方程，并采用微分求积法进行离散，进而得出了较为精确的分析变截面构件的单元力学模型。该模型的刚度方程不仅反映了单元的刚度性质，而且反映了单元的实际荷载作用，可较为精确地分析变截面门式刚架结构在分布载荷作用下的受力性能。通过与有限元法计算结果的比较，表明了微分求积单元法在变截面刚架的力学分析中的正确性和优越性。微分求积单元法可用于任意形状的刚架结构的静力分析。     

等截面梁有限变形的传递函数增量算法

本文介绍了一种计算等截面梁有限变形的新方法－传递函数增量算法，它是一种半解析数值计算方法。此算法充分利用增量失空法Ｇａｕｓｓ求积公式计算非线性有限变形的特点，并将这些特点与传递函数方法，有效地结合起来，既避免了数值方法计算量大的困难，又使得求解高阶非线性微分方程的解析解成为可能，算例分析表明，这是一种易编程，计算量小，收敛快，求解精度高的行之有效的计算方法。     

A Cross-Entropy Approach to the Estimation of Generalized Linear Multilevel Models

In this article, we use the cross-entropy method for noisy optimization for fitting generalized linear multilevel models through maximum likelihood. We propose specifications of the instrumental distributions for positive and bounded parameters that improve the computational performance. We also introduce a new stopping criterion, which has the advantage of being problem-independent. In a second step we find, by means of extensive Monte Carlo experiments, the most suitable values of the input parameters of the algorithm. Finally, we compare the method to the benchmark estimation technique based on numerical integration. The cross-entropy approach turns out to be preferable from both the statistical and the computational point of view. In the last part of the article, the method is used to model the probability of firm exits in the healthcare industry in Italy. Supplemental materials are available online.     

Generalized Bessel functions: Theory and their applications

This paper presents 2 new classes of the Bessel functions on a compact domain [0,T] as generalized‐tempered Bessel functions of the first‐ and second‐kind which are denoted by GTBFs‐1 and GTBFs‐2. Two special cases corresponding to the GTBFs‐1 and GTBFs‐2 are considered. We first prove that these functions are as the solutions of 2 linear differential operators and then show that these operators are self‐adjoint on suitable domains. Some interesting properties of these sets of functions such as orthogonality, completeness, fractional derivatives and integrals, recursive relations, asymptotic formulas, and so on are proved in detail. Finally, these functions are performed to approximate some functions and also to solve 3 practical differential equations of fractionalorders.