function [v]=hw2p3(N,v_1)

% for k=1:N
%     v_1(k)=sin(2*pi*(k-1)/N);
% end; %vector initialized

for k=1:N/2
    ve(k)=v_1(2*k);
end;

for k=1:N/2
    vo(k)=v_1(2*(k-1)+1);
end; %split into even/odd parts

F2=zeros(N/2,N/2);

for m=1:N/2
    for n=1:N/2
        F2(m,n)=exp(-2*pi*i*(m-1)*(n-1)/(N/2));
    end; %N/2 dimensional transformation
end;
v1=F2*ve';
v2=F2*vo';

for k=1:N/2
    vf(k,1)=v2(k,1)*exp(-2*pi*i*k/N);
end; 

v=v1+vf;
