A modified Verlet method which involves a kind of mid-point rule is constructed and applied to the one-dimensional motion of elastic balls of finite size, falling under constant gravity in space and then under the chemical potential in the interface region of phase separation within a two-liquid film. When applied to the simulation of two balls falling under constant gravity in space, the new method is found to be computationally superior to the usual Verlet method and to Runge–Kutta methods, as it allows a larger time step for comparable accuracy. The main purpose of this paper is to develop an efficient numerical method to simulate balls in the interface region of phase separation within the two-liquid film, where the ball motion is coupled with two-phase flow. The two-phase flow in the film is described via shallow water equations, using an invariant finite difference scheme that accurately resolves the interface region. A larger time step in computing the ball motion, more comparable with the time step in computing the two-phase flow, is a significant advantage. The computational efficiency of the new method in the coupled problem is demonstrated for the case of four elastic balls in the two-liquid film.

This paper concerns the Monte Carlo method in pricing American-style options under the general class of exponential Lévy models. Traditionally, one must store all the intermediate asset prices so that they can be used for the backward pricing in the least squares algorithm. Therefore the storage requirement grows like , where m is the number of time steps and n is the number of simulated paths. In this paper, we propose a simulation method where the storage requirement is only . The total computational cost is less than twice that of the traditional method. For machines with limited memory, one can now enlarge m and n to improve the accuracy in pricing the options. In numerical experiments, we illustrate the efficiency and accuracy of our method by pricing American options where the log-prices of the underlying assets follow typical Lévy processes such as Brownian motion, lognormal jump-diffusion process, and variance gamma process.

We present a technique to find threshold values that allows the user to separate signal from noise in fluorescence grey-level images. It can be classified as a purely comparative method based upon the amount of “Mutual Information” between two or more florescence images, and we apply it to stacks of such images produced using the newly-developed MELK technology. Our results are compared to results obtained by another research group using a quite different (completely independent and more technology-based) approach; and also to results obtained using Otsu's Thresholding Method, yet another completely distinct approach invented to separate foreground and background in a grey-level image, based on minimising “intra-class variance” [9,10]. The remarkably good agreement found suggests that our proposed comparative information based method not only accounts for the biological mechanisms governing cellular protein networks very well, but also (and probably much more importantly) shows that cells actually organise the spatial structure of their protein networks in a highly non-random fashion as might be expected – and thereby try to optimise their “mutual information content”, and thus most probably their efficiency.

In this paper, we prove existence and uniqueness results for the ground states of the coupled Gross-Pitaevskii equations for describing two-component Bose-Einstein condensates with an internal atomic Josephson junction, and obtain the limiting behavior of the ground states with large parameters. Efficient and accurate numerical methods based on continuous normalized gradient flow and gradient flow with discrete normalization are presented, for computing the ground states numerically. A modified backward Euler finite difference scheme is proposed to discretize the gradient flows. Numerical results are reported, to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomena of the ground sates in the problem.

We consider a nonsymmetric Toeplitz system which arises in the discretization of a partial integro-differential equation in option pricing problems. The preconditioned conjugate gradient method with a tri-diagonal preconditioner is used to solve this system. Theoretical analysis shows that under certain conditions the tri-diagonal preconditioner leads to a superlinear convergence rate. Numerical results exemplify our theoretical analysis.

The paper studies the production inventory problem of minimizing the expected discounted present value of production cost control in manufacturing systems with degenerate stochastic demand. We have developed the optimal inventory production control problem by deriving the dynamics of the inventory-demand ratio that evolves according to a stochastic neoclassical differential equation through Ito's Lemma. We have also established the Riccati based solution of the reduced (one- dimensional) HJB equation corresponding to production inventory control problem through the technique of dynamic programming principle. Finally, the optimal control is shown to exist from the optimality conditions in the HJB equation.