We improve the calculation of the equation of state for two flavor QCD by simulating on $N_t=6$ lattices at appropriate values of the couplings for the deconfinement/chiral symmetry restoration crossover. For $am_q=0.0125$ the energy density rises rapidly to approximately 1 ${\rm GeV/fm^3}$ just after the crossover($m_\pi/m_\rho\approx 0.4$ at this point). Comparing with our previous result for $N_t=4$~\cite{eos}, we find large finite $N_t$ corrections as expected from free field theory on finite lattices. We also provide formulae for extracting the speed of sound from the measured quantities.