判断程序是否算稳的小程序

function [ ] = wentai( para )
%wentai 判断是否算稳
%   画径向中心处的电子密度和电子温度
path0 = pwd;
listing = dir([path0]); %获取路径下的所有文件（夹）
length_listing = length(listing); %获取文件（夹）数目
count_num = 0;
for k = 1:1:length_listing  %对所有文件（夹）进行循环
    if ~isempty(strfind(listing(k).name,para)) %对文件（夹）名字中含有MHz的进行操作  
        count_num = count_num +1 ;        
        para1(count_num) = str2num( listing(k).name(1:end - length(para)));
        %获取放电参数值                      
    end
end
para1 = sort(para1);
for k = 1:1:count_num 
   path1 = [path0,'\',num2str(para1(k)),'.0',para]; 
   str_loadrhoe =strcat(path1,'\','rhoem','.dat'); 
   rhoe0 = load(str_loadrhoe); 
   rhoe1 = rhoe0(:,1); %改
   str_loadte =strcat(path1,'\','Te','.dat'); %  
   Te0 = load(str_loadte); 
   Te1 = Te0(:,1); %改
   figure(k)
   subplot(2,1)
   plot(rhoe1);title('径向中心处的N_{e}');
   subplot(2,2)
   plot(Te1);title('径向中心处的T_{e}')
end
end