Multisymplectic integrators like Preissman and six-point schemes and a semi-explicit symplectic method are applied to the coupled nonlinear Schr&ouml / dinger equations (CNLSE). Energy, momentum and additional conserved quantities are preserved by the multisymplectic integrators, which are shown using modified equations. The multisymplectic schemes are backward stable and non-dissipative. A semi-explicit method which is symplectic in the space variable and based on linear-nonlinear, even-odd splitting in time is derived. These methods are applied to the CNLSE with plane wave and soliton solutions for various combinations of the parameters of the equation. The numerical results confirm the excellent long time behavior of the conserved quantities and preservation of the shape of the soliton solutions in space and time.

In this thesis, the main perturbation problems encountered in quantum mechanics have been studied.Since the special functions and orthogonal polynomials appear very extensively in such problems, we emphasize on those topics as well. In this context, the classical quantum mechanical anharmonic oscillators described mathematically by the one-dimensional Schr&uml / odinger equation have been treated perturbatively in both finite and infinite intervals, corresponding to confined and non-confined systems, respectively.

In this thesis, a general pseudospectral formulation for a class of Sturm-Liouville eigenvalue problems is consructed. It is shown that almost all, regular or singular, Sturm-Liouville eigenvalue problems in the Schr&ouml / dinger form may be transformed into a more tractable form. This tractable form will be called here a weighted equation of hypergeometric type with a perturbation (WEHTP) since the non-weighted and unperturbed part of it is known as the equation of hypergeometric type (EHT). It is well known that the EHT has polynomial solutions which form a basis for the Hilbert space of square integrable functions. Pseudospectral methods based on this natural expansion basis are constructed to approximate the eigenvalues of WEHTP, and hence the energy eigenvalues of the Schr&ouml / dinger equation. Exemplary computations are performed to support the convergence numerically.

In this thesis, a survey on pseudospectral methods for differential equations is presented. Properties of the classical orthogonal polynomials required in this context are reviewed. Differentiation matrices corresponding to Jacobi, Laguerre,and Hermite cases are constructed. A fairly detailed investigation is made for the Hermite spectral methods, which is applied to the Schr&ouml / dinger eigenvalue equation with several potentials. A discussion of the numerical results and comparison with other methods are then introduced to deduce the effciency of the method.

In this thesis, the properties of the Bose-Einstein condensation (BEC) in low dimensions are reviewed. Three dimensional weakly interacting Bose systems are examined by the variational method. The effects of both the attractive and the repulsive interatomic forces are studied. Thomas-Fermi approximation is applied to find the ground state energy and the chemical potential. The occurrence of the BEC in low dimensional systems, is studied for ideal gases confined by both harmonic and power-law potentials. The properties of BEC in highly anisotropic trap are investigated and the conditions for reduced dimensionality are derived.

The one dimensional position dependent effective mass problem is studied by solving the Schr&ouml / dinger equation for some well known potentials, such as the deformed Hulthen, the Mie, the Kratzer, the pseudoharmonic, and the Morse potentials. Nikiforov-Uvarov method is used in the calculations to get energy eigenvalues and the corresponding wave functions exactly. By introducing a free parameter in the transformation of the wave function, the position dependent effective mass problem is reduced to the solution of the Schr&ouml / dinger equation for the constant mass case. At the same time, the deformed Hulthen potential is solved for the position dependent effective mass case by applying the method directly. The Morse potential is also solved for a mass distribution function, such that the solution can be reduced to the constant mass case.