function [ kernel ] = gaussian( n, m )
%GAUSSIAN Returns n x m gauss kernel.
%   Detailed explanation goes here

row = pascal( n, 2 );
row = abs( row( : , 1 ));

col = pascal( m, 2 );
col = abs( col( : , 1 ));

% Preallocation
kernel = zeros( n, m );

% Calculate kernel
for i = 1:n
    for j = 1:m
        kernel( i, j ) = row( i ) * col( j );
    end
end % End of for

end % End of function