Equations of state for neutron-star matter at densities above nuclear-matter density are usually joined onto an equation of state for matter at lower densities calculated using a different physical treatment. Since the way of making the join is not clearly defined, we have tested the sensitivity of derived neutron-star models to the way in which the matching is made. We consider the joining between Skyrme equations of state for the higher densities and the Baym--Bethe--Pethick equation of state for densities between the neutron-drip point and nuclear-matter density.