This paper deals with the numerical solution of the Lax system L' = [B(L), L], L(0) = L-0 (*), where L-0 is a constant symmetric matrix, B(-) maps symmetric matrices into skew-symmetric matrices, and [B(L), L] is the commutator of B(L) and L. Here two different procedures, based on the approach recently proposed by Calvo, Iserles and Zanna (the MGLRK methods), are suggested. Such an approach is a computational form for the Flaschka formulation of (*). Our numerical procedures consist in solving (*) by a Runge-Kutta method, then, a single step of a Gauss-Legendre Runge-Kutta (GLRK) method may be applied to the Flaschka formulation of (*). In the first procedure we compute the approximation of the Lax system by a continuous explicit RK method, instead, the second procedure computes the approximation of the Lax system by a GLRK method (the same method used for the Flaschka system). The computational costs have been derived and compared with the ones of the MGLRK methods. Finally, several numerical tests and computational comparisons will be shown.
|Titolo:||Numerical Procedures Based on Runge-Kutta Methods for Solving Isospectral Flows|
|Data di pubblicazione:||1997|
|Digital Object Identifier (DOI):||10.1016/S0168-9274(97)00051-2|
|Appare nelle tipologie:||1.1 Articolo in rivista|