Solitons in Optical Fibers: Fundamentals and Applications

The following pages are a printout of a program for efficient numerical computation of pulse behavior in a dispersion map by the ODE method described in Chapter 2. Section 2.1.4. The program was written by Bell Labs colleague J rgen Gripp in the late 1990s and is based on the Maple mathematics software package. It has since proved to be an extremely handy tool for exploration of the behavior of dispersion-managed solitons, and for map design. (Many of the graphs presented in Chapter 2 were generated with this program.)
Following is an outline of the program sections:
Section 1 (Initialization) simply sets up a few general computing parameters and defines a few physical constants.
Section 2 (Fundamental Input Parameters) allows the user to specify details of the dispersion map, viz., the number of different fiber segments in the map, and for each segment, its length, its D value, effective core area, loss coefficient, and Raman gain coefficient. Usually, only the length of the + (or ?) D fiber is specified, and the program then calculates the length of the other type to provide a user-specified D. In the convention used here, the fiber segments are numbered 1, 2, 3, , from left to right, i.e., in order of increasing distance z. Provision is also made for specifying splice losses, and the free spectral range of an etalon filter if one is used. Here also the user specifies an unchirped pulse width.
Section 3 (Calculate Derivative...