Numerical simulations of Nonlinear Schrodinger Equation are studied using differential quadrature method based on cosine expansion. Propogation of a soliton, interaction of two solitons, birth of standing and mobile solitons and bound state solutions are simulated. The accuracy of the method (DQ) is measured using maximum error norm. The results are compared with some earlier works. The lowest two conserved quantities are computed numerically for all cases. (C) 2008 Elsevier Ltd. All rights reserved.