
temp1=read_file2mat(105,1033,'log.lammps');
pe1=temp1(:,2);
ke1=temp1(:,3);
te1=pe1+ke1;
dpe1=pe1(end)-pe1(1);dke1=ke1(end)-ke1(1);
dte1=te1(end)-te1(1);
% plot(temp1(:,1),pe1,'.');hold on;
plot(temp1(:,1),te1,'.');

hold on;
%-------------------------------------%

temp2=read_file2mat(1074,2075,'log.lammps');
pe2=temp2(:,2);
ke2=temp2(:,3);
te2=pe2+ke2;
dpe2=pe2(end)-pe2(1);dke2=ke2(end)-ke2(1);
dte2=te2(end)-te2(1);
% plot(temp1(:,1),pe1,'.');hold on;
plot(temp2(:,1),te2,'.');

dte1
