We will transform Numerov’s method into a representation of matrix form to solve Schrödinger equation. The validity of the method is tested by applying it to calculate spectra of bottomonium. We compare our results with the experimentally observed masses and theoretically predicted results. The obtained results are found to be in good agreement with the experimental results.