The Zeta Functions of Picard Modular Surfaces