In this paper, a new approach to a characterization of solvability of a nonlinear nonsmooth multiobjective programming problem with inequality constraints is introduced. A family of η-approximated vector optimization problems is constructed by a modification of the objective and the constraint functions in the original nonsmooth multiobjective programming problem. The connection between (weak) efficient points in the original nonsmooth multiobjective programming problem and its equivalent η-approximated vector optimization problems is established under V-invexity. It turns out that, in most cases, solvability of a nonlinear nonsmooth multiobjective programming problem can be characterized by solvability of differentiable vector optimization problems.

A boundary-value problem for cell growth leads to an eigenvalue problem. In this paper some properties of the eigenfunctions are studied. The first eigenfunction is a probability density function and is of importance in the cell growth model. We sharpen an earlier uniqueness result and show that the distribution is unimodal. We then show that the higher eigenfunctions have nested zeros. We show that the eigenfunctions are not mutually orthogonal, but that there are certain orthogonality relations that effectively partition the set of eigenfunctions into two sets.

In this paper we shall describe a numerical method for the solution of curve flow problems in which the normal velocity of the curve depends locally on the position, normal and curvature of the curve. The method involves approximating the curve by a number of line elements (segments) which are only allowed to move in a direction normal to the element. Hence the normal of each line element remains constant throughout the evolution. In regions of high curvature elements naturally tend to accumulate. The method easily deals with the formation of cusps as found in flame propagation problems and is computationally comparable to a naive marker particle method. As a test of the method we present a number of numerical experiments related to mean curvature flow and flows associated with flame propagation and bushfires.

The allometric hypothesis which relates the shape (y) of biological organs to the size of the plant or animal (x), as a function of the relative growth rates, is ubiquitous in biology. This concept has been especially useful in studies of carcass composition of farm animals, and is the basis for the definition of maintenance requirements in animal nutrition.

When the size variable is random the differential equation describing the relative growth rates of organs becomes a stochastic differential equation, with a solution different from that of the deterministic equation normally used to describe allometry. This is important in studies of carcass composition where animals are slaughtered in different sizes and ages, introducing variance between animals into the size variable.

This paper derives an equation that relates values of the shape variable to the expected values of the size variable at any point. This is the most easily interpreted relationship in many applications of the allometric hypothesis such as the study of the development of carcass composition in domestic animals by serial slaughter. The change in the estimates of the coefficients of the allometric equation found through the usual deterministc equation is demonstrated under additive and multiplicative errors. The inclusion of a factor based on the reciprocal of the size variable to the usual log - log regression equation is shown to produce unbiased estimates of the parameters when the errors can be assumed to be multiplicative.

The consequences of stochastic size variables in the study of carcass composition are discussed.

Both Evans et al. and Li et al. have presented preconditioned methods for linear systems to improve the convergence rates of AOR-type iterative schemes. In this paper, we present a new preconditioner. Some comparison theorems on preconditioned iterative methods for solving L-matrix linear systems are presented. Comparison results and a numerical example show that convergence of the preconditioned Gauss-Seidel method is faster than that of the preconditioned AOR iterative method.

This paper deals with a class of network optimization problems in which the flow is a function of time rather than static as in the classical network flow problem, and storage is permitted at the nodes. A solution method involving discretization will be presented as an application of the ASG algorithm. We furnish a proof that the discretized solution converges to the exact continuous solution. We also apply the method to a water distribution network where we minimize the cost of pumping water to meet supply and demand, subject to both linear and nonlinear constraints.

A delayed predator-prey system with Holling type III functional response is investigated. It is proved that the system is uniformly persistent under some appropriate conditions. By means of suitable Lyapunov functionals, sufficient conditions are derived for the local and global asymptotic stability of a positive equilibrium of the system. Numerical simulations are presented to illustrate the feasibility of our main results.

Mathematical modelling has been used to support the response to the COVID-19 pandemic in countries around the world including Australia and New Zealand. Both these countries have followed similar pandemic response strategies, using a combination of strict border measures and community interventions to minimize infection rates until high vaccine coverage was achieved. This required a different set of modelling tools to those used in countries that experienced much higher levels of prevalence throughout the pandemic.

In this article, we provide an overview of some of the mathematical modelling and data analytics work that has helped to inform the policy response to the pandemic in Australia and New Zealand. This is a reflection on our experiences working at the modelling–policy interface and the impact this has had on the pandemic response. We outline the various types of model outputs, from short-term forecasts to longer-term scenario models, that have been used in different contexts. We discuss issues relating to communication between mathematical modellers and stakeholders such as health officials and policymakers. We conclude with some future challenges and opportunities in this area.

Sufficient conditions are derived for the relative controllability of nonlinear neutral Volterra integrodifferential systems with distributed delays in the control variables. The results are a generalization of previous results and are obtained by using Schauder's fixed-point theorem.

In a paper published in 1949, E. R. Love [10] found an integral relation between a separated solution (in oblate spheroidal coordinates) to a particular mixed boundary-value problem and a solution to the same problem using an integral representation. This note examines further representations of the same type which occur in some simple two and three-dimensional potential problems.

In this paper we consider an infinite horizon, continuous time model of economic growth. We prove two theorems; one on the existence of optimal paths of capital accumulation and the other on the dependence of the set of optimal paths on the initial capital stock (sensitivity analysis). In the existence result the underlying technology set is nonconvex and only its “investment’ slices are convex. The proof is direct, without any use of necessary conditions. In the sensitivity analysis, the technology set is convex and so we have that the value function is concave. Then having that, we show that the set of optimal paths is an upper semicaontinuous multifunction of the initial capital stock.

A mathematical model is developed for the process of gas exchange in lung capillaries, taking into account the transport mechanisms of molecular diffusion and the facilitated diffusion of the species due to haemoglobin. We have assumed here equilibrium conditions which enable us to neglect advection effects. The nth order one-step kinetics of oxygen uptake by haemoglobin, proposed by Sharan and Singh [8], have been incorporated. The solution of this coupled nonlinear facilitated diffusion-reaction problem together with the physiologically-relevant boundary conditions is obtained in the closed form.

It is found that about 97.15% of total haemoglobin has combined with oxygen and 2.85% free pigment is left, which is present as carbaminohaemoglobin, met haemoglobin, carboxy haemoglobin etc. It is also shown that the percentage of free haemoglobin at a given PO2 and PCO2 is independent of total haemoglobin content present in the blood.

The well-known Hill's empirical relation is deduced from our solution. The results obtained from our model, based on physical formulation, are in good agreement with the documented data [6] and those computed from the Kelman [3] empirical relation.

The probability that an interval I is free of eigenvalues in a matrix ensemble with unitary symmetry is given by a Fredholm determinant. When the weight function in the matrix ensemble is a classical weight function, and the interval I includes an endpoint of the support, Tracy and Widom have given a formalism which gives coupled differential equations for the required probability and some auxiliary quantities. We summarize and extend earlier work by expressing the probability and some of the auxiliary quantities in terms of Painlevé transcendents.

For Gauss–Turán quadrature formulae with an even weight function on the interval [−1, 1] and functions analytic in regions of the complex plane which contain in their interiors a circle of radius greater than I, the error term is investigated. In some particular cases we prove that the error decreases monotonically to zero. Also, for certain more general cases, we illustrate how to check numerically if this property holds. Some ℓ2-error estimates are considered.

In this paper we study the best constant in the Sobolev trace embedding H1 (Ω) →Lq(∂Ω) in a bounded smooth domain for 1 < q < 2+ = 2(N - 1)/(N - 2), that is, critical or subcritical q. First, we consider a domain with periodically distributed holes inside which we impose that the involved functions vanish. There exists a critical size of the holes for which the limit problem has an extra term. For sizes larger than critical the best trace constant diverges to infinity and for sizes smaller than critical it converges to the best constant in the domain without holes. Also, we study the problem with the holes located on the boundary of the domain. In this case another critical exists and its extra term appears on the boundary.

Solutions of a homogeneous (r + 1)-term linear difference equation are given in two different forms. One involves the elements of a certain matrix, while the other is in terms of certain lower Hessenberg determinants. The results generalize some earlier results of Brown [1] for the solution of a 3-term linear difference equation.

We consider the effect the competing mechanisms of buoyancy-driven acceleration (arising from heating a surface) and streamline curvature (due to curvature of a surface) have on the stability of boundary-layer flows. We confine our attention to vortex type instabilities (commonly referred to as Görtler vortices) which have been identified as one of the dominant mechanisms of instability in both centrifugally and buoyancy driven boundary layers. The particular model we consider consists of the boundary-layer flow over a heated (or cooled) curved rigid body. In the absence of buoyancy forcing the flow is centrifugally unstable to counter-rotating vortices aligned with the direction of the flow when the curvature is concave (in the fluid domain) and stable otherwise. Heating the rigid plate to a level sufficiently above the fluid's ambient (free-stream) temperature can also serve to render the flow unstable. We determine the level of heating required to render an otherwise centrifugally stable flow unstable and likewise, the level of body cooling that is required to render a centrifugally unstableflow stable.

A numerical comparison of the Monte Carlo (MC) simulation and the finite-difference method for pricing European options under a regime-switching framework is presented in this paper. We consider pricing options on stocks having two to four volatility regimes. Numerical results show that the MC simulation outperforms the Crank–Nicolson (CN) finite-difference method in both the low-frequency case and the high-frequency case. Even though both methods have linear growth, as the number of regimes increases, the computational time of CN grows much faster than that of MC. In addition, for the two-state case, we propose a much faster simulation algorithm whose computational time is almost independent of the switching frequency. We also investigate the performances of two variance-reduction techniques: antithetic variates and control variates, to further improve the efficiency of the simulation.

A mathematical model is presented in which the long Jump is treated as the motion of a projectile under gravity with slight drag. The first two terms of a perturbation solution are obtained and are shown to be more accurate than earlier approximate analytical solutions. Results from the perturbation analysis are just as accurate as results from various numerical schemes, and require far less computer time.

The model is modified to include the observation that a long-jumper's centre of mass is forward of his feet at take-off and behind his feet on landing.

The modified model is used to determine the take-off angle for the current world long jump record, resulting in several interesting observations for athletic coaches.

An improvement is made to two results of Brown and Shepp which are useful in calculations with fractal sets.