Simulations of Quantum Chromodynamics, the theory of strong interactions, currently cannot be performed at finite density. Nevertheless, extrapolations and other expansion methods allow us to access a portion of the finite-density region. I will review the state-of-the-art equation of state for strongly interacting matter from first principles. I will then discuss how to extend these calculations to larger density, approaching the relevant regime for neutron star mergers.