img = imread('cameraman.bmp');
img = rgb2gray(img);

[r, c] = size(img);
fimg = zeros(r,c);

part1 = ( (2/r)^0.5 ) * ( (2/c)^0.5 );

for i=1:r
    for j=1:c

        if(i==0)
            k1 = 0.5^0.5;
        else
            k1 = 1;
        end

        if(j==0)
            k2 = 0.5^0.5;
        else
            k2 = 1;
        end

        part2 = k1 * k2 * img(i, j);

        cos1 = cos()