function z=zmn1(b,c,d1,d2,e,f,g)
n=9;    
a=0.005;         
k=2*pi;        
w=3*10^8*2*pi;      
E=8.854*10^(-12);    
u=4*pi*10^(-7);      
for p=1:n
    for q=1:n
        s(p,q)=exp((-i*k*e((p+b*n),(q+c*n))))/(4*pi*e((p+b*n),(q+c*n)));
    end
end
for p=1:n
    for q=1:n
         s1(p,q)=exp((-i*k*f((p+b*n),(q+c*n))))/(4*pi*f((p+b*n),(q+c*n)));
     end
 end
 for p=1:n
    for q=1:n
         s2(p,q)=exp((-i*k*g((p+b*n),(q+c*n))))/(4*pi*g((p+b*n),(q+c*n)));
     end
 end
for p=1:n
    for q=1:n
        z(p,q)=i*w*u*d1*d2*s(p,q)+(2*s(p,q)-s1(p,q)-s2(p,q))/(i*w*E);
    end
end