EDIT: I used an IFIT (inverse fast idiot transform) to change the xaxis of the plot from Hz to m. I think xlabel('Frequency [Hz]') is in my muscle memory now..
I have redone the beam fit, this time omitting the M^{2}, which I believe was superfluous. I have made the requested changes to the plot, save for the error analysis, which I am still trying to work out (the function I used for the least squares fit does not work out standard error in fit parameters). I will figure out a way to do this and amend the plot to have error bars.
