Harmonic-balance simulation of strongly nonlinear very large-size microwave circuits by inexact Newton methods