for i := 0 to n do
           begin
             SUM[i] := 0;
             begin
               x := (n-1) + j;
               SUM[i] := SUM[i] + x;
               a[i,j] := x;
             end
           end
           for i:= 0 to n do 
             PROD := PROD * SUM[i]