This paper investigates applying the approximate method; Graphical Solution Method (GSM), to theoretically solve the Time Independent Schrödinger Equation (TISE) in one dimension for a finite square well using MATLAB. With just few lines of MATLAB coding, calculating and plotting accurate eigenvalues (energy), eigenvectors (wave functions) and the bound eigenstates are possible for the finite square well of a negative potential (depth of the well) of -400 eV and a well width of 0.1 nm for an electron confined to this quantum well. These eigenvalues, eigenvectors and eigenstates are obtained and discussed. The found energy eigenvalues and states are discrete and yield physical acceptable solutions. The even and odd solutions of the TISE are also considered. The graphical solutions for the finite potential well are shown. The locations of discrete eigenvalues for even and odd solutions are also presented. These eigenvalues are tested confirming the correct eigenfunctions. The precision of these solutions depend on well width L and on the interval dx used to integrate the equation. Exact analytical solutions for this case are obtained and compared with results from the GSM. The accuracy and the convergence of the numerical results are easily checked. The results showed that the GSM can be considered as a suitable mean for determining the one dimensional solutions for the finite square well.