neutron star mass vary between 11 . 3 − 14 . 8km. It is interesting to note that the values of the maximum neutron star mass for both the APR and TM1 EOSs are equal, but, the radius at the canonical mass of the neutron star is reasonably smaller for the APR EOS. This is due to the fact that the APR EOS is softer relative to the TM1 at intermediate densities and it becomes sti ff er at high densities as can be seen from Fig. 1. Similar is the case with SLY4 and BSR10 EOSs. We shall see in next section that these pairs of nuclear matter EOSs, for which the maximum neutron star masses are the same, yield significantly di ff erent structure for the hybrid stars. The EOSs for the CSQM in the 2SC or CFL phase are obtained within the NJL model. The input variables of the NJL model are the chemical potentials for all the quark flavours and colours in the chemical equilibrium which is given by the matrix: where, µ is the quark chemical potential, µ$_{Q}$ is the chemical potential of the electric charge equal to minus the electron chemical potential µ$_{e}$ and µ$_{3}$ and µ$_{8}$ are the colour chemical potentials asso- ciated with the two mutually commuting colour charges of the S U (3)$_{c}$ gauge group. The explicit form of the electric charge matrix Q$_{f}$ = diag$_{f}$ ( 2 $^{$_{3}$,}$− 1 $^{$_{3}$,}$− 1 $_{3}$), and for the colour charge matrices T$_{3}$ = diag$_{c}$ ( 1 $^{$_{2}$,}$− 1 $_{2}$, 0), and √$_{3}$$_{T}$$_{8}$ = diag$_{c}$ ( 1 $_{2}$, 1 $^{$_{2}$,}$− 1). In the mean-field approximation, the pressure at vanishing temperature reads as, where, σ$_{u}$$_{,}$$_{d}$$_{,}$$_{s}$ are the quark-antiquark condensates and Δ$_{c}$ are the three diquark condensates. The values of σ$_{i}$ and Δ$_{c}$ are determined using, In Eq. (2) ω$_{0}$ is the mean field expectation value for isoscalar vector likemeson ω given as [32] ω$_{0}$ = 2 $^{G$_{V}$}$〈QM | $_{ψ}$† $_{u}$ψ$_{u}$ + $_{ψ}$† $_{d}$ψ$_{d}$ + $_{ψ}$† $_{s}$ψ$_{s}$ | $^{QM}$〉 . This field modifies also the chemical potentials: µ$_{u}$$_{,}$$_{d}$$_{,}$$_{s}$ → µ$_{u}$$_{,}$$_{d}$$_{,}$$_{s}$ − ω$_{0}$ . The ϵ$_{i}$ are the dispersion relations computed by following the Ref. [3]. The ϵ$_{i}$ depend explicitly on the values of current quark masses, quark-antiquark and diquark condensates and various chemical potentials appearing in Eq. (1). The P$_{e}$ = µ 4 $_{e}$/ (12 π $^{2}$) is the contribution to 4 µ αβ ab = ( µδ αβ + µ$_{Q}$ Q αβ f ) δ$_{ab}$ + [ µ$_{3}$ ( T$_{3}$ )$_{ab}$ + µ$_{8}$ ( T$_{8}$ )$_{ab}$ ] δ$^{αβ}$, (1) p = 4 K σ$_{u}$σ$_{d}$σ$_{s}$ − 1 4 G$_{D}$ 3 ∑ c = 1 | $^{Δ$_{c}$}$| 2 − 2 G$_{S}$ 3 ∑ α = 1 σ 2 α + ω 2 0 4 G$_{V}$ + 1 2 π 2 18 ∑ i = 1 ∫ Λ 0 dkk $^{2}$| $^{ϵ$_{i}$}$| + P$_{e}$ − B (2) ∂ p ∂σ$_{i}$ = 0 (3) ∂ p ∂ Δ$_{c}$ = 0 . (4)