To simulate plasma shape and position during tokamak discharges a new plasma evolution code was developed in MATLAB. To obtain plasma position and coil currents the plasma motion equation, Faraday’s law equations for PF coils, vacuum vessel and plasma circuits are numerically solved with input voltages. Forces acting on plasma and inductance matrices are obtained from a plasma current density distribution, which is calculated to satisfy Grad-Shafranov equation. To design feedback controllers, linear plasma models were obtained on the Globus-M (Ioffe Inst.) experimental data. Then these controllers were applied on the developed plasma evolution code. The hierarchical multi-loop control system includes cascades for controlling the plasma position and current with robust SISO PID controllers tuned by QFT. The outer cascade of plasma shape control incorporates a plasma reconstruction code, which is created on moving filaments and has sufficient accuracy and speed of response. A plasma shape MIMO controller was designed by the H∞ robust loop-shaping approach for magnetic field at X-point Br, Bz and poloidal fluxes on the plasma separatrix Ψ1, Ψ2. Full thyristor current inverter models operating in self-oscillation mode up to 3 kHz were used as original fast actuators for plasma position control. The plasma shape controller is turned on at 0.18 s of the discharge and transfers the plasma boundary to the desired position in less than 5 ms during the divertor phase. That is critically important due to comparatively short plasma discharges on the Globus-M tokamak of about 200-250 ms.