``atoms8.m''

%# Ge cubic lattice constant (first in Angstrom, converted to atomic units)
a=5.66/0.52917721;

%# Standard cubic cell crystalline coordinates (as columns) for Ge
%#   ("diamond" lattice)
dl=[0.00 0.25 0.00 0.25 0.50 0.75 0.50 0.75
    0.00 0.25 0.50 0.75 0.00 0.25 0.50 0.75
    0.00 0.25 0.50 0.75 0.50 0.75 0.00 0.25];

%# Scale atom locations to actual cell size and form lattice vectors for cell
X=dl*a;
R=diag([a,a,a]);