{
	a = vectorsmall(10 000);
	u = 1;

	for (n = 1, oo,
		h = u;	\\ initial term
		k = 0;
		forstep (v = h, #a, h,
			if (a[v]==0,
				k++;
				a[v] = n*(n-1)/2 + k;
				while (a[u],
					print (u " " a[u]);
					if (u++ > #a,
						break (3);
					);
				);

				if (k==n,
					break;	\\ row complete
				);
			);
		);
	);
}

quit