Calculation of magnetic fields using surface currents and law of full current conservation

A method of magnetic fields calculation using the boundary integral equations has been realized. The total field of magnetic induction B in an arbitrary point is presented as the sum of two components. The current field Bc is defined by Biot-Savart law. The contribution Bs due to the iron surface currents is calculated by the boundary integral equations method involving the law of full current conservation to take into account the volume magnetic flux distribution.