我的第一个for循环不循环

zeng119

第一个for循环不循环，我没法把图绘制在一起，诸君有什么办法嘛


InE = 0.01*10^(-3);  %mA
DnB = 2;  %cm^2/s
WB = 50*10^(-7);  %cm
q = 1.6e-19;  %元电荷C
NB0 = 1.5e16;  % /cm^3                    参量定义
Wb=linspace(0,50*10^(-7),50);            %为了提高竞速，将基区宽度划分50等分
for n = 0: 2: 8                              %主循环，n分别循环0、2、4、6、8

end
for value = 1:50                         %循环求积分
    syms x;
    f1 = @ (x) NB0.*(1-x./WB);
    c1(value) = integral(f1,Wb(value),WB);
end

for value = 1:50                          %循环求积分
    syms x;
    f2 = @ (x) NB0.*exp(-n.*x./WB);
    c2(value) = integral(f2,Wb(value),WB);
end

z1=vpa(c1,5);                             %5位精度
z2=vpa(c2,5);                             %5位精度
x=linspace(0,50*10^(-7),50);              %X分割，函数定义
majority1=NB0.*(1-x./WB);                 %计算线性区杂质浓度
majority2=NB0.*exp(-n.*x./WB);            %计算指数区杂质浓度
nB1=InE.*z1./(q*DnB*majority1);           %计算线性区少子浓度
nB2=InE.*z2./(q*DnB*majority2);           %计算指数区少子浓度

%绘制图形，4张子图，上边为第一种分布，下边为第二种分布

%***********线性基区少子浓度分布图********************
subplot(1,2,1);
plot(x,nB1);
title('线性基区少子浓度分布图');
xlabel('基区宽度/cm');
ylabel('少子浓度/cm^3');
hold on;

%*************缓变基区少子浓度分布图****************************
subplot(1,2,2);
plot(x,nB2);
title('缓变基区少子浓度分布图');
text(WB./10,nB2(5),['n=' num2str(n) ]);
xlabel('基区宽度/cm');
ylabel('少子浓度/cm^3');
hold on;