This paper examines and develops matrix methods to approximate the eigenvalues of a fourth order Sturm-Liouville problem subjected to a kind of fixed boundary conditions, furthermore, it extends the matrix methods for a kind of general boundary conditions. The idea of the methods comes from finite difference and Numerov's method as well as boundary value methods for second order regular Sturm-Liouville problems. Moreover, the determination of the correction term formulas of the matrix methods are investigated in order to obtain better approximations of the problem with fixed boundary conditions since the exact eigenvalues for q = 0 are known in this case. Finally, some numerical examples are illustrated.

This article emphasizes on fitting GARCH option pricing model with Taiwan stock market. Duan¡¦s(1995) NGARCH option pricing model is adopted. Duan solved the European option by simulation, this article follow the method and extents to pricing American option. In general, simulation approach is not convenient to solve American options as well as European options. However, the least-squares method proposed by Longstaff and Schwartz is a simple and powerful tool, so this article tests the method. The NGARCH model has parameters, and base on loglikelihood function, we fit the model with empirical observations to obtain parameters. Then we can simulate the stock prices, once stock prices are simulated, the option value can be priced. Since the article simulates the option, there should be the antithetic approaches instead of simulation. In practice, the Black-Schoels model is the benchmark for pricing European option, so this article compares the simulated European options with Black-Scholes. For American option, this article compares the simulated American options which are priced by least-squares method with trinomial tree (finite difference method).

One of the goals of Electroencephalography (EEG) is to correctly localize brain activities by the help of voltage measurements taken on scalp. However, due to computational difficulties of the problem and technological limitations, the accuracy level of the activity localization is not perfect and should be improved. To increase accuracy level of the solution, realistic, i.e. patient dependent, head models should be created. Such head models are created via assigning realistic conductivity values of head tissues onto realistic tissue positions. This study initially focuses on obtaining patient dependent spatial information from T1-weighted Magnetic Resonance (MR) head images. Existing segmentation algorithms are modified according to our needs for classifying eye tissues, white matter, gray matter, cerebrospinal fluid, skull and scalp from volumetric MR head images. Determination of patient dependent conductivity values, on the other hand, is not considered as a part of this study, and isotropic conductivity values anticipated in literature are assigned to each segmented MR-voxel accordingly. Upon completion of the tissue classification, forward problem of EEG is solved using the Finite Difference (FD) method employing a realistic head model. Utilization of the FD method aims to lower computational complexity and to simplify the process of mesh creation for brain, which has a very complex boundary. Accuracy of the employed numerical method is investigated both on Electrical Impedance Tomography (EIT) and EEG forward problems, for which analytical solutions are available. The purpose of EIT forward problem integration into this study is to evaluate reciprocal solution of the EEG forward problem.

The mortgage system has been used for many years in many countries of the world. Although the system has undergone many changes over the passing years, the basics remain the same. So, it can be thought that the earlier systems form the basis of today&rsquo / s mortgage system even though it represents some differences in practice among the countries. However, this system is very new for Turkish financial market as compared with developed countries. The aim of this study is estimating the default and prepayment risk of mortgage contract and pricing the contract in emerging markets like Turkey. In this study, a classical option pricing technique based on Cox, Ingersoll and Ross [8] is used in order to evaluate Turkish fixed-rate mortgages. In this methodology, the spot interest rate and the house price are used as state variables and it is assumed that the termination decision of mortgage is driven by a economic rationale. Under this framework, the model evaluates the embedded options, namely prepayment and default options, and the future payments which corresponds to the mortgage monthly payments. Another aim of this study is the pricing of mortgage insurance policy which has not been used yet in Turkish mortgage market but thought as potential derivative in this market. Therefore, the model used in the study also provides values for mortgage insurance policy. The partial differential equation which is derived for the mortgage, its components and mortgage insurance policy does not have closed form solutions. To cope with this problem, an explicit finite difference method is used to solve the partial differential equation. Numerical results for the value of mortgage-related assets are determined under different economic scenarios. Results obtained in the basic economic scenario show that Turkish banks apply lower contract rates as compared with the optimal ones. This observation indicates that the primary mortgage market in Turkey is still in its infancy stage. Numerical results also suggest that it is beneficial for the lenders to have mortgage default insurance, especially for the high LTV ratio mortgages.

Estimating the temperature distribution in transformer components in the design stage and during the operation is crucial since temperatures above the thermal limits of these components might seriously damage them. Thermal models are used to predict this vital information prior to actual operations. In this study, a two-dimensional, steady-state model based on the Finite Difference Method (FDM) is proposed to estimate the temperature distribution in the three-phase, SF6 gas insulatedcooled power transformer. The model can predict the temperature distribution at the specific discredited locations in the transformer successfully. This study also compares predicted temperatures of the model proposed in this study with the results of the previous study which is based on Finite Element Method (FEM) and the results of the research performed by the designers of the transformer. The results show that FDM model proposed in this study can be used to estimate the temperature distribution in the transformer with an acceptable accuracy and can be an alternative of the previous study which is based on FEM.

Optical whispering-gallery mode (WGM) cavities, which exhibit extraordinary spatial and temporal confinement of light, are one of the leading transducers for examining molecular recognition at low particle counts. With the advent of hybrid photonic-plasmonic and increasingly sophisticated forms of these resonators, the importance of supporting numerical methods has correspondingly become evident. In response, we adopt a full-vector finite difference approximation in order to solve for WGM's in terms of their field distributions, resonant wavelengths, and quality factors in the context of naturally discontinuous permittivity structure. A segmented Taylor series and alignment/rotation operator are utilized at such singularities in conjunction with arbitrarily spaced grid points. Simulations for microtoroids, with and without dielectric nanobeads, and plasmonic microdisks are demonstrated for short computation times and shown to be in agreement with data in the literature. Constricted surface plasmon polariton (SPP) WGM's are also featured within this document. The module of this thesis is devised as a keystone for composite WGM models that may guide experiments in the field. / Graduate

For a uniquely defined subset of phase space, solutions of non-linear, coupled reaction-diffusion equations may converge to heterogeneous steady states, organic in appearance. Hence, many theoretical models for pattern formation, as in the theory of morphogenesis, include the mechanics of reaction-diffusion equations. The standard method of simulation for such pattern formation models does not facilitate reproducibility of results, or the verification of convergence to a solution of the problem via the method of mesh refinement. In this thesis we explore a new methodology circumventing the aforementioned issues, which is independent of the choice of programming language. While the new method allows more control over solutions, the user is required to make more choices, which may or may not have a determining effect on the nature of resulting patterns. In an attempt to quantify the extent of the possible effects, we study heterogeneous steady states for two well known reaction-diffusion models, the Gierer-Meinhardt model and the Schnakenberg model. / Alexander Graham Bell Canada Graduate Scholarship provides financial support to high calibre scholars who are engaged in master's or doctoral programs in the natural sciences or engineering. / Natural Sciences and Engineering Research Council of Canada

Scientific computing and computer simulation play an increasingly important role in scientific investigation and engineering designs, supplementing traditional experiments, such as in automotive crash studies, global climate change, ocean modeling, medical imaging, and nuclear weapons. The numerical simulation is much cheaper than experimentation for these application areas and it can be used as the third way of science discovery beyond the experimental and theoretical analysis. However, the increasing demand of high resolution solutions of the Partial Differential Equations (PDEs) with less computational time has increased the importance for researchers and engineers to come up with efficient and scalable computational techniques that can solve very large-scale problems. In this dissertation, we build an efficient and highly accurate computational framework to solve PDEs using high order discretization schemes and multiscale multigrid method. Since there is no existing explicit sixth order compact finite difference schemes on a single scale grids, we used Gupta and Zhang’s fourth order compact (FOC) schemes on different scale grids combined with Richardson extrapolation schemes to compute the sixth order solutions on coarse grid. Then we developed an operator based interpolation scheme to approximate the sixth order solutions for every find grid point. We tested our method for 1D/2D/3D Poisson and convection-diffusion equations. We developed a multiscale multigrid method to efficiently solve the linear systems arising from FOC discretizations. It is similar to the full multigrid method, but it does not start from the coarsest level. The major advantage of the multiscale multigrid method is that it has an optimal computational cost similar to that of a full multigrid method and can bring us the converged fourth order solutions on two grids with different scales. In order to keep grid independent convergence for the multiscale multigrid method, line relaxation and plane relaxation are used for 2D and 3D convection diffusion equations with high Reynolds number, respectively. In addition, the residual scaling technique is also applied for high Reynolds number problems. To further optimize the multiscale computation procedure, we developed two new methods. The first method is developed to solve the FOC solutions on two grids using standardW-cycle structure. The novelty of this strategy is that we use the coarse level grid that will be generated in the standard geometric multigrid to solve the discretized equations and achieve higher order accuracy solution. It is more efficient and costs less CPU and memory compared with the V-cycle based multiscale multigrid method. The second method is called the multiple coarse grid computation. It is first proposed in superconvergent multigrid method to speed up the convergence. The basic idea of multigrid superconvergent method is to use multiple coarse grids to generate better correction for the fine grid solution than that from the single coarse grid. However, as far as we know, it has never been used to increase the order of solution accuracy for the fine grid. In this dissertation, we use the idea of multiple coarse grid computation to approximate the fourth order solutions on every coarse grid and fine grid. Then we apply the Richardson extrapolation for every fine grid point to get the sixth order solutions. For parallel implementation, we studied the parallelization and vectorization potential of the Gauss-Seidel relaxation by partitioning the grid space with four colors for solving 3D convection-diffusion equations. We used OpenMP to parallelize the loops in relaxation and residual computation. The numerical results show that the parallelized and the sequential implementation have the same convergence rate and the accuracy of the computed solutions.

Barrier options are becoming more popular, mainly due to the reduced cost to hold a barrier option when compared to holding a standard call/put options, but exotic options are difficult to price since the payoff functions depend on the whole path of the underlying process, rather than on its value at a specific time instant. It is a path dependent option, which implies that the payoff depends on the path followed by the price of the underlying asset, meaning that barrier options prices are especially sensitive to volatility. For basic exchange traded options, analytical prices, based on the Black-Scholes formula, can be computed. These prices are influenced by supply and demand. There is not always an analytical solution for an exotic option. Hence it is advantageous to have methods that efficiently provide accurate numerical solutions. This study gives a literature overview and compares implementation of some available numerical methods applied to barrier options. The three numerical methods that will be adapted and compared for the pricing of barrier options are: • Binomial Tree Methods • Monte-Carlo Methods • Finite Difference Methods / Thesis (MSc (Applied Mathematics))--North-West University, Potchefstroom Campus, 2013

Barrier options are becoming more popular, mainly due to the reduced cost to hold a barrier option when compared to holding a standard call/put options, but exotic options are difficult to price since the payoff functions depend on the whole path of the underlying process, rather than on its value at a specific time instant. It is a path dependent option, which implies that the payoff depends on the path followed by the price of the underlying asset, meaning that barrier options prices are especially sensitive to volatility. For basic exchange traded options, analytical prices, based on the Black-Scholes formula, can be computed. These prices are influenced by supply and demand. There is not always an analytical solution for an exotic option. Hence it is advantageous to have methods that efficiently provide accurate numerical solutions. This study gives a literature overview and compares implementation of some available numerical methods applied to barrier options. The three numerical methods that will be adapted and compared for the pricing of barrier options are: • Binomial Tree Methods • Monte-Carlo Methods • Finite Difference Methods / Thesis (MSc (Applied Mathematics))--North-West University, Potchefstroom Campus, 2013