Recent Monte Carlo calculations from Paris and from Livermore for dense one and two component plasmas have led to systematic and accurate results for the thermodynamic properties of dense coulombic fluids. This talk will summarize the results of these numerical experiments, and the simple analytic expressions for the equation of state and other thermodynamic functions that have been obtained. The thermal energy for the one component plasma has a simple power law dependence on temperature that is identical to Monte Carlo results on strongly coupled fluids governed by l/rn potentials. A universal model for fluids governed by simple repulsive forces is suggested. For two component plasmas the ion-sphere model is shown to accurately reproduce the Monte Carlo data for the static portion of the energy. Electron screening is included using the Lindhard electric function and linear response theory. Free energy expressions have been constructed for one and two component plasmas that allow easy computation of all thermodynamic functions.