In this paper we deal with construction and analysis of a multiwavelet spectral element scheme for a generalized Cauchy type problem with Caputo fractional derivative. Numerical schemes for this type of problems, often suffer from the draw-back of spurious oscillations. A common remedy is to render the problem to an equivalent integral equation. For the generalized Cauchy type problem, a corresponding integral equation is of nonlinear Volterra type. In this paper we investigate wellposedness and convergence of a stabilizing multiwavelet scheme for a, one-dimensional case (in [a, b] or [0, 1]), of this problem. Based on multiwavelets, we construct an approximation procedure for the fractional integral operator that yields a linear system of equations with sparse coefficient matrix. In this setting, choosing an appropriate threshold, the number of non-zero coefficients in the system is substantially reduced. A severe obstacle in the convergence analysis is the lack of continuous derivatives in the vicinity of the inflow/ starting boundary point. We overcome this issue through separating a J (mesh)-dependent, small, neighborhood of a (or origin) from the interval, where we only take $$L_2$$L 2 -norm. The estimate in this part relies on Chebyshev polynomials, viz. As reported by Richardson( Chebyshev interpolation for functions with endpoint singularities via exponential and double-exponential transforms, Oxford University, UK, 2012) and decreases, almost, exponentially by raising J. At the remaining part of the domain the solution is sufficiently regular to derive the desired optimal error bound. We construct such a modified scheme and analyze its wellposedness, efficiency and accuracy. The robustness of the proposed scheme is confirmed implementing numerical examples.