“…These harmonics can be substantial and require a more complex computation. There have been a series of attempts to model this effect, ranging from the purely analytical methods of [2,3,4] through to hybrid finite element/analytical methods employed in [5,6,7]. Full transient solutions will provide the most accurate model, but have hitherto been avoided because, even in two dimensions, they are very slow and cumbersome.…”