Function y(r) characterising the tidal deformability of a 1.4 M_⊙ NS for values of the radial coordinate r corresponding to the outer region of the star. For accreted NS, two approaches are presented: (1) treating y(r) as a continuous function (dotted line) and (2) taking into account discontinuities (solid line). Results are compared to those obtained for catalyzed NS (dashed line). Calculations were performed for the EoS based on the BSk21 EDF.