Modified homotopy analysis method for solving systems of second-order BVPs