A simple error control strategy using MATLAB BVP solvers for Yb3+-doped fiber lasers