The function implements an horizontal magnetic dipole (along the x-axis) and calculates the x and z component of the magnetic field for y=0. The scr ipt allows to set a wanted value of multiple layers, with their own relative electric permittivity, r