program-id.y.data division.working-storage section.
1 i pic 9999.
1 j pic 9.
1 k pic 9999.
1 n pic 9(7).
procedure division.a.add 1 to j
compute k=k/j
compute n=n+function mod(k,1+j)*10^j
if k>0 go a.display n(7- function log10(n):)if i<1000
add 1 to i
move i to k
move 0 to n j
go a.