function y = dct_2(x)
%
% Homebrew 2-D DCT Program Based on fft2
%
[n1 n2]=size(x);
%
% Make 2n1 x 2n2 function
%
intermediate=zeros(2*n1,2*n2);

intermediate=??

%
% Take 2n1 x 2n2 fft2
% And Define Twiddle Factor Matrix
%

intermediate=fft2(intermediate);

twiddle=??;

%
% Compute n1 x n2 DCT Coefficients
%
y=??;

y=real(y);
%
% Finish the Function
%