请大家帮忙看看这个框图的的程序怎么编，我自己是初学者，刚编了一个程序，只能运行一次，请大家帮帮忙，谢谢了

[hitzhjg]

p=8407.484658*10^(-6)
ht=297.648033*10^(-9)
E=9.03*10^10;
e_sheding=0.1;
D=2*10^(-6);
n=3.39*ones(1,50);
sigema0=930*10^6*ones(1,50);
i=1;
e=ones(1,50);
while(e(i)>e_sheding)
c(i)=sqrt(0.5*((2*n(i)+1)/(2*n(i)))^(2*n(i)-2));
C=(c).^2;
d(i)=2*sqrt(C(i)*ht*D-(C(i)*ht^2));
kesi_p(i)=0.2*d(i)/D;
fai(i)=(6*n(i)/(2*n(i)+1))*(40/(9*pi))^(1/n(i));
sigema(i)=4*p/(pi*(d(i))^2*fai(i));
K(i)= E*(sigema0(i)/E)^(1-1/(n(i)));
m(i)=log(kesi_p)/log(sigema(i)/K(i));
sigema_y(i)=E*(K(i)/E)^(m(i)/(m(i)-1));
kesi_0(i)=sigema_y(i)/E;
e(i+1)=max([abs((sigema_y(i)-sigema0(i))/sigema0(i)) abs((m(i)-n(i))/n(i))]);
sigema0(i+1)=sigema_y(i);
n(i+1)=m(i);
end
disp('适合数组为：'),disp(i-1),disp('sigema='),disp(sigema_y(i-1)),disp('n='),disp(m(i-1)),disp('误差='),disp(e(i-1))

框图在附件里面，谢谢各位MATLAB精通者帮忙，谢谢。

[barcelona]

我没有仔细看你的代码内容，只是针对你说的只能运行一次简单看了一下;
只能运行一次是因为在while循环语句中i的值一直没有变，只是1；
你可以在end之前添加一个关于i的赋值语句i=i+1;