Numerical simulation of a Gaussian beam in a circular optical waveguide

[Aniket1]

I am trying to simulate a Gaussian beam through an optical waveguide having a circular cross-section in matlab. I am familiar with the theory of modes in an optical fiber and can analytically calculate the evolution of the beam by breaking down the beam into a sum of infinite modes.
However, I wish to know if there is any numerical technique to simulate the same process numerically.
Thank you.

 

[UltrafastPED]

Is group velocity dispersion considered? Bending radius? Fiber cladding? What are your objectives?

 

[Aniket1]

Is group velocity dispersion considered? Bending radius? Fiber cladding? What are your objectives?

Yes. Group velocity dispersion is considered. The fiber has a step index profile (constant refractive index inside the core and cladding with a small difference between their values). I'm not considering the effect of bending right now. I only want to get an idea of how to start with the simulation in terms of mathematics.
Thank you

 

[UltrafastPED]

Here is someone who has done it for you: http://www.rp-photonics.com/fiberpower_sm_launch.html

Here is the model: http://www.rp-photonics.com/fiberpower_model.html
And here: http://www.rp-photonics.com/spotlight_2008_11_08.html

You can write them for references to papers/texts that were used in the design of this software.