Investigation of mode coupling in strained and unstrained multimode step-index POFs using the Langevin equation