The aim of this thesis is to address some mathematical issues on the application of pseudo spectral method on several important models in financial engineering, such as scope of application, speed of convergence, numerical stability, etc. This method has been proven superior with respect to the traditional ones especially for problems with special boundary conditions.