Added support for the double rotor model to the frame order matrix element simulation script.
The double rotation is constructed in the new rotation_double_xy_axes() method, and the checks for
the violation of the two torsion angles in the inside_double_rotor() method. In the main loop, the
theta, phi and sigma angles correspond to sigma1, sigma2, and nothing.