function u = U0(x,y)
% U0 : initial values for diffusion equation
%
% u = U0(x,y), for 0 <= x,y <= 1

% something hot in a corner

u = 0;

if norm([x y] - [0.0 1.0]) < .3
    u = 150;
end;