Analysis of numerical methods efficiency for EDFA modelling

A general numerical technique for EDFA modelling is presented. An original and efficient algorithm is used to compute the light intensity distribution in the active fibre section. A fast discretization method is applied reducing the time for ASE spectra computation. Different numerical methods have been investigated for solving ordinary differential equations (ODEs): it appears that their efficiencies depend on the system size. Runge-Kutta and Bulirsch-Stoer methods seem to be most suitable for EDFA modelling. Moreover, the Dormand and Prince method allows to save 25% of computation time.<<ETX>>