%***************start*******************%
clear;
clc;
lamda=1;%波長歸一化為1
l=0.5*lamda;%半波長偶極子天線的長度，通過更改l的值，可
????????????%以隨意求出其他長度的偶極子天線的電流分布圖
k=2*pi/lamda;%k為波數
a=0.001*lamda;%偶極子天線半徑
omega=2*log（l/a）;
V0=1;%V0為常數，這里取為1
step1=0.01;
step2=0.001;
x=1;
Iz=zeros（l/step1+11）;%預先為Iz矩陣分配空間

%****************求Iz*********************************************%
for?z=-0.5*l:step1:0.5*l

%***************求積分****************************%
???M=0;I=0;%賦初值???
???for?z1=-0.5*l:step2:0.5*l???
??????r=sqrt（a^2+（z-z1）^2）;
??????r1=sqrt（a^2+（0.5*l-z1）^2）;
??????F0z=cos（k*z）-cos（k*l*0.5）;
??????F0z1=cos（k*z1）-cos（k*l*0.5）;
??????I=I+（（F0z1*exp（-1i*k*r）-F0z）/r）*step2;?%化積分為累加求和
??????M=M+（（F0z1*exp（-1i*k*r1））/r1）*step2;???%化積分為累加求和???
???end
?%****************求Iz**********************