Eddy currents play a significant role in the evolution of tokamak plasmas and must therefore be correctly taken into account in time-dependent simulations. In this paper, a computational method for solving the evolution of tokamak plasma considering eddy currents utilising VMEC (Hirshman & Whitson, Phys. Fluids, vol. 26, 1983, pp. 3553–3568), a commonly used static magnetohydrodynamic equilibrium solver, is proposed. This method is convenient since it does not modify the equilibrium solver internally and achieves convergence calculation through external processing. By allowing the components of the magnetic field to be treated separately, this method provides convergence for cases with displacements in arbitrary directions, which has been difficult to achieve with the previous methods.