%%%%%%%%
% WC.m %
%%%%%%%%

c = 5.361 / 5.492;

x = [
    0   0          0
    1   0          0
    0.5 sqrt(3)/2  0
    0   0          c
    1   0          c
    0.5 sqrt(3)/2  c
    ];

y = [0.5 sqrt(3)/2/3 c/2];

dx = x - [1 1 1 1 1 1]' * y;

for i = 1 : 5,
  for j = i+1 : 6,
    fprintf(1, '%f %f %f\n', norm(dx(i,:)), norm(dx(j,:)), ...
            acos(dx(i,:)*dx(j,:)'/norm(dx(i,:))/norm(dx(j,:)))*180/pi);
  end;
end;
