Lorenzo Luini

@Greg: The problem arises due to an error in the inputs. If you open the Help of the GUI, you will read for Electric Permittivity:

"Electric permittivity

...Complex values of the electric permittivity (indicating power losses associated to the material) are allowed (the imaginary part has to be negative, according to the temporal convention exp(j*omega*t))..."

This means that the right input in your case should be:
eps = [1 -100-100*j 1]

@Alexander: The primary field is the radiation emitted by the EM source, in this case, the magnetic dipole; the secondary field originates from the interation of the primary field with the multilayered structure, i.e. consists in the reflected and transmitted components of the field. Coming to your first question (if I have understood it correctly), the method can be applied to most of the radiation patterns that are known in a closed-form solution (both concerning E and H, therefore it is possible to calculate the full EM verctor) . The tough part may be to mathematically decompose the primary field into several elementary waves. Afterwards, the propagation of simpler waves is relatively easy.

@Jveer: by the way, just comment all the "impixelinfo" commands in the script to make the GUI work correctly: you won't be able to inspect single pixel values by moving the mouse on the plots, but the GUI will work correctly.

Lorenzo, thanks, looks like the right answer now!

Hi Lorenzo, Ive noticed what I believe is a mistake:

frequency = 3.23e13
eps = [1 -100+100i 1]
zcoor = [0 -1000E-9]

returns G_ = -1.05+.165i, abs(G_) = 1.066>1. Reflection must be less than 1, so how can this be?