Solving a kind of BVP for ODEs on heterogeneous CPU + CUDA-enabled GPU systems