Summary: We build a stable scheme for the solution of a deconvolution problem of the Abel integral equation type. This scheme is obtained by further developing the orthogonal polynomial-based techniques for solving the Abel integral equation of A. Ammari and A. Karoui [Inverse Probl. 26, No. 10, Article ID 105005 (2010; Zbl 1200.47017)]. More precisely, this method is based on the simultaneous use of the two families of orthogonal polynomials of the Legendre and Jacobi types. In particular, we provide an explicit formula for the computation of the Legendre expansion coefficients of the solution. This explicit formula is based on some known formulae for the exact computation of the integrals of the product of some Jacobi polynomials with the derivatives of the Legendre polynomials. Besides the explicit and the exact computation of the expansion coefficients of the solution, our proposed method has the advantage of ensuring the stability of the solution under a fairly weak condition on the functional space to which the data function belongs. Finally, we provide the reader with some numerical examples that illustrate the results of this work.