This version of the code contains the new implementations of the function for teh second and third order Magnus corrections.