n=16;
harmonic1d(n);
system('lmp_daily -in dis.lammps');

% system('bash rowChain.sh'); %and edit rowChain.sh

%findallf;F{n-1,:}=ff