Test of Grad-Shafranov condition for the experiment with h = 0.7 and p₀ = 0.375. Current density (black, y-axis on the left) and plasma pressure (blue, y-axis on the right) are plotted against the flux function A_z, for every single point in the numerical domain. Positive values of A_z refer to inside of the cusp, while negative A_z are outside the cusp. The functional form j_z(A_z) = j_max(1 + mA_z) ^− n is overplotted here.