[ÇóÖú] ¹ØÓÚÎÞÏß´«¸ÐÆ÷ÍøÂç·ÓÉË㷨ʣÓàÄÜÁ¿µÄ·ÂÕæÎÊÌâ

etotal=0;
for r=0:1:rmax
    r
        %Çóµ±Ç°ÍøÂçµÄ×ÜÄÜÁ¿
        %Etotal(r+1)=0;
    for i=1:1:n;
        %Etotal(r+1) = S(i).E + Etotal(r+1);
        etotal = etotal + S(i).E;
    end     
ÒÔÉϵĴúÂëÒâ˼ÊÇÍøÂçÿһ´ÎÑ­»·£¬Í³¼Æ±¾ÂÖÑ­»·ÖÐÍøÂçµÄÊ£ÓàÄÜÁ¿
ÍøÂçÓÉһЩ¸ö½Úµã×é³É S(i).E±íʾ½ÚµãiµÄÄÜÁ¿

plot(1:r+1,etotal,'b-');
axis([0,4999,0,65]);
title('½ÚµãÊ£ÓàÄÜÁ¿  n=100ʱ');
xlabel('ÂÖÊý  0\leq\itX\rm\leq4999');
ylabel('½¹¶ú  0\leq\itY\rm\leq65');