We formulate the fourth order correction to a paraxial Gaussian beam propagated along the axis of symmetry of a parabolic index lens. First we examine the evolution of a complex-source-point spherical wave (equivalent paraxially to a Gaussian beam) through the lens in a two-dimensional xz plane. Taking into account the terms of up to fourth order in aperture variables, we find a ray-optical solution to the exit beam that is represented in terms of aberration function. We also analyze the effect of the lens aberration exerted on the degradation in the quality of a Gaussian beam. The fourth order-corrected wave function derived here may be used to evaluate the quality of a Gaussian beam focused with a parabolic index lens. Further it may be applied to the case of an orthogonal system in which the index variations are different in the xz and yz planes.