function main(n,p)

global source;

source = [5000,5000,4000];

%figure;
%hold on;

global next_mark;
global sites;

nsites = n;

%[sites,nsites,xlimts,ylimts] = makesites(nsites,p);

get_sites;
%plot_sites(sites);
nsites = length(sites);


%axis([xlimts(1) xlimts(2) ylimts(1) ylimts(2)]);

[e1,e2] = delaunay_build (sites, nsites);

next_mark = 1;

%quad_enum(e1,'visit_proc');
%quad_enum(e2,'visit_proc');

%int main(void)
%  {
%    edge_ref e;
%
%    makesites();
%
%    e = delaunay_build (sites, NSITES);
%
%    printf("Delaunay returned %x:%d \n", (e & 0xfffffffcu), (e&3));
%
%    plot_delaunay(e);
%
%    return(0);
%  }

