Multivariate macromodeling with stability and passivity constraints