%Buzási Bence
%TXEVVQ

function [kernel_matrix] = TXEVVQ_kernel_creator()

n1 = 8;
n2 = 8;
x = 1:8;
y = 1:8;
[X,Y] = meshgrid(x,y);
n = @(x) ((x~=1)*1) + ((x==1)*(1/sqrt(2)));
kernel = @(u,v) 1/4 .* n(u) .* n(v) .* cos(((2.*(X-1)+1).*(u-1).*pi) ./ ...
(2* n1)) .* cos(((2.*(Y-1)+1).*(v-1).*pi) ./ (2* n2));

figure(2)
kernel_matrix = zeros(8,8,64);
for u = 1 : n1
    for v = 1 : n2
       kernel_matrix(:,:,(u-1)*8+v) = kernel(u,v); 
    end
end

end