Multigrid method for numerical solution of ordinary differential equations

.

We consider the initial value problem for systems of ordinary differential equations such that the solution vector can be split into subvectors and each subvector represented as a product of a scalar amplitude and a shape vector which changes slowly with time. The equations for the shape vectors can be solved with much larger time steps than those required for the original equations. The numerical results show that a substantial reduction in the computing time may be achieved