Weibull renewal equation solution

The existing solution methods for the Weibull Renewal Equation suffer from a lack of sufficient accuracy due to the singularity at the origin for some parameter values of the weibull density. The proposed method of solution provides accuracy to any desired degree of precision for all parameter values particularly in the singular range. The method utilizes a cubic spline approximation of the unknown renewal function and applies the Galerkin technique of integral equation solution. Gaussian quadratures are used to evaluate integrals. The singular nature of the integrand is handled by the Gauss-Jacobi quadrature. Results are compared with those obtained by simulation.