##method
https://github.com/ntpl/ntpl/wiki/Methods-Dynamic-Structure-Factor
##code
https://raw.github.com/jasonlarkin/disorder/master/matlab/m_dsf_long.m
##result
#c=0.5
#amorphous
choosing max(S(k)) for k in the [100] direction, up to k = [pi/a 0 0] where a=lattice constant of the LJ crystal conventional cubic unit cell. So, for the following systems: 4x, 8x, 10x have box lengths of L = 4a, 8a, 10*a:
