% SYMBOLS UTILIZED IN THIS CODE
%-------------------------------------------------------------------------%
% W = Group number (Assigned by TA)
% eo = free space permittivity
% er = Relative permittivity
% uo = free space permeability
% ur = Relative permeability
% Eo = Electric field magnitude
% wavelength_o = distance between two consecutive valleys in free space
% wavelength = distance between two consecutive valleys in dielectric media
% RadianFreq = 2*pi*frequency
% Ko = propagation constant in free space
% K = propagation constant for dielectric media
% no = intrinsic impedance of free space
% n = intrinsic impedance of the dielectric media
% Tcoef = Transmission coefficient
% Rcoef = Reflection coefficient
% Nx = Number of grid points in X-Direction
% Ny = Number of grid points in Y-Direction
% Nz = Number of grid points in Z-Direction
% L = Distance coverture
% X,Y,Z = coordinate system
% R = distance between a selected point and the location of charge
% Ex = X-Component of Electric-Field
% Ey = Y-Component of Electric-Field
% Ez = Z-Component of Electric-Field
%-------------------------------------------------------------------------%
% INITIALIZATION OF INPUT PARAMETERS
% Here, all the input constants are defined
%-------------------------------------------------------------------------%
w = 1; %Use your group number here. Remember to modify before running the code!
eo = 8.85*10^-12;
er = (w);
uo = 4*pi*10^-7;
ur = 1;
wavelength_o = w*10^-2;
wavelength = wavelength_o / sqrt(ur*er);%complete the equation for wavelength in dielectric media
Ko = (2*pi)/wavelength_o
K = (2*pi)/wavelength %complete the equation for propagation constant in dielectric media
RadianFreq = Ko*(1/(sqrt(uo*eo)))
Eo = 5
no = 377;
n = no * (sqrt(ur/er)); % complete the equation for intrinsic impedance as a function of no
Rcoef = (n-no)/(n+no)
Tcoef = (2*n)/(n+no)

% X Y Z SYSTEM CALCULATION
%-------------------------------------------------------------------------%
L = 20; Nx = 5; Ny = 5; Nz = 20;
xinc=-L:2*L/(Nx-1):L; yinc=-L:2*L/(Ny-1):L; zinc=-wavelength_o: wavelength_o/(Nz-1):0;
[Xinc,Yinc,Zinc]= meshgrid(xinc,yinc,zinc);
xtra=-L:2*L/(Nx-1):L; ytra=-L:2*L/(Ny-1):L; ztra=0: wavelength_o/(Nz-1): wavelength_o;
[Xtra,Ytra,Ztra]= meshgrid(xtra,ytra,ztra);
xref=-L+1:2*L/(Nx-1):L-1;
yref=-L+1:2*L/(Ny-1):L-1;
zref=-wavelength_o+(wavelength_o/(Nz-1)/2): wavelength_o/(Nz-1):-(wavelength_o/(Nz-1)/2);
[Xref,Yref,Zref]= meshgrid(xref,yref,zref);

%-------------------------------------------------------------------------%
% PLOTTING THE ELECTRIC FIELD
%-------------------------------------------------------------------------%
for k = 1:1:21
%-------------------------------------------------------------------------%
% CALCULATION OF THE x AND y COMPONENTS OF THE ELECTRIC FIELD
%-------------------------------------------------------------------------%
set(gcf,'units','normalized','outerposition',[0.2 0.2 0.6 0.6])
t = ((k-1)/20)*4*wavelength_o;
Exinc =Eo*cos(RadianFreq*t- Ko*Zinc); % complete the equation
Eyinc =zeros(size(Zinc)); % complete the equation
Ezinc =zeros(size(Zinc)); % complete the equation
h=quiver3(Xinc,Yinc,Zinc,Exinc,Eyinc,Ezinc, 'color','b','AutoScale','off');
hold on
Extra = Tcoef*Eo*cos(RadianFreq*t- K*Ztra); % complete the equation
Eytra = zeros(size(Ztra));% complete the equation
Eztra = zeros(size(Ztra)); % complete the equation
h=quiver3(Xtra,Ytra,Ztra,Extra,Eytra,Eztra, 'color','g','AutoScale','off');
Exref = Rcoef*Eo*cos(RadianFreq*t+ Ko*Zref) ; % complete the equation
Eyref = zeros(size(Zref)); % complete the equation
Ezref = zeros(size(Zref)); % complete the

1) The solution has to be done first by hand using only a calculato...