A numerical study of the free-boundary axisymmetric (n = 0) ideal magnetohydrodynamical (MHD) motions of the EXTRAP device is presented. The dependence of stability on current profiles in the plasma and currents in the external conductors is investigated. Results are shown for linear growth-rates and nonlinear saturation amplitudes and their dependence on plasma radius as well as on the conducting shell radius. The interior region of the plasma is simulated by means of a Lagrangian Finite Element Method (FEM) for ideal magnetohydrodynamics. The method is based on a nonlinear variational principle for the Lagrangian description of ideal MHD. The Boundary Element Method (BEM) is used together with the Lagrangian FEM to simulate nonlinear motion of an ideal MHD plasma behaviour in a vacuum region under the influence of external magnetic fields.