Fast algorithm for broadband macromodeling of high-order frequency responses