function patterns
    
close all;
clear all;
clc;

%  ------------- Wave Vectors  ------------

% --------  Rolls -------
Nvec   = 1;
k(1,:) = [1.0 0.0];

% % --------  Squares -------
% Nvec   = 2;
% k(1,:) = [1.0 0.0];
% k(2,:) = [0.0 1.0];

% --------  Hexagons -------
% Nvec   = 3;
% k(1,:) = [1.0 0.0];
% k(2,:) = [-0.5 sqrt(3)/2];
% k(3,:) = [-0.5 -sqrt(3)/2];
% 
[x,y] = meshgrid(-10:0.2:10,-10:0.2:10);

%  ------------- Plotting  ------------
u = 0.0;
for j=1:Nvec
    u = u + exp( (k(j,1)*x+k(j,2)*y)*i) ...
          + exp(-(k(j,1)*x+k(j,2)*y)*i);
end

figure(1);
imagesc(u);

figure(2);
surface(x,y,u);

figure(3);
contour(x,y,u);