A method of full wave analysis with improved stability