Gibbs phenomenon is the particular manner how a global spectral approximation of a piecewise analytic function behaves at the jump discontinuity. The truncated spectral series has large oscillations near the jump, and the overshoot does not decay as the number of terms in the truncated series increases. There is therefore no convergence in the maximum norm, and convergence in smooth regions away from the discontinuity is also slow. In earlier work, a methodology is proposed to completely overcome this difficulty in the context of spectral collocation methods, resulting in the recovery of exponential accuracy from collocation point values of a piecewise analytic function. In this paper, we extend this methodology to handle spectral collocation methods for functions which are analytic in the open interval but have singularities at end-points. With this extension, we are able to obtain exponential accuracy from collocation point values of such functions. Similar to earlier work, the proof is constructive and uses the Gegenbauer polynomials $C^\lambda_n(x)$. The result implies that the Gibbs phenomenon can be overcome for smooth functions with endpoint singularities.