Mathematical modeling of physical processes often leads to systems of differential and algebraic equations involving quantities of interest. A computer model created based on these equations can be numerically integrated to predict future states of the system and its evolution in time. This thesis investigates current methods in numerical time-stepping schemes, identifying a number of important features needed to speed up and increase the accuracy of the solutions. The focus is on developing new methods suitable for large-scale applications with multiple physical processes, potentially with significant differences in their time-scales. Various families of new methods are introduced with special attention to multirating, low computational cost implicitness, high order of convergence, and robustness. For each family, the order condition theory is discussed and a number of examples are derived. The accuracy and stability of the methods are investigated using standard analysis techniques and numerical experiments are performed to verify the abilities of the new methods. / Doctor of Philosophy / Mathematical descriptions of physical processes are often in the form of systems of differential equations describing the time-evolution of a phenomenon. Computer simulations are realizations of these equations using well-known discretization schemes. Numerical time-stepping methods allow us to advance the state of a computer model using a sequence of time-steps. This thesis investigates current methods in time-stepping schemes, identifying a number of additional features needed to improve the speed and accuracy of simulations, and devises new methods suitable for large-scale applications where multiple processes of different physical nature drive the equations, potentially with significant differences in their time-scales. Various families of new methods are introduced with proper mathematical formulations provided for creating new ones on demand. The accuracy and stability of the methods are investigated using standard analysis techniques. These methods are then used in numerical experiments to investigate their abilities.

Doctor Scientiae - DSc / Robust Spectral Methods for Solving Option Pricing Problems by Edson Pindza PhD thesis, Department of Mathematics and Applied Mathematics, Faculty of Natural Sciences, University of the Western Cape Ever since the invention of the classical Black-Scholes formula to price the financial derivatives, a number of mathematical models have been proposed by numerous researchers in this direction. Many of these models are in general very complex, thus closed form analytical solutions are rarely obtainable. In view of this, we present a class of efficient spectral methods to numerically solve several mathematical models of pricing options. We begin with solving European options. Then we move to solve their American counterparts which involve a free boundary and therefore normally difficult to price by other conventional numerical methods. We obtain very promising results for the above two types of options and therefore we extend this approach to solve some more difficult problems for pricing options, viz., jump-diffusion models and local volatility models. The numerical methods involve solving partial differential equations, partial integro-differential equations and associated complementary problems which are used to model the financial derivatives. In order to retain their exponential accuracy, we discuss the necessary modification of the spectral methods. Finally, we present several comparative numerical results showing the superiority of our spectral methods.