Quarkonium in particle physics refers to meson whose constituents are a quark and its own antiquark. The famous quarkonium system is charmonium and bottomonium. Bottomonium meson has discovered recently with the ATLAS detector at the Large Hadron Collider (LHC) [1] . Bottomonium family is the set of particles that contain both a bottom quark and an anti-bottom quark but are bound together with different energies. A number of botommonium properties are well described by the quark model [2] - [8] where mesons have quantum numbers, and; and are the quantum numbers for the quark-anti- quark orbital angular momentum and their net spin angular momentum respectively [9] . The main aim of our work is to study the spectra of heavy mesons and the corresponding wave functions. Bottomonium spectra, as an example of heavy meson, are investigated by using matrix Numerov’s method [10] [11] via non-relativistic potential model [12] - [14] . However, a vast majority of numerical methods have been used to solve the Schrödinger equation (SE) numerically, for instance, Runge-Kutta method [15] , Shooting method [16] , Numerov’s method [17] , four-step exponentially fitted method [18] and the factorization method [19] . But, here we show that the matrix Numerov’s algorithm is a more efficient and fast one to achieve our goal; we hope this approximation gives the reliability features of heavy meson investigation. Moreover, the heavy-meson wave functions determined in this work can be employed to make predictions of other properties. On the other hand, the main motivation is to calculate the root mean square radius of different states for bottomonium and the numerical values of coefficient, which can be used to calculate the decay widths [20] , and differential cross sections [21] for quarkonium states. Besides, an additional aim of our work is to investigate the mass-radius dependence for states of bottomonium. The remainder of this paper is organized as follows. In Section 2, we present some characteristics properties of bottomonium mesons which in turn depend on the potential model. In Section 3, we present our main problem and its analytic solution. In Section 4, results and discussion are given. Finally in the last section, we summarize our main results and conclusions.

A non-relativistic potential model is used to study some properties of bottomonium meson by using the matrix Numerov’s method. The eigenvalues and the corresponding wave functions are found by using the same method. Then we normalized the wave functions and found the root mean square radius of bottomonium mesons by using Equation (3). Moreover, we can obtain computational values of coefficient by using Equation (4). The normalized radial wave functions for bottomonium mesons are graphically represented in Figures 1-3 respectively. For bottomonium mesons, our calculated masses and root mean square radius are reported in Table 1 in case of S, P and D States respectively. We observe that our results are in good agreement with the experiment [27] and existing theoretically predicted values [24] , which shows the validity of the used method. Some of our calculated root mean square radii are found to be in good agreement with the published one of ref. [28] . The values could be used to calculate the decay constants [29] , decay widths [29] , and differential cross sections [30] for quarkonium states with high accuracy as we used complicated potential model. The predictions about these quantities are also reported in Table 1 for bottomonium S, P, and D States respectively. Finally, we investigated the mass-radius dependence for states of bottomonium. We confirmed a leading linear relation between masses

and the radius for bottomonium and found that, with the exception of the 1S-state, the linear relation is also a good approximation for bottomonium. The relation between mass and radius in case of S-state, P-state and D-state are shown in Figure 4. Moreover, the mass-radius relation for bottomonium in case of P-state and D-state are shown in Figure 5. It is easily to behold that the S-state seems to be special as it is, in contrast to the other states. It worth noting that, a linear dependence of mass and radius might only be a good guess for states other than the S-state.

Bottomium is predicted by the standard model of particle physics but other models and techniques are required to calculate the particles’ properties. The bottomium meson spectroscopy is studied experimentally

according to PDG [27] . In this work we use the matrix Numerov’s method to obtain the radial wave functions of bottomonium meson to calculate the bottomonium root mean square and coefficient. As a remarkable result, we can point out that it is recommended to use the obtained values of and coefficient to calculate the decay widths and differential cross sections for bottomonium system. We indicate a high level of accuracy by comparing the results of decay width with available published results of decay width. Moreover, the matrix Numerov’s method [11] is tested again to obtain some mesons properties. The obtained results are in good agreement with the published data [28] . Then, the method could be safely used to solve SE. In addition, the relation between theoretical masses we obtained and values have been drawn. As pointed out previously, the fitted straight line reflects perfectly the dependence of mass and radius. Eventually, we may notice that the calculated values of and other parameters are the newer outputs where we didn’t find others for comparison. So, we are looking forward to take these data in consideration by other experimental or theoretical researchers.