求助

[cbj21448]

希望各位高手帮看一下程序的问题。我是个菜鸟。
mu0=4*pi*10^(-7);N=1000;I=1;a=0.2; %设定电流线圈的各参量
l1=a;l2=0.9*a;l3=1.1*a;
c=4*mu0*N*I*a^2/pi;
dx=0.001; %设定计算位置步长
x=[0:dx:a]; %设定计算位置的范围
B1=c*(1/(4*x.^2+a^2)*sqrt(4*x.^2+2*a^2)+1/(4*(l1-x).^2+a^2)*sqrt(4*(l1-x).^2+2*a^2))
B2=c*(1/(4*x.^2+a^2)*sqrt(4*x.^2+2*a^2)+1/(4*(l2-x).^2+a^2)*sqrt(4*(l2-x).^2+2*a^2))
B3=c*(1/(4*x.^2+a^2)*sqrt(4*x.^2+2*a^2)+1/(4*(l3-x).^2+a^2)*sqrt(4*(l3-x).^2+2*a^2))
plot(x,B1,'b',x,B2,'r',x,B3,'g:')
xlabel('距离/m');ylabel('磁感应强度/T')
legend('l1=a','l2=0.9a','l3=1.1a')


公式在附件