[en] A numerical solution of the linear and nonlinear time-dependent Schrödinger equation is obtained, using the strong form MLPG Collocation method. Schrödinger equation is replaced by a system of coupled partial differential equa tions in terms of particle density and velocity potential, by separating the real and imaginary parts of a general solution, called a quantum hydrodynamic (QHD) equa tion, which is formally analogous to the equations of irrotational motion in a classical fluid. The approximation of the field variables is obtained with the Moving Least Squares (MLS) approximation and the implicit Crank-Nicolson scheme is used for time discretization. For the two-dimensional nonlinear Schrödinger equation, the lagging of coefficients method has been utilized to eliminate the non-linearity of the corresponding examined problem. A Type-I nodal distribution is used in order to provide convergence for the discrete Laplacian operator used at the governing equation. Numerical results are validated, comparing them with analyti cal and numerical solutions.