function T = temperature2 ( X, Y )
%
%% TEMPERATURE2 defines a temperature field T(X,Y).
%
%  We now assume that T=0 on the boundary.
%
% X, Y, and T are arrays so use dot operators!
%
  arg1 = -0.4  * (      (X-1.0) .^ 2 + 0.7 * (Y-3.0) .^ 2 );
  arg2 = -0.2 * ( 2.0 * (X-5.0) .^ 2 + 1.5 * (Y-1.0) .^ 2 );
  T = 100.0 * exp ( arg1 ) + 80.0 * exp ( arg2 );
%
%  Force T = 0 on boundary.
%
  T(1,:) = 0.0;
  T(end,:) = 0.0;
  T(:,1) = 0.0;
  T(:,end) = 0.0;

  return
end
