Matlab Practica matlab (urgente)

[necoleto]

Hola buenas, les cuento somos estudiantes de ingeniería de teleco de segundo curso y nos han mandado realizar una práctica de ondas electromagneticas la cual nadie está capacitado para realizar.
Les dejamos la práctica y algo que tenemos echo para ver si nos pueden ayudar. la fecha limite seria el día 7.
.

El enunciado de la práctica está adjuntado



f=(428000000000000:545000000000000:4);


t1=(0:10);
L = 7*10^-6;
z = (0:L);
c = 3*10^8; %velocidad de la luz en el vacio
u1 = 1;
u3 = 1;
sigma1 = 0; %conductividades
sigma2 = 0;
sigma3 = 0;
n1 = 1.54; %indice de refraccion medio 1
n3 = 1.54; %indice de refraccion medio 3
A = 0.2054*10^-6;
ep = 2.3716;
em = 0.0924;
e1 = n1.^2; %epsilon medio 1
e3 = n3.^2; %epsilon medio 3
e2 = ep+em*cos((2*pi*z)/A); %epsilon medio 2
n2 = e2.^2; %indice de refraccion medio 2
nu0= 120*pi; %impedancia en el vacio






w=2*pi.*f;
Er = exp(z- w);
Et = exp(z -w);
Ei=1;
r = Er/Ei; %reflectancia
t = Et/Ei; %transmitancia






E1 = Ei*(exp((1i*w*n1.*z)/c)+r*exp((-1i*w*n1.*z)/c));
H1 = (n1/nu0)*Ei*(exp((1i*w*n1.*z)/c)-r*exp((-1i*w*n1.*z)/c));
E3 = t*Ei*exp((t*w*n3.*(z-L)/c));
H3 = (n3/nu0)*t*Ei*exp((1i*w*n3.*(z-L)/c)); %z>L
solinit = bvpinit(linspace(0,7e-6,1000),[0,1]);
sol = bvp4c(@eode,@ebc,solinit);
plot(f,sol);








% Método para resolver las ecuaciones diferenciales


function res = ebc (ya,yb)


res = [ya(1)+ya(2)*(c/li*w*n1)-2*Ei
yb(1)-yb(2) * (c/1*w*n1)
]
end


function dxdy = eode(x,y)




dxy = [y(2)
-(ep+em*cos((2*pi.*z)/A) .* y(1)
]


end