Simultaneous inference of neutron star equation of state and the Hubble constant with a population of merging neutron stars