螺线管磁场

R=24/2/1000; %螺线管直径，单位：mm
L=50/100; %螺线管长度，单位：cm
Sn=30; %每米线圈匝数
I=0.2; %电流大小，单位:A;
n=1;
for i=-1:0.00001:1
x(n)=i;
n=n+1;
end
k=1;
for k=1:(n-1)
B(k)=4pi0.0000001SnI10000/2((x(k)+L/2)/sqrt(R^2+(x(k)+L/2)2)-((x(k)-L/2)/sqrt(R^2+(x(k)-L/2)2))); %磁场强度，单位为高斯(Gs)
k=k+1;
end
hold on
plot(x,B,‘k-’)
% set(gca,‘xtick’,[-0.2:0.1:0.2])
% axis([inf,inf,0,inf])
grid on