Understanding the modal response of an atomic force microscope is important for the identification of image artefacts captured using contact-mode atomic force microscopy (AFM). As the scan rate of high speed AFM increases, these modes present themselves as ever clearer noise patterns as the frequency of cantilever vibration falls under the frequency of pixel collection. An Euler-Bernoulli beam equation is used to simulate the flexural modes of the cantilever of an atomic force microscope as it images a hard surface in contact mode. Theoretical results are compared with experimental recordings taken in the high speed regime, as well as previous analytical results. It is shown that the model can capture the mode shapes and resonance properties of the first four eigenmodes.