Numerical simulation for vertical cavity surface emitting lasers with external optical feedback