A Method for Solving Fredholm Integral Equations of the First Kind Based on Chebyshev Wavelets

In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.     

Necessary Maximum Principle of Stochastic Optimal Control with Delay and Jump Diffusion

In this paper, we have studied the necessary maximum principle of stochastic optimal control problem with delay and jump diffusion.     

A Remark on Pal Type Interpolation on Non-Uniformly Distributed Nodes on the Unit Circle

In this paper we study the problem of explicit representation and convergence of Pal type （0;1） interpolation and its converse, with some additional conditions, on the non-uniformly distributed nodes on the unit circle obtaIned by projecting the interlaced zeros of Pn （x） and Pn′ （x） on the unit circle. The motivation to this problem can be traced to the recent studies on the regularity of Birkhoff interpolation and Pal type interpolations on non-uniformly distributed zeros on the unit circle.     

A MIXED FINITE ELEMENT METHOD ON A STAGGERED MESH FOR NAVIER-STOKES EQUATIONS

In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.     

A Complement to the Valiron-Titchmarsh Theorem for Subharmonic Functions

The Valiron-Titchmarsh theorem on asymptotic behavior of entire functions with negative zeros has been recently generalized onto subharmonic functions with the Riesz measure on a half-line in Rn, n ≥ 3. Here we extend the Drasin complement to the Valiron-Titchmarsh theorem and show that if u is a subharmonic function of this class and of order 0 ρ 1, then the existence of the limit limr →∞ logu(r)/N(r),where N(r) is the integrated counting function of the masses of u, implies the regular asymptotic behavior for both u and its associated measure.     

BIOLUMINESCENCE TOMOGRAPHY： BIOMEDICAL BACKGROUND,MATHEMATICAL THEORY,AND NUMERICAL APPROXIMATION

Over the last couple of years molecular imaging has been rapidly developed to study physiological and pathological processes in vivo at the cellular and molecular levels. Among molecular imaging modalities, optical imaging stands out for its unique advantages, especially performance and cost-effectiveness. Bioluminescence tomography （BLT） is an emerging optical imaging mode with promising biomedical advantages. In this survey paper, we explain the biomedical significance of BLT, summarize theoretical results on the analysis and numerical solution of a diffusion based BLT model, and comment on a few extensions for the study of BLT.     

A PRIORI ERROR ESTIMATES OF A FINITE ELEMENT METHOD FOR DISTRIBUTED FLUX RECONSTRUCTION＊

This paper is concerned with a priori error estimates of a finite element method for numerical reconstruction of some unknown distributed flux in an inverse heat conduction problem. More precisely, some unknown distributed Neumann data are to be recovered on the interior inaccessible boundary using Dirichlet measurement data on the outer ac- cessible boundary. The main contribution in this work is to establish the some a priori error estimates in terms of the mesh size in the domain and on the accessible/inaccessible boundaries, respectively, for both the temperature u and the adjoint state p under the lowest regularity assumption. It is revealed that the lower bounds of the convergence rates depend on the geometry of the domain. These a priori error estimates are of immense interest by themselves and pave the way for proving the convergence analysis of adaptive techniques applied to a general classes of inverse heat conduction problems. Numerical experiments are presented to verify our theoretical prediction.     

LOCALLY STABILIZED FINITE ELEMENT METHOD FOR STOKES PROBLEM WITH NONLINEAR SLIP BOUNDARY CONDITIONS

Based on the low-order conforming finite element subspace （Vh, Mh） such as the P1-P0 triangle element or the Q1-P0 quadrilateral element, the locally stabilized finite element method for the Stokes problem with nonlinear slip boundary conditions is investigated in this paper. For this class of nonlinear slip boundary conditions including the subdifferential property, the weak variational formulation associated with the Stokes problem is an variational inequality. Since （Vh, Mh） does not satisfy the discrete inf-sup conditions, a macroelement condition is introduced for constructing the locally stabilized formulation such that the stability of （Vh, Mh） is established. Under these conditions, we obtain the H1 and L2 error estimates for the numerical solutions.     

A quadrature tau method for fractional differential equations with variable coefficients

In this article, we develop a direct solution technique for solving multi-order fractional differential equations (FDEs) with variable coefficients using a quadrature shifted Legendre tau (Q-SLT) method. The spatial approximation is based on shifted Legendre polynomials. A new formula expressing explicitly any fractional-order derivatives of shifted Legendre polynomials of any degree in terms of shifted Legendre polynomials themselves is proved. Extension of the tau method for FDEs with variable coefficients is treated using the shifted Legendre–Gauss–Lobatto quadrature. Numerical results are given to confirm the reliability of the proposed method for some FDEs with variable coefficients.     

FINITE ELEMENT APPROXIMATIONS FOR SCHROEDINGER EQUATIONS WITH APPLICATIONS TO ELECTRONIC STRUCTURE COMPUTATIONS

In this paper, both the standard finite element discretization and a two-scale finite element discretization for SchrSdinger equations are studied. The numerical analysis is based on the regularity that is also obtained in this paper for the Schroedinger equations. Very satisfying applications to electronic structure computations are provided, too.     

CONVERGENCE RATES FOR DIFFERENCE SCHEMES FOR POLYHEDRAL NONLINEAR PARABOLIC EQUATIONS

We build finite difference schemes for a class of fully nonlinear parabolic equations. The schemes are polyhedral and grid aligned. While this is a restrictive class of schemes, a wide class of equations are well approximated by equations from this class. For regular （C2,α） solutions of uniformly parabolic equations, we also establish of convergence rate of O（α）. A case study along with supporting numerical results is included.     

ERROR ESTIMATES FOR THE TIME DISCRETIZATION FOR NONLINEAR MAXWELL＇S EQUATIONS

This paper is devoted to the study of a nonlinear evolution eddy current model of the type δtB（H） ＋ △↓× （ △↓ × H） =0 subject to homogeneous Dirichlet boundary conditions H × v = 0 and a given initial datum. Here, the magnetic properties of a soft ferromagnet are linked by a nonlinear material law described by B（H）. We apply the backward Euler method for the time discretization and we derive the error estimates in suitable function spaces. The results depend on the nonlinearity of B（H）.     

THE FINITE DIFFERENCE METHOD FOR DISSIPATIVE KLEIN-GORDON-SCHRODINGER EQUATIONS IN THREE SPACE DIMENSIONS

A fully discrete finite difference scheme for dissipative Klein-Gordon-SchrSdinger equations in three space dimensions is analyzed. On the basis of a series of the time-uniform priori estimates of the difference solutions and discrete version of Sobolev embedding the- orems, the stability of the difference scheme and the error bounds of optimal order for the difference solutions are obtained in H2 × H2 ×H1 over a finite time interval. Moreover, the existence of a maximal attractor is proved for a discrete dynamical system associated with the fully discrete finite difference scheme.     

插值型求积公式的代数精度

借助于勒让德多项式的零点性质,证明了N阶插值型求积公式的代数精度可取N到2 N+1之间的任意整数值,计算得到了两点插值型求积公式的代数精度与求积节点位置的关系.简化了[1]中关于3次代数精度的条件的讨论.     

WEAK APPROXIMATION OF OBLIQUELY REFLECTED DIFFUSIONS IN TIME-DEPENDENT DOMAINS

In an earlier paper, we proved the existence of solutions to the Skorohod problem with oblique reflection in time-dependent domains and, subsequently, applied this result to the problem of constructing solutions, in time-dependent domains, to stochastic differential equations with oblique reflection. In this paper we use these results to construct weak approximations of solutions to stochastic differential equations with oblique reflection, in time-dependent domains in Rd, by means of a projected Euler scheme. We prove that the constructed method has, as is the case for normal reflection and time-independent domains, an order of convergence equal to 1/2 and we evaluate the method empirically by means of two numerical examples. Furthermore, using a well-known extension of the Feynman-Kac formula, to stochastic differential equations with reflection, our method gives, in addition, a Monte Carlo method for solving second order parabolic partial differential equations with Robin boundary conditions in time-dependent domains.     

A New Approach for Solving a Class of Delay Fractional Partial Differential Equations

In this article, a new numerical approach has been proposed for solving a class of delay time-fractional partial differential equations. The approximate solutions of these equations are considered as linear combinations of Müntz–Legendre polynomials with unknown coefficients. Operational matrix of fractional differentiation is provided to accelerate computations of the proposed method. Using Padé approximation and two-sided Laplace transformations, the mentioned delay fractional partial differential equations will be transformed to a sequence of fractional partial differential equations without delay. The localization process is based on the space-time collocation in some appropriate points to reduce the fractional partial differential equations into the associated system of algebraic equations which can be solved by some robust iterative solvers. Some numerical examples are also given to confirm the accuracy of the presented numerical scheme. Our results approved decisive preference of the Müntz–Legendre polynomials with respect to the Legendre polynomials.     

Markov调制的随机微分延迟方程的函数稳定性

随机微分延迟方程的指数稳定性被人们广泛研究,但讨论带Markov调制的随机微分延迟方程的函数稳定性的不多.本文主要研究了两种类型的函数稳定性.我们采用了一例特定的Lyapunov函数,来研究带Markov调制的随机微分延迟方程的p阶矩ψα-函数稳定性,并对其几乎必然ψβ/p-函数稳定性也进行了探讨.     

Solutions to Boundary Value Problem of Nonlinear Impulsive Differential Equation of Fractional Order

This paper is devoted to study the existence and uniqueness of solutions to a boundary value problem of nonlinear fractional differential equation with impulsive effects.The arguments are based upon Schauder and Banach fixed-point theorems. We improve and generalize the results presented in[B.Ahmad,S.Sivasundaram, Existence results for nonlinear impulsive hybrid boundary value problems involving fractional differential equations,Nonlinear Analysis：Hybrid Systems,3（2009）,251- 258].     

Error bounds for quadrature methods involving lower order derivatives

Quadrature methods for approximating the definite integral of a function f(t) over an interval [a,b] are in common use. Examples of such methods are the Newton–Cotes formulas (midpoint, trapezoidal and Simpson methods etc.) and the Gauss–Legendre quadrature rules, to name two types of quadrature. Error bounds for these approximations involve higher order derivatives. For the Simpson method, in particular, the error bound involves a fourth-order derivative. Discounting the fact that calculating a fourth-order derivative requires a lot of differentiation, the main concern is that an error bound for the Simpson method, for example, is only relevant for a function that is four times differentiable, a rather stringent condition. This paper caters for functions for which derivatives exist only of order lower than normally required. A number of quadrature methods are considered and error bounds derived involving only lower order derivatives that can be used depending on the smoothness of the function.     

Estimates of the error in Gauss–Legendre quadrature for double integrals

Error estimates are a very important aspect of numerical integration. It is desirable to know what level of truncation error might be expected for a given number of integration points. Here, we determine estimates for the truncation error when Gauss–Legendre quadrature is applied to the numerical evaluation of two dimensional integrals which arise in the boundary element method. Two examples are considered; one where the integrand contains poles, when its definition is extended into the complex plane, and another which contains branch points. In both cases we obtain error estimates which agree with the actual error to at least one significant digit.