In general relativity, the motion of expanding shear-free perfect fluids is governed by the ordinary differential equation , where F is an arbitrary function from which the equation of state can be computed. A complete symmetry analysis of this differential equation is given; its solutions are classified according to this scheme, and in particular the relation to Wyman's Painlevé analysis is clarified.