A numerical method for pricing European options with proportional transaction costs

In the paper, we propose a numerical technique based on a finite difference scheme in space and an implicit time-stepping scheme for solving the Hamilton–Jacobi–Bellman (HJB) equation arising from the penalty formulation of the valuation of European options with proportional transaction costs. We show that the approximate solution from the numerical scheme converges to the viscosity solution of the HJB equation as the mesh sizes in space and time approach zero. We also propose an iterative scheme for solving the nonlinear algebraic system arising from the discretization and establish a convergence theory for the iterative scheme. Numerical experiments are presented to demonstrate the robustness and accuracy of the method.     

Symmetry of the Barochronous Motions of a Gas

We study the system of nonlinear differential equations which expresses the constancy of the algebraic invariants of the Jacobian matrix for smooth vector fields in three-dimensional space. This system is equivalent to the equations of gas dynamics which describe the barochronous motions of a gas (the pressure and density depend only on the time). We present the results of computation of the admissible local Lie group and construction of the general solution of the system. We mention a few new problems that arise here.     

Two-Grid Algorithm of $H^{1}$-Galerkin Mixed Finite Element Methods for Semilinear Parabolic Integro-Differential Equations

In this paper, we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by $H^{1}$-Galerkin mixed finite element methods. We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization, and backward Euler scheme for temporal discretization. Firstly, a priori error estimates and some superclose properties are derived. Secondly, a two-grid scheme is presented and its convergence is discussed. In the proposed two-grid scheme, the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy. Finally, a numerical experiment is implemented to verify theoretical results of the proposed scheme. The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice $h=H^2$.     

A linearized finite-difference scheme for the numerical solution of the nonlinear cubic Schrödinger equation

A linearized finite-difference scheme is used to transform the initial/boundary-value problem associated with the nonlinear Schrödinger equation into a linear algebraic system. This method is developed by re placing the time and the space partial derivatives by parametric finite-difference re placements and the nonlinear term by an appropriate parametric linearized scheme based on Taylor’s expansion. The resulting finite-difference method is analysed for stability and convergence. The results of a number of numerical experiments for the single-soliton wave are given.     

A numerical scheme for multi-point special boundary-value problems and application to fluid flow through porous layers

In this paper we propose a numerical scheme based on finite differences for the numerical solution of nonlinear multi-point boundary-value problems over adjacent domains. In each subdomain the solution is governed by a different equation. The solutions are required to be smooth across the interface nodes. The approach is based on using finite difference approximation of the derivatives at the interface nodes. Smoothness across the interface nodes is imposed to produce an algebraic system of nonlinear equations. A modified multi-dimensional Newton’s method is proposed for solving the nonlinear system. The accuracy of the proposed scheme is validated by examples whose exact solutions are known. The proposed scheme is applied to solve for the velocity profile of fluid flow through multilayer porous media.     

A Linearized Spectral-Galerkin Method for Three-Dimensional Riesz-Like Space Fractional Nonlinear Coupled Reaction-Diffusion Equations

In this paper, we establish a novel fractional model arising in the chemical reaction and develop an efficient spectral method for the three-dimensional Riesz-like space fractional nonlinear coupled reaction-diffusion equations. Based on the backward difference method for time stepping and the Legendre-Galerkin spectral method for space discretization, we construct a fully discrete numerical scheme which leads to a linear algebraic system. Then a direct method based on the matrix diagonalization approach is proposed to solve the linear algebraic system, where the cost of the algorithm is of a small multiple of $N^4$ ($N$ is the polynomial degree in each spatial coordinate) flops for each time level. In addition, the stability and convergence analysis are rigorously established. We obtain the optimal error estimate in space, and the results also show that the fully discrete scheme is unconditionally stable and convergent of order one in time. Furthermore, numerical experiments are presented to confirm the theoretical claims. As the applications of the proposed method, the fractional Gray-Scott model is solved to capture the pattern formation with an analysis of the properties of the fractional powers.     

A spectral element method using the modal basis and its application in solving second‐order nonlinear partial differential equations

We present a high‐order spectral element method (SEM) using modal (or hierarchical) basis for modeling of some nonlinear second‐order partial differential equations in two‐dimensional spatial space. The discretization is based on the conforming spectral element technique in space and the semi‐implicit or the explicit finite difference formula in time. Unlike the nodal SEM, which is based on the Lagrange polynomials associated with the Gauss–Lobatto–Legendre or Chebyshev quadrature nodes, the Lobatto polynomials are used in this paper as modal basis. Using modal bases due to their orthogonal properties enables us to exactly obtain the elemental matrices provided that the element‐wise mapping has the constant Jacobian. The difficulty of implementation of modal approximations for nonlinear problems is treated in this paper by expanding the nonlinear terms in the weak form of differential equations in terms of the Lobatto polynomials on each element using the fast Fourier transform (FFT). Utilization of the Fourier interpolation on equidistant points in the FFT algorithm and the enough polynomial order of approximation of the nonlinear terms can lead to minimize the aliasing error. Also, this approach leads to finding numerical solution of a nonlinear differential equation through solving a system of linear algebraic equations. Numerical results for some famous nonlinear equations illustrate efficiency, stability and convergence properties of the approximation scheme, which is exponential in space and up to third‐order in time. Copyright © 2014 John Wiley & Sons, Ltd.     

Solving the 2D and 3D nonlinear inverse source problems of elliptic type partial differential equations by a homogenization function method

In this article, the inverse source problems of 2D and 3D elliptic type nonlinear partial differential equations are resolved. For this purpose, a family of single-parameter homogenization functions that automatically meet the given boundary conditions are deduced and employed as the bases to expand the solution. We solve a linear algebraic equations system which satisfies the over-specified Neumann boundary condition to obtain the unspecified coefficients, and then the solution in the entire domain is permitted. Taking the solution into the governing equation, the unknown source function can be determined quickly. The present novel method is verified to be an accurate, effective, and robust scheme which is without solving nonlinear equations and iterations, and the additional data used are quite economical.     

Analyzing the identifiability of linear dynamic models with parameter space separators

We take a new approach to studying the structural global identifiability of linear dynamic models in the state space using the concept of separators of the parameter space. We offer some criteria for the truth of separators which are based on specially constructed matrices, thereby avoiding the laborious analytic solution of a systemof nonlinear algebraic equations. Some examples are given that illustrate applications of the proposed approach.     

On the best parametrization

The numerical solution to a system of nonlinear algebraic or transcendental equations with several parameters is examined in the framework of the parametric continuation method. Necessary and sufficient conditions are proved for choosing the best parameters, which provide the best condition number for the system of linear continuation equations. Such parameters have to be sought in the subspace tangent to the solution space of the system of nonlinear equations. This subspace is obtained if the original system of nonlinear equations is solved at the various parameter values from a given set. The parametric approximation of curves and surfaces is considered.     

A numerical algorithm for simulation of the Q-switched fiber laser using the travelling wave model

We develop a finite-difference scheme for approximation of a system of nonlinear PDEs describing the Q-switching process. We construct it by using staggered grids. The transport equations are approximated along characteristics, and quadratic nonlinear functions are linearized using a special selection of staggered grids. The stability analysis proves that a connection between time and space steps arises only due to approximation requirements in order to follow exactly the directions of characteristics. The convergence analysis of this scheme is done in two steps. First, some estimates of the uniform boundedness of the discrete solution are proved. This part of the analysis is done locally, in some neighborhood of the exact solution. Second, on the basis of the obtained estimates, the main stability inequality is proved. The second-order convergence rate with respect to the space and time coordinates follows from this stability estimate. Using the obtained convergence result, we prove that the local stability analysis in the selected neighborhood of the exact solution is sufficient.     

微分变换法求解二维非线性Volterra积分微分方程

为了解决二维非线性Volterra积分微分方程的求解问题,本文给出微分变换法.利用该方法将方程中的微分部分和积分部分进行变换,这样简化了原方程,进而得到非线性代数方程组,从而将原问题转换为求解非线性代数方程组的解,使得计算更简便.文中最后数值算例说明了该方法的可行性和有效性.     

Rational General Solutions of Trivariate Rational Differential Systems

We generalize the method of Ngô and Winkler (J Symbolic Comput 46:1173–1186, 2011) for finding rational general solutions of a plane rational differential system to the case of a trivariate rational differential system. We give necessary and sufficient conditions for the trivariate rational differential system to have a rational solution based on proper reparametrization of invariant algebraic space curves. In fact, the problem for computing a rational solution of the trivariate rational differential system can be reduced to finding a linear rational solution of an autonomous differential equation. We prove that the linear rational solvability of this autonomous differential equation does not depend on the choice of proper parametrizations of invariant algebraic space curves. In addition, two different rational solutions corresponding to the same invariant algebraic space curve are related by a shifting of the variable. We also present a criterion for a rational solution to be a rational general solution.     

Solution of Volterra operator-integral equations in the nonregular case by the successive approximation method

The branches of solutions of a nonlinear integral equation of Volterra type in a Banach space are constructed by the successive approximation method. We consider the case in which a solution may have an algebraic branching point. We reduce the equation to a system regular in a neighborhood of the branching point. Continuous and generalized solutions are considered. General existence theorems are used to study an initial-boundary value problem with degeneration in the leading part.     

Be careful with the Exp-function method – Additional remarks

It is shown that the solution produced by the Exp-function method may not hold for all initial conditions. Riccati and Maccari nonlinear differential equations are used to illustrate that fact. Conditions of existence for the produced solution in the space of initial conditions and in the space of system’s parameters are derived using the operator method based on the generalized operator of differentiation. The concept of the expansion of an ordinary differential equation is introduced and it is shown that the algebraic–analytical solution of Maccari equation can be produced by solving Riccati equation.     

A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations

In this paper we propose a new modified recursion scheme for the resolution of multi-order and multi-point boundary value problems for nonlinear ordinary and partial differential equations by the Adomian decomposition method (ADM). Our new approach, including Duan’s convergence parameter, provides a significant computational advantage by allowing for the acceleration of convergence and expansion of the interval of convergence during calculations of the solution components for nonlinear boundary value problems, in particular for such cases when one of the boundary points lies outside the interval of convergence of the usual decomposition series. We utilize the boundary conditions to derive an integral equation before establishing the recursion scheme for the solution components. Thus we can derive a modified recursion scheme without any undetermined coefficients when computing successive solution components, whereas several prior recursion schemes have done so. This modification also avoids solving a sequence of nonlinear algebraic equations for the undetermined coefficients fraught with multiple roots, which is required to complete calculation of the solution by several prior modified recursion schemes using the ADM.     

Nonlinear Stability of General Linear Methods

This paper is devoted to a study of stability of general linear methods for the numerical solution of nonlinear stiff initial value problems in a Hilbert space. New stability concepts are introduced. A criterion of weak algebraic stability is established, which is an improvement and extension of the existing criteria of algebraic stability.     

Algebraic approach for the exploration of the onset of chaos in discrete nonlinear dynamical systems

An algebraic approach based on the rank of a sequence is proposed for the exploration of the onset of chaos in discrete nonlinear dynamical systems. The rank of the partial solution is identified and a special technique based on Hankel matrices is used to decompose the solution into algebraic primitives comprising roots of the modified characteristic equation. The distribution of roots describes the dynamical complexity of a solution and is used to explore properties of the nonlinear system and the onset of chaos.     

A method of solving a nonlinear boundary value problem with a parameter for a loaded differential equation

A nonlinear loaded differential equation with a parameter on a finite interval is studied. The interval is partitioned by the load points, at which the values of the solution to the equation are set as additional parameters. A nonlinear boundary value problem for the considered equation is reduced to a nonlinear multipoint boundary value problem for the system of nonlinear ordinary differential equations with parameters. For fixed parameters, we obtain the Cauchy problems for ordinary differential equations on the subintervals. Substituting the values of the solutions to these problems into the boundary condition and continuity conditions at the partition points, we compose a system of nonlinear algebraic equations in parameters. A method of solving the boundary value problem with a parameter is proposed. The method is based on finding the solution to the system of nonlinear algebraic equations composed.