An approximate variational method for improved thermodynamics of molecular fluids