. . "The Dormand-Prince method (commonly used by commercial solving software such as MATLAB) has been found to be susceptible to numerical error, which manifests itself in a 3.8% increase in the total power over 200 m of nonlinear interaction length." . .