Запропоновано багаторазове використання принципу безперервності магнітного потоку в інтегральній формі як додаткових співвідношень до інтегральних рівнянь (ІР) для поверхневої густини простого шару магнітних зарядів та подвійного шару струмів з метою збільшення точності та зменшення часу розрахунку магнітного поля (МП) в кусково-однорідному середовищі (КОС). Наведено нові ефективні математичні моделі розрахунку МП постійних струмів (магнітів) у КОС для розрахунку та оптимізації індуктивності трифазної магнітної системи з лінійними магнітними характеристиками та силової характеристики втяжного електромагніта з урахуванням нелінійної характеристики магнітом'яких матеріалів, що входять до його складу. Обгрунтовано переваги (за точністю та часом розрахунку поля) використання поверхневої густини подвійного шару магнітних зарядів як шуканої функції в ІР, що описує МП у КОС.