A modified formulation of the volume integral equations method for 3D magnetostatics