Fig. 4. Field solution based on the eigenfunction expansion Eq. (42a,etc), Eq. (43a,etc) which can be seen to correspond closely to the benchmark field of Fig. 1, with the same values of , , . The images here only include the first 7 modes in the radial series.