function [ X, Y ] = temperature_grid ( k )

%% temperature_grid defines arrays X and Y for the temperature problem.
%
%  K controls the fineness of the grid.
%
  xvec = linspace ( 0.0, 6.0, k*6+1 );
  yvec = linspace ( 0.0, 4.0, k*4+1 );

  [ X, Y ] = meshgrid ( xvec, yvec );

  return
end
