In this project, we need access to parallel computing with sufficient capacity. This is necessary because, to simulate the statistical properties of random coupling, a large number of random unitary matrices must be generated and concatenated to model the real fiber. Using this method, the transfer function of the fiber is determined as a function of frequency and time, which contains all the important properties of the fiber under test. To analyze long-time behavior, many realizations (e.g., 25,000) of the transfer function are needed. we must then collect all realizations to obtain important statistical properties such as the probability distribution function, autocorrelation, cross-correlation, and other relevant metrics.