The Method of Lines (MOL) has been one of the simplest but effective technique for solving partial differential equations (PDEs) in which all but one dimension is discretized. This work describes numerical solution of 1-dimensional Schrödinger equation using the method of lines approach (MOL) where spatial dimensions were discretized using some finite difference approximation leaving the time dimension to be the only independent variable in the resulting system of initial value problems. The effect of changing the discretization size on the accuracy of the solution procedure versus changing the step size in the integration of the resulting differential equation was also studied with the incorporation of Simpson s rule function in MATLAB.