One of the challenges of using Raman spectroscopy for biological applications is the inherent fluorescence generated by many biological molecules that underlies the measured spectra. This fluorescence can sometimes be several orders of magnitude more intense than the weak Raman scatter, and its presence must be minimized in order to resolve and analyze the Raman spectrum. Several techniques involving hardware and software have been devised for this purpose; these include the use of wavelength shifting, time gating, frequency-domain filtering, first- and second-order derivatives, and simple curve fitting of the broadband variation with a high-order polynomial. Of these, polynomial fitting has been found to be a simple but effective method. However, this technique typically requires user intervention and thus is time consuming and prone to variability. An automated method for fluorescence subtraction, based on a modification to least-squares polynomial curve fitting, is described. Results indicate that the presented automated method is proficient in fluorescence subtraction, repeatability, and in retention of Raman spectral lineshapes.