1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

for i := 0 to 1440 do
begin
    for j := 0 to 1440 do
    begin
    ninside := 0;
      for k := 0 to noofevents do
      begin
        dx := i - xvalues[k];
        dy := j - yvalues[k];
          if (((dx*dx/(a*a))+(dy*dy/(b*b))) <= 1) then
          begin
            ninside := ninside + 1;
          end;
      end;
    grid[i,j] := ninside;
    if ninside > max then
    max := ninside;
    end;
end;