Magnus–Moldavansky method