% make Deluanay triangulation of random point set
% chosen uniformly for area measure in a disk

function [x,y,tri] = make_rand_tri_disk(n)
count=0;
while count <n
p=2*rand-1;
q=2*rand-1;
if p^2+q^2 < 1
count=count+1;
x(count)=p;
y(count)=q;
end
end % while
tri=delaunay(x,y);
return