Improved secant method for getting proper initial magnetization in transformer DC bias simulation