We study the steady, three-dimensional creeping and viscoelastic flow around a rigid sphere subject to steady uniaxial extensional flow imposed at infinity. The viscoelastic response of the ambient fluid to the flow deformation is modeled using the second-order-fluid model, the Upper Convected Maxwell, the exponential affine Phan-Thien and Tanner and the Giesekus constitutive equations. A spherical coordinate system with origin at the center of the sphere is used to describe the flow field and the solution of the governing equations is expanded as a series in the Deborah number. The resulting sequence of differential equations is solved analytically up to second order and numerically up to fourth order in Deborah number by employing fully spectral representations for all the primary variables. In particular, Chebyshev polynomials are utilized in the radial coordinate and the Double Fourier Series in the longitudinal and latitudinal coordinates. The numerical results up to second-order agree within machine accuracy with the available analytical solutions clearly indicating the correctness and accuracy of the numerical method used here. (C) 2011 Elsevier B.V. All rights reserved.