Simulations of the absorption efficiency using the discrete dipole approximation (DDA) method and taking into account the real shape of gold nanorods are reported. A dominant surface plasma band corresponding to the longitudinal resonance is observed. Its maximum position lambda(max), shifts to the red as the aspect ratio increases. The transversal dipolar and multipolar mode wavelength positions are also discussed. These data are in good agreement with previous theoretical work based on classical electrostatic predictions and assuming that gold nanorods behave as ellipsoidal particles. From the experimental point of view, good agreement with the published data for gold nanorods is obtained.