cos_f.m

利用matlab解海伦方程并画出电流分布图

a=-1/4;b=1/4;
n=1000;
h=(b-a)/n;
xx=1/4;
x=a:h:b;
r=sqrt((xx-x).^2+(7.022/1000*1)^2);
f=cos(2*pi/1.*x).*(exp(-i*2*pi/1.*r)./r);
z=0;
for m=1:n
    z=z+(f(m)+f(m+1))*h/2;
end
z