  In this work, we apply Adomian decomposition method for solving nonlinear
derivative-dependent doubly singular boundary value problems: $(py')'=
qf(x,y,y')$. This method is based on the modification of ADM and new two-fold
integral operator. The approximate solution is obtained in the form of series
with easily determinable components. The effectiveness of the proposed approach
is examined by considering three examples and numerical results are compared
with known results.
