function y=fft2_4_8(x);
% 8 point input only
% y=zeros(8,1);
% y(1:4)=x(1:4)+x(5:end);
% y(5:8)=x(1:4)-x(5:end);
% z=zeros(8,1);
% z(1:2)=y(1:2)+y(3:4);
% z(3:4)=y(1:2)-y(3:4);
% z(5:6)=y(5:6)-j*y(7:8);
% z(7:8)=y(5:6)+j*y(7:8);
% y(1)=z(1)+z(2);
% y(2)=z(1)-z(2);
% y(3)=z(3)-j*z(4);
% y(4)=z(3)+j*z(4);
% y(5)=z(5)+exp(-j*2*pi/8)*z(6);
% y(6)=z(5)-exp(-j*2*pi/8)*z(6);
% y(7)=z(7)+exp(-3*j*2*pi/8)*z(8);
% y(8)=z(7)-exp(-3*j*2*pi/8)*z(8);
% y=y([0,4,2,6,1,5,3,7]+1);
% 
y=fft64_8(x);
z=fft64_4(y);
y=fft64_2(z);
% y=y([0,4,2,6,1,5,3,7]'+1);