This paper addresses a computational technique for the numerical solutions of nonlinear distributed-order fractional boundary value equations with singular coefficients. The proposed strategy is based on the shifted Legendre-series expansion and the composite midpoint quadrature rule. Moreover, a collocation technique is utilized to reduce the understudy equations to a system of nonlinear algebraic equations solved by Newton’s iteration formula. The l2 and l∞-norm errors and experimental convergence order are selected as criteria to analyze the accuracy and precision of the proposed strategy. The results of the performed numerical simulations illustrate the reliability and validity of the proposed approach.
