Content deleted Content added
corrected the kn terms (sinc() instead or sin()) + added indent to show Kz for loop more clearly |
|||
Line 91:
Y=y-Mmax;
kn=(1 -1/eps_r)*d_over_a.*sinc((X-Y).*d_over_a) + ((X-Y)==0)*1/eps_r;
Q(x,y)=(2*pi*(Y-1)/a + Kz).^2*kn;% -Mmax<=(Y-1)<=Mmax
end
end
| |||