Philosophiae Doctor - PhD / Many problems in financial world are being modeled by means of differential equation. These problems are time dependent, highly nonlinear, stochastic and heavily depend on the previous history of time. A variety of financial products exists in the market, such as forwards, futures, swaps and options. Our main focus in this thesis is to use the numerical analysis tools to solve some option pricing problems. Depending upon the inter-relationship of the financial derivatives, the dimension of the associated problem increases drastically and hence conventional methods (for example, the finite difference methods or finite element methods) for solving them do not provide satisfactory results. To resolve this issue, we use a special class of numerical methods, namely, the mesh free methods. These methods are often better suited to cope with changes in the geometry of the domain of interest than classical discretization techniques. In this thesis, we apply these methods to solve problems that price standard and non-standard options. We then extend the proposed approach to solve Heston' volatility model. The methods in each of these cases are analyzed for stability and thorough comparative numerical results are provided. / South Africa

Radial basis function methods exhibit several very attractive properties such as a high order convergence of the approximated solution and flexibility to the domain geometry. However the method in its classical formulation becomes impractical for problems with relatively large numbers of degrees of freedom due to the ill-conditioning and dense structure of coefficient matrix. To overcome the latter issue we employ a localisation technique, namely a partition of unity method, while the former issue was previously addressed by several authors and was of less concern in this thesis. In this thesis we develop radial basis function partition of unity methods for partial differential equations arising in financial mathematics and glaciology. In the applications of financial mathematics we focus on pricing multi-asset equity and credit derivatives whose models involve several stochastic factors. We demonstrate that localised radial basis function methods are very effective and well-suited for financial applications thanks to the high order approximation properties that allow for the reduction of storage and computational requirements, which is crucial in multi-dimensional problems to cope with the curse of dimensionality. In the glaciology application we in the first place make use of the meshfree nature of the methods and their flexibility with respect to the irregular geometries of ice sheets and glaciers. Also, we exploit the fact that radial basis function methods are stated in strong form, which is advantageous for approximating velocity fields of non-Newtonian viscous liquids such as ice, since it allows to avoid a full coefficient matrix reassembly within the nonlinear iteration. In addition to the applied problems we develop a least squares radial basis function partition of unity method that is robust with respect to the node layout. The method allows for scaling to problem sizes of a few hundred thousand nodes without encountering the issue of large condition numbers of the coefficient matrix. This property is enabled by the possibility to control the coefficient matrix condition number by the rate of oversampling and the mode of refinement.

The aim of this diploma thesis is to analyze and implement selected option pricing models using statistical software. The first chapter introduces theoretical basics of options as financial instruments ideal for hedging and speculation. The second chapter constitutes the core part of this thesis since it unveils theoretical concepts of risk-neutral pricing and at the same time analyze some basic, as well as highly sophisticated option pricing models. In addition, each model is accompanied by a practical example of their effective implementation. The final chapter characterize the most widely used option trading strategies and defines the ideal expected market development linked to each strategy.

An option is a contract between two parties where the holder has the option to buy or sell some underlying asset after a predefined exercise time. Options where the holder only has the right to buy or sell at the exercise time is said to be of European-style, while options that can be exercised any time before the exercise time is said to be of American-style. Asian options are options where the payoff is determined by some average value of the underlying asset, e.g., the arithmetic or the geometric average. For arithmetic Asian options, there are no closed-form pricing formulas, and one must apply numerical methods. Several methods have been proposed and tested for Asian options. For example, the Monte Carlo method isslowforEuropean-styleAsianoptionsandnotapplicableforAmerican-styleAsian options. In contrast, the finite difference method have successfully been applied to price both European- and American-style Asian options. But from a financial point of view, one is also interested in different measures of sensitivity, called the Greeks, which are hard approximate with the finite difference method. For more accurate approximations of the Greeks, researchers have turned to the finite element method with promising results for European-style Asian options. However, the finite element method has never been applied to American-style Asian options, which still lack accurate approximations of the Greeks. Here we present a study of pricing European- and American-style Asian options using the finite element method. For European-style options, we consider two different pricing PDEs. The first equation we consider is a convection-dominated problem, which we solve by applying the so-called streamline-diffusion method. The second equation comes from modelling Asian options as options on a traded account, which we solve by using the so-called cG(1)cG(1) method. For American-style options, the model based on options on a traded account is not applicable. Therefore, we must consider the first convection-dominated problem. To handle American-style options, we study two different methods, a penalty method and the projected successive over-relaxation method. For European-style Asian options, both approaches give good results, but the model based on options on a traded account show more accurate results. For American-style Asian options, the penalty method give accurate results. Meanwhile, the projected successive over-relaxation method does not converge properly for the tested parameters. Our result is a first step towards an accurate and fast method to calculate the price and the Greeks of both European- and American-style Asian options. Because good estimations of the Greeks are crucial when hedging and trading of options, we anticipate that the ideas presented in this work can lead to new ways of trading with Asian options.

The purpose of this thesis is to compare option pricing models. We have investigated the constant volatility models Black-Scholes-Merton (BSM) and Merton’s Jump Diffusion (MJD) as well as the stochastic volatility models Heston and Bates. The data used were option prices from Microsoft, Advanced Micro Devices Inc, Walt Disney Company, and the S&P 500 index. The data was then divided into training and testing sets, where the training data was used for parameter calibration for each model, and the testing data was used for testing the model prices against prices observed on the market. Calibration of the parameters for each model were carried out using the nonlinear least-squares method. By using the calibrated parameters the price was calculated using the method of Carr and Madan. Generally it was found that the stochastic volatility models, Heston and Bates, replicated the market option prices better than both the constant volatility models, MJD and BSM for most data sets. The mean average relative percentage error for Heston and Bates was found to be 2.26% and 2.17%, respectively. Merton and BSM had a mean average relative percentage error of 6.90% and 5.45%, respectively. We therefore suggest that a stochastic volatility model is to be preferred over a constant volatility model for pricing options. / Syftet med denna tes är att jämföra prissättningsmodeller för optioner. Vi har undersökt de konstanta volatilitetsmodellerna Black-Scholes-Merton (BSM) och Merton’s Jump Diffusion (MJD) samt de stokastiska volatilitetsmodellerna Heston och Bates. Datat vi använt är optionspriser från Microsoft, Advanced Micro Devices Inc, Walt Disney Company och S&P 500 indexet. Datat delades upp i en träningsmängd och en test- mängd. Träningsdatat användes för parameterkalibrering med hänsyn till varje modell. Testdatat användes för att jämföra modellpriser med priser som observerats på mark- naden. Parameterkalibreringen för varje modell utfördes genom att använda den icke- linjära minsta-kvadratmetoden. Med hjälp av de kalibrerade parametrarna kunde priset räknas ut genom att använda Carr och Madan-metoden. Vi kunde se att de stokastiska volatilitetsmodellerna, Heston och Bates, replikerade marknadens optionspriser bättre än båda de konstanta volatilitetsmodellerna, MJD och BSM för de flesta dataseten. Medelvärdet av det relativa medelvärdesfelet i procent för Heston och Bates beräknades till 2.26% respektive 2.17%. För Merton och BSM beräknades medelvärdet av det relativa medelvärdesfelet i procent till 6.90% respektive 5.45%. Vi anser därför att en stokastisk volatilitetsmodell är att föredra framför en konstant volatilitetsmodell för att prissätta optioner.

In this paper, we study mathematical properties of a generalized bivariate Ornstein-Uhlenbeck model for financial assets. Originally introduced by Lo and Wang, this model possesses a stochastic drift term which influences the statistical properties of the asset in the real (observable) world. Furthermore, we generali- ze the model with respect to a time-dependent (but still non-random) volatility function. Although it is well-known, that drift terms - under weak regularity conditions - do not affect the behaviour of the asset in the risk-neutral world and consequently the Black-Scholes option pricing formula holds true, it makes sense to point out that these regularity conditions are fulfilled in the present model and that option pricing can be treated in analogy to the Black-Scholes case.

Companiesissuing insurance cover, in return for insurance premiums, face the payments ofclaims occurring according to a loss distribution. Hence, capital must be heldby the companies so that they can guarantee the fulfilment of the claims ofeach line of insurance. The increased incidence of insurance insolvencymotivates the birth of new legislations as the European Solvency II Directive.Companies have to determine the required amount of capital and the optimalcapital allocation across the different lines of insurance in order to keep therisk of insolvency at an adequate level. The capital allocation problem may betreated in different ways, starting from the insurance company balance sheet.Here, the running process and efficiency of four methods are evaluated andcompared so as to point out the characteristics of each of the methods. TheValue-at-Risk technique is straightforward and can be easily generated for anyloss distribution. The insolvency put option principle is easily implementableand is sensitive to the degree of default. The capital asset pricing model isone of the oldest reliable methods and still provides very helpful intermediateresults. The Myers and Read marginal capital allocation approach encouragesdiversification and introduces the concept of default value. Applications ofthe four methods to some fictive and real insurance companies are provided. Thethesis further analyses the sensitivity of those methods to changes in the economiccontext and comments how insurance companies can anticipate those changes.

>Magister Scientiae - MSc / Warrant pricing has become very crucial in the present market scenario. See, for example, M. Hanke and K. Potzelberger, Consistent pricing of warrants and traded options, Review Financial Economics 11(1) (2002) 63-77 where the authors indicate that warrants issuance affects the stock price process of the issuing company. This change in the stock price process leads to subsequent changes in the prices of options written on the issuing company's stocks. Another notable work is W.G. Zhang, W.L. Xiao and C.X. He, Equity warrant pricing model under Fractional Brownian motion and an empirical study, Expert System with Applications 36(2) (2009) 3056-3065 where the authors construct equity warrants pricing model under Fractional Brownian motion and deduce the European options pricing formula with a simple method. We study this paper in details in this mini-thesis. We also study some of the mathematical models on warrant pricing using the Black-Scholes framework. The relationship between the price of the warrants and the price of the call accounts for the dilution effect is also studied mathematically. Finally we do some numerical simulations to derive the value of warrants.

Philosophiae Doctor - PhD / Many problems in financial world are being modeled by means of differential equation. These problems are time dependent, highly nonlinear, stochastic and heavily depend on the previous history of time. A variety of financial products exists in the market, such as forwards, futures, swaps and options. Our main focus in this thesis is to use the numerical analysis tools to solve some option pricing problems. Depending upon the inter-relationship of the financial derivatives, the dimension of the associated problem increases drastically and hence conventional methods (for example, the finite difference methods or finite element methods) for solving them do not provide satisfactory results. To resolve this issue, we use a special class of numerical methods, namely, the mesh free methods. These methods are often better suited to cope with changes in the geometry of the domain of interest than classical discretization techniques. In this thesis, we apply these methods to solve problems that price standard and non-standard options. We then extend the proposed approach to solve Heston's volatility model. The methods in each of these cases are analyzed for stability and thorough comparative numerical results are provided.

Many phenomena in our world can be described as differential equations in high dimensions. However, they are notoriously challenging to solve numerically due to the exponential growth in computational cost with increasing dimensions. This thesis explores an algorithm, known as deep BSDE, for solving high-dimensional partial differential equations and applies it to finance, namely European option pricing. In addition, an implementation of the method is provided that seemingly shortens the runtime by a factor of two, compared with the results in previous studies. From the results, we can conclude that the deep BSDE method does handle high-dimensional problems well. Lastly, the thesis gives the relevant prerequisites required to be able to digest the theory from an undergraduate level.