function f = myfunction(epsr,~,~,~,~,~,~,~,~,~,~,~,~,~,~,~,~,M,L,~,A6)
eps0=1;
mu0=1;mu1=1;mu2=1;mu3=1;mu4=1;mu5=1;
f=3e8/(500e-9);
w=2*pi*f;
r1=0.015*(500e-9);r2=0.030*(500e-9);r3=0.070*(500e-9);r4=0.085*(500e-9);r5=0.100*(500e-9);
rs=0.105*(500e-9);
k0=w*sqrt(8.854187817*(1e-12))*sqrt(mu0*4*pi*(1e-7));
k1=w*sqrt(epsr(1)*8.854187817*(1e-12))*sqrt(mu1*4*pi*(1e-7));
k2=w*sqrt(epsr(2)*8.854187817*(1e-12))*sqrt(mu2*4*pi*(1e-7));
k3=w*sqrt(epsr(3)*8.854187817*(1e-12))*sqrt(mu3*4*pi*(1e-7));
k4=w*sqrt(epsr(4)*8.854187817*(1e-12))*sqrt(mu4*4*pi*(1e-7));
k5=w*sqrt(epsr(5)*8.854187817*(1e-12))*sqrt(mu5*4*pi*(1e-7));
for x=0:4
M(1,1,x+1)=epsr(1)*besselj(x,k1*r1);M(1,2,x+1)=-epsr(2)*besselj(x,k2*r1);M(1,3,x+1)=-epsr(3)*bess