The numerical solution requires that the Eq. (REF ) must be integrated self consistently with the TOV equations using the boundary conditions \(y(0)=2\) , \(P(0)=P_c\)  and \(M(0)=0\)  [1]}, [2]}. From the solution of TOV equations the mass \(M\)  and radius \(R\)  of the neutron star can be extracted, while the corresponding solution of the differential Eq. (REF ) provides the value of \(y_R=y(R)\) . This parameter along with the quantity \(\beta \)  are the basic ingredients of the tidal Love number \(k_2\) .
