Cubic Hermite Finite Element Method for Nonlinear Black-Scholes Equation Governing European Options

A numerical algorithm for solving a generalized Black-Scholes partial differential equation, which arises in European option pricing considering transaction costs is developed. The Crank-Nicolson method is used to discretize in the temporal direction and the Hermite cubic interpolation method to discretize in the spatial direction. The efficiency and accuracy of the proposed method are tested numerically, and the results confirm the theoretical behaviour of the solutions, which is also found to be in good agreement with the exact solution.

European Option Pricing Formula in Risk-Aversive Markets

In this study, using the method of discounting the terminal expectation value into its initial value, the pricing formulas for European options are obtained under the assumptions that the financial market is risk-aversive, the risk measure is standard deviation, and the price process of underlying asset follows a geometric Brownian motion. In particular, assuming the option writer does not need the risk compensation in a risk-neutral market, then the obtained results are degenerated into the famous Black–Scholes model (1973); furthermore, the obtained results need much weaker conditions than those of the Black–Scholes model. As a by-product, the obtained results show that the value of European option depends on the drift coefficient μ of its underlying asset, which does not display in the Black–Scholes model only because μ = r in a risk-neutral market according to the no-arbitrage opportunity principle. At last, empirical analyses on Shanghai 50 ETF options and S&P 500 options show that the fitting effect of obtained pricing formulas is superior to that of the Black–Scholes model.

European Option Pricing under Wishart Processes

This study deals with a single risky asset pricing model whose volatility is described by Wishart affine processes. This multifactor model with two dependency matrices describing the correlation between the asset dynamic and Wishart processes makes it more flexible enough to fit the market data for short or long maturities. The aim of the study is to derive and solve the call option pricing problem under the double Wishart stochastic volatility model. The Fourier transform techniques combined with perturbation methods are employed in order to price the European call options. The numerical illustrations on pricing predictions show similar behavior of price movements under the double Wishart model with respect to the market price.

Optimal exercise of American put options near maturity: A new economic perspective

The critical price $$S^{*}\left( t\right) $$ S ∗ t of an American put option is the underlying stock price level that triggers its immediate optimal exercise. We provide a new perspective on the determination of the critical price near the option maturity T when the jump-adjusted dividend yield of the underlying stock is either greater than or weakly smaller than the riskfree rate. Firstly, we prove that $$S^{*}\left( t\right) $$ S ∗ t coincides with the critical price of the covered American put (a portfolio that is long in the put as well as in the stock). Secondly, we show that the stock price that represents the indifference point between exercising the covered put and waiting until T is the European-put critical price, at which the European put is worth its intrinsic value. Finally, we prove that the indifference point’s behavior at T equals $$S^{*}\left( t\right) $$ S ∗ t ’s behavior at T when the stock price is either a geometric Brownian motion or a jump-diffusion. Our results provide a thorough economic analysis of $$S^{*}\left( t\right) $$ S ∗ t and rigorously show the correspondence of an American option problem to an easier European option problem at maturity .

Guaranteed Deterministic Approach to Superhedging: Case of Binary European Option

For the superreplication problem with discrete time, a guaranteed deterministic formulation is considered: the problem is to guarantee coverage of the contingent liability on sold option under all admissible scenarios. These scenarios are defined by means of a priori defined compacts dependent on price prehistory: the price increments at each point in time must lie in the corresponding compacts. In a general case, we consider a market with trading constraints and assume the absence of transaction costs. The formulation of the problem is game theoretic and leads to the Bellman–Isaacs equations. This paper analyses the solution to these equations for a specific pricing problem, i.e., for a binary option of the European type, within a multiplicative market model, with no trading constraints. A number of solution properties and an algorithm for the numerical solution of the Bellman equations are derived. The interest in this problem, from a mathematical prospective, is related to the discontinuity of the option payoff function.

Bond and Option Prices under Skew Vasicek Model with Transaction Cost

This paper studies the European option pricing on the zero-coupon bond in which the Skew Vasicek model uses to predict the interest rate amount. To do this, we apply the skew Brownian motion as the random part of the model and show that results of the model predictions are better than other types of the model. Besides, we obtain an analytical formula for pricing the zero-coupon bond and find the European option price by constructing a portfolio that contains the option and a share of the bond. Since the skew Brownian motion is not a martingale, thus we add transaction costs to the portfolio, where the time between trades follows the exponential distribution. Finally, some numerical results are presented to show the efficiency of the proposed model.

REPLICATION SCHEME FOR THE PRICING OF EUROPEAN OPTIONS

This paper proposes an efficient method for calculating European option prices under local, stochastic, and fractional volatility models. Instead of directly calculating the density function of a target underlying asset, we replicate it from a simpler diffusion process with a known analytical solution for the European option. For this purpose, we derive six functions that characterize the density function of a diffusion process, for both the original and simpler processes and match these functions so that the latter mimics the former. Using the analytical formula, we then approximate the option price of the target asset. By comparison with previous works and numerical experiments, we show that the accuracy of our approximation is high, and the calculation is fast enough for practical purposes; hence, it is suitable for calibration purposes.