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

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)));        %获取放电参数值                          endendpara1 = 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}')endend