Accurate and rapid simulation of the free boundary tokamak plasma equilibrium evolution is essential for modern plasma control, stability analysis, and scenario development. This paper presents the Free-Boundary Grad–Shafranov evolutive (FGE) code, a highly flexible and control-oriented solver designed to address the challenges posed by advanced plasma configurations across a range of devices. FGE evolved from the FBT (Hofmann 1988 Comput. Phys. Commun.48 207–21) and LIUQE (Moret et al 2015 Fusion Eng. Des.91 1–15) codes and is part of the Matlab EQuilibrium suite (MEQ), sharing many of the low-level optimized functions. It self-consistently solves the free-boundary Grad–Shafranov equation coupled with circuit equations for external conductors and models for plasma profile evolution. The code implements a fully non-linear, Newton-based framework with multiple, highly optimized solver options, state representations, and residual formulations that enable rapid computation across different simulation setups. A key capability is the self-consistent integration of various 0D and 1D current diffusion equations to model the resistive evolution of the plasma current profile as well as the ability to model plasmas with multiple magnetic axes (Doublets). Furthermore, FGE allows to linearize the plasma dynamics around a given equilibrium to generate a state-space model suitable for controller design and analysis. Numerical studies are presented demonstrating the code’s speed of convergence and validating its performance against experimental data. Furthermore, benchmarks against the RAPTOR and KINX codes for profile evolution and vertical growth rate estimates are presented highlighting FGE’s capabilities for both prediction and analysis.