The vibrational ground state and the fundamental excited states of (Ar)(13) were studied by vibrational self-consistent field (VSCF) calculations. These calculations treat the interaction between different modes through a mean potential approximation, and incorporate anharmonicity in full. The good accuracy of VSCF for such systems was demonstrated by test calculations for (AT)(3) and other clusters. The study of (Ar)(13) focused on the properties of the wave functions and the excitation energies, on the role of the coupling between the modes and on the deviation from the harmonic approximation. It was found that SCF excitation energies for the fundamental transitions differ from the harmonic values by about 25% for the softest modes, and by about 10% for the stiffest modes. Coupling between the modes, treated by SCF; was found to be much more important than the intrinsic anharmonicity of the individual modes. For the ground state, the harmonic wave function compares well with VSCF, but for the fundamental excited states appreciable differences were found. The results for a potential field expanded to fourth-order polynomial in the normal mode displacements are found to be valid, almost indentical with those for a more elaborate sixth-order polynomial expansion. The fundamental excitation frequencies computed using the Aziz-Slaman Ar-Ar pair potential are very similar, with some quantitative deviations, to the values obtained with a Lennard-Jones potential. The differences are larger for certain specific modes, and very small for the others. These calculations demonstrate the computational power of VSCF as a tool for quantum-mechanical calculations for large clusters, at the level of specific wave functions.