m = 15;
n = 10;
x = 2*(1:n);
xx = ones(m, 1)*x;
xxx = [ ones(m-1, 1); 1.5 ]*x;
y = 0.5*(1:m);
yy = y'*ones(1, n);
yyy = y'*[ 0.5 ones(1, n-1) ];
c = 2*xx.^2 + yy.^2;

figure(1)
pcolor(x, y, c);
figure(2)
flat_pcolor(x, y, c);

figure(3)
pcolor(x', y', c);
figure(4)
flat_pcolor(x', y', c);

figure(5)
pcolor(xxx, yyy, c);
figure(6)
flat_pcolor(xxx, yyy, c);

figure(7)
flat_pcolor(c);
