We calculate fluctuations and correlations of conserved charges in finite temperature QCD. We
also present the equation up to order $(mu_B/T)^6$.
This allows reliable calculations of thermodynamic
quantities up to $\mu_B/T \approx 2$, which covers most
of the Beam Energy Scan program at RHIC. Our simulations
use staggered quarks with physical quark masses. All of
our results are extrapolated to the continuum limit.