  We propose and analyze a symmetric version of the Zassenhaus formula for
disentangling the exponential of two non-commuting operators. A recursive
procedure for generating the expansion up to any order is presented which also
allows one to get an enlarged domain of convergence when it is formulated for
matrices. It is shown that the approximations obtained by truncating the
infinite expansion considerably improve those arising from the standard
Zassenhaus formula.
