tic
i = ( 1 : m );
j = ( 1 : n );
a = sin ( i * pi / m )' * cos ( j * pi / n );
toc
size ( a )