The article presents a method for calculating the magnetic systems of volume models by the finite difference method, taking into account the separation of dissimilar media for an axially symmetric system. A special feature of the presented calculation method is the automation of determining the boundaries of the media of volumetric models by taking into account the magnetic resistance of adjacent nodes relative to the calculated node. The results of the program for calculating the magnetic field pattern are presented on the example of a threedimensional model of a synchronous generator with a two-core magnetic system used as part of a wind generator. Three-dimensional graphical dependences of the change in the magnetic induction value inside the stator winding coil of a synchronous generator are constructed for a rotor with an external magnetic core, for a rotor with an internal magnetic core, and for a rotor with two magnetic cores simultaneously. The calculation of the magnetic system is presented in article, the method of the example synchronous generator with dual magnetic system, determining the maximum value of magnetic induction, penetrating the winding of the synchronous generator in the direction of the r axis. Based on the picture of the three-dimensional distribution of the magnetic induction vector and the constructed three-dimensional graphs of the changes of the values of the magnetic induction vector along the r axis for an external magnetic core, the magnetic core and the double-core magnetic system is determined by the efficiency of the use of the dual magnetic system, where the magnetic induction inside the stator winding increases in average 2 times, which entails an increase in the power of the synchronous generator 4 times.