function array1 ( m, n )

if ( nargin < 1 )
    m = 1000;
end

if ( nargin < 2 )
    n = m;
end

tic
for i = 1 : m
  for j = 1 : n
    a(i,j) = sin ( i * pi / m ) * cos ( j * pi /n );
  end
end
toc
fprintf ( 1, '  M = %d, N = %d\n', m, n );
return
end
