In this paper we present a finite volume method for a two-dimensional Black-Scholes equation with stochastic volatility governing European option pricing. In this work, we first formulate the Black-Scholes equation with a tensor (or matrix) diffusion coefficient into a conversative form. We then present a finite volume method for the resulting equation, based on a fitting technique proposed for a one-dimensional Black-Scholes equation. We show that the method is monotone by proving that the system matrix of the discretized equation is an M-matrix. Numerical experiments, performed to demonstrate the usefulness of the method, will be presented.

There are some nonlinear models for pricing financial derivatives which can improve the linear Black-Scholes model introduced by Black, Scholes and Merton. In these models volatility is not constant anymore, but depends on some extra variables. It can be, for example, transaction costs, a risk from a portfolio, preferences of a large trader, etc. In this thesis we focus on these models. In the first chapter we introduce some important theory of financial derivatives. The second chapter is devoted to the volatility models. We derive three models concerning transaction costs (RAPM, Leland's  and Barles-Soner's model) and Frey's model which assumes a large (dominant) trader on the market. In the third and in the forth chapter we derive portfolio and make numerical experiments with a free boundary. We use the first order additive and the second order Strang splitting methods. We also use approximations of Barles-Soner's model using the identity function and introduce an approximation with the logarithm function of Barles-Soner's model. These models we finally compare with models where the volatility includes constant transaction costs.

Radial basis function (RBF) approximation, is a new extremely powerful tool that is promising for high-dimensional problems, such as those arising from pricing of basket options using the Black-Scholes partial differential equation. The main problem for RBF methods have been ill-conditioning as the RBF shape parameter becomes small, corresponding to flat RBFs. This thesis employs a recently developed method called the RBF-QR method to reduce computational cost by improving the conditioning, thereby allowing for the use of a wider range of shape parameter values. Numerical experiments for the one-dimensional case are presented  and a MATLAB implementation is provided. In our thesis, the RBF-QR method performs better  than the RBF-Direct method for small shape parameters. Using Chebyshev points, instead of a standard uniform distribution, can increase the accuracy through clustering of the nodes towards the boundary. The least squares formulation for RBF methods is preferable to the collocation approach because it can result in smaller errors  for the same number of basis functions.

In this dissertation we first formulate the Black-Scholes equation with a tensor (or matrix) diffusion coefficient into a conservative form and present a convergence analysis for the two-dimensional Black-Scholes equation arising in the Hull-White model for pricing European options with stochastic volatility. We formulate a non-conforming Petrov-Galerkin finite element method with each basis function of the trial space being determined by a set of two-point boundary value problems defined on element edges. We show that the bilinear form of the finite element method is coercive and continuous and establish an upper bound of order O(h) on the discretization error of method, where h denotes the mesh parameter of the discretization. We then present a finite volume method for the resulting equation, based on a fitting technique proposed for a one-dimensional Black-Scholes equation. We show that the method is monotone by proving that the system matrix of the discretized equation is an M-matrix. Numerical experiments, performed to demonstrate the usefulness of the method, will be presentd.

<p>In this work we consider the nonlinear case of Black-Scholes equation and apply it to American options. Also, method of Liao and Khaliq of high order was applied to nonlinear Black-Scholes equation in case of American options. Here, we use this method oh fourth order in time and space to raise American option price accuracy.</p>

In this work we consider the nonlinear case of Black-Scholes equation and apply it to American options. Also, method of Liao and Khaliq of high order was applied to nonlinear Black-Scholes equation in case of American options. Here, we use this method oh fourth order in time and space to raise American option price accuracy.

Impulsos são perturbações abruptas que ocorrem em curto espaço de tempo e podem ser consideradas instantâneas. E os mercados financeiros estão sujeitos a choques bruscos como mudanças de governos, quebra de empresas, entre outros. Assim, é natural considerarmos a ação de tais eventos na precificação de ativos financeiros. Nosso objetivo neste trabalho é obtermos uma formulação para a equação diferencial parcial de Black-Scholes com ação impulsiva de modo que os impulsos representem estes choques. Utilizaremos a teoria de integração não-absoluta em espaço de funções para obtenção desta formulação / Impulses describe the evolution of systems where the continuous development of a process is interrupted by abrupt changes of state. Financial markets are subject to extreme events or shocks as government changes, companies colapse, etc. Thus it seems natural to consider the action of these events in the valuation of derivative securities. The aim of this work is to obtain a formulation for the Black-Scholes equation with impulse action where the impulses can represent these shocks. We use the non-absolute integration theory in functional spaces to obtain such formulation

Options are a special type of derivative securities because their values are derived from the value of some underlying security. Most options can be grouped into either of the two categories: European options which can be exercised only on the expiration date, and American options which can be exercised on or before the expiration date. American options are much harder to deal with than European ones. The reason being the optimal exercise policy of these options which led to free boundary problems. Ever since the seminal work of Black and Scholes [J. Pol. Econ. 81(3) (1973), 637-659], the differential equation approach in pricing options has attracted many researchers. Recently, numerical singular perturbation techniques have been used extensively for solving many differential equation models of sciences and engineering. In this thesis, we explore some of those methods which are based on spline approximations to solve the option pricing problems. We show a systematic construction and analysis of these methods to solve some European option problems and then extend the approach to solve problems of pricing American options as well as some exotic options. Proposed methods are analyzed for stability and convergence. Thorough numerical results are presented and compared with those seen in the literature.

Options are a special type of derivative securities because their values are derived from the value of some underlying security. Most options can be grouped into either of the two categories: European options which can be exercised only on the expiration date, and American options which can be exercised on or before the expiration date. American options are much harder to deal with than European ones. The reason being the optimal exercise policy of these options which led to free boundary problems. Ever since the seminal work of Black and Scholes [J. Pol. Econ. 81(3) (1973), 637-659], the differential equation approach in pricing options has attracted many researchers. Recently, numerical singular perturbation techniques have been used extensively for solving many differential equation models of sciences and engineering. In this thesis, we explore some of those methods which are based on spline approximations to solve the option pricing problems. We show a systematic construction and analysis of these methods to solve some European option problems and then extend the approach to solve problems of pricing American options as well as some exotic options. Proposed methods are analyzed for stability and convergence. Thorough numerical results are presented and compared with those seen in the literature.

Impulsos são perturbações abruptas que ocorrem em curto espaço de tempo e podem ser consideradas instantâneas. E os mercados financeiros estão sujeitos a choques bruscos como mudanças de governos, quebra de empresas, entre outros. Assim, é natural considerarmos a ação de tais eventos na precificação de ativos financeiros. Nosso objetivo neste trabalho é obtermos uma formulação para a equação diferencial parcial de Black-Scholes com ação impulsiva de modo que os impulsos representem estes choques. Utilizaremos a teoria de integração não-absoluta em espaço de funções para obtenção desta formulação / Impulses describe the evolution of systems where the continuous development of a process is interrupted by abrupt changes of state. Financial markets are subject to extreme events or shocks as government changes, companies colapse, etc. Thus it seems natural to consider the action of these events in the valuation of derivative securities. The aim of this work is to obtain a formulation for the Black-Scholes equation with impulse action where the impulses can represent these shocks. We use the non-absolute integration theory in functional spaces to obtain such formulation