A numerical study for multiple solutions of a singular boundary value problem arising from laminar flow in a porous pipe with moving wall