• The Pure and Applied Mathematics
• /
• v.22 no.2
• /
• pp.159-168
• /
• 2015

Abstract. We propose a fast and robust finite difference method for Merton's jump diffusion model, which is a partial integro-differential equation. To speed up a computational time, we compute a matrix so that we can calculate the non-local integral term fast by a simple matrix-vector operation. Also, we use non-uniform grids to increase efficiency. We present numerical experiments such as evaluation of the option prices and Greeks to demonstrate a performance of the proposed numerical method. The computational results are in good agreements with the exact solutions of the jump-diffusion model.

• Bulletin of the Korean Mathematical Society
• /
• v.35 no.2
• /
• pp.269-278
• /
• 1998

Let ($G,\;{\upsilon}$) be a bounded smooth domain and reflection vector field on $\partial$G, which points uniformly into G. Under the condition that locally for some coordinate system, ${\mid}{\upsilon^i}{\mid}\;i\;=\;1,{\cdot},{\cdot}$,d - 1, where is constant depending on the Lipschitz constant of G, we have tightness for reflected diffusion with jump on G with reflection $\upsilon$ depending only on c. From this, we obtain some properties of L-harmonic function where L is a sum of Laplacian and integro one.

• Journal of applied mathematics & informatics
• /
• v.17 no.1_2_3
• /
• pp.519-536
• /
• 2005

The purpose of this paper is to propose several approximating methods to obtain the American option prices under jump-diffusion processes. The first method is to extend an approximating method to the optimal exercise boundary by a multipiece exponential function suggested by Ju [17]. The second approach is to modify the analytical methods of MacMillan [20] and Zhang [25] in a discrete time space. The third approach is to apply the simulation technique of Ibanez and Zapareto [14] to the problem of American option pricing when the jumps are allowed. Finally, we compare the numerical performance of each suggesting method with those of the previous numerical approaches.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.1
• /
• pp.23-45
• /
• 2015

We consider counterparty risk in CDS rates. To do so, we use a multivariate jump diffusion process for obligors' default intensity, where jumps (i.e. magnitude of contribution of primary events to default intensities) occur simultaneously and their sizes are dependent. For these simultaneous jumps and their sizes, a homogeneous Poisson process. We apply copula-dependent default intensities of multivariate Cox process to derive the joint Laplace transform that provides us with joint survival/default probability and other relevant joint probabilities. For that purpose, the piecewise deterministic Markov process (PDMP) theory developed in [7] and the martingale methodology in [6] are used. We compute survival/default probability using three copulas, which are Farlie-Gumbel-Morgenstern (FGM), Gaussian and Student-t copulas, with exponential marginal distributions. We then apply the results to calculate CDS rates assuming deterministic rate of interest and recovery rate. We also conduct sensitivity analysis for the CDS rates by changing the relevant parameters and provide their figures.

• Communications for Statistical Applications and Methods
• /
• v.10 no.3
• /
• pp.655-663
• /
• 2003

The fundamental tenn structure model is based on the modelling of the short rate. It is well-known that the short rate depends on the interest rate policy of monetary authorities, especially on the official rate. Babbs and Webber(1994) modelled the tenn structure of interest rates using the official rate. They assume that the official rate follows a jump process. This reflects that the official rate infrequently changes. In this paper, we test this official tenn structure model and compare the jump-diffusion model with the pure diffusion model.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.6
• /
• pp.1105-1119
• /
• 2010

We consider jump processes which has only downward jumps with size a fixed fraction of the current process. The jumps of the pro cesses are interpreted as crashes and we assume that the jump intensity is a nondecreasing function of the current process say $\lambda$(X) (X = X(t) process). For the case of $\lambda$(X) = $X^{\alpha}$, $\alpha$ > 0, we show that the process X shold explode in finite time, say $t_e$, conditional on no crash For the case of $\lambda$(X) = (lnX)$^{\alpha}$, we show that $\alpha$ = 1 is the borderline of two different classes of processes. We generalize the model by adding a Brownian noise and examine the blow up properties of the sample paths.

• Proceedings of the Membrane Society of Korea Conference
• /
• 1997.10a
• /
• pp.27-30
• /
• 1997

1. Introduction : Molecular diffusion of small molecules in polymers plays an important role in many areas where polymers are acting as barriers, and in separation processes, such as selective diffusion. Different applications of polymers have different requirements on their transport properties. Therefore, reliable predictions of diffusion coefficients for small molecules in polymeric materials could be a useful tool to design appropriate materials. For many years, the theories based on free-volume concepts have been widely used to correlate and predict diffusion behavior in polymer/solvent systems. In the theory derived by Vrentas and Duda, the empty space between molecules that is available for molecular transport, referred to as hole free-volume, is being redistributed. Molecular transport will occur only when a free-volume of sufficient size appears adjacent to a molecule and the molecule has enough energy to jump into this void. The diffusive jump is considered complete when the void left behind is closed before the molecule returns to its original position. In this paper, the Vrentas-Duda free-volume theory is presented and the methods to estimate free-volume parameters for predicting polymer/ solvent diffusion coefficients are described in detail.

• Journal of the Korean Mathematical Society
• /
• v.43 no.2
• /
• pp.383-398
• /
• 2006

We consider a geometric Levy process for an underlying asset. We prove first that the option price is the unique solution of certain integro-differential equation without assuming differentiability and boundedness of derivatives of the payoff function. Second result is to provide convergence rate for option prices when the small jumps are removed from the Levy process.

• The Korean Journal of Applied Statistics
• /
• v.29 no.1
• /
• pp.99-112
• /
• 2016

Asymmetric jump-diffusion models are effectively used to model the dynamic behavior of asset prices with abrupt asymmetric upward and downward changes. However, the estimation of their extension to the multivariate asymmetric jump-diffusion model has been hampered by the analytically intractable likelihood function. This article confronts the problem using a data augmentation method and proposes a new Bayesian method for a multivariate asymmetric Laplace jump-diffusion model. Unlike the previous models, the proposed model is rich enough to incorporate all possible correlated jumps as well as mention individual and common jumps. The proposed model and methodology are illustrated with a simulation study and applied to daily returns for the KOSPI, S&P500, and Nikkei225 indices data from January 2005 to September 2015.

• The Pure and Applied Mathematics
• /
• v.27 no.4
• /
• pp.231-249
• /
• 2020

We present an efficient and robust finite difference method for a two-asset jump diffusion model, which is a partial integro-differential equation (PIDE). To speed up a computational time, we compute a matrix so that we can calculate the non-local integral term fast by a simple matrix-vector operation. In addition, we use bilinear interpolation to solve integral term of PIDE. We can obtain more stable value by using the payoff-consistent extrapolation. We provide numerical experiments to demonstrate a performance of the proposed numerical method. The numerical results show the robustness and accuracy of the proposed method.