tic
for i = 1 : m
  for j = 1 : n
    a(i,j) = sin ( i * pi / m ) * cos ( j * pi /n );
  end
end
toc
