The theory of moving frames has been used successfully to solve one dimensional (1D) variational problems invariant under a Lie group symmetry. In the one dimensional case, Noether's laws give first integrals of the Euler–Lagrange equations. In higher dimensional problems, the conservation laws do not enable the exact integration of the Euler–Lagrange system. In this paper we use the theory of moving frames to help solve, numerically, some higher dimensional variational problems, which are invariant under a Lie group action. In order to find a solution to the variational problem, we need first to solve the Euler Lagrange equations for the relevant differential invariants, and then solve a system of linear, first order, compatible, coupled partial differential equations for a moving frame, evolving on the Lie group. We demonstrate that Lie group integrators may be used in this context. We show first that the Magnus expansions on which one dimensional Lie group integrators are based, may be taken sequentially in a well defined way, at least to order 5; that is, the exact result is independent of the order of integration. We then show that efficient implementations of these integrators give a numerical solution of the equations for the frame, which is independent of the order of integration, to high order, in a range of examples. Our running example is a variational problem invariant under a linear action of \begin{document}$ SU(2) $\end{document} . We then consider variational problems for evolving curves which are invariant under the projective action of \begin{document}$ SL(2) $\end{document} and finally the standard affine action of \begin{document}$ SE(2) $\end{document} .