Побудовано математичну модель і проведено розрахунок розподілу магнітного поля у магнітотвердому шарі, що знаходиться за умов одночасної дії гармонійного за часом та постійного магнітного полів. Задача розв'язується двома способами: за допомогою методу малого параметру, а також з використанням нелінійних обчислювальних формул. Запропоновано розрахункові формули, що дозволяють у кожній вузловій точці лише за чотири звертання до правої частини диференційного рівняння отримати метод четвертого порядку точності, а також два двосторонні наближення третього порядку точності. Ключові слова. Магнітне поле, магнітотвердий шар, рівняння електродинаміки, індукція, напруженість, задача Коші, неперервні дроби, нелінійні чисельні методи, двосторонні наближення. Вступ. Математичне моделювання є одним із сучасних наукових методів дослідження прикладних задач і у багатьох важливих випадках дозволяє замінити реальний процес, а також отримати як якісну, так і кількісну характеристику модельованого процесу. При розрахунку напружено-деформівного стану тонкостінних елементів конструкцій (стержнів, пластин, оболонок), а також задач гідроакустики, кінетики, електроніки, і т.д., виникає потреба знаходити не тільки наближені розв'язки досліджуваних математичних моделей, але й гарантовану оцінку похибки. Одним з ефективних способів побудови таких наближень є ланцюгові (неперервні) дроби. Процес їх обчислень є циклічним і легко програмується на ПК. Формулювання задачі. Розглянемо електропровідний шар із магнітотверлого матеріалу товщиною l. Шар знаходиться під впливом зовнішнього магнітного поля, коли на поверхнях. (14) Тут z z z z z z z z z z z z z z z z z z z z z z z z sh z z z z z z z z z z z z z z z z , (15) де z z z z z b b z z z b z z z z COMPUTER TECHNOLOGIES OF PRINTING 2019/1 (41) 6 9 z z z z z z z z z b z z z z z z z z z z diff erential equation, allow us to obtain not only a one-step method of the fourth order of accuracy, but also two-sided calculation formulas of the third order of accuracy. Note that when using the above calculation formulas, at each node we get several approximations to the exact solution. Comparison of these approximations provides additional information for the choice of the integration step or for evaluation of the accuracy of the result.