function s = S(x,y)
% S : source term for diffusion equation
%
% s = S(x,y), for 0 < x,y < 1

% a source of heat in a corner

s = 0;

if norm([x y] - [0.9 0.9]) < .05
    s = 750;
end;