%%% begin octave code %%%
 pkg load io
 w = exp(-1i*pi/3);
 P_6x6 = [ 1, -1,    1,    1,    1,    -1; ...
           1, -w,    w^2,  w^3,  w^4,  -w^5; ...
           1, -w^2,  w^4,  w^6,  w^8,  -w^10; ...
           1, -w^3,  w^6,  w^9,  w^12, -w^15; ...
           1, -w^4,  w^8,  w^12, w^16, -w^20; ...
           1, -w^5,  w^10, w^15, w^20, -w^25];
 Pi_6x6 = inv(P_6x6); % causes segfault
 D_6x6  = det(P_6x6); % causes segfault
%%% end octave code %%%
