function ret = p_checker(fname, rows) fid = fopen(fname); sprintf('target file => %s\n', fname) a = fscanf(fid, '%e ',[((rows * rows) + 1) inf]); fclose(fid); pack; [row_column, items] = size(a); items = items - 1; sprintf('total data => %d\n', items) t = []; %zeros(items); ret = []; for loop = 1:items v = zeros(rows, rows); num = 1; ret(loop).ts = a(num, loop); t(loop) = ret(loop).ts; for i = 1:rows for j = 1:rows num = num + 1; v(i,j) = a(num, loop); end end ret(loop).p = v; end clear a; pack; for target_index = 1:rows y = []; %zeros(items); for i = 1:items y(i) = ret(i).p(target_index, target_index); end figure; semilogy(t, y); end