The magnetic field distribution in saturable iron part of electromagnetic energy conversion divices is defined by the nonlinear quasi-Poisson enquation that is described the electromagnetic field characteristics and satisfied the natural boundary condition. The solution of this equation is obtained by minimizing an energy functional by means of trial function that defined in triangular subregion of two-dimensional field region. As a result, the accuracy of the machine design is increased by use of its solution. In this respect, this study is developed the basic theory to analyze the magnetic flux distribution in saturable iron part and air gap of induction motor that its secondary part is short circuit by the variational principle, the minimized theory of energy functional, the application of F.E.M., and treatment of computer. As theoritical data compared with the practics, the validity of the theory in this study is supported by experimental findings.