This paper presents a short survey of convergence results and properties of the Lebesgue function λm,n(x) for(0, 1, …, m)Hermite-Fejér interpolation based on the zeros of the nth Chebyshev polynomial of the first kind. The limiting behaviour as n → ∞ of the Lebesgue constant Λm,n = max{λm,n(x): −1 ≤ x ≤ 1} for even m is then studied, and new results are obtained for the asymptotic expansion of Λm,n. Finally, graphical evidence is provided of an interesting and unexpected pattern in the distribution of the local maximum values of λm,n(x) if m ≥ 2 is even.

A Fredholm operator exists which maps the solutions of a system of linear partial differential equations of the form ∂u/∂t = DLu + Au coupled by a matrix A onto those solutions of a similar system coupled by a matrix B which have the same initial values. The kernels of this operator satisfy a hyperbolic system of equations. Since these equations are independent of the linear partial differential operator L, the same operator serves as a mapping for a large class of equations. If B is chosen diagonal, the solutions of a coupled system with matrix A may be obtained from the uncoupled system with matrix B.

Consider the prototype ill-posed problem of a first kind integral equation ℛ with discrete noisy data di, = f(xi) + εi, i = 1, …, n. Let u0 be the true solution and unα a regularised solution with regularisation parameter α. Under certain assumptions, it is known that if α → 0 but not too quickly as n → ∞, then unα converges to u0. We examine the dependence of the optimal sequence of α and resulting optimal convergence rate on the smoothness of f or u0, the kernel K, the order of regularisation m and the error norm used. Some important implications are made, including the fact that m must be sufficiently high relative to the smoothness of u0 in order to ensure optimal convergence. An optimal filtering criterion is used to determine the order where is the maximum smoothness of u0. Two practical methods for estimating the optimal α, the unbiased risk estimate and generalised cross validation, are also discussed.

The purpose of this work is to begin the development of a theory of generating functions that will not only include the generating functions which are partly bilateral and partly unilateral but also provide a set of expansions, by taking successive partial derivatives with respect to one of the variables of the generating relations. Our starting point is a result of Exton [4] on associated Laguerre polynomials whose application gives certain generating functions of the polynomials of Jacobi and Appell, and functions of n variables of Lauricella.

Let L, T, S, and R be closed densely defined linear operators from a Hubert space X into X where L can be factored as L = TS + R. The equation Lu = f is equivalent to the linear system Tv + Ru = f and Su = v. If Lu = f is a two-point boundary value problem, numerical solution of the split system admits cruder approximations than the unsplit equations. This paper develops the theory of such splittings together with the theory of the Methods of Least Squares and of Collocation for the split system. Error estimates in both L2 and L∞ norms are obtained for both methods.

In this paper, we present a method for the construction of a robust observer-based H∞ controller for an uncertain time-delay system. Cases of both single and multiple delays are considered. The parameter uncertainties are time-varying and norm-bounded. Observer and controller are designed to be such that the uncertain system is stable and a disturbance attenuation is guaranteed, regardless of the uncertainties. It has been shown that the above problem can be solved in terms of two linear matrix inequalities (LMIs). Finally, an illustrative example is given to show the effectiveness of the proposed techniques.

We propose the second-order time discretization scheme with the finite-element approximation for the partial integro-differential equations with a weakly singular kernel. The space discretization is based on the finite element method and the time discretization is based on the Crank-Nicolson scheme with a graded mesh. We show the stability of the scheme and obtain the second-order convergence result for the fully discretized scheme.

We consider the boundary-value problems corresponding to the scattering of a time-harmonic acoustic plane wave by a multi-layered obstacle with a sound-soft, hard or penetrable core. Firstly, we construct in closed forms the normalized scattering amplitudes and prove the classical reciprocity and scattering theorems for these problems. These results are then used to study the spectrum of the scattering amplitude operator. The scattering cross-section is expressed in terms of the forward value of the corresponding normalized scattering amplitude. Finally, we develop a more general theory for scattering relations.

In this paper we propose a systematic theoretical procedure for the constructive approximation of non-linear operators and show how this procedure can be applied to the modelling of dynamical systems. We extend previous work to show that the model is stable to small disturbances in the input signal and we pay special attention to the role of real number parameters in the modelling process. The implications of computability are also discussed. A number of specific examples are presented for the particular purpose of illustrating the theoretical procedure.

The maximum principle for optimal control problems of stochastic systems consisting of forward and backward state variables is proved, under the assumption that the diffusion coefficient does not contain the control variable, but the control domain need not be convex.

In previous papers, three terms have been included in Ogden's stress-deformation function for incompressible isotropic elastic materials. The material constants have been calculated by elementary methods and the resulting fits to sets of experimental data have been moderately good.

The purpose of the present paper is to improve upon established correlation between theory and experiment by means of a systematic optimization procedure for calculating material constants. For purposes of illustration the Levenberg-Marquardt non-linear least squares optimization algorithm is adapted to determine the material constants in Ogden's stress-deformation function.

The use of this algorithm for three-term stress-deformation functions improves somewhat on previous results. Calculations are also carried out in respect of a four-term stress-deformation function and further improvement in the fit is achieved over a large range of deformation.

The problem of determining a square integrable function from its modulus and that of its Fourier transform has been considered in an article by Corbett and Hurst ([1]). In this work we point out an error in one of the main results of the cited article concerning the Pauli uniqueness of real states, provide a proof of uniqueness for non-negative states, and present various related examples and discussion of the problem.

We improve some results of [17], which relate to key tools given in [7] for establishing canonical inequalities used in the analysis of sum sets and fractals.

The Shannon system is generalized and the expansion of a function in the generalized Shannon system is considered. No study of a wavelet expansion exists without the assumption of ‘fast’ decay of wavelets. The wavelet ψ which is associated with the generalized Shannon system has a ‘slow’ decay. The expansion of a function in the system is shown to converge at a point which satisfies the Lipschitz condition of order α > 0. On the other hand, there is a continuous function whose wavelet expansion in the generalized Shannon system diverges. An observation of Gibbs' phenomenon is also given.

The storage of bagasse, which is principally cellulose, presents many problems for the sugar industry, one of which is bagasse loss due to spontaneous combustion. This is an expensive problem for the industry as bagasse is used as a fuel by sugar mills, and for cogeneration of electricity. Self-heating occurs in the pile through an oxidation mechanism as well as a moisture dependent reaction. The latter reaction is now known to exhibit a local maximum, similar to the heat release curves found in cool-flame problems. Bagasse typically contains 45–55% by weight of water when milling is completed and the question of how to reduce the moisture content is important for two reasons. Firstly, wet bagasse does not burn nearly as efficiently as dry bagasse, and secondly, self-heating is greatly enhanced in the presence of water, for temperatures less than 60–70°C.

An existing mathematical model is used, but modified to take into account the newly observed peak in the moisture dependent reaction. Most of the previously reported complex bifurcation behaviour possible in this model is not realized when physically realistic parameter values are used. The bifurcation diagram describing the long-time steady-state solution is the familiar S-shaped hysteresis curve. In the presence of the new form of the moisture dependent reaction, an intermediate state can be found which is not a true steady-state of the system as, in reality, the characteristics of the pile slowly change as water is lost. This state corresponds to observations of an elevated temperature (around 60–70°C) which persists for long periods of time. Approximate equations can then be defined which predict this intermediate state, and hence a different hysteresis curve is found. A simple explanation for the process by which water is lost from the pile is obtained from these equations and an analytical expression is given for the exponential decay of water levels in the pile.

The nonlinear convection-diffusion equation has been studied for 40 years in the context of nonhysteretic water movement in unsaturated soil. We establish new similarity solutions for instantaneous sources of finite strength redistributed by nonlinear convection-diffusion obeying the dimensionless equation ∂θ/∂t = ∂ (θm ∂θ/∂z) − θm+1 ∂θ/∂z (m ≥ 0). For m = 0 (Burgers’ equation) solutions involve the error function, and for m = 1 Airy functions. Problems 1, 2, and 3 relate, respectively, to the regions 0 ≤ z ≤ ∞, –∞ ≤ z ≤ ∞, and –∞ ≤ z ≤ 0. Solutions for m = 0 have infinite tails, but for m > 0 and finite t, θ > 0 inside, and θ = 0 outside, a finite interval in z. At the slug boundary, θ(z) is tangential to the z-axis for 0 < m < 1; and it meets the axis obliquely for m = 1 and normally for m > 1. Illustrative results are presented. For Problems 1 and 2 (but not 3) finiteness of source strength sets an upper bound on Θ0, the similarity “concentration” at z = 0. The magnitude of convection relative to diffusion increases with Θ0; and apparently the dynamic equilibrium between the two processes, implied by the similarity solutions, ceases to be possible when Θ0 is large enough.

For the cubic Schrödinger equation iut = uxx + k|u|2u, 0 ≤ x, t < ∞, initial data u(x, 0) = u0(x) ∈ H2[0, ∞), and Robin boundary data ux(0, t) + αu(0, t) = R(t) ∈ C2[0, ∞) (where α is real), we show that the solution u depends continuously on u0 and R.

Critical point behaviour of the diffusion length γ for the solutions of the radiative transfer equation deep in a homogenous medium is studied. The Legendre expansion of the medium's phase function P(cos ψ) is taken to be an infinite series and is characterized by the parameters h0, h1h2,…. A characteristic equation for γ is given in terms of an infinite continued fraction. From this equation it is shown that as any one of the hn, say hp, approaches zero, the others being held constant, γ behaves as , where the critical exponent is found to be vp = ½ for all p = 0, 1, 2,….

A large number of mathematical models are expressed as differential equations. Such models are often derived through a slowly-varying approximation under the assumption that the domain of interest is arbitrarily large; however, typical solutions and the physical problem of interest possess finite domains. The issue is: what are the correct boundary conditions to be used at the edge of the domain for such model equations? Centre manifold theory [24] and its generalisations may be used to derive these sorts of approximations, and higher-order refinements, in an appealing and systematic fashion. Furthermore, the centre manifold approach permits the derivation of appropriate initial conditions and forcing for the models [25, 7]. Here I show how to derive asymptotically-correct boundary conditions for models which are based on the slowly-varying approximation. The dominant terms in the boundary conditions typically agree with those obtained through physical arguments. However, refined models of higher order require subtle corrections to the previously-deduced boundary conditions, and also require the provision of additional boundary conditions to form a complete model.