Numerical solution of isospectral flows

In this paper we are concerned with the problem of solving numerically isospec L' = [B(L),L], L(0) = L-0, where L-0 is a d x d symmetric matrix, B(L) is a skew-symmetric matrix function of L and [B,L] is the Lie bracket operator. We show that standard Runge-Kutta schemes fail in recovering the main qualitative feature of these flows, that is isospectrality, since they cannot recover arbitrary cubic conservation laws. This failure motivates us to introduce an alternative approach and establish a framework for generation of isospectral methods of arbitrarily high order.