n = 49;
m = 20/n^2;
k = 4.0e5;
a = 2.e-2;
p = zeros(n,n);
v = p;
q1 = 45;
q2 = 10;
for i=-4:4
     for j=-4:4
        v(q1+i,q2+j) = 1.e4*exp(-0.07*(i^2+j^2));
     end
end
dom = ones(n,n);
dom(38:39,1:28)=0;
dom(20:39,29:30)=0;
dom(1:12,29:30)=0;
for i=30:n
	for j=30:n
		if i+j>=30+n
			dom(i,j)=0;
		end
	end
end
s = 2; xi = 5; xj = 5;
