% a wall of tiles
%
% returns a gridwork
%
% USE:  contour(tiles(4,4,100,100,pi/8, 1))
%               tiles(u,v, M , N ,angle,f1,f2,X0,Y0)
%
% f2 defaults to f1     X0,Y0 default to 0


function G = dummy(u,v,M,N,angle,f1,f2,X0,Y0);

if nargin<9
  Y0=0;
end%if
if nargin<8
  X0=0;
end%if
if nargin<7
  f2 = f1;
end%if

%x = cos(2*pi*u*((0:(N-1))/N));
%y = cos(2*pi*v*((0:(M-1))/M));

[X, Y] = meshdom(u*((0:(N-1))/N),v*((0:(M-1))/M));
[X2, Y2] = project2(X-X0,Y-Y0,angle,f1,f2);

%A = cos(2*pi*X);  % non chirped
%B = cos(2*pi*Y);  % non chirped
A = cos(2*pi*X2); % chirped
B = cos(2*pi*Y2);

G = abs(A).*abs(B);