We report on a continuum extrapolated result [arXiv:1309.5258] for the equation of state (EoS) of QCD with $N_f=2+1$ dynamical quark flavors. In this study, all systematics are controlled, quark masses are set to their physical values, and the continuum limit is taken using at least three lattice spacings corresponding to temporal extents up to $N_t=16$. A Symanzik improved gauge and stout-link improved staggered fermion action is used. Our results are available online [ancillary file to arXiv:1309.5258].