A numerical method for calculating plasma equilibrium in a tokamak with an iron core transformer is presented. The influence of the iron core on plasma equilibrium is considered as an action of the currents in imaginary external poloidal coils. Coupling coefficients allow precomputation of iron core effects on equilibrium. The problem is then converted into solving the MHD equation for an air core tokamak. The calculation region is reduced to a minimum